Reverse genetic studies of Benyvirus - Polymyxa betae molecular interaction: Role of the RNA4-encoded protein in virus transmission

Beet necrotic yellow vein virus (BNYVV), the leading infectious agent that affects sugar beet, is included within viruses transmitted through the soil from plasmodiophorid as Polymyxa betae. BNYVV is the causal agent of Rhizomania, which induces abnormal rootlet proliferation and is widespread in the sugar beet growing areas in Europe, Asia and America; for review see (Peltier et al., 2008). In this latter continent, Beet soil-borne mosaic virus (BSBMV) has been identified (Lee et al., 2001) and belongs to the benyvirus genus together with BNYVV, both vectored by P. betae. BSBMV is widely distributed only in the United States and it has not been reported yet in others countries. It was first identified in Texas as a sugar beet virus morphologically similar but serologically distinct to BNYVV. Subsequent sequence analysis of BSBMV RNAs evidenced similar genomic organization to that of BNYVV but sufficient molecular differences to distinct BSBMV and BNYVV in two different species (Rush et al., 2003). Benyviruses field isolates usually consist of four RNA species but some BNYVV isolates contain a fifth RNA. RNAs -1 contains a single long ORF encoding polypeptide that shares amino acid homology with known viral RNA-dependent RNA polymerases (RdRp) and helicases. RNAs -2 contains six ORFs: capsid protein (CP), one readthrough protein, triple gene block proteins (TGB) that are required for cell-to-cell virus movement and the sixth 14 kDa ORF is a post-translation gene silencing suppressor. RNAs -3 is involved on disease symptoms and is essential for virus systemic movement. BSBMV RNA-3 can be trans-replicated, trans-encapsidated by the BNYVV helper strain (RNA-1 and -2) (Ratti et al., 2009). BNYVV RNA-4 encoded one 31 kDa protein and is essential for vector interactions and virus transmission by P. betae (Rahim et al., 2007). BNYVV RNA-5 encoded 26 kDa protein that improve virus infections and accumulation in the hosts. We are interest on BSBMV effect on Rhizomania studies using powerful tools as full-length infectious cDNA clones. B-type full-length infectious cDNA clones are available (Quillet et al., 1989) as well as A/P-type RNA-3, -4 and -5 from BNYVV (unpublished). A-type BNYVV full-length clones are also available, but RNA-1 cDNA clone still need to be modified. During the PhD program, we start production of BSBMV full-length cDNA clones and we investigate molecular interactions between plant and Benyviruses exploiting biological, epidemiological and molecular similarities/divergences between BSBMV and BNYVV. During my PhD researchrs we obtained full length infectious cDNA clones of BSBMV RNA-1 and -2 and we demonstrate that they transcripts are replicated and packaged in planta and able to substitute BNYVV RNA-1 or RNA-2 in a chimeric viral progeny (BSBMV RNA-1 + BNYVV RNA-2 or BNYVV RNA-1 + BSBMV RNA-2). During BSBMV full-length cDNA clones production, unexpected 1,730 nts long form of BSBMV RNA-4 has been detected from sugar beet roots grown on BSBMV infected soil. Sequence analysis of the new BSBMV RNA-4 form revealed high identity (~100%) with published version of BSBMV RNA-4 sequence (NC_003508) between nucleotides 1-608 and 1,138-1,730, however the new form shows 528 additionally nucleotides between positions 608-1,138 (FJ424610). Two putative ORFs has been identified, the first one (nucleotides 383 to 1,234), encode a protein with predicted mass of 32 kDa (p32) and the second one (nucleotides 885 to 1,244) express an expected product of 13 kDa (p13). As for BSBMV RNA-3 (Ratti et al., 2009), full-length BSBMV RNA-4 cDNA clone permitted to obtain infectious transcripts that BNYVV viral machinery (Stras12) is able to replicate and to encapsidate in planta. Moreover, we demonstrated that BSBMV RNA-4 can substitute BNYVV RNA-4 for an efficient transmission through the vector P. betae in Beta vulgaris plants, demonstrating a very high correlation between BNYVV and BSBMV. At the same time, using BNYVV helper strain, we studied BSBMV RNA-4’s protein expression in planta. We associated a local necrotic lesions phenotype to the p32 protein expression onto mechanically inoculated C. quinoa. Flag or GFP-tagged sequences of p32 and p13 have been expressed in viral context, using Rep3 replicons, based on BNYVV RNA-3. Western blot analyses of local lesions contents, using FLAG-specific antibody, revealed a high molecular weight protein, which suggest either a strong interaction of BSBMV RNA4’s protein with host protein(s) or post translational modifications. GFP-fusion sequences permitted the subcellular localization of BSBMV RNA4’s proteins. Moreover we demonstrated the absence of self-activation domains on p32 by yeast two hybrid system approaches. We also confirmed that p32 protein is essential for virus transmission by P. betae using BNYVV helper strain and BNYVV RNA-3 and we investigated its role by the use of different deleted forms of p32 protein. Serial mechanical inoculation of wild-type BSBMV on C. quinoa plants were performed every 7 days. Deleted form of BSBMV RNA-4 (1298 bp) appeared after 14 passages and its sequence analysis shows deletion of 433 nucleotides between positions 611 and 1044 of RNA-4 new form. We demonstrated that this deleted form can’t support transmission by P. betae using BNYVV helper strain and BNYVV RNA-3, moreover we confirmed our hypothesis that BSBMV RNA-4 described by Lee et al. (2001) is a deleted form. Interesting after 21 passages we identifed one chimeric form of BSBMV RNA-4 and BSBMV RNA-3 (1146 bp). Two putative ORFs has been identified on its sequence, the first one (nucleotides 383 to 562), encode a protein with predicted mass of 7 kDa (p7), corresponding to the N-terminal of p32 protein encoded by BSBMV RNA-4; the second one (nucleotides 562 to 789) express an expected product of 9 kDa (p9) corresponding to the C-terminal of p29 encoded by BSBMV RNA-3. Results obtained by our research in this topic opened new research lines that our laboratories will develop in a closely future. In particular BSBMV p32 and its mutated forms will be used to identify factors, as host or vector protein(s), involved in the virus transmission through P. betae. The new results could allow selection or production of sugar beet plants able to prevent virus transmission then able to reduce viral inoculum in the soil.

