Naoyuki Miyashita

Structures of β-Amyloid Peptide 1−40, 1−42, and 1−55—the 672−726 Fragment of APP—in a Membrane Environment with Implications for Interactions with γ-Secretase

Aggregation of Amyloid beta (Abeta) peptide has been linked to the neurodegenerative Alzheimers Disease and implicated in other amyloid diseases including cerebral amyloid angiopathy. Abeta peptide is generated by cleavage of the amyloid precursor protein (APP) by transmembrane proteases. It is crucial to determine the structures of beta-amyloid peptides in a membrane to provide a molecular basis for the cleavage mechanism. We report the structures of amyloid beta peptide (Abeta(1-40) and Abeta(1-42)) as well as the 672-726 fragment of APP (referred to as Abeta(1-55)) in a membrane environment determined by replica-exchange molecular dynamics simulation. Abeta(1-40) is found to have two helical domains A (13-22) and B(30-35) and a type I beta-turn at 23-27. The peptide is localized at the interface between membrane and solvent. Substantial fluctuations in domain A are observed. The dominant simulated tertiary structure of Abeta(1-40) is observed to be similar to the simulated Abeta(1-42) structure. However, there are differences observed in the overall conformational ensemble, as characterized by the two-dimensional free energy surfaces. The fragment of APP (Abeta(1-55)) is observed to have a long transmembrane helix. The position of the transmembrane region and ensemble of membrane structures are elucidated. The conformational transition between the transmembrane Abeta(1-55) structure, prior to cleavage, and the Abeta(1-40) structure, following cleavage, is proposed.

Transmembrane Structures of Amyloid Precursor Protein Dimer Predicted by Replica-Exchange Molecular Dynamics Simulations

Abeta peptide is an essential protein in the pathogenesis of Alzheimers disease and is derived from amyloid precursor protein (APP) in the membrane by beta- and gamma-secretase cleavage. An experimental study has shown that a pairwise replacement of Gly with Leu in APP enhances homodimerization but leads to a drastic reduction of Abeta secretion. To resolve this apparent discrepancy, we predicted the wild-type (WT) and mutant APP dimer conformations by replica-exchange molecular dynamics simulations using the implicit membrane model IMM1. The simulations illustrate large conformational differences between the WT and mutant APP fragments in the membrane. Dimerization of the WT is due to two C(alpha)-H...O hydrogen bonds between two APP fragments, whereas dimerization of the mutant is due to hydrophobic interactions. In the mutant, each APP fragment is more tilted, and the gamma-cleavage site is shifted toward the center of the membrane. This position produces a mismatch between the active site of gamma-secretase and the gamma-cleavage site of APP that might prohibit Abeta production.

Effect of bisecting GlcNAc and core fucosylation on conformational properties of biantennary complex-type N-glycans in solution.

The introduction of bisecting GlcNAc and core fucosylation in N-glycans is essential for fine functional regulation of glycoproteins. In this paper, the effect of these modifications on the conformational properties of N-glycans is examined at the atomic level by performing replica-exchange molecular dynamics (REMD) simulations. We simulate four biantennary complex-type N-glycans, namely, unmodified, two single-substituted with either bisecting GlcNAc or core fucose, and disubstituted forms. By using REMD as an enhanced sampling technique, five distinct conformers in solution, each of which is characterized by its local orientation of the Manα1-6Man glycosidic linkage, are observed for all four N-glycans. The chemical modifications significantly change their conformational equilibria. The number of major conformers is reduced from five to two and from five to four upon the introduction of bisecting GlcNAc and core fucosylation, respectively. The population change is attributed to specific inter-residue hydrogen bonds, including water-mediated ones. The experimental NMR data, including nuclear Overhauser enhancement and scalar J-coupling constants, are well reproduced taking the multiple conformers into account. Our structural model supports the concept of conformer selection, which emphasizes the conformational flexibility of N-glycans in protein-glycan interactions.

Molecular dynamics simulations of biological membranes and membrane proteins using enhanced conformational sampling algorithms

