Terbium-mediated footprinting probes a catalytic conformational switch in the antigenomic hepatitis delta virus ribozyme

The two forms of the hepatitis delta virus ribozyme are derived from the genomic and antigenomic RNA strands of the human hepatitis delta virus (HDV), where they serve a crucial role in pathogen replication by catalyzing site-specific self-cleavage reactions. The HDV ribozyme requires divalent metal ions for formation of its tertiary structure, consisting of a tight double-nested pseudoknot, and for efficient self- (or cis-) cleavage. Comparison of recently solved crystal structures of the cleavage precursor and 3' product indicates that a significant conformational switch is required for catalysis by the genomic HDV ribozyme. Here, we have used the lanthanide metal ion terbium(III) to footprint the precursor and product solution structures of the cis-acting antigenomic HDV ribozyme. Inhibitory Tb(3+) binds with high affinity to similar sites on RNA as Mg(2+) and subsequently promotes slow backbone scission. We find subtle, yet significant differences in the terbium(III) footprinting pattern between the precursor and product forms of the antigenomic HDV ribozyme, consistent with differences in conformation as observed in the crystal structures of the genomic ribozyme. In addition, UV melting profiles provide evidence for a less tight tertiary structure in the precursor. In both the precursor and product we observe high-affinity terbium(III) binding sites in joining sequence J4/2 (Tb(1/2) approximately 4 microM) and loop L3, which are key structural components forming the catalytic core of the HDV ribozyme, as well as in several single-stranded regions such as J1/2 and the L4 tetraloop (Tb(1/2) approximately 50 microM). Sensitized luminescence spectroscopy confirms that there are at least two affinity classes of Tb(3+) binding sites. Our results thus demonstrate that a significant conformational change accompanies catalysis in the antigenomic HDV ribozyme in solution, similar to the catalytic conformational switch observed in crystals of the genomic form, and that structural and perhaps catalytic metal ions bind close to the catalytic core.     

Local helix content and RNA-binding activity of the N-terminal leucine-repeat region of hepatitis delta antigen

Hepatitis delta virus (HDV) is a satellite virus of the hepatitis B virus (HBV) which provides the surface antigen for the viral coat. Our results show that the N-terminal leucine-repeat region of hepatitis delta antigen (HDAg), encompassing residues 24–50, binds to the autolytic domain of HDV genomic RNA and attenuates its autolytic activity. The solution conformation of a synthetic peptide corresponding to residues 24–50 of HDAg as determined by two-dimensional ¹H NMR and circular dichroism techniques is found to be an -helix. The local helix content of this peptide was analyzed by NOEs and coupling constants. Mutagenesis studies indicate that Lys³⁸, Lys³⁹, and Lys⁴⁰ within this -helical peptide may be directly involved in RNA binding. A structural knowledge of the N-terminal leucine-repeat region of HDAg thus provides a molecular basis for understanding its role in the interaction with RNA.     

Structure and interactions of RecA: plasticity revealed by molecular dynamics simulations

Eleven independent simulations, each involving three consecutive molecules in the RecA filament, carried out on the protein from Mycobacterium tuberculosis, Mycobacterium smegmatis and Escherichia coli and their Adenosine triphosphate (ATP) complexes, provide valuable information which is complementary to that obtained from crystal structures, in addition to confirming the robust common structural framework within which RecA molecules from different eubacteria function. Functionally important loops, which are largely disordered in crystal structures, appear to adopt in each simulation subsets of conformations from larger ensembles. The simulations indicate the possibility of additional interactions involving the P-loop which remains largely invariant. The phosphate tail of the ATP is firmly anchored on the loop while the nucleoside moiety exhibits substantial structural variability. The most important consequence of ATP binding is the movement of the ‘switch’ residue. The relevant simulations indicate the feasibility of a second nucleotide binding site, but the pathway between adjacent molecules in the filament involving the two nucleotide binding sites appears to be possible only in the mycobacterial proteins.     

Mutant bovine odorant-binding protein: Temperature affects the protein stability and dynamics as revealed by infrared spectroscopy and molecular dynamics simulations

