A Novel Dimeric Inhibitor Targeting Beta2GPI in Beta2GPI/Antibody Complexes Implicated in Antiphospholipid Syndrome

Background: b2GPI is a major antigen for autoantibodies associated with antiphospholipid syndrome (APS), an autoimmune disease characterized by thrombosis and recurrent pregnancy loss. Only the dimeric form of b2GPI generated by anti-b2GPI antibodies is pathologically important, in contrast to monomeric b2GPI which is abundant in plasma. Principal Findings: We created a dimeric inhibitor, A1-A1, to selectively target b2GPI in b2GPI/antibody complexes. To make this inhibitor, we isolated the first ligand-binding module from ApoER2 (A1) and connected two A1 modules with a flexible linker. A1-A1 interferes with two pathologically important interactions in APS, the binding of b2GPI/antibody complexes with anionic phospholipids and ApoER2. We compared the efficiency of A1-A1 to monomeric A1 for inhibition of the binding of b2GPI/antibody complexes to anionic phospholipids. We tested the inhibition of b2GPI present in human serum, b2GPI purified from human plasma and the individual domain V of b2GPI. We demonstrated that when b2GPI/antibody complexes are formed, A1-A1 is much more effective than A1 in inhibition of the binding of b2GPI to cardiolipin, regardless of the source of b2GPI. Similarly, A1-A1 strongly inhibits the binding of dimerized domain V of b2GPI to cardiolipin compared to the monomeric A1 inhibitor. In the absence of anti-b2GPI antibodies, both A1-A1 and A1 only weakly inhibit the binding of pathologically inactive monomeric b2GPI to cardiolipin. Conclusions: Our results suggest that the approach of using a dimeric inhibitor to block b2GPI in the pathologica

Free Energy Simulations of a GTPase: GTP and GDP Binding to Archaeal Initiation Factor 2

International audienceArchaeal initiation factor 2 (aIF2) is a protein involved in the initiation of protein biosynthesis. In its GTP-bound, "ON" conformation, aIF2 binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To understand the specific binding of GTP and GDP and its dependence on the ON or OFF conformational state of aIF2, molecular dynamics free energy simulations (MDFE) are a tool of choice. However, the validity of the computed free energies depends on the simulation model, including the force field and the boundary conditions, and on the extent of conformational sampling in the simulations. aIF2 and other GTPases present specific difficulties; in particular, the nucleotide ligand coordinates a divalent Mg(2+) ion, which can polarize the electronic distribution of its environment. Thus, a force field with an explicit treatment of electronic polarizability could be necessary, rather than a simpler, fixed charge force field. Here, we begin by comparing a fixed charge force field to quantum chemical calculations and experiment for Mg(2+):phosphate binding in solution, with the force field giving large errors. Next, we consider GTP and GDP bound to aIF2 and we compare two fixed charge force fields to the recent, polarizable, AMOEBA force field, extended here in a simple, approximate manner to include GTP. We focus on a quantity that approximates the free energy to change GTP into GDP. Despite the errors seen for Mg(2+):phosphate binding in solution, we observe a substantial cancellation of errors when we compare the free energy change in the protein to that in solution, or when we compare the protein ON and OFF states. Finally, we have used the fixed charge force field to perform MDFE simulations and alchemically transform GTP into GDP in the protein and in solution. With a total of about 200 ns of molecular dynamics, we obtain good convergence and a reasonable statistical uncertainty, comparable to the force field uncertainty, and somewhat lower than the predicted GTP/GDP binding free energy differences. The sign and magnitudes of the differences can thus be interpreted at a semiquantitative level, and are found to be consistent with the experimental binding preferences of ON- and OFF-aIF2

Fast calculation of thermodynamic and structural parameters of solutions using the 3DRISM model and the multi-grid method

In the paper a new method to solve the tree-dimensional reference interaction site model (3DRISM) integral equations is proposed. The algorithm uses the multi-grid technique which allows to decrease the computational expanses. 3DRISM calculations for aqueous solutions of four compounds (argon, water, methane, methanol) on the different grids are performed in order to determine a dependence of the computational error on the parameters of the grid. It is shown that calculations on the grid with the step 0.05\Angstr and buffer 8\Angstr give the error of solvation free energy calculations less than 0.3 kcal/mol which is comparable to the accuracy of the experimental measurements. The performance of the algorithm is tested. It is shown that the proposed algorithm is in average more than 12 times faster than the standard Picard direct iteration method.Comment: the information in this preprint is not up to date. Since the first publication of the preprint (9 Nov 2011) the algorithm was modified which allowed to achieve better results. For the new algorithm see the JCTC paper: DOI: 10.1021/ct200815v, http://pubs.acs.org/doi/abs/10.1021/ct200815

