Block of NMDA receptor channels by endogenous neurosteroids: implications for the agonist induced conformational states of the channel vestibule

N-methyl-D-aspartate receptors (NMDARs) mediate synaptic plasticity, and their dysfunction is implicated in multiple brain disorders. NMDARs can be allosterically modulated by numerous compounds, including endogenous neurosteroid pregnanolone sulfate. Here, we identify the molecular basis of the use-dependent and voltage-independent inhibitory effect of neurosteroids on NMDAR responses. The site of action is located at the extracellular vestibule of the receptor's ion channel pore and is accessible after receptor activation. Mutations in the extracellular vestibule in the SYTANLAAF motif disrupt the inhibitory effect of negatively charged steroids. In contrast, positively charged steroids inhibit mutated NMDAR responses in a voltage-dependent manner. These results, in combination with molecular modeling, characterize structure details of the open configuration of the NMDAR channel. Our results provide a unique opportunity for the development of new therapeutic neurosteroid-based ligands to treat diseases associated with dysfunction of the glutamate system

Use of molecular modelling to probe the mechanism of the nucleoside transporter NupG.

Nucleosides play key roles in biology as precursors for salvage pathways of nucleotide synthesis. Prokaryotes import nucleosides across the cytoplasmic membrane by proton- or sodium-driven transporters belonging to the Concentrative Nucleoside Transporter (CNT) family or the Nucleoside:H(+) Symporter (NHS) family of the Major Facilitator Superfamily. The high resolution structure of a CNT from Vibrio cholerae has recently been determined, but no similar structural information is available for the NHS family. To gain a better understanding of the molecular mechanism of nucleoside transport, in the present study the structures of two conformations of the archetypical NHS transporter NupG from Escherichia coli were modelled on the inward- and outward-facing conformations of the lactose transporter LacY from E. coli, a member of the Oligosaccharide:H(+) Symporter (OHS) family. Sequence alignment of these distantly related proteins (∼ 10% sequence identity), was facilitated by comparison of the patterns of residue conservation within the NHS and OHS families. Despite the low sequence similarity, the accessibilities of endogenous and introduced cysteine residues to thiol reagents were found to be consistent with the predictions of the models, supporting their validity. For example C358, located within the predicted nucleoside binding site, was shown to be responsible for the sensitivity of NupG to inhibition by p-chloromercuribenzene sulphonate. Functional analysis of mutants in residues predicted by the models to be involved in the translocation mechanism, including Q261, E264 and N228, supported the hypothesis that they play important roles, and suggested that the transport mechanisms of NupG and LacY, while different, share common features

Simultaneous alignment and folding of protein sequences

Accurate comparative analysis tools for low-homology proteins remains a difficult challenge in computational biology, especially sequence alignment and consensus folding problems. We presentpartiFold-Align, the first algorithm for simultaneous alignment and consensus folding of unaligned protein sequences; the algorithm’s complexity is polynomial in time and space. Algorithmically,partiFold-Align exploits sparsity in the set of super-secondary structure pairings and alignment candidates to achieve an effectively cubic running time for simultaneous pairwise alignment and folding. We demonstrate the efficacy of these techniques on transmembrane β-barrel proteins, an important yet difficult class of proteins with few known three-dimensional structures. Testing against structurally derived sequence alignments,partiFold-Align significantly outperforms state-of-the-art pairwise sequence alignment tools in the most difficult low sequence homology case and improves secondary structure prediction where current approaches fail. Importantly, partiFold-Align requires no prior training. These general techniques are widely applicable to many more protein families. partiFold-Align is available at http://partiFold.csail.mit.edu

Structure-based statistical analysis of transmembrane helices

Recent advances in determination of the high-resolution structure of membrane proteins now enable analysis of the main features of amino acids in transmembrane (TM) segments in comparison with amino acids in water-soluble helices. In this work, we conducted a large-scale analysis of the prevalent locations of amino acids by using a data set of 170 structures of integral membrane proteins obtained from the MPtopo database and 930 structures of water-soluble helical proteins obtained from the protein data bank. Large hydrophobic amino acids (Leu, Val, Ile, and Phe) plus Gly were clearly prevalent in TM helices whereas polar amino acids (Glu, Lys, Asp, Arg, and Gln) were less frequent in this type of helix. The distribution of amino acids along TM helices was also examined. As expected, hydrophobic and slightly polar amino acids are commonly found in the hydrophobic core of the membrane whereas aromatic (Trp and Tyr), Pro, and the hydrophilic amino acids (Asn, His, and Gln) occur more frequently in the interface regions. Charged amino acids are also statistically prevalent outside the hydrophobic core of the membrane, and whereas acidic amino acids are frequently found at both cytoplasmic and extra-cytoplasmic interfaces, basic amino acids cluster at the cytoplasmic interface. These results strongly support the experimentally demonstrated biased distribution of positively charged amino acids (that is, the so-called the positive-inside rule) with structural data

A candidate ion-retaining state in the inward-facing conformation of sodium/galactose symporter: Clues from atomistic simulations