Bovine odorant-binding protein (bOBP), a member of the lipocalin family, presents the so-called 3D "domain-swapped" protein structure. In fact, in solution, it appears as a dimer in which each monomer is composed by the classical lipocalin fold, with a central beta-barrel followed by a stretch of residues and the alpha-helix domain protruding out of the barrel and crossing the dimer interface. Recently, a deswapped mutant form of bOBP was obtained, in which a Gly residue was inserted after position 121 and the two residues in position 64 and 156 were replaced by Cys residues for restoring the disulfide bridge common to the lipocalin family. In this work, we used Fourier transform infrared spectroscopy and molecular dynamics simulations to investigate the effect of temperature on the structural stability and conformational dynamics of the mutant bOBP. The spectroscopic and molecular simulation data pointed out that the hydrophobic regions of the protein matrix appear to be an important factor for the protein stability and integrity. In addition, it was also found that the mutant bOBP is significantly stabilized by the binding of the ligand, which may have an impact on the biological function of bOBP. The obtained results will allow for a better use of this protein as probe for the design of advanced protein-based biosensors for the detection of compounds used in the fabrication of explosive powders.     

Dissociation mechanism of GDP from Cdc42 via DOCK9 revealed by molecular dynamics simulations

Cell division control protein 42 homolog (Cdc42) influences a variety of cellular responses such as cell migration and polarity. Deregulation of Cdc42 has been associated with several human diseases and developmental disorders. Over-activation of Cdc42 through guanine nucleotide exchange factor (GEF) is a critical event for Cdc42 involved cancer metastasis. Members of DOCK family of GEF are important activators of Cdc42. However, this activation mechanism is still unknown. Molecular dynamics (MD) simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) calculations were employed to investigate the central step of the activation of Cdc42: the dissociation mechanism of GDP from Cdc42 via DOCK9. Simulation results show that Mg²⁺ ion has a remarkable influence on the conformational change of switch I of Cdc42 through residue Pro34 which functions as a “clasp” to control the flexibility of switch I. In the GDP dissociation process, the Mg²⁺ ion leave first to result in a suitable conformation of Cdc42 for following DOCK9 binding to. When DOCK9 binds to Cdc42, it changes the orientations of residues Lys16, Thr17, Cys18 and Phe28 of Cdc42 to weaken the interactions between Cdc42 and GDP to release GDP. This study first elucidates the dissociation mechanism of GDP from Cdc42 via DOCK9 and identifies the essential residues of Cdc42 in this process. These simulation results are consistent with the recent findings of biochemical and amino acid mutational studies, and the observations are beneficial to understand the activation mechanism of Cdc42 and to provide insights for designing compounds targeting on Cdc42 related cancer metastasis.     

Cations and hydration in catalytic RNA: molecular dynamics of the hepatitis delta virus ribozyme

The hepatitis delta virus (HDV) ribozyme is an RNA enzyme from the human pathogenic HDV. Cations play a crucial role in self-cleavage of the HDV ribozyme, by promoting both folding and chemistry. Experimental studies have revealed limited but intriguing details on the location and structural and catalytic functions of metal ions. Here, we analyze a total of approximately 200 ns of explicit-solvent molecular dynamics simulations to provide a complementary atomistic view of the binding of monovalent and divalent cations as well as water molecules to reaction precursor and product forms of the HDV ribozyme. Our simulations find that an Mg2+ cation binds stably, by both inner- and outer-sphere contacts, to the electronegative catalytic pocket of the reaction precursor, in a position to potentially support chemistry. In contrast, protonation of the catalytically involved C75 in the precursor or artificial placement of this Mg2+ into the product structure result in its swift expulsion from the active site. These findings are consistent with a concerted reaction mechanism in which C75 and hydrated Mg2+ act as general base and acid, respectively. Monovalent cations bind to the active site and elsewhere assisted by structurally bridging long-residency water molecules, but are generally delocalized.     

RNA simulations: probing hairpin unfolding and the dynamics of a GNRA tetraloop