This paper reviews various enhanced conformational sampling methods and explicit/implicit solvent/membrane models, as well as their recent applications to the exploration of the structure and dynamics of membranes and membrane proteins. Molecular dynamics simulations have become an essential tool to investigate biological problems, and their success relies on proper molecular models together with efficient conformational sampling methods. The implicit representation of solvent/membrane environments is reasonable approximation to the explicit all-atom models, considering the balance between computational cost and simulation accuracy. Implicit models can be easily combined with replica-exchange molecular dynamics methods to explore a wider conformational space of a protein. Other molecular models and enhanced conformational sampling methods are also briefly discussed. As application examples, we introduce recent simulation studies of glycophorin A, phospholamban, amyloid precursor protein, and mixed lipid bilayers and discuss the accuracy and efficiency of each simulation model and method. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Functionality mapping on internal surfaces of multidrug transporter AcrB based on molecular theory of solvation: implications for drug efflux pathway.

AcrB is a membrane protein and a multidrug efflux transporter. Although the recently solved X-ray crystal structures of AcrB provide a rough sketch of how drugs efflux, the pathway and mechanism have not been completely elucidated. In this study, a ligand-mapping method based on the 3D-RISM molecular theory of solvation, which we recently developed, is applied to AcrB in order to identify the drug efflux pathway. We use a fragment-based approach as a strategy to map chemical functionality on the internal surfaces. A few multifunctional ligand-binding sites, which recognize various types of functional groups, are detected inside the porter domain. Spatial links between the multifunctional sites indicate a probable multidrug efflux pathway. The frustrated environment of the protein cavity constructed of weak interactions between ligand and protein may be a mechanism for allowing smooth transportation through the protein. Guided diffusion appears to be the main mechanism for efflux.

Structural Diversity and Changes in Conformational Equilibria of Biantennary Complex-Type N-Glycans in Water Revealed by Replica-Exchange Molecular Dynamics Simulation

Structural diversity of N-glycans is essential for specific binding to their receptor proteins. To gain insights into structural and dynamic aspects in atomic detail not normally accessible by experiment, we here perform extensive molecular-dynamics simulations of N-glycans in solution using the replica-exchange method. The simulations show that five distinct conformers exist in solution for the N-glycans with and without bisecting GlcNAc. Importantly, the population sizes of three of the conformers are drastically reduced upon the introduction of bisecting GlcNAc. This is caused by a local hydrogen-bond rearrangement proximal to the bisecting GlcNAc. These simulations show that an N-glycan modification like the bisecting GlcNAc selects a certain key (or group of keys) within the framework of the bunch of keys mechanism. Hence, the range of specific glycan-protein interactions and affinity changes need to be understood in terms of the structural diversity of glycans and the alteration of conformational equilibria by core modification.

Multidimensional umbrella sampling and replica-exchange molecular dynamics simulations for structure prediction of transmembrane helix dimers.

Structural information of a transmembrane (TM) helix dimer is useful in understanding molecular mechanisms of important biological phenomena such as signal transduction across the cell membrane. Here, we describe an umbrella sampling (US) scheme for predicting the structure of a TM helix dimer in implicit membrane using the interhelical crossing angle and the TM–TM relative rotation angles as the reaction coordinates. This scheme conducts an efficient conformational search on TM–TM contact interfaces, and its robustness is tested by predicting the structures of glycophorin A (GpA) and receptor tyrosine kinase EphA1 (EphA1) TM dimers. The nuclear magnetic resonance (NMR) structures of both proteins correspond to the global free‐energy minimum states in their free‐energy landscapes. In addition, using the landscape of GpA as a reference, we also examine the protocols of temperature replica‐exchange molecular dynamics (REMD) simulations for structure prediction of TM helix dimers in implicit membrane. A wide temperature range in REMD simulations, for example, 250–1000 K, is required to efficiently obtain a free‐energy landscape consistent with the US simulations. The interhelical crossing angle and the TM–TM relative rotation angles can be used as reaction coordinates in multidimensional US and be good measures for conformational sampling of REMD simulations.

Conformational flexibility of N-glycans in solution studied by REMD simulations