Molecular mechanisms involved in the pathogenesis of beet soil-borne viruses.

The genus Benyvirus includes the most important and widespread sugar beet viruses transmitted through the soil by the plasmodiophorid Polymyxa betae. In particular Beet necrotic yellow vein virus (BNYVV), the leading infectious agent that affects sugar beet, causes an abnormal rootlet proliferation known as rhizomania. Beet soil-borne mosaic virus (BSBMV) is widely distributed in the United States and, up to date has not been reported in others countries. My PhD project aims to investigate molecular interactions between BNYVV and BSBMV and the mechanisms involved in the pathogenesis of these viruses. BNYVV full-length infectious cDNA clones were available as well as full-length cDNA clones of BSBMV RNA-1, -2, -3 and -4. Handling of these cDNA clones in order to produce in vitro infectious transcripts need sensitive and expensive steps, so I developed agroclones of BNYVV and BSBMV RNAs, as well as viral replicons allowing the expression of different proteins. Chenopodium quinoa and Nicotiana benthamiana plants have been infected with in vitro transcripts and agroclones to investigate the interaction between BNYVV and BSBMV RNA-1 and -2 and the behavior of artificial viral chimeras. Simultaneously I characterized BSBMV p14 and demonstrated that it is a suppressor of post-transcriptional gene silencing sharing common features with BNYVV p14.

Molecular interaction between perthiolated [beta]-cyclodextrin (CD) and the guests molecules adamantaneacetic acid (AD) and ferroceneacetic acid (FC); and the effect of the interaction on the electron transition of CD anchored particles

