Proton-Binding Sites of Acid-Sensing Ion Channel 1

Acid-sensing ion channels (ASICs) are proton-gated cation channels that exist throughout the mammalian central and peripheral nervous systems. ASIC1 is the most abundant of all the ASICs and is likely to modulate synaptic transmission. Identifying the proton-binding sites of ASCI1 is required to elucidate its pH-sensing mechanism. By using the crystal structure of ASIC1, the protonation states of each titratable site of ASIC1 were calculated by solving the Poisson-Boltzmann equation under conditions wherein the protonation states of all these sites are simultaneously in equilibrium. Four acidic-acidic residue pairs—Asp238-Asp350, Glu220-Asp408, Glu239-Asp346, and Glu80-Glu417—were found to be highly protonated. In particular, the Glu80-Glu417 pair in the inner pore was completely protonated and possessed 2 H+, implying its possible importance as a proton-binding site. The pKa of Glu239, which forms a pair with a possible pH-sensing site Asp346, differs among each homo-trimer subunit due to the different H-bond pattern of Thr237 in the different protein conformations of the subunits. His74 possessed a pKa of ≈6–7. Conservation of His74 in the proton-sensitive ASIC3 that lacks a residue corresponding to Asp346 may suggest its possible pH-sensing role in proton-sensitive ASICs

Order through Disorder: Hyper-Mobile C-Terminal Residues Stabilize the Folded State of a Helical Peptide. A Molecular Dynamics Study

Conventional wisdom has it that the presence of disordered regions in the three-dimensional structures of polypeptides not only does not contribute significantly to the thermodynamic stability of their folded state, but, on the contrary, that the presence of disorder leads to a decrease of the corresponding proteins' stability. We have performed extensive 3.4 µs long folding simulations (in explicit solvent and with full electrostatics) of an undecamer peptide of experimentally known helical structure, both with and without its disordered (four residue long) C-terminal tail. Our simulations clearly indicate that the presence of the apparently disordered (in structural terms) C-terminal tail, increases the thermodynamic stability of the peptide's folded (helical) state. These results show that at least for the case of relatively short peptides, the interplay between thermodynamic stability and the apparent structural stability can be rather subtle, with even disordered regions contributing significantly to the stability of the folded state. Our results have clear implications for the understanding of peptide energetics and the design of foldable peptides

A mathematical and computational review of Hartree-Fock SCF methods in Quantum Chemistry

We present here a review of the fundamental topics of Hartree-Fock theory in Quantum Chemistry. From the molecular Hamiltonian, using and discussing the Born-Oppenheimer approximation, we arrive to the Hartree and Hartree-Fock equations for the electronic problem. Special emphasis is placed in the most relevant mathematical aspects of the theoretical derivation of the final equations, as well as in the results regarding the existence and uniqueness of their solutions. All Hartree-Fock versions with different spin restrictions are systematically extracted from the general case, thus providing a unifying framework. Then, the discretization of the one-electron orbitals space is reviewed and the Roothaan-Hall formalism introduced. This leads to a exposition of the basic underlying concepts related to the construction and selection of Gaussian basis sets, focusing in algorithmic efficiency issues. Finally, we close the review with a section in which the most relevant modern developments (specially those related to the design of linear-scaling methods) are commented and linked to the issues discussed. The whole work is intentionally introductory and rather self-contained, so that it may be useful for non experts that aim to use quantum chemical methods in interdisciplinary applications. Moreover, much material that is found scattered in the literature has been put together here to facilitate comprehension and to serve as a handy reference.Comment: 64 pages, 3 figures, tMPH2e.cls style file, doublesp, mathbbol and subeqn package

Molecular dynamics simulations and in silico peptide ligand screening of the Elk-1 ETS domain

Background: The Elk-1 transcription factor is a member of a group of proteins called ternary complex factors, which serve as a paradigm for gene regulation in response to extracellular signals. Its deregulation has been linked to multiple human diseases including the development of tumours. The work herein aims to inform the design of potential peptidomimetic compounds that can inhibit the formation of the Elk-1 dimer, which is key to Elk-1 stability. We have conducted molecular dynamics simulations of the Elk-1 ETS domain followed by virtual screening. Results: We show the ETS dimerisation site undergoes conformational reorganisation at the a1b1 loop. Through exhaustive screening of di- and tri-peptide libraries against a collection of ETS domain conformations representing the dynamics of the loop, we identified a series of potential binders for the Elk-1 dimer interface. The di-peptides showed no particular preference toward the binding site; however, the tri-peptides made specific interactions with residues: Glu17, Gln18 and Arg49 that are pivotal to the dimer interface. Conclusions: We have shown molecular dynamics simulations can be combined with virtual peptide screening to obtain an exhaustive docking protocol that incorporates dynamic fluctuations in a receptor. Based on our findings, we suggest experimental binding studies to be performed on the 12 SILE ranked tri-peptides as possible compounds for the design of inhibitors of Elk-1 dimerisation. It would also be reasonable to consider the score ranked tri-peptides as a comparative test to establish whether peptide size is a determinant factor of binding to the ETS domain

