Computation of Hemagglutinin Free Energy Difference by the Confinement Method

Hemagglutinin (HA) mediates membrane fusion, a crucial step during influenza virus cell entry. How many HAs are needed for this process is still subject to debate. To aid in this discussion, the confinement free energy method was used to calculate the conformational free energy difference between the extended intermediate and postfusion state of HA. Special care was taken to comply with the general guidelines for free energy calculations, thereby obtaining convergence and demonstrating reliability of the results. The energy that one HA trimer contributes to fusion was found to be 34.2 ± 3.4<i>k</i>_B<i>T</i>, similar to the known contributions from other fusion proteins. Although computationally expensive, the technique used is a promising tool for the further energetic characterization of fusion protein mechanisms. Knowledge of the energetic contributions per protein, and of conserved residues that are crucial for fusion, aids in the development of fusion inhibitors for antiviral drugs

CD44 modulates Smad1 activation in the BMP-7 signaling pathway

Bone morphogenetic protein 7 (BMP-7) regulates cellular metabolism in embryonic and adult tissues. Signal transduction occurs through the activation of intracellular Smad proteins. In this paper, using a yeast two-hybrid screen, Smad1 was found to interact with the cytoplasmic domain of CD44, a receptor for the extracellular matrix macromolecule hyaluronan. Coimmunoprecipitation experiments confirmed the interaction of Smad1 with full-length CD44—interactions that did not occur when CD44 receptors truncated within the cytoplasmic domain were tested. Chondrocytes overexpressing a truncated CD44 on a background of endogenous full-length CD44 no longer exhibited Smad1 nuclear translocation upon BMP-7 stimulation. Further, pretreatment of chondrocytes with Streptomyces hyaluronidase to disrupt extracellular hyaluronan–cell interactions inhibited BMP-7–mediated Smad1 phosphorylation, nuclear translocation of Smad1 or Smad4, and SBE4–luciferase reporter activation. These results support a functional link between the BMP signaling cascade and CD44. Thus, changes in hyaluronan–cell interactions may serve as a means to modulate cellular responsiveness to BMP

Reaction profiling of a set of acrylamide-based human tissue transglutaminase inhibitors

The major function of the enzyme human tissue transglutaminase (TG2) is the crosslinking of proteins via a transamidation between the γ-carboxamide of a glutamine and the ε-amino group of a lysine. Overexpression of TG2 can lead to undesirable outcomes and has been linked to conditions such as fibrosis, celiac disease and neurodegenerative diseases. Accordingly, TG2 is a tempting drug target. The most effective TG2 inhibitors to date are small-molecule peptidomimetics featuring electrophilic warheads that irreversibly modify the active site catalytic cysteine (CYS277). In an effort to facilitate the design of such TG2 inhibitors, we undertook a quantum mechanical reaction profiling of the Michael reaction between a set of six acrylamide-based known TG2 inhibitors and the TG2 CYS277. The inhibitors were docked into the active site and the coordinates were refined by MD simulations prior to modelling the covalent modification of the CYS277 thiolate. The results of QM/MM MD umbrella sampling applied to reaction coordinates driving the Michael reaction are presented for two approximations of the Michael reaction: a concerted reaction (simultaneous thiolate attack onto the acrylamide warhead and pronation from the adjacent HIS335) and a two-stage reaction (consecutive thiolate attack and protonation). The two-stage approximation of the Michael reaction gave the better results for the evaluation of acrylamide-based potential TG2 inhibitors in silico. Good correlations were observed between the experimental TG2 IC50 data and the calculated activation energies over the range 0.0061 – 6.3 µM (three orders of magnitude) and we propose that this approach may be used to evaluate acrylamide-based potential TG2 inhibitors

Multiple Binding Sites for Fatty Acids on the Potassium Channel KcsA

