RNA force field with accuracy comparable to state-of-the-art protein force fields.

Molecular dynamics (MD) simulation has become a powerful tool for characterizing at an atomic level of detail the conformational changes undergone by proteins. The application of such simulations to RNA structures, however, has proven more challenging, due in large part to the fact that the physical models ("force fields") available for MD simulations of RNA molecules are substantially less accurate in many respects than those currently available for proteins. Here, we introduce an extensive revision of a widely used RNA force field in which the parameters have been modified, based on quantum mechanical calculations and existing experimental information, to more accurately reflect the fundamental forces that stabilize RNA structures. We evaluate these revised parameters through long-timescale MD simulations of a set of RNA molecules that covers a wide range of structural complexity, including single-stranded RNAs, RNA duplexes, RNA hairpins, and riboswitches. The structural and thermodynamic properties measured in these simulations exhibited dramatically improved agreement with experimentally determined values. Based on the comparisons we performed, this RNA force field appears to achieve a level of accuracy comparable to that of state-of-the-art protein force fields, thus significantly advancing the utility of MD simulation as a tool for elucidating the structural dynamics and function of RNA molecules and RNA-containing biological assemblies.

Molecular mechanisms for the regulation of histone mRNA stem-loop-binding protein by phosphorylation.

Replication-dependent histone mRNAs end with a conserved stem loop that is recognized by stem-loop-binding protein (SLBP). The minimal RNA-processing domain of SLBP is phosphorylated at an internal threonine, and Drosophila SLBP (dSLBP) also is phosphorylated at four serines in its 18-aa C-terminal tail. We show that phosphorylation of dSLBP increases RNA-binding affinity dramatically, and we use structural and biophysical analyses of dSLBP and a crystal structure of human SLBP phosphorylated on the internal threonine to understand the striking improvement in RNA binding. Together these results suggest that, although the C-terminal tail of dSLBP does not contact the RNA, phosphorylation of the tail promotes SLBP conformations competent for RNA binding and thereby appears to reduce the entropic penalty for the association. Increased negative charge in this C-terminal tail balances positively charged residues, allowing a more compact ensemble of structures in the absence of RNA.

Structure and function of pre-mRNA 5'-end capping quality control and 3'-end processing.

Messenger RNA precursors (pre-mRNAs) are produced as the nascent transcripts of RNA polymerase II (Pol II) in eukaryotes and must undergo extensive maturational processing, including 5'-end capping, splicing, and 3'-end cleavage and polyadenylation. This review will summarize the structural and functional information reported over the past few years on the large machinery required for the 3'-end processing of most pre-mRNAs, as well as the distinct machinery for the 3'-end processing of replication-dependent histone pre-mRNAs, which have provided great insights into the proteins and their subcomplexes in these machineries. Structural and biochemical studies have also led to the identification of a new class of enzymes (the DXO family enzymes) with activity toward intermediates of the 5'-end capping pathway. Functional studies demonstrate that these enzymes are part of a novel quality surveillance mechanism for pre-mRNA 5'-end capping. Incompletely capped pre-mRNAs are produced in yeast and human cells, in contrast to the general belief in the field that capping always proceeds to completion, and incomplete capping leads to defects in splicing and 3'-end cleavage in human cells. The DXO family enzymes are required for the detection and degradation of these defective RNAs.

Novel conjugated microporous polymers for efficient tetracycline adsorption: insights from theoretical investigations.

This paper presents a detailed theoretical understanding of the noncovalent interactions between antibiotics tetracycline and conjugated microporous polymer (CMP), which is important to understand the recent experimental finding of efficient removal of antibiotics by CMP materials. We show that the co-work of π-π and H-π interactions determines the final equilibrium structures, when a tetracycline molecule spontaneously adsorbs to the surface or within the pores of the CMP network at physisorption distances. The binding energies for tetracycline/CMP systems are calculated to be -0.31 âˆ¼ -1.15 eV, demonstrating the reliability of the adsorption. The electronic structures of CMP nanostructures remain basically undamaged upon the tetracycline adsorption. The replacement of benzothiadiazole unit with S and N heteroatoms to the phenyl moiety in the linker effectively enhanced the molecular polarity of CMP molecule and increases the interaction area between tetracycline and CMP network, consequently enhancing the average binding energies notably. Our calculations provide useful theoretical guidance for design of novel carbon-based porous adsorbents with good adsorption performance to remove residual tetracycline and other antibiotics in water.

Acene-modified triphenylamine dyes for dye-sensitized solar cells: a computational study.