Visualization of Allostery in P-Selectin Lectin Domain Using MD Simulations

Allostery of P-selectin lectin (Lec) domain followed by an epithelial growth factor (EGF)-like domain is essential for its biological functionality, but the underlying pathways have not been well understood. Here the molecular dynamics simulations were performed on the crystallized structures to visualize the dynamic conformational change for state 1 (S1) or state 2 (S2) Lec domain with respective bent (B) or extended (E) EGF orientation. Simulations illustrated that both S1 and S2 conformations were unable to switch from one to another directly. Instead, a novel S1' conformation was observed from S1 when crystallized B-S1 or reconstructed “E-S1” structure was employed, which was superposed well with that of equilibrated S1 Lec domain alone. It was also indicated that the corresponding allosteric pathway from S1 to S1' conformation started with the separation between residues Q30 and K67 and terminated with the release of residue N87 from residue C109. These results provided an insight into understanding the structural transition and the structure-function relationship of P-selectin allostery

Unfolding of the Amyloid β-Peptide Central Helix: Mechanistic Insights from Molecular Dynamics Simulations

Alzheimer's disease (AD) pathogenesis is associated with formation of amyloid fibrils caused by polymerization of the amyloid β-peptide (Aβ), which is a process that requires unfolding of the native helical structure of Aβ. According to recent experimental studies, stabilization of the Aβ central helix is effective in preventing Aβ polymerization into toxic assemblies. To uncover the fundamental mechanism of unfolding of the Aβ central helix, we performed molecular dynamics simulations for wild-type (WT), V18A/F19A/F20A mutant (MA), and V18L/F19L/F20L mutant (ML) models of the Aβ central helix. It was quantitatively demonstrated that the stability of the α-helical conformation of both MA and ML is higher than that of WT, indicating that the α-helical propensity of the three nonpolar residues (18, 19, and 20) is the main factor for the stability of the whole Aβ central helix and that their hydrophobicity plays a secondary role. WT was found to completely unfold by a three-step mechanism: 1) loss of α-helical backbone hydrogen bonds, 2) strong interactions between nonpolar sidechains, and 3) strong interactions between polar sidechains. WT did not completely unfold in cases when any of the three steps was omitted. MA and ML did not completely unfold mainly due to the lack of the first step. This suggests that disturbances in any of the three steps would be effective in inhibiting the unfolding of the Aβ central helix. Our findings would pave the way for design of new drugs to prevent or retard AD

Homology modeling and molecular dynamics simulations of MUC1-9/H-2Kb complex suggest novel binding interactions

International audienceHuman MUC1 is over-expressed in human adenocarcinomas and has been used as a target for immunotherapy studies. The 9-mer MUC1-9 peptide has been identified as one of the peptides which binds to murine MHC class I H-2K. The structure of MUC1-9 in complex with H-2K has been modeled and simulated with classical molecular dynamics, based on the x-ray structure of the SEV9 peptide/H-2K complex. Two independent trajectories with the solvated complex (10 ns in length) were produced. Approximately 12 hydrogen bonds were identified during both trajectories to contribute to peptide/MHC complex, as well as 1-2 water mediated hydrogen bonds. Stability of the complex was also confirmed by buried surface area analysis, although the corresponding values were about 20% lower than those of the original x-ray structure. Interestingly, a bulged conformation of the peptide's central region, partially characterized as a -turn, was found exposed form the binding groove. In addition, P1 and P9 residues remained bound in the A and F binding pockets, even though there was a suggestion that P9 was more flexible. The complex lacked numerous water mediated hydrogen bonds that were present in the reference peptide x-ray structure. Moreover, local displacements of residues Asp4, Thr5 and Pro9 resulted in loss of some key interactions with the MHC molecule. This might explain the reduced affinity of the MUC1-9 peptide, relatively to SEV9, for the MHC class I H-2K

Mechanism of eIF6 release from the nascent 60S ribosomal subunit.