The molecular interactions between the host molecule, perthiolated beta-cyclodextrin (CD), and the guest molecules, adamantaneacetic acid (AD) and ferroceneacetic acid (FC), have been inestigated theoretically in both the gas and aqueous phases. The major computations have been carried out at the theoretical levels, RHF/6-31G and B3LYP/6- 31G. MP2 electronic energies were also computed based at the geometries optimized by both the RHF and B3LYP methods in the gas phase to establish a better estimate of the correlation effect. The solvent phase computations were completed at the RHF/6-31G and B3LYP/6-31G levels using the PCM model. The most stable structures optimized in gas phase by both the RHF and B3LYP methods were used for the computations in solution. A method to systematically manipulate the relative position and orientation between the interacting molecules is proposed. In the gas phase, six trials with different host-guest relative positions and orientations were completed successfully with the B3LYP method for both the CD-AD and CD-FC complexes. Only four trials were completed with RHF method. In the gas phase, the best results from the RHF method gives for the association Gibbs free energy (ΔG°) values equal to -32.21kj/mol for CD-AD and -25.73kj/mol for CD-FC. And the best results from the B3LYP method have ΔG° equal to -47.57kj/mol for CD-AD and -41.09kj/mol for CD-FC. The MP2 correction significantly lowers ΔG° based on the geometries from both methods. For the RHF structure, the MP2 computations lowered ΔG° to -60.64kj/mol for CD-AD and -54.10 for CD-FC. For the structure from the B3LYP method, it was reduced to -59.87 kj/mol for CD-AD and -54.84 kj/mol for CDFC. The RHF solvent phase calculations yielded following results: ΔG°(aq) equals 107.2kj/mol for CD-AD and 111.4kj/mol for CD-FC. Compared with the results from the RHF method, the B3LYP method provided clearly better solvent phase results with ΔG° (aq) equal to 38.64kj/mol for CD-AD and 39.61kj/mol for CD-FC. These results qualitatively explain the experimental observations. However quantitatively they are in poor agreement with the experimental values available in the literature and those recently published by Liu et al. And the reason is believed to be omission of hydrophobic contribution to the association. Determining the global geometrical minima for these very large systems was very difficult and computationally time consuming, but after a very thorough search, these were identified. A relevant result of this search is that when the complexes, CD-AD and CD-FC, are formed, the AD and FC molecules are only partially embedded inside the CD cavity. The totally embedded complexes were found to have significantly higher energies. The semiempirical method, ZINDO, was employed to investigate the effect of complexation on the first electronic excitation of CD anchored to a metal nano-particle. The computational results revealed that after complexation to FC, the transition intensity declines to about 25% of the original value, and after complexation with AD, the intensity drops almost 50%. The tighter binding and transition intensity of CD-AD qualitatively agrees with the experimental result that the addition of AD to a solution of CD and FC restores the fluorescence of CD that was quenched by the addition of FC. A method to evaluate the “hydrophobic force” effect is proposed for future work.

T Cell Motility as Modulator of Interactions with Dendritic Cells

It is well established that the balance of costimulatory and inhibitory signals during interactions with dendritic cells (DCs) determines T cell transition from a naïve to an activated or tolerant/anergic status. Although many of these molecular interactions are well reproduced in reductionist in vitro assays, the highly dynamic motility of naïve T cells in lymphoid tissue acts as an additional lever to fine-tune their activation threshold. T cell detachment from DCs providing suboptimal stimulation allows them to search for DCs with higher levels of stimulatory signals, while storing a transient memory of short encounters. In turn, adhesion of weakly reactive T cells to DCs presenting peptides presented on major histocompatibility complex with low affinity is prevented by lipid mediators. Finally, controlled recruitment of CD8(+) T cells to cognate DC-CD4(+) T cell clusters shapes memory T cell formation and the quality of the immune response. Dynamic physiological lymphocyte motility therefore constitutes a mechanism to mitigate low avidity T cell activation and to improve the search for "optimal" DCs, while contributing to peripheral tolerance induction in the absence of inflammation.

Pathway Semantics: An Algebraic Data Driven Algorithm to Generate Hypotheses about Molecular Patterns Underlying Disease Progression