The recent Vibrio parahaemolyticus sodium/galactose (vSGLT) symporter crystal structure captures the protein in an inward-facing substrate-bound conformation, with the sodium ion placed, by structural alignment, in a site equivalent to the Na2 site of the leucine transporter (LeuT). A recent study, based on molecular dynamics simulations, showed that the sodium ion spontaneously leaves its initial position diffusing outside vSGLT, toward the intracellular space. This suggested that the crystal structure corresponds to an ion-releasing state of the transporter. Here, using metadynamics, we identified a more stable Na+ binding site corresponding to a putative ion-retaining state of the transporter. In addition, our simulations, consistently with mutagenesis studies, highlight the importance of D189 that, without being one of the NA(+)-coordinating residues, regulates its binding/release

Determining the Orientation of Protegrin-1 in DLPC Bilayers Using an Implicit Solvent-Membrane Model

Continuum models that describe the effects of solvent and biological membrane molecules on the structure and behavior of antimicrobial peptides, holds a promise to improve our understanding of the mechanisms of antimicrobial action of these peptides. In such methods, a lipid bilayer model membrane is implicitly represented by multiple layers of relatively low dielectric constant embedded in a high dielectric aqueous solvent, while an antimicrobial peptide is accounted for by a dielectric cavity with fixed partial charge at the center of each one of its atoms. In the present work, we investigate the ability of continuum approaches to predict the most probable orientation of the β-hairpin antimicrobial peptide Protegrin-1 (PG-1) in DLPC lipid bilayers by calculating the difference in the transfer free energy from an aqueous environment to a membrane-water environment for multiple orientations. The transfer free energy is computed as a sum of two terms; polar/electrostatic and non-polar. They both include energetic and entropic contributions to the free energy. We numerically solve the Poisson-Boltzmann equation to calculate the electrostatic contribution to the transfer free energy, while the non-polar contribution to the free energy is approximated using a linear solvent accessible surface area relationships. The most probable orientation of PG-1 is that with the lowest relative transfer free energy. Our simulation results indicate that PG-1 assumes an oblique orientation in DLPC lipid bilayers. The predicted most favorable orientation was with a tilt angle of 19°, which is in qualitative agreement with the experimentally observed orientations derived from solid-state NMR data

Structure and dynamics of the active Gs-coupled human secretin receptor

The class B secretin GPCR (SecR) has broad physiological effects, with target potential for treatment of metabolic and cardiovascular disease. Molecular understanding of SecR binding and activation is important for its therapeutic exploitation. We combined cryo-electron microscopy, molecular dynamics, and biochemical cross-linking to determine a 2.3 Å structure, and interrogate dynamics, of secretin bound to the SecR:Gs complex. SecR exhibited a unique organization of its extracellular domain (ECD) relative to its 7-transmembrane (TM) core, forming more extended interactions than other family members. Numerous polar interactions formed between secretin and the receptor extracellular loops (ECLs) and TM helices. Cysteine-cross-linking, cryo-electron microscopy multivariate analysis and molecular dynamics simulations revealed that interactions between peptide and receptor were dynamic, and suggested a model for initial peptide engagement where early interactions between the far N-terminus of the peptide and SecR ECL2 likely occur following initial binding of the peptide C-terminus to the ECD

APBSmem: A Graphical Interface for Electrostatic Calculations at the Membrane

Electrostatic forces are one of the primary determinants of molecular interactions. They help guide the folding of proteins, increase the binding of one protein to another and facilitate protein-DNA and protein-ligand binding. A popular method for computing the electrostatic properties of biological systems is to numerically solve the Poisson-Boltzmann (PB) equation, and there are several easy-to-use software packages available that solve the PB equation for soluble proteins. Here we present a freely available program, called APBSmem, for carrying out these calculations in the presence of a membrane. The Adaptive Poisson-Boltzmann Solver (APBS) is used as a back-end for solving the PB equation, and a Java-based graphical user interface (GUI) coordinates a set of routines that introduce the influence of the membrane, determine its placement relative to the protein, and set the membrane potential. The software Jmol is embedded in the GUI to visualize the protein inserted in the membrane before the calculation and the electrostatic potential after completing the computation. We expect that the ease with which the GUI allows one to carry out these calculations will make this software a useful resource for experimenters and computational researchers alike. Three examples of membrane protein electrostatic calculations are carried out to illustrate how to use APBSmem and to highlight the different quantities of interest that can be calculated

Extracellular loops 2 and 3 of the calcitonin receptor selectively modify agonist binding and efficacy.

Class B peptide hormone GPCRs are targets for the treatment of major chronic disease. Peptide ligands of these receptors display biased agonism and this may provide future therapeutic advantage. Recent active structures of the calcitonin (CT) and glucagon-like peptide-1 (GLP-1) receptors reveal distinct engagement of peptides with extracellular loops (ECLs) 2 and 3, and mutagenesis of the GLP-1R has implicated these loops in dynamics of receptor activation. In the current study, we have mutated ECLs 2 and 3 of the human CT receptor (CTR), to interrogate receptor expression, peptide affinity and efficacy. Integration of these data with insights from the CTR and GLP-1R active structures, revealed marked diversity in mechanisms of peptide engagement and receptor activation between the CTR and GLP-1R. While the CTR ECL2 played a key role in conformational propagation linked to Gs/cAMP signalling this was mechanistically distinct from that of GLP-1R ECL2. Moreover, ECL3 was a hotspot for distinct ligand- and pathway- specific effects, and this has implications for the future design of biased agonists of class B GPCRs