Protein–glycan recognition regulates a wide range of biological and pathogenic processes. Conformational diversity of glycans in solution is apparently incompatible with specific binding to their receptor proteins. One possibility is that among the different conformational states of a glycan, only one conformer is utilized for specific binding to a protein. However, the labile nature of glycans makes characterizing their conformational states a challenging issue. All-atom molecular dynamics (MD) simulations provide the atomic details of glycan structures in solution, but fairly extensive sampling is required for simulating the transitions between rotameric states. This difficulty limits application of conventional MD simulations to small fragments like di- and tri-saccharides. Replica-exchange molecular dynamics (REMD) simulation, with extensive sampling of structures in solution, provides a valuable way to identify a family of glycan conformers. This article reviews recent REMD simulations of glycans carried out by us or other research groups and provides new insights into the conformational equilibria of N-glycans and their alteration by chemical modification. We also emphasize the importance of statistical averaging over the multiple conformers of glycans for comparing simulation results with experimental observables. The results support the concept of “conformer selection” in protein–glycan recognition.

Energetics of the Presequence-Binding Poses in Mitochondrial Protein Import Through Tom20

Tom20 is located at the outer membrane of mitochondria and functions as a receptor for the N-terminal presequence of mitochondrial-precursor proteins. Recently, three atomic structures of the Tom20-presequence complex were determined using X-ray crystallography and classified into A-, M-, and Y-poses in terms of their presequence-binding modes. Combined with biochemical and NMR data, a dynamic-equilibrium model between the three poses has been proposed. To investigate this mechanism in further detail, we performed all-atom molecular dynamics (MD) simulations and replica-exchange MD (REMD) simulations of the Tom20-presequence complex in explicit water. In the REMD simulations, one major distribution and another minor one were observed in the converged free-energy landscape at 300 K. In the major distribution, structures similar to A- and M-poses exist, whereas those similar to Y-pose are located in the minor one, suggesting that A-pose in solution is more stable than Y-pose. A k-means clustering algorithm revealed a new pose not yet obtained by X-ray crystallography. This structure has double salt bridges between Arg14 in the presequence and Glu78 or Glu79 in Tom20 and can explain the binding affinity of the complex in previous pull-down assay experiments. Structural clustering and analyses of contacts between Tom20 and the presequence suggest smooth conformational changes from Y- to A-poses through low activation barriers. M-pose lies between Y- and A-poses as a metastable state. The REMD simulations thus provide insights into the energetics of the multiple-binding forms and help to detail the progressive conformational states in the dynamic-equilibrium model based on the experimental data.

The Intertransmembrane Region of Kaposi's Sarcoma-Associated Herpesvirus Modulator of Immune Recognition 2 Contributes to B7-2 Downregulation

ABSTRACT Kaposis sarcoma-associated herpesvirus (KSHV), a human tumor virus, encodes two homologous membrane-associated E3 ubiquitin ligases, modulator of immune recognition 1 (MIR1) and MIR2, to evade host immunity. Both MIR1 and MIR2 downregulate the surface expression of major histocompatibility complex class I (MHC I) molecules through ubiquitin-mediated endocytosis followed by lysosomal degradation. Since MIR2 additionally downregulates a costimulatory molecule (B7-2) and an integrin ligand (intercellular adhesion molecule 1 [ICAM-1]), MIR2 is thought to be a more important molecule for immune evasion than MIR1; however, the molecular basis of the MIR2 substrate specificity remains unclear. To address this issue, we determined which regions of B7-2 and MIR2 are required for MIR2-mediated B7-2 downregulation. Experiments with chimeras made by swapping domains between human B7-2 and CD8α, a non-MIR2 substrate, and between MIR1 and MIR2 demonstrated a significant contribution of the juxtamembrane (JM) region of B7-2 and the intertransmembrane (ITM) region of MIR2 to MIR2-mediated downregulation. Structure prediction and mutagenesis analyses indicate that Phe119 and Ser120 in the MIR2 ITM region and Asp244 in the B7-2 JM region contribute to the recognition of B7-2 by MIR2. This finding provides new insight into the molecular basis of substrate recognition by MIR family members.