The overarching goal of the Pathway Semantics Algorithm (PSA) is to improve the in silico identification of clinically useful hypotheses about molecular patterns in disease progression. By framing biomedical questions within a variety of matrix representations, PSA has the flexibility to analyze combined quantitative and qualitative data over a wide range of stratifications. The resulting hypothetical answers can then move to in vitro and in vivo verification, research assay optimization, clinical validation, and commercialization. Herein PSA is shown to generate novel hypotheses about the significant biological pathways in two disease domains: shock / trauma and hemophilia A, and validated experimentally in the latter. The PSA matrix algebra approach identified differential molecular patterns in biological networks over time and outcome that would not be easily found through direct assays, literature or database searches. In this dissertation, Chapter 1 provides a broad overview of the background and motivation for the study, followed by Chapter 2 with a literature review of relevant computational methods. Chapters 3 and 4 describe PSA for node and edge analysis respectively, and apply the method to disease progression in shock / trauma. Chapter 5 demonstrates the application of PSA to hemophilia A and the validation with experimental results. The work is summarized in Chapter 6, followed by extensive references and an Appendix with additional material.

Molecular Dynamics Simulations of Couette flow

The first steps towards developing a continuum-molecular coupled simulations techniques are presented, for the purpose of computing macroscopic systems of confined fluids. The idea is to compute the interface wall-fluid by Molecular Dynamics simulations, where Lennard-Jones potential (and others) have been employed for the molecular interactions, so the usual non slip boundary condition is not specified. Instead, a shear rate can be imposed at the wall, which allows to obtain the properties of the wall material by means of an iterative method. The remaining fluid region will be computed by a spectral hp method. We present MD simulations of a Couette flow, and the results of the developed boundary conditions from the wall fluid interaction.

Quantitative modeling of stochastic systems in molecular biology by using stochastic Petri nets

An integrated understanding of molecular and developmental biology must consider the large number of molecular species involved and the low concentrations of many species in vivo. Quantitative stochastic models of molecular interaction networks can be expressed as stochastic Petri nets (SPNs), a mathematical formalism developed in computer science. Existing software can be used to define molecular interaction networks as SPNs and solve such models for the probability distributions of molecular species. This approach allows biologists to focus on the content of models and their interpretation, rather than their implementation. The standardized format of SPNs also facilitates the replication, extension, and transfer of models between researchers. A simple chemical system is presented to demonstrate the link between stochastic models of molecular interactions and SPNs. The approach is illustrated with examples of models of genetic and biochemical phenomena where the UltraSAN package is used to present results from numerical analysis and the outcome of simulations.

Determination of the lifetime and force dependence of interactions of single bonds between surface-attached CD2 and CD48 adhesion molecules

We studied single molecular interactions between surface-attached rat CD2, a T-lymphocyte adhesion receptor, and CD48, a CD2 ligand found on antigen-presenting cells. Spherical particles were coated with decreasing densities of CD48–CD4 chimeric molecules then driven along CD2-derivatized glass surfaces under a low hydrodynamic shear rate. Particles exhibited multiple arrests of varying duration. By analyzing the dependence of arrest frequency and duration on the surface density of CD48 sites, it was concluded that (i) arrests were generated by single molecular bonds and (ii) the initial bond dissociation rate was about 7.8 s−1. The force exerted on bonds was increased from about 11 to 22 pN; the detachment rate exhibited a twofold increase. These results agree with and extend studies on the CD2–CD48 interaction by surface plasmon resonance technology, which yielded an affinity constant of ≈104 M−1 and a dissociation rate of ≥6 s−1. It is concluded that the flow chamber technology can be an useful complement to atomic force microscopy for studying interactions between isolated biomolecules, with a resolution of about 20 ms and sensitivity of a few piconewtons. Further, this technology might be extended to actual cells.

A yeast genetic system for selecting small molecule inhibitors of protein–protein interactions in nanodroplets