Simulations of an RNA hairpin containing a GNRA tetraloop were conducted to allow the characterization of its secondary structure formation and dynamics. Ten 10 ns trajectories of the folded hairpin 5'-GGGC[GCAA]GCCU-3' were generated using stochastic dynamics and the GB/SA implicit solvent model at 300 K. Overall, we find the stem to be a very stable subunit of this molecule, whereas multiple loop conformations and transitions between them were observed. These trajectories strongly suggest that extension of the C6 base away from the loop occurs cooperatively with an N-type-->S-type sugar pucker conversion in that residue and that similar pucker transitions are necessary to stabilize other looped-out bases. In addition, a short-lived conformer with an extended fourth loop residue (A8) lacking this stabilizing 2'-endo pucker mode was observed. Results of thermal perturbation at 400 K support this model of loop dynamics. Unfolding trajectories were produced using this same methodology at temperatures of 500 to 700 K. The observed unfolding events display three-state behavior kinetically (including native, globular, and unfolded populations) and, based on these observations, we propose a folding mechanism that consists of three distinct events: (i) collapse of the random unfolded structure and sampling of the globular state; (ii) passage into the folded region of configurational space as stem base-pairs form and gain helicity; and (iii) attainment of proper loop geometry and organization of loop pairing and stacking interactions. These results are considered in the context of current experimental knowledge of this and similar nucleic acid hairpins.     

The folding pathway of the genomic hepatitis delta virus ribozyme is dominated by slow folding of the pseudoknots

Hepatitis delta virus (HDV) replicates by a double rolling-circle mechanism that requires self-cleavage by closely related genomic and antigenomic versions of a ribozyme. We have previously shown that the uncleaved genomic ribozyme is subject to a variety of alternative (Alt) pairings. Sequence upstream of the ribozyme can regulate self-cleavage activity by formation of an Alt 1 ribozyme-containing structure that severely inhibits self-cleavage, or a P(-1) self-structure that permits rapid self-cleavage. Here, we test three other alternative pairings: Alt P1, Alt 2, and Alt 3. Alt P1 and Alt 3 contain primarily ribozyme-ribozyme interactions, while Alt 2 involves ribozyme-flanking sequence interaction. A number of single and double mutant ribozymes were prepared to increase or decrease the stability of the alternative pairings, and rates of self-cleavage were determined. Results of these experiments were consistent with the existence of the proposed alternative pairings and their ability to inhibit the overall rate of native ribozyme folding. Local misfolds are treated as internal equilibrium constants in a binding polynomial that modulates the intrinsic rate of cleavage. This model of equilibrium effects of misfolds should be general and apply to other ribozymes. All of the alternative pairings sequester a pseudoknot-forming strand. Folding of ribozymes containing Alt P1 and Alt 2 was accelerated by urea as long as the native ribozyme fold was sufficiently stable, while folding of mutants in which both of these alternative pairings had been removed were not stimulated by urea. This behavior suggests that the pseudoknots form by capture of an unfolded or appropriately rearranged alternative pairing by its complementary native strand. Fast-folding mutants were prepared by either weakening alternative pairings or by strengthening native pairings. A kinetic model was developed that accommodates these features and explains the position of the rate-limiting step for the G11C mutant. Implications of these results for structural and dynamic studies of the uncleaved HDV ribozyme are discussed.     

Identification and functional characterization of the nascent RNA contacting residues of the hepatitis C virus RNA-dependent RNA polymerase

Understanding how the hepatitis C virus (HCV) RNA-dependent RNA polymerase (RdRp) interacts with nascent RNA would provide valuable insight into the virus's mechanism for RNA synthesis. Using a peptide mass fingerprinting method and affinity capture of peptides reversibly cross-linked to an alkyn-labeled nascent RNA, we identified a region below the Δ1 loop in the fingers domain of the HCV RdRp that contacts the nascent RNA. A modification protection assay was used to confirm the assignment. Several mutations within the putative nascent RNA binding region were generated and analyzed for RNA synthesis in vitro and in the HCV subgenomic replicon. All mutations tested within this region showed a decrease in primer-dependent RNA synthesis and decreased stabilization of the ternary complex. The results from this study advance our understanding of the structure and function of the HCV RdRp and the requirements for HCV RNA synthesis. In addition, a model of nascent RNA interaction is compared with results from structural studies.     