A series of metal-free acene-modified triphenylamine dyes (benzene to pentacene, denoted as TPA-AC1 to TPA-AC5) are investigated as organic sensitizers for application in dye-sensitized solar cells (DSSCs). A combination of density functional theory (DFT), density functional tight-binding (DFTB), and time-dependent DFT (TDDFT) approaches is employed. The effects of acene units on the spectra and electrochemical properties of the acene-modified TPA organic dyes are demonstrated. The dye/(TiO(2))(46) anatase nanoparticle systems are also simulated to show the electronic structures at the interface. The results show that from TPA-AC1 to TPA-AC5 with increasing sizes of the acenes, the absorption and fluorescence spectra are systematically broadened and red-shifted, but the oscillator strength and electron injection properties are reduced. The molecular orbital contributions show increasing localization on the bridging acene units from TPA-AC1 to TPA-AC5. From the theoretical examination of some key parameters including free enthalpy related to the electron injection, light-harvesting efficiency, and the shift of semiconductor conduction band, TPA-AC3 with an anthracene moiety demonstrates a balance of the above crucial factors. TPA-AC3 is expected to be a promising dye with desirable energetic and spectroscopic parameters in the DSSC field, which is consistent with recent experimental work. This study is expected to deepen our understanding of TPA-based organic dyes and assist the molecular design of new metal-free dyes for the further optimization of DSSCs.

Development of Force Field Parameters for the Simulation of Single- and Double-Stranded DNA Molecules and DNA-Protein Complexes.

Although molecular dynamics (MD) simulations have been used extensively to study the structural dynamics of proteins, the role of MD simulation in studies of nucleic acid based systems has been more limited. One contributing factor to this disparity is the historically lower level of accuracy of the physical models used in such simulations to describe interactions involving nucleic acids. By modifying nonbonded and torsion parameters of a force field from the Amber family of models, we recently developed force field parameters for RNA that achieve a level of accuracy comparable to that of state-of-the-art protein force fields. Here we report force field parameters for DNA, which we developed by transferring nonbonded parameters from our recently reported RNA force field and making subsequent adjustments to torsion parameters. We have also modified the backbone charges in both the RNA and DNA parameter sets to make the treatment of electrostatics compatible with our recently developed variant of the Amber protein and ion force field. We name the force field resulting from the union of these three parameter sets (the new DNA parameters, the revised RNA parameters, and the existing protein and ion parameters) DES-Amber. Extensive testing of DES-Amber indicates that it can describe the thermal stability and conformational flexibility of single- and double-stranded DNA systems with a level of accuracy comparable to or, especially for disordered systems, exceeding that of state-of-the-art nucleic acid force fields. Finally, we show that, in certain favorable cases, DES-Amber can be used for long-timescale simulations of protein-nucleic acid complexes.

Modification of benzoindenothiophene-based organic dye with fused thiophenes for efficient dye-sensitized solar cells.

Recently low-energy-gap benzoindenothiophene (BIT)-based organic dyes have been experimentally sensitized to dye-sensitized solar cells (DSSCs) with impressive 10.9% power conversion efficiency. This paper presents a computational study of the modification of BIT-based dyes with fused thiophene moieties to build novel low band gap sensitizers. Density functional theory (DFT), tight-binding DFT, and time dependent DFT (TDDFT) approaches are used to demonstrate the electronic and optical properties of the BIT dyes and dye/(TiO2)46 complexes. Our calculations show that the structural modification by using fused thiophenes can effectively lower the band gap of the BIT dyes by 0.07-0.12 eV and affect the optical properties of BIT dyes. Enlarging the thiophene unit in BIT with thienothiophene and dithienothiophene improves the oscillator strength by 14%-25%, while the lowest-energy absorption peak basically remains at 559 nm. The incorporation of cyclopentadithiophene unit leads to a significant 47 nm red-shift of absorption peak and a 25% enhanced oscillator strength, compared to the original BIT dye. Those fused thiophenes modified BIT dyes also demonstrate ideal molecular orbital distribution patterns and ultra-fast injection time at the dye/(TiO2)46 interface. Our calculations provide useful guidance for the molecular design of novel naphthalene-based dyes for DSSC optimizations.

Theoretical investigation of triphenylamine dye/titanium dioxide interface for dye-sensitized solar cells.