Cellular processes are mediated by complex networks of molecular interactions. Dissection of their role most commonly is achieved by using genetic mutations that alter, for example, protein–protein interactions. Small molecules that accomplish the same result would provide a powerful complement to the genetic approach, but it generally is believed that such molecules are rare. There are several natural products, however, that illustrate the feasibility of this approach. Split-pool synthesis now provides a simple mechanical means to prepare vast numbers of complex, even natural product-like, molecules individually attached to cell-sized polymer beads. Here, we describe a genetic system compatible with split-pool synthesis that allows the detection of cell-permeable, small molecule inhibitors of protein–protein interactions in 100- to 200-nl cell culture droplets, prepared by a recently described technique that arrays large numbers of such droplets. These “nanodroplets” contain defined media, cells, and one or more beads containing ≈100 pmol of a photoreleasable small molecule and a controlled number of cells. The engineered Saccharomyces cerevisiae cells used in this study express two interacting proteins after induction with galactose whose interaction results in cell death in the presence of 5-fluoroorotic acid (inducible reverse two-hybrid assay). Disruption of the interaction by a small molecule allows growth, and the small molecule can be introduced into the system hours before induction of the toxic interaction. We demonstrate that the interaction between the activin receptor R1 and the immunophilin protein FKBP12 can be disrupted by the small molecule FK506 at nanomolar concentrations in nanodroplets. This system should provide a general method for selecting cell-permeable ligands that can be used to study the relevance of protein–protein interactions in living cells or organisms.

Protein-protein interactions in the rigor actomyosin complex.

Since it has not been possible to crystallize the actomyosin complex, the x-ray structures of the individual proteins together with data obtained by fiber diffraction and electron microscopy have been used to build detailed models of filamentous actin (f-actin) and the actomyosin rigor complex. In the f-actin model, a single monomer uses 10 surface loops and two alpha-helices to make sometimes complicated interactions with its four neighbors. In the myosin molecule, both the essential and regulatory light chains show considerable structural homology to calmodulin. General principles are evident in their mode of attachment to the target alpha-helix of the myosin heavy chain. The essential light chain also makes contacts with other parts of the heavy chain and with the regulatory light chain. The actomyosin rigor interface is extensive, involving interaction of a single myosin head with regions on two adjacent actin monomers. A number of hydrophobic residues on the apposing faces of actin and myosin contribute to the main binding site. This site is flanked on three sides by charged myosin surface loops that form predominantly ionic interactions with adjacent regions of actin. Hydrogen bonding is likely to play a significant role in actin-actin and actin-myosin interactions since many of the contacts involve loops. The model building approach used with actomyosin is applicable to other multicomponent assemblies of biological interest and is a powerful method for revealing molecular interactions and providing insights into the mode of action of the assemblies.

Functional implications of protein-protein interactions in icosahedral viruses.

Biological processes often require that a single gene product participate in multiple types of molecular interactions. Viruses with quasiequivalent capsids provide an excellent paradigm for studying such phenomena because identical protein subunits are found in different structural environments. Differences in subunit joints may be controlled by protein segments, duplex or single-stranded RNA, metal ions, or some combination of these. Each of the virus groups examined display a distinctive mechanism for switching interface interactions, illustrating the magnitude of options that are likely to be found in other biological systems. In addition to determining capsid morphology, assembly controls the timing of autocatalytic maturation cleavage of the viral subunits that is required for infectivity in picorna-, noda-, and tetraviruses. The mechanism of assembly-dependent cleavage is conserved in noda- and tetraviruses, although the quaternary structures of the capsids are different as are the molecular switches that control subunit interfaces. The function of the cleavage in picorna-, noda-, and tetraviruses is probably to release polypeptides that participate in membrane translocation of RNA.

Analysis of RNA-protein interactions by flow cytometry

Flow cytometry, in combination with advances in bead coding technologies, is maturing as a powerful high-throughput approach for analyzing molecular interactions. Applications of this technology include antibody assays and single nucleotide polymorphism mapping. This review describes the recent development of a microbead flow cytometric approach to analyze RNA-protein interactions and discusses emerging bead coding strategies that together will allow genome-wide identification of RNA-protein complexes. The microbead flow cytometric approach is flexible and provides new opportunities for functional genomic studies and small-molecule screening.

Identifying the molecular phenotype of renal progenitor cells

