Study of Human RIG-I Polymorphisms Identifies Two Variants with an Opposite Impact on the Antiviral Immune Response

Background

RIG-I is a pivotal receptor that detects numerous RNA and DNA viruses. Thus, its defectiveness may strongly impair the host antiviral immunity. Remarkably, very little information is available on RIG-I single-nucleotide polymorphisms (SNPs) presenting a functional impact on the host response.

Methodology/Principal Findings

Here, we studied all non-synonymous SNPs of RIG-I using biochemical and structural modeling approaches. We identified two important variants: (i) a frameshift mutation (P₂₂₉fs) that generates a truncated, constitutively active receptor and (ii) a serine to isoleucine mutation (S₁₈₃I), which drastically inhibits antiviral signaling and exerts a down-regulatory effect, due to unintended stable complexes of RIG-I with itself and with MAVS, a key downstream adapter protein.

Conclusions/Significance

Hence, this study characterized P₂₂₉fs and S₁₈₃I SNPs as major functional RIG-I variants and potential genetic determinants of viral susceptibility. This work also demonstrated that serine 183 is a residue that critically regulates RIG-I-induced antiviral signaling.     

A novel approach to fold recognition using sequence-derived properties from sets of structurally similar local fragments of proteins

Comparative modeling methods can consistently produce reliable structural models for protein sequences with more than 25% sequence identity to proteins with known structure. However, there is a good chance that also sequences with lower sequence identity have their structural components represented in structural databases. To this end, we present a novel fragment-based method using sets of structurally similar local fragments of proteins. The approach differs from other fragment-based methods that use only single backbone fragments. Instead, we use a library of groups containing sets of sequence fragments with geometrically similar local structures and extract sequence related properties to assign these specific geometrical conformations to target sequences. We test the ability of the approach to recognize correct SCOP folds for 273 sequences from the 49 most popular folds. 49% of these sequences have the correct fold as their top prediction, while 82% have the correct fold in one of the top five predictions. Moreover, the approach shows no performance reduction on a subset of sequence targets with less than 10% sequence identity to any protein used to build the library.     

Peptide nucleic acid-metal complex conjugates: facile modulation of PNA-DNA duplex stability

Conjugates of peptide nucleic acids (PNA) and metal binding ligands were prepared using solid-phase synthesis. Stability of duplexes of bis-picolylamine-PNA conjugates and DNA was found to be modulated by equimolar concentrations of bioavailable metal ions: Ni(2+), Zn(2+)>Cu(2+). Sequence specificity of PNA was not compromised in the presence of these metal ions.     

Gain-of-function mutations reveal expanded intermediate states and a sequential action of two gates in MscL

The tension-driven gating transition in the large mechanosensitive channel MscL proceeds through detectable states of intermediate conductance. Gain-of-function (GOF) mutants with polar or charged substitutions in the main hydrophobic gate display altered patterns of subconducting states, providing valuable information about gating intermediates. Here we present thermodynamic analysis of several GOF mutants to clarify the nature and position of low-conducting conformations in the transition pathway. Unlike wild-type (WT) MscL, which predominantly occupies the closed and fully open states with very brief substates, the mild V23T GOF mutant frequently visits a multitude of short-lived subconducting states. Severe mutants V23D and G22N open in sequence: closed (C) --> low-conducting substate (S) --> open (O), with the first subtransition occurring at lower tensions. Analyses of equilibrium state occupancies as functions of membrane tension show that the C-->S subtransition in WT MscL is associated with only a minor conductance increment, but the largest in-plane expansion and free energy change. The GOF substitutions strongly affect the first subtransition by reducing area ((Delta)A) and energy ((Delta)E) changes between C and S states commensurably with the severity of mutation. GOF mutants also exhibited a considerably larger (Delta)E associated with the second (S-->O) subtransition, but a (Delta)A similar to WT. The area changes indicate that closed conformations of GOF mutants are physically preexpanded. The tension dependencies of rate constants for channel closure (k(off)) predict different positions of rate-limiting barriers on the energy-area profiles for WT and GOF MscL. The data support the two-gate mechanism in which the first subtransition (C-->S) can be viewed as opening of the central (M1) gate, resulting in an expanded water-filled "leaky" conformation. Strong facilitation of this step by polar GOF substitutions suggests that separation of M1 helices associated with hydration of the pore in WT MscL is the major energetic barrier for opening. Mutants with a stabilized S1 gate demonstrate impeded transitions from low-conducting substates to the fully open state, whereas extensions of S1-M1 linkers result in a much higher probability of reverse O-->S transitions. These data strongly suggest that the bulk of conductance gain in the second subtransition (S-->O) occurs through the opening of the NH2-terminal (S1) gate and the linkers are coupling elements between the M1 and S1 gates.     

Phosphodiesterase 4D deficiency in the ryanodine-receptor complex promotes heart failure and arrhythmias

Phosphodiesterases (PDEs) regulate the local concentration of 3',5' cyclic adenosine monophosphate (cAMP) within cells. cAMP activates the cAMP-dependent protein kinase (PKA). In patients, PDE inhibitors have been linked to heart failure and cardiac arrhythmias, although the mechanisms are not understood. We show that PDE4D gene inactivation in mice results in a progressive cardiomyopathy, accelerated heart failure after myocardial infarction, and cardiac arrhythmias. The phosphodiesterase 4D3 (PDE4D3) was found in the cardiac ryanodine receptor (RyR2)/calcium-release-channel complex (required for excitation-contraction [EC] coupling in heart muscle). PDE4D3 levels in the RyR2 complex were reduced in failing human hearts, contributing to PKA-hyperphosphorylated, "leaky" RyR2 channels that promote cardiac dysfunction and arrhythmias. Cardiac arrhythmias and dysfunction associated with PDE4 inhibition or deficiency were suppressed in mice harboring RyR2 that cannot be PKA phosphorylated. These data suggest that reduced PDE4D activity causes defective RyR2-channel function associated with heart failure and arrhythmias.     