Interactions of fatty acids with the potassium channel KcsA were studied using Trp fluorescence quenching and electron paramagnetic resonance (EPR) techniques. The brominated analogue of oleic acid was shown to bind to annular sites on KcsA and to the nonannular sites at each protein-protein interface in the homotetrameric structure with binding constants relative to dioleoylphosphatidylcholine of 0.67 ± 0.04 and 0.87 ± 0.08, respectively. Mutation of the two Arg residues close to the nonannular binding sites had no effect on fatty acid binding. EPR studies with a spin-labeled analogue of stearic acid detected a high-affinity binding site for the fatty acid with strong immobilization. Fluorescence quenching studies with the spin-labeled analogue showed that the binding site detected in the EPR experiments could not be one of the annular or nonannular binding sites. Instead, it is proposed that the EPR studies detect binding to the central hydrophobic cavity of the channel, with a binding constant in the range of ~0.1-1 ?M

PLUMED: a portable plugin for free-energy calculations with molecular dynamics

Here we present a program aimed at free-energy calculations in molecular systems. It consists of a series of routines that can be interfaced with the most popular classical molecular dynamics (MD) codes through a simple patching procedure. This leaves the possibility for the user to exploit many different MD engines depending on the system simulated and on the computational resources available. Free-energy calculations can be performed as a function of many collective variables, with a particular focus on biological problems, and using state-of-the-art methods such as metadynamics, umbrella sampling and Jarzynski-equation based steered MD. The present software, written in ANSI-C language, can be easily interfaced with both fortran and C/C++ codes.Comment: to be submitted to Computer Physics Communication

Interaction of short modified peptides deriving from glycoprotein gp36 of feline immunodeficiency virus with phospholipid membranes

A tryptophan-rich octapeptide, C8 (Ac-Trp-Glu-Asp-Trp-Val-Gly-Trp-Ile-NH2), modelled on the membrane-proximal external region of the feline immunodeficiency virus (FIV) gp36 glycoprotein ectodomain, exhibits potent antiviral activity against FIV. A mechanism has been proposed by which the peptide, being positioned on the surface of the cell membrane, inhibits its fusion with the virus. In the present work, peptide–lipid interactions of C8 with dimyristoyl phosphatidylcholine liposomes are investigated using electron spin resonance spectroscopy of spin-labelled lipids. Three other peptides, obtained from modifications of C8, have also been investigated, in an attempt to clarify the essential molecular features of the interactions involving the tryptophan residues. The results show that C8 adsorbs strongly on the bilayer surface. Membrane binding requires not only the presence of the Trp residues in the sequence, but also their common orientation on one side of the peptide that is engendered by the WX2 WX2 W motif. Membrane interaction correlates closely with peptide antiviral activity, indicating that the membrane is essential in stabilizing the peptide conformation that will be able to inhibit viral infection

The Roles of Transmembrane Domain Helix-III during Rhodopsin Photoactivation