Protonation States of Remote Residues Affect Binding-Release Dynamics of the Ligand but not the Conformation of apo Ferric Binding Protein

We have studied the apo (Fe3+ free) form of periplasmic ferric binding protein (FbpA) under different conditions and we have monitored the changes in the binding and release dynamics of H2PO4- that acts as a synergistic anion in the presence of Fe3+. Our simulations predict a dissociation constant of 2.2$\pm$0.2 mM which is in remarkable agreement with the experimentally measured value of 2.3$\pm$0.3 mM under the same ionization strength and pH conditions. We apply perturbations relevant for changes in environmental conditions as (i) different values of ionic strength (IS), and (ii) protonation of a group of residues to mimic a different pH environment. Local perturbations are also studied by protonation or mutation of a site distal to the binding region that is known to mechanically manipulate the hinge-like motions of FbpA. We find that while the average conformation of the protein is intact in all simulations, the H2PO4- dynamics may be substantially altered by the changing conditions. In particular, the bound fraction which is 20$\%$ for the wild type system is increased to 50$\%$ with a D52A mutation/protonation and further to over 90$\%$ at the protonation conditions mimicking those at pH 5.5. The change in the dynamics is traced to the altered electrostatic distribution on the surface of the protein which in turn affects hydrogen bonding patterns at the active site. The observations are quantified by rigorous free energy calculations. Our results lend clues as to how the environment versus single residue perturbations may be utilized for regulation of binding modes in hFbpA systems in the absence of conformational changes.Comment: 26 pages, 4 figure

All-d-Enantiomer of β-Amyloid Peptide Forms Ion Channels in Lipid Bilayers

Alzheimer’s disease (AD) is the most common type of senile dementia in aging populations. Amyloid β (Aβ)-mediated dysregulation of ionic homeostasis is the prevailing underlying mechanism leading to synaptic degeneration and neuronal death. Aβ-dependent ionic dysregulation most likely occurs either directly via unregulated ionic transport through the membrane or indirectly via Aβ binding to cell membrane receptors and subsequent opening of existing ion channels or transporters. Receptor binding is expected to involve a high degree of stereospecificity. Here, we investigated whether an Aβ peptide enantiomer, whose entire sequence consists of d-amino acids, can form ion-conducting channels; these channels can directly mediate Aβ effects even in the absence of receptor–peptide interactions. Using complementary approaches of planar lipid bilayer (PLB) electrophysiological recordings and molecular dynamics (MD) simulations, we show that the d-Aβ isomer exhibits ion conductance behavior in the bilayer indistinguishable from that described earlier for the l-Aβ isomer. The d isomer forms channel-like pores with heterogeneous ionic conductance similar to the l-Aβ isomer channels, and the d-isomer channel conductance is blocked by Zn2+, a known blocker of l-Aβ isomer channels. MD simulations further verify formation of β-barrel-like Aβ channels with d- and l-isomers, illustrating that both d- and l-Aβ barrels can conduct cations. The calculated values of the single-channel conductance are approximately in the range of the experimental values. These findings are in agreement with amyloids forming Ca2+ leaking, unregulated channels in AD, and suggest that Aβ toxicity is mediated through a receptor-independent, nonstereoselective mechanism

Ionic liquids at electrified interfaces

Until recently, “room-temperature” (<100–150 °C) liquid-state electrochemistry was mostly electrochemistry of diluted electrolytes(1)–(4) where dissolved salt ions were surrounded by a considerable amount of solvent molecules. Highly concentrated liquid electrolytes were mostly considered in the narrow (albeit important) niche of high-temperature electrochemistry of molten inorganic salts(5-9) and in the even narrower niche of “first-generation” room temperature ionic liquids, RTILs (such as chloro-aluminates and alkylammonium nitrates).(10-14) The situation has changed dramatically in the 2000s after the discovery of new moisture- and temperature-stable RTILs.(15, 16) These days, the “later generation” RTILs attracted wide attention within the electrochemical community.(17-31) Indeed, RTILs, as a class of compounds, possess a unique combination of properties (high charge density, electrochemical stability, low/negligible volatility, tunable polarity, etc.) that make them very attractive substances from fundamental and application points of view.(32-38) Most importantly, they can mix with each other in “cocktails” of one’s choice to acquire the desired properties (e.g., wider temperature range of the liquid phase(39, 40)) and can serve as almost “universal” solvents.(37, 41, 42) It is worth noting here one of the advantages of RTILs as compared to their high-temperature molten salt (HTMS)(43) “sister-systems”.(44) In RTILs the dissolved molecules are not imbedded in a harsh high temperature environment which could be destructive for many classes of fragile (organic) molecules

Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations

Alkali (Li+, Na+, K+, Rb+, and Cs+) and halide (F−, Cl−, Br−, and I−) ions play an important role in many biological phenomena, roles that range from stabilization of biomolecular structure, to influence on biomolecular dynamics, to key physiological influence on homeostasis and signaling. To properly model ionic interaction and stability in atomistic simulations of biomolecular structure, dynamics, folding, catalysis, and function, an accurate model or representation of the monovalent ions is critically necessary. A good model needs to simultaneously reproduce many properties of ions, including their structure, dynamics, solvation, and moreover both the interactions of these ions with each other in the crystal and in solution and the interactions of ions with other molecules. At present, the best force fields for biomolecules employ a simple additive, nonpolarizable, and pairwise potential for atomic interaction. In this work, we describe our efforts to build better models of the monovalent ions within the pairwise Coulombic and 6-12 Lennard-Jones framework, where the models are tuned to balance crystal and solution properties in Ewald simulations with specific choices of well-known water models. Although it has been clearly demonstrated that truly accurate treatments of ions will require inclusion of nonadditivity and polarizability (particularly with the anions) and ultimately even a quantum mechanical treatment, our goal was to simply push the limits of the additive treatments to see if a balanced model could be created. The applied methodology is general and can be extended to other ions and to polarizable force-field models. Our starting point centered on observations from long simulations of biomolecules in salt solution with the AMBER force fields where salt crystals formed well below their solubility limit. The likely cause of the artifact in the AMBER parameters relates to the naive mixing of the Smith and Dang chloride parameters with AMBER-adapted Åqvist cation parameters. To provide a more appropriate balance, we reoptimized the parameters of the Lennard-Jones potential for the ions and specific choices of water models. To validate and optimize the parameters, we calculated hydration free energies of the solvated ions and also lattice energies (LE) and lattice constants (LC) of alkali halide salt crystals. This is the first effort that systematically scans across the Lennard-Jones space (well depth and radius) while balancing ion properties like LE and LC across all pair combinations of the alkali ions and halide ions. The optimization across the entire monovalent series avoids systematic deviations. The ion parameters developed, optimized, and characterized were targeted for use with some of the most commonly used rigid and nonpolarizable water models, specifically TIP3P, TIP4PEW, and SPC/E. In addition to well reproducing the solution and crystal properties, the new ion parameters well reproduce binding energies of the ions to water and the radii of the first hydration shells

Identification and evaluation of novel protective antigens for the development of a candidate TB subunit vaccine

The development of a vaccine against tuberculosis (TB), a disease caused by Mycobacterium tuberculosis (M.tb), is urgently needed. The only currently available vaccine, BCG, has variable efficacy. One approach in the global vaccine development effort is focused on boosting BCG using subunit vaccines. The identification of novel antigens for inclusion in subunit vaccines is a critical step in the TB vaccine development pathway. We selected four novel mycobacterial antigens recognised during the course of human infection. A replication deficient chimpanzee adenovirus (ChAdOx1) was constructed to express each antigen individually and these vectors were evaluated for protective efficacy in murine M.tb challenge experiments. One antigen, PPE15 (Rv1039c), conferred significant and reproducible protection when administered alone and as a boost to BCG vaccination. We identified immunodominant epitopes to define the protective immune responses using tetramers and intravascular staining. Lung parenchymal CD4+ and CD8+ CXCR3+ KLRG1- T cells, previously associated with M.tb protection, were enriched in vaccinated compared to control groups. Further work to evaluate the protective efficacy of PPE15 in more stringent preclinical animal models is now merited, together with the identification of further novel protective antigens using this selection strategy

Identification and evaluation of novel protective antigens for the development of a candidate TB subunit vaccine

The development of a vaccine against tuberculosis (TB), a disease caused by Mycobacterium tuberculosis (M.tb), is urgently needed. The only currently available vaccine, BCG, has variable efficacy. One approach in the global vaccine development effort is focused on boosting BCG using subunit vaccines. The identification of novel antigens for inclusion in subunit vaccines is a critical step in the TB vaccine development pathway. We selected four novel mycobacterial antigens recognised during the course of human infection. A replication deficient chimpanzee adenovirus (ChAdOx1) was constructed to express each antigen individually and these vectors were evaluated for protective efficacy in murine M.tb challenge experiments. One antigen, PPE15 (Rv1039c), conferred significant and reproducible protection when administered alone and as a boost to BCG vaccination. We identified immunodominant epitopes to define the protective immune responses using tetramers and intravascular staining. Lung parenchymal CD4+ and CD8+ CXCR3+ KLRG1- T cells, previously associated with M.tb protection, were enriched in vaccinated compared to control groups. Further work to evaluate the protective efficacy of PPE15 in more stringent preclinical animal models is now merited, together with the identification of further novel protective antigens using this selection strategy