SBDS protein (deficient in the inherited leukemia-predisposition disorder Shwachman-Diamond syndrome) and the GTPase EFL1 (an EF-G homolog) activate nascent 60S ribosomal subunits for translation by catalyzing eviction of the antiassociation factor eIF6 from nascent 60S ribosomal subunits. However, the mechanism is completely unknown. Here, we present cryo-EM structures of human SBDS and SBDS-EFL1 bound to Dictyostelium discoideum 60S ribosomal subunits with and without endogenous eIF6. SBDS assesses the integrity of the peptidyl (P) site, bridging uL16 (mutated in T-cell acute lymphoblastic leukemia) with uL11 at the P-stalk base and the sarcin-ricin loop. Upon EFL1 binding, SBDS is repositioned around helix 69, thus facilitating a conformational switch in EFL1 that displaces eIF6 by competing for an overlapping binding site on the 60S ribosomal subunit. Our data reveal the conserved mechanism of eIF6 release, which is corrupted in both inherited and sporadic leukemias.Supported by a Federation of European Biochemical Societies Long term Fellowship (to FW), Specialist Programme from Bloodwise [12048] (AJW), the Medical Research Council [MC_U105161083] (AJW) and [U105115237] (RRK), Wellcome Trust strategic award to the Cambridge Institute for Medal Research [100140], Tesni Parry Trust (AJW), Ted’s Gang (AJW) and the Cambridge NIHR Biomedical Research Centre.This is the author accepted manuscript. The final version is available from Nature Publishing Group via http://dx.doi.org/10.1038/nsmb.311

Modeling a New Water Channel That Allows SET9 to Dimethylate p53

SET9, a protein lysine methyltransferase, has been thought to be capable of transferring only one methyl group to target lysine residues. However, some reports have pointed out that SET9 can dimethylate Lys372 of p53 (p53-K372) and Lys4 of histone H3 (H3-K4). In order to understand how p53 can be dimethylated by SET9, we measured the radius of the channel that surrounds p53-K372, first on the basis of the crystal structure of SET9, and we show that the channel is not suitable for water movement. Second, molecular dynamic (MD) simulations were carried out for 204 ns on the crystal structure of SET9. The results show that water leaves the active site of SET9 through a new channel, which is made of G292, A295, Y305 and Y335. In addition, the results of molecular docking and MD simulations indicate that the new water channel continues to remain open when S-adenosyl-L-methionine (AdoMet) or S-adenosyl-L-homocysteine (AdoHcy) is bound to SET9. The changes in the radii of these two channels were measured in the equilibrium phase at the constant temperature of 300 K. The results indicate that the first channel still does not allow water to get into or out of the active site, but the new channel is large enough to allow this water to circulate. Our results indicate that water can be removed from the active site, an essential process for allowing the dimethylation reaction to occur

Quantifying Intramolecular Binding in Multivalent Interactions: A Structure-Based Synergistic Study on Grb2-Sos1 Complex

Numerous signaling proteins use multivalent binding to increase the specificity and affinity of their interactions within the cell. Enhancement arises because the effective binding constant for multivalent binding is larger than the binding constants for each individual interaction. We seek to gain both qualitative and quantitative understanding of the multivalent interactions of an adaptor protein, growth factor receptor bound protein-2 (Grb2), containing two SH3 domains interacting with the nucleotide exchange factor son-of-sevenless 1 (Sos1) containing multiple polyproline motifs separated by flexible unstructured regions. Grb2 mediates the recruitment of Sos1 from the cytosol to the plasma membrane where it activates Ras by inducing the exchange of GDP for GTP. First, using a combination of evolutionary information and binding energy calculations, we predict an additional polyproline motif in Sos1 that binds to the SH3 domains of Grb2. This gives rise to a total of five polyproline motifs in Sos1 that are capable of binding to the two SH3 domains of Grb2. Then, using a hybrid method combining molecular dynamics simulations and polymer models, we estimate the enhancement in local concentration of a polyproline motif on Sos1 near an unbound SH3 domain of Grb2 when its other SH3 domain is bound to a different polyproline motif on Sos1. We show that the local concentration of the Sos1 motifs that a Grb2 SH3 domain experiences is approximately 1000 times greater than the cellular concentration of Sos1. Finally, we calculate the intramolecular equilibrium constants for the crosslinking of Grb2 on Sos1 and use thermodynamic modeling to calculate the stoichiometry. With these equilibrium constants, we are able to predict the distribution of complexes that form at physiological concentrations. We believe this is the first systematic analysis that combines sequence, structure, and thermodynamic analyses to determine the stoichiometry of the complexes that are dominant in the cellular environment