Although many of the molecular interactions in kidney development are now well understood, the molecules involved in the specification of the metanephric mesenchyme from surrounding intermediate mesoderm and, hence, the formation of the renal progenitor population are poorly characterized. In this study, cDNA microarrays were used to identify genes enriched in the murine embryonic day 10.5 (E10.5) uninduced metanephric mesenchyme, the renal progenitor population, in comparison with more rostral derivatives of the intermediate mesoderm. Microarray data were analyzed using R statistical software to determine accurately genes differentially expressed between these populations. Microarray outliers were biologically verified, and the spatial expression pattern of these genes at E10.5 and subsequent stages of early kidney development was determined by RNA in situ hybridization. This approach identified 21 genes preferentially expressed by the E10.5 metanephric mesenchyme, including Ewing sarcoma homolog, 14-3-3 theta, retinoic acid receptor-alpha, stearoyl-CoA desaturase 2, CD24, and cadherin-11, that may be important in formation of renal progenitor cells. Cell surface proteins such as CD24 and cadherin-11 that were strongly and specifically expressed in the uninduced metanephric mesenchyme and mark the renal progenitor population may prove useful in the purification of renal progenitor cells by FACS. These findings may assist in the isolation and characterization of potential renal stem cells for use in cellular therapies for kidney disease.

Molecular determinants of ginkgolide binding in the glycine receptor pore

Ginkgolides are potent blockers of the glycine receptor Cl- channel (GlyR) pore. We sought to identify their binding sites by comparing the effects of ginkgolides A, B and C and bilobalide on alpha 1, alpha 2, alpha 1 beta and alpha 2 beta GlyRs. Bilobalide sensitivity was drastically reduced by incorporation of the beta subunit. In contrast, the sensitivities to ginkgolides B and C were enhanced by beta subunit expression. However, ginkgolide A sensitivity was increased in the alpha 2 beta GlyR relative to the alpha 2 GlyR but not in the alpha 1 beta GlyR relative to the alpha 1 GlyR. We hypothesised that the subunit-specific differences were mediated by residue differences at the second transmembrane domain 2' and 6' pore-lining positions. The increased ginkgolide A sensitivity of the alpha 2 beta GlyR was transferred to the alpha 1 beta GlyR by the G2'A (alpha 1 to alpha 2 subunit) substitution. In addition, the alpha 1 subunit T6'F mutation abolished inhibition by all ginkgolides. As the ginkgolides share closely related structures, their molecular interactions with pore-lining residues were amenable to mutant cycle analysis. This identified an interaction between the variable R2 position of the ginkgolides and the 2' residues of both alpha 1 and beta subunits. These findings provide strong evidence for ginkgolides binding at the 2' pore-lining position.

Molecular dissection of the Munc18c/syntaxin4 interaction: Implications for regulation of membrane trafficking

Sec1p/Munc18 (SM) proteins are believed to play an integral role in vesicle transport through their interaction with SNAREs. Different SM proteins have been shown to interact with SNAREs via different mechanisms, leading to the conclusion that their function has diverged. To further explore this notion, in this study, we have examined the molecular interactions between Munc18c and its cognate SNAREs as these molecules are ubiquitously expressed in mammals and likely regulate a universal plasma membrane trafficking step. Thus, Munc18c binds to monomeric syntaxin4 and the N-terminal 29 amino acids of syntaxin4 are necessary for this interaction. We identified key residues in Munc18c and syntaxin4 that determine the N-terminal interaction and that are consistent with the N-terminal binding mode of yeast proteins Sly1p and Sed5p. In addition, Munc18c binds to the syntaxin4/SNAP23/VAMP2 SNARE complex. Pre-assembly of the syntaxin4/Munc18c dimer accelerates the formation of SNARE complex compared to assembly with syntaxin4 alone. These data suggest that Munc18c interacts with its cognate SNAREs in a manner that resembles the yeast proteins Sly1p and Sed5p rather than the mammalian neuronal proteins Munc18a and syntaxin1a. The Munc18c-SNARE interactions described here imply that Munc18c could play a positive regulatory role in SNARE assembly.