A detailed unfolding pathway of a (beta/alpha)8-barrel protein as studied by molecular dynamics simulations

The (beta/alpha)(8)-barrel is the most common protein fold. Similar structural properties for folding intermediates of (beta/alpha)(8)-barrel proteins involved in tryptophan biosynthesis have been reported in a number of experimental studies; these intermediates have the last two beta-strands and three alpha-helices partially unfolded, with other regions of the polypeptide chain native-like in conformation. To investigate the detailed folding/unfolding pathways of these (beta/alpha)(8)-barrel proteins, temperature-induced unfolding simulations of N-(5'-phosphoribosyl)anthranilate isomerase from Escherichia coli were carried out using a special-purpose parallel computer system. Unfolding simulations at five different temperatures showed a sequential unfolding pathway comprised of several events. Early events in unfolding involved disruption of the last two strands and three helices, producing an intermediate ensemble similar to those detected in experimental studies. Then, denaturation of the first two betaalpha units and separation of the sixth strand from the fifth took place independently. The remaining central betaalphabetaalphabeta module persisted the longest during all simulations, suggesting an important role for this module as the incipient folding scaffold. Our simulations also predicted the presence of a nucleation site, onto which several hydrophobic residues condensed forming the foundation for the central betaalphabetaalphabeta module.     

Structural differences between allelic variants of the ovine prion protein revealed by molecular dynamics simulations

We have modeled ovine prion protein (residues 119-233) based on NMR structures of PrP from other mammalian species. Modeling of the C-terminal domain of ovine PrP predicts three helices: helix-1 (residues 147-155), flanked by two short beta-strands; helix-2 (residues 176-197), and helix-3 (residues 203-229). Molecular dynamics simulations on this model of ovine PrP have determined structural differences between allelic variants. At neutral pH, limited root mean-squared (RMS) fluctuations were seen in the region of helix-1; between beta-strand-2 and residue 171, and the loop connecting helix-2 and helix-3. At low pH, these RMS fluctuations increased and showed allelic variation. The extent of RMS fluctuation between beta-strand 2 and residue 171 was ARR > ARQ > VRQ. This order was reversed for the loop region connecting helix-2 and helix-3. Although all three variants have the potential to display an extended helix at the C-terminal region of helix-1, the major influence of the VRQ allele was to restrict the conformations of the Asn162 and Arg139 side-chains. Variations observed in the simulations in the vicinity of helix-1 correlated with reactivity of C-terminal specific anti-PrP monoclonal antibodies with peripheral blood cells from scrapie-susceptible and -resistant genotypes of sheep: cells from VRQ homozygous sheep showed uniform reactivity, while cells from ARQ and ARR homozygous sheep showed variable binding. Our data show that molecular dynamics simulations can be used to determine structural differences between allelic variants of ovine PrP. The binding of anti-PrP monoclonal antibodies to ovine blood cells may validate these structural predictions.     

Insight into the early stages of thermal unfolding of peanut agglutinin by molecular dynamics simulations

Peanut agglutinin is a homotetrameric nonglycosylated protein. The protein has a unique open quaternary structure. Molecular dynamics simulations have been employed to follow the atomistic details of its unfolding at different temperatures. The early events of the deoligomerization of the protein have been elucidated in the present study. Simulation trajectories of the monomer as well as those of the tetramer have been compared and the tetramer is found to be substantially more stable than its monomeric counterpart. The tetramer shows retention of most of its secondary structure but considerable loss of the tertiary structure at high temperature. This observation implies the generation of a molten globule-like intermediate in the later stages of deoligomerization. The quaternary structure of the protein has weakened to a large extent, but none of the subunits are separated. In addition, the importance of the metal-binding to the stability of the protein structure has also been investigated. Binding of the metal ions not only enhances the local stability of the metal-ion binding loop, but also imparts a global stability to the overall structure. The dynamics of different interfaces vary significantly as probed through interface clusters. The differences are substantially enhanced at higher temperatures. The dynamics and the stability of the interfaces have been captured mainly by cluster analysis, which has provided detailed information on the thermal deoligomerization of the protein.     

Assessment of metal-assisted nucleophile activation in the hepatitis delta virus ribozyme from molecular simulation and 3D-RISM

The hepatitis delta virus ribozyme is an efficient catalyst of RNA 2′-O-transphosphorylation and has emerged as a key experimental system for identifying and characterizing fundamental features of RNA catalysis. Recent structural and biochemical data have led to a proposed mechanistic model whereby an active site Mg²⁺ ion facilitates deprotonation of the O2′ nucleophile, and a protonated cytosine residue (C75) acts as an acid to donate a proton to the O5′ leaving group as noted in a previous study. This model assumes that the active site Mg²⁺ ion forms an inner-sphere coordination with the O2′ nucleophile and a nonbridging oxygen of the scissile phosphate. These contacts, however, are not fully resolved in the crystal structure, and biochemical data are not able to unambiguously exclude other mechanistic models. In order to explore the feasibility of this model, we exhaustively mapped the free energy surfaces with different active site ion occupancies via quantum mechanical/molecular mechanical (QM/MM) simulations. We further incorporate a three-dimensional reference interaction site model for the solvated ion atmosphere that allows these calculations to consider not only the rate associated with the chemical steps, but also the probability of observing the system in the presumed active state with the Mg²⁺ ion bound. The QM/MM results predict that a pathway involving metal-assisted nucleophile activation is feasible based on the rate-controlling transition state barrier departing from the presumed metal-bound active state. However, QM/MM results for a similar pathway in the absence of Mg²⁺ are not consistent with experimental data, suggesting that a structural model in which the crystallographically determined Mg²⁺ is simply replaced with Na⁺ is likely incorrect. It should be emphasized, however, that these results hinge upon the assumption of the validity of the presumed Mg²⁺-bound starting state, which has not yet been definitively verified experimentally, nor explored in depth computationally. Thus, further experimental and theoretical study is needed such that a consensus view of the catalytic mechanism emerges.     

Dissolution of cellulose in ionic liquid and water mixtures as revealed by molecular dynamics simulations

Abstract

Increasing population growth and industrialization are continuously oppressing the existing energy resources, elevating the pollution and global fuel demand. Various alternate energy resources can be utilized to cope with these problems in an environment-friendly fashion. Currently, bioethanol (sugarcane, corn-derived) is one of the most widely consumed biofuels in the world. Lignocellulosic biomass is yet another attractive resource for sustainable bioethanol production. Pretreatment step plays a crucial role in the lignocellulose to bioethanol conversion by enhancing cellulose susceptibility to enzymatic hydrolysis. However, economical lignocellulose pretreatment still remains a challenging job. Ionic liquids (ILs), especially 1-ethyl-3-methylimidazolium acetate (EmimAc), is an efficient solvent for cellulose dissolution with improved enzymatic saccharification kinetics. To increase the process efficiency as well as recyclability of IL, water is shown as a compatible cosolvent for lignocellulosic pretreatment. The performance analysis of IL–water mixture based on the molecular level understanding may help to design effective pretreatment solvents. In this study, all-atom molecular dynamics simulation has been performed using EmimAc–water mixtures to understand the behavior of cellulose microcrystal containing eight glucose octamers at room and pretreatment temperatures. High-temperature simulation results show effective cellulose chain separation where cellulose–acetate interaction is found to be the driving force behind dissolution. It is also observed that pretreatment with 50 and 80% IL mixture is efficient in decreasing cellulose crystallinity. At a high IL concentration, water exists in a clustered network which gradually spans into the medium with increasing water fraction leading to loss of its cosolvation activity.

Communicated by Ramaswamy H. Sarma     

Mechanisms for allosteric activation of protease DegS by ligand binding and oligomerization as revealed from molecular dynamics simulations

A local perturbation of a protein may lead to functional changes at some distal site, a phenomenon denoted as allostery. Here, we study the allosteric control of a protease using molecular dynamics simulations. The system considered is the bacterial protein DegS which includes a protease domain activated on ligand binding to an adjacent PDZ domain. Starting from crystallographic structures of DegS homo‐trimers, we perform simulations of the ligand‐free and ‐bound state of DegS at equilibrium. Considering a single protomer only, the trimeric state was mimicked by applying restraints on the residues in contact with other protomers in the DegS trimer. In addition, the bound state was also simulated without any restraints to mimic the monomer. Our results suggest that not only ligand release but also disassembly of a DegS trimer inhibits proteolytic activity. Considering various observables for structural changes, we infer allosteric pathways from the interface with other protomers to the active site. Moreover, we study how ligand release leads to (i) catalytically relevant changes involving residues 199–201 and (ii) a transition from a stretched to a bent conformation for residues 217–219 (which prohibits proper substrate binding). Finally, based on ligand‐induced C_α shifts we identify residues in contact with other protomers in the DegS trimer that likely transduce the perturbation from ligand release from a given protomer to adjacent protomers. These residues likely play a key role in the experimentally known effect of ligand release from a protomer on the proteolytic activity of the other protomers. Proteins 2016; 84:1690–1705. © 2016 Wiley Periodicals, Inc.     

Evaluating the suitability of RNA intervention mechanism exerted by some flavonoid molecules against dengue virus MTase RNA capping site: a molecular docking,molecular dynamics simulation,and binding free energy study

Abstract

Dengue virus (DENV) is one of the most dangerous mosquito-borne human pathogens known to the mankind. Currently, no vaccines or standard therapy is avaliable to treate DENV infection. This makes the drug development against DENV more significant and challenging. The MTase domain of DENV RNA RdRp NS5 is a promising drug target, because this domain hosts the RNA capping process of DENV RNA to escape from human immune system. In the present study, we have analysed the RNA intervention mechanism exerted by flavoniod molecules against NS5 MTase RNA capping site by using molecular docking, molecular dynamics simulation and the binding free energy calculations. The results from the docking analysis confirmed that the RNA intervention mecanism is exerted by the quercetagetin (QGN) molecule with all necessary intermolecular interactions and high binding affinity. Notably, QGN forms strong hydrogen bonding interactions with Asn18, Leu20 and Ser150 residues and π???π stacking interaction with Phe25 residue. The apo and QGN bound NS5 MTase and QGN-NS5 MTase complex were used for MD simulation. The results of MD simulation reveal that the RMSD and RMSF values of QGN-MTase complex have increased on comparing the apo protein due to the effect of ligand binding. The binding free energy calulation includes prediction of total binding free energy of ligand-protein complex and per-residue free energy decomposition. The QGN binding to NS5 MTase affects it’s native motion, this result is found from Principal component analysis.

Communicated by Ramaswamy H. Sarma     

Structural characterization of the 5' untranslated RNA of hepatitis C virus by vibrational spectroscopy

Raman and FTIR spectroscopy have been used to characterize the structure of 5'untranslated region (5'UTR, 342-mer RNA) of the HCV genome. The study of the 750-850 cm(-1) Raman spectral domain of the ribose-phosphate backbone reveals that the percentage of nucleobases involved in double helix-loop junctions is 19+/-1%, which is very close to that of a theoretical secondary structure model (18.7%) proposed on the basis of comparative sequence analysis and thermodynamic modelling. In addition, about 68+/-2% of the bases are helically ordered having C(3')-endo ribofuranose pucker. FTIR-monitored H/D exchange provides the following results: (a) base-paired guanine and cytosine nucleobases show the lowest rate of isotopic exchange, and some synchronous intensity changes of marker bands of A.U pair and single stranded adenine are consistent with the presence of A(*)A.U triplets; (b) the vibrational coupling between the ribose ether C-O stretching and 2'OH bending motions reveals that helical regions of 5'UTR RNA are characterized by hydrogen bonding between the 2'OH ribose groups and the ether oxygen atoms of neighbouring ribose residues.