Background: Rhodopsin, the prototypic member of G protein-coupled receptors (GPCRs), undergoes isomerization of 11- cis-retinal to all-trans-retinal upon photoactivation. Although the basic mechanism by which rhodopsin is activated is well understood, the roles of whole transmembrane (TM) helix-III during rhodopsin photoactivation in detail are not completely clear. Principal Findings: We herein use single-cysteine mutagenesis technique to investigate conformational changes in TM helices of rhodopsin upon photoactivation. Specifically, we study changes in accessibility and reactivity of cysteine residues introduced into the TM helix-III of rhodopsin. Twenty-eight single-cysteine mutants of rhodopsin (P107C-R135C) were prepared after substitution of all natural cysteine residues (C140/C167/C185/C222/C264/C316) by alanine. The cysteine mutants were expressed in COS-1 cells and rhodopsin was purified after regeneration with 11-cis-retinal. Cysteine accessibility in these mutants was monitored by reaction with 4, 49-dithiodipyridine (4-PDS) in the dark and after illumination. Most of the mutants except for T108C, G109C, E113C, I133C, and R135C showed no reaction in the dark. Wide variation in reactivity was observed among cysteines at different positions in the sequence 108–135 after photoactivation. In particular, cysteines at position 115, 119, 121, 129, 131, 132, and 135, facing 11-cis-retinal, reacted with 4-PDS faster than neighboring amino acids. The different reaction rates of mutants with 4-PDS after photoactivation suggest that the amino acids in different positions in helix-III are exposed to aqueous environment to varying degrees. Significance: Accessibility data indicate that an aqueous/hydrophobic boundary in helix-III is near G109 and I133. The lack of reactivity in the dark and the accessibility of cysteine after photoactivation indicate an increase of water/4-PDS accessibility for certain cysteine-mutants at Helix-III during formation of Meta II. We conclude that photoactivation resulted in water-accessible at the chromophore-facing residues of Helix-III.National Institutes of Health (U.S.) (grant GM28289)National Eye Institute (Grant Grant EY11716)National Science Foundation (U.S.) (grant EIA-0225609

Molecular fingerprinting of radiation resistant tumors: Can we apprehend and rehabilitate the suspects?

Radiation therapy continues to be one of the more popular treatment options for localized prostate cancer. One major obstacle to radiation therapy is that there is a limit to the amount of radiation that can be safely delivered to the target organ. Emerging evidence suggests that therapeutic agents targeting specific molecules might be combined with radiation therapy for more effective treatment of tumors. Recent studies suggest that modulation of these molecules by a variety of mechanisms (e.g., gene therapy, antisense oligonucleotides, small interfering RNA) may enhance the efficacy of radiation therapy by modifying the activity of key cell proliferation and survival pathways such as those controlled by Bcl-2, p53, Akt/PTEN and cyclooxygenase-2. In this article, we summarize the findings of recent investigations of radiosensitizing agents in the treatment of prostate cancer

Effects of Restrained Sampling Space and Nonplanar Amino Groups on Free-Energy Predictions for RNA with Imino and Sheared Tandem GA Base Pairs Flanked by GC, CG, iGiC or iCiG Base Pairs

Guanine-adenine (GA) base pairs play important roles in determining the structure, dynamics, and stability of RNA. In RNA internal loops, GA base pairs often occur in tandem arrangements and their structure is context and sequence dependent. Calculations reported here test the thermodynamic integration (TI) approach with the amber99 force field by comparing computational predictions of free energy differences with the free energy differences expected on the basis of NMR determined structures of the RNA motifs (5′-GCGGACGC-3′)2, (5′-GCiGGAiCGC-3′)2, (5′-GGCGAGCC-3′)2, and (5′-GGiCGAiGCC-3′)2. Here, iG and iC denote isoguanosine and isocytidine, which have amino and carbonyl groups transposed relative to guanosine and cytidine. The NMR structures show that the GA base pairs adopt either imino (cis Watson−Crick/Watson−Crick A-G) or sheared (trans Hoogsteen/Sugar edge A-G) conformations depending on the identity and orientation of the adjacent base pair. A new mixing function for the TI method is developed that allows alchemical transitions in which atoms can disappear in both the initial and final states. Unrestrained calculations gave ΔG° values 2−4 kcal/mol different from expectations based on NMR data. Restraining the structures with hydrogen bond restraints did not improve the predictions. Agreement with NMR data was improved by 0.7 to 1.5 kcal/mol, however, when structures were restrained with weak positional restraints to sample around the experimentally determined NMR structures. The amber99 force field was modified to partially include pyramidalization effects of the unpaired amino group of guanosine in imino GA base pairs. This provided little or no improvement in comparisons with experiment. The marginal improvement is observed when the structure has potential cross-strand out-of-plane hydrogen bonding with the G amino group. The calculations using positional restraints and a nonplanar amino group reproduce the signs of ΔG° from the experimental results and are, thus, capable of providing useful qualitative insights complementing the NMR experiments. Decomposition of the terms in the calculations reveals that the dominant terms are from electrostatic and interstrand interactions other than hydrogen bonds in the base pairs. The results suggest that a better description of the backbone is key to reproducing the experimental free energy results with computational free energy predictions