The structural, electronic and optical features of two metal-free triphenylamine (TPA) organic dyes (namely C206 and C217) before and after binding to a TiO(2) anatase nanoparticle have been investigated in detail, as a model for the corresponding dye-sensitized solar cells (DSSCs). The combination of density functional tight-binding (DFTB), density functional theory (DFT), and time-dependent DFT (TDDFT) approaches are employed. To understand the effects of the linker part in the TPA organic dyes on the energy conversion efficiency of the DSSCs, C217 and C206, which share the same donor and anchor parts but different linker parts, are theoretically evaluated. Our results show that compared with C206 containing just one thienothiophene unit as the linker, for C217 the introduction of one electron-rich 3,4-ethylenedioxythiophene group to the linker part results in stronger couplings with the TiO(2) conduction band and more efficient electron transfer. This difference contributes to the higher efficiency of C217 in DSSCs experiments. This study is expected to assist the molecular design of new and more efficient TPA-based organic dyes for the optimization of the DSSCs.

Design of novel conjugated microporous polymers for efficient adsorptive desulfurization of small aromatic sulfur molecules.

This paper presents a computational study of the adsorptive desulfurization of small aromatic sulfur compounds by conjugated microporous polymers (CMPs). The density-functional tight-binding method augmented with an R-6 dispersion correction is employed to investigate the physisorption binding mechanism and electronic properties of the CMP-aromatic sulfur complexes. We show that the widely extended π conjugation in the CMP skeletons is favorable for the non-covalent adsorption of aromatic thiophene and dibenzothiophene via π-π, H-π, and S-π interactions. The average binding energies are calculated to be -6.2 âˆ¼ -15.2 kcal/mol for CMP- thiophene/dibenzothiophene systems. For the dibenzothiophene molecule with larger size and more extended conjugation, it binds more than twice stronger to CMP than the thiophene molecule. We show that the replacement of quinoline unit to the phenylene group in the network linker effectively enhances the average binding capacities by around 0.8-1.8 kcal/mol. Our calculations theoretically demonstrate that CMPs materials are kind of promising candidates for the adsorptive desulfurization of small aromatic sulfur compounds. This paper provides useful theoretical guidance for design of novel carbon-based adsorbents for adsorptive desulfurization.

Development of a Force Field for the Simulation of Single-Chain Proteins and Protein-Protein Complexes.

The accuracy of atomistic physics-based force fields for the simulation of biological macromolecules has typically been benchmarked experimentally using biophysical data from simple, often single-chain systems. In the case of proteins, the careful refinement of force field parameters associated with torsion-angle potentials and the use of improved water models have enabled a great deal of progress toward the highly accurate simulation of such monomeric systems in both folded and, more recently, disordered states. In living organisms, however, proteins constantly interact with other macromolecules, such as proteins and nucleic acids, and these interactions are often essential for proper biological function. Here, we show that state-of-the-art force fields tuned to provide an accurate description of both ordered and disordered proteins can be limited in their ability to accurately describe protein-protein complexes. This observation prompted us to perform an extensive reparameterization of one variant of the Amber protein force field. Our objective involved refitting not only the parameters associated with torsion-angle potentials but also the parameters used to model nonbonded interactions, the specification of which is expected to be central to the accurate description of multicomponent systems. The resulting force field, which we call DES-Amber, allows for more accurate simulations of protein-protein complexes, while still providing a state-of-the-art description of both ordered and disordered single-chain proteins. Despite the improvements, calculated protein-protein association free energies still appear to deviate substantially from experiment, a result suggesting that more fundamental changes to the force field, such as the explicit treatment of polarization effects, may simultaneously further improve the modeling of single-chain proteins and protein-protein complexes.

The N-terminal domains of FLASH and Lsm11 form a 2:1 heterotrimer for histone pre-mRNA 3'-end processing.

Unlike canonical pre-mRNAs, animal replication-dependent histone pre-mRNAs lack introns and are processed at the 3'-end by a mechanism distinct from cleavage and polyadenylation. They have a 3' stem loop and histone downstream element (HDE) that are recognized by stem-loop binding protein (SLBP) and U7 snRNP, respectively. The N-terminal domain (NTD) of Lsm11, a component of U7 snRNP, interacts with FLASH NTD and these two proteins recruit the histone cleavage complex containing the CPSF-73 endonuclease for the cleavage reaction. Here, we determined crystal structures of FLASH NTD and found that it forms a coiled-coil dimer. Using solution light scattering, we characterized the stoichiometry of the FLASH NTD-Lsm11 NTD complex and found that it is a 2:1 heterotrimer, which is supported by observations from analytical ultracentrifugation and crosslinking.

Computational study of diketopyrrolopyrrole-based organic dyes for dye sensitized solar cell applications.