Optimal Therapy Scheduling Based on a Pair of Collaterally Sensitive Drugs

Despite major strides in the treatment of cancer, the development of drug resistance remains a major hurdle. One strategy which has been proposed to address this is the sequential application of drug therapies where resistance to one drug induces sensitivity to another drug, a concept called collateral sensitivity. The optimal timing of drug switching in these situations, however, remains unknown. To study this, we developed a dynamical model of sequential therapy on heterogeneous tumors comprised of resistant and sensitive cells. A pair of drugs (DrugA, DrugB) are utilized and are periodically switched during therapy. Assuming resistant cells to one drug are collaterally sensitive to the opposing drug, we classified cancer cells into two groups, \(A_\mathrm{R}\) and \(B_\mathrm{R}\), each of which is a subpopulation of cells resistant to the indicated drug and concurrently sensitive to the other, and we subsequently explored the resulting population dynamics. Specifically, based on a system of ordinary differential equations for \(A_\mathrm{R}\) and \(B_\mathrm{R}\), we determined that the optimal treatment strategy consists of two stages: an initial stage in which a chosen effective drug is utilized until a specific time point, T, and a second stage in which drugs are switched repeatedly, during which each drug is used for a relative duration (i.e., \(f \Delta t\)-long for DrugA and \((1-f) \Delta t\)-long for DrugB with \(0 \le f \le 1\) and \(\Delta t \ge 0\)). We prove that the optimal duration of the initial stage, in which the first drug is administered, T, is shorter than the period in which it remains effective in decreasing the total population, contrary to current clinical intuition. We further analyzed the relationship between population makeup, \(\mathcal {A/B} = A_\mathrm{R}/B_\mathrm{R}\), and the effect of each drug. We determine a critical ratio, which we term \(\mathcal {(A/B)}^{*}\), at which the two drugs are equally effective. As the first stage of the optimal strategy is applied, \(\mathcal {A/B}\) changes monotonically to \(\mathcal {(A/B)}^{*}\) and then, during the second stage, remains at \(\mathcal {(A/B)}^{*}\) thereafter. Beyond our analytic results, we explored an individual-based stochastic model and presented the distribution of extinction times for the classes of solutions found. Taken together, our results suggest opportunities to improve therapy scheduling in clinical oncology.     

The study of complex-formation of DNA with the antimicrobial drug decamethoxine

The interaction of effective antibacterial drug decametoxyn with natural DNA was studied by UV-spectroscopy. Decametoxyn shows a specificity to nucleotides: it decreases the cooperativity of melting and the thermal stability of DNA parts enriched by AT pairs. The characteristics of the helix-coil transition on the DNA parts enriched by GC-pairs are invariable. Interaction with AT-pairs results in their partial or complete melting at room temperature, followed by intermolecule aggregation. Interacting with phosphates decametoxyn manifests itself not as a dication but as two single-charged ions.     

A strategy for cloning glycosyltransferase genes involved in natural product biosynthesis

The soil-borne and marine gram-positive Actinomycetes are a particularly rich source of carbohydrate-containing metabolites. With the advent of molecular tools and recombinant methods applicable to Actinomycetes, it has become feasible to investigate the biosynthesis of glycosylated compounds at genetic and biochemical levels, which has finally set the basis for engineering novel natural product derivatives. Glycosyltransferases (GT) are key enzymes for the biosynthesis of many valuable natural products that contain sugar moieties and they are most important for drug engineering. So far, the direct cloning of unknown glycosyltransferase genes by polymerase chain reaction (PCR) has not been described because glycosyltransferases do not share strongly conserved amino acid regions. In this study, we report a method for cloning of novel so far unidentified glycosyltransferase genes from different Actinomycetes strain. This was achieved by designing primers after a strategy named consensus-degenerate hybrid oligonucleotide primer (CODEHOP). Using this approach, 22 novel glycosyltransferase encoding genes putatively involved in the decoration of polyketides were cloned from the genomes of 10 Actinomycetes. In addition, a phylogenetic analysis of glycosyltransferases from Actinomycetes is shown in this paper.     

Locating missing water molecules in protein cavities by the three-dimensional reference interaction site model theory of molecular solvation

Water molecules confined in protein cavities are of great importance in understanding the protein structure and functions. However, it is a nontrivial task to locate such water molecules in protein by the ordinary molecular simulation and modeling techniques as well as experimental methods. The present study proves that the three-dimensional reference interaction site model (3D-RISM) theory, a recently developed statistical-mechanical theory of molecular solvation, has an outstanding advantage in locating such water molecules. In this paper, we demonstrate that the 3D-RISM theory is able to reproduce the structure and the number of water molecules in cavities of hen egg-white lysozyme observed commonly in the X-ray structures of different resolutions and conditions. Furthermore, we show that the theory successfully identified a water molecule in a cavity, the existence of which has been ambiguous even from the X-ray results. In contrast, we confirmed that molecular dynamics simulation is helpless at present to find such water molecules because the results substantially depend on the initial coordinates of water molecules. Possible applications of the theory to problems in the fields of biochemistry and biophysics are also discussed.