Four diketopyrrolopyrrole (DPP)-based organic dyes utilizing the donor-π-acceptor motif were investigated by density functional theory (DFT) and time-dependent DFT (TDDFT) approaches. The four dyes were composed of different donor groups, i.e. indoline, carbazole, triphenylamine, and coumarin. We investigated the effects of the DPP unit and different donors on the spectra and electrochemical properties of the dyes, respectively. In comparison with the model dye which adopts a phenylene unit as the π-spacer, the DPP dyes all display remarkably enhanced spectral responses in the visible region of the solar spectrum. The key to this increase was the incorporation of electron-deficient DPP moieties to the molecular core, which significantly lowers LUMO levels and therefore reduces the band gap. The dye/(TiO2)46 anatase nanoparticle systems were also simulated to show the electronic structures at the interface. We studied some key properties including absorption spectra, light-harvesting efficiency, molecular orbital distributions, and injection time of electrons from the excited state of dye to the conduction band of TiO2. The dye DPP-I with indoline moiety as the electron donor demonstrates desirable energetic, electronic, and spectroscopic parameters for dye sensitized solar cells (DSSCs) applications. Our theoretical study is expected to provide valuable insights into the molecular design of novel DPP-based organic dyes for the optimizations of DSSCs.

The ROQ domain of Roquin recognizes mRNA constitutive-decay element and double-stranded RNA.

A conserved stem-loop motif of the constitutive decay element (CDE) in the 3' UTR of mRNAs is recognized by the ROQ domain of Roquin, which mediates mRNA degradation. Here we report two crystal structures of the Homo sapiens ROQ domain in complex with CDE RNA. The ROQ domain has an elongated shape with three subdomains. The 19-nt Hmgxb3 CDE is bound as a stem-loop to domain III. The 23-nt TNF RNA is bound as a duplex to a separate site at the interface between domains I and II. Mutagenesis studies confirm that the ROQ domain has two separate RNA-binding sites, one for stem-loop RNA (A site) and the other for double-stranded RNA (B site). Mutation in either site perturbs the Roquin-mediated degradation of HMGXB3 and IL6 mRNAs in human cells, demonstrating the importance of both sites for mRNA decay.

Structure of histone mRNA stem-loop, human stem-loop binding protein, and 3'hExo ternary complex.

Metazoan replication-dependent histone messenger RNAs (mRNAs) have a conserved stem-loop (SL) at their 3'-end. The stem-loop binding protein (SLBP) specifically recognizes the SL to regulate histone mRNA metabolism, and the 3'-5' exonuclease 3'hExo trims its 3'-end after processing. We report the crystal structure of a ternary complex of human SLBP RNA binding domain, human 3'hExo, and a 26-nucleotide SL RNA. Only one base of the SL is recognized specifically by SLBP, and the two proteins primarily recognize the shape of the RNA. SLBP and 3'hExo have no direct contact with each other, and induced structural changes in the loop of the SL mediate their cooperative binding. The 3' flanking sequence is positioned in the 3'hExo active site, but the ternary complex limits the extent of trimming.

Crystal structure of DmdD, a crotonase superfamily enzyme that catalyzes the hydration and hydrolysis of methylthioacryloyl-CoA.

Dimethyl-sulphoniopropionate (DMSP) is produced in abundance by marine phytoplankton, and the catabolism of this compound is an important source of carbon and reduced sulfur for marine bacteria and other organisms. The enzyme DmdD catalyzes the last step in the methanethiol (MeSH) pathway of DMSP catabolism. DmdD is a member of the crotonase superfamily of enzymes, and it catalyzes both the hydration and the hydrolysis of methylthioacryloyl-CoA (MTA-CoA), converting it to acetaldehyde, CO2, MeSH, and CoA. We report here the crystal structure of Ruegeria pomeroyi DmdD free enzyme at 1.5 Å resolution and the structures of the E121A mutant in complex with MTA-CoA and 3-methylmercaptopropionate-CoA (MMPA-CoA) at 1.8 Å resolution. DmdD is a hexamer, composed of a dimer of trimers where the three monomers of each trimer are related by a crystallographic 3-fold axis. The overall structure of this hexamer is similar to those of canonical crotonases. However, the C-terminal loops of DmdD in one of the trimers assume a different conformation and contribute to CoA binding in the active site of a neighboring monomer of the trimer, while these loops in the second trimer are disordered. MTA-CoA is bound deep in the active site in the first trimer, but shows a 1.5 Å shift in its position in the second trimer. MMPA-CoA has a similar binding mode to MTA-CoA in the first trimer. MMPA-CoA cannot be hydrated and is only hydrolyzed slowly by DmdD. Replacement of the sulfur atom in MMPA-CoA with a methylene group abolishes hydrolysis, suggesting that the unique property of the substrate is a major determinant of the hydrolysis activity of DmdD.