RESUMEN

Interoperability in computational chemistry is elusive, impeded by the independent development of software packages and idiosyncratic nature of their output files. The cclib library was introduced in 2006 as an attempt to improve this situation by providing a consistent interface to the results of various quantum chemistry programs. The shared API across programs enabled by cclib has allowed users to focus on results as opposed to output and to combine data from multiple programs or develop generic downstream tools. Initial development, however, did not anticipate the rapid progress of computational capabilities, novel methods, and new programs; nor did it foresee the growing need for customizability. Here, we recount this history and present cclib 2, focused on extensibility and modularity. We also introduce recent design pivots-the formalization of cclib's intermediate data representation as a tree-based structure, a new combinator-based parser organization, and parsed chemical properties as extensible objects.

RESUMEN

The bond-valence model is a reliable way to validate assumed oxidation states based on structural data. It has successfully been employed for analyzing metal-binding sites in macromolecule structures. However, inconsistent results for heme-based structures suggest that some widely used bond-valence R0 parameters may need to be adjusted in certain cases. Given the large number of experimental crystal structures gathered since these initial parameters were determined and the similarity of binding sites in organic compounds and macromolecules, the Cambridge Structural Database (CSD) is a valuable resource for refining metal-organic bond-valence parameters. R0 bond-valence parameters for iron(II), iron(III) and other metals have been optimized based on an automated processing of all CSD crystal structures. Almost all R0 bond-valence parameters were reproduced, except for iron-nitrogen bonds, for which distinct R0 parameters were defined for two observed subpopulations, corresponding to low-spin and high-spin states, of iron in both oxidation states. The significance of this data-driven method for parameter discovery, and how the spin state affects the interpretation of heme-containing proteins and iron-binding sites in macromolecular structures, are discussed.

Asunto(s)

RESUMEN

The low reproducibility of published experimental results in many scientific disciplines has recently garnered negative attention in scientific journals and the general media. Public transparency, including the availability of `raw' experimental data, will help to address growing concerns regarding scientific integrity. Macromolecular X-ray crystallography has led the way in requiring the public dissemination of atomic coordinates and a wealth of experimental data, making the field one of the most reproducible in the biological sciences. However, there remains no mandate for public disclosure of the original diffraction data. The Integrated Resource for Reproducibility in Macromolecular Crystallography (IRRMC) has been developed to archive raw data from diffraction experiments and, equally importantly, to provide related metadata. Currently, the database of our resource contains data from 2920 macromolecular diffraction experiments (5767 data sets), accounting for around 3% of all depositions in the Protein Data Bank (PDB), with their corresponding partially curated metadata. IRRMC utilizes distributed storage implemented using a federated architecture of many independent storage servers, which provides both scalability and sustainability. The resource, which is accessible via the web portal at http://www.proteindiffraction.org, can be searched using various criteria. All data are available for unrestricted access and download. The resource serves as a proof of concept and demonstrates the feasibility of archiving raw diffraction data and associated metadata from X-ray crystallographic studies of biological macromolecules. The goal is to expand this resource and include data sets that failed to yield X-ray structures in order to facilitate collaborative efforts that will improve protein structure-determination methods and to ensure the availability of `orphan' data left behind for various reasons by individual investigators and/or extinct structural genomics projects.

Asunto(s)

RESUMEN

Bacteria are increasingly resistant to existing antibiotics, which target a narrow range of pathways. New methods are needed to identify targets, including repositioning targets among distantly related species. We developed a novel combination of systems and structural modeling and bioinformatics to reposition known antibiotics and targets to new species. We applied this approach to Mycoplasma genitalium, a common cause of urethritis. First, we used quantitative metabolic modeling to identify enzymes whose expression affects the cellular growth rate. Second, we searched the literature for inhibitors of homologs of the most fragile enzymes. Next, we used sequence alignment to assess that the binding site is shared by M. genitalium, but not by humans. Lastly, we used molecular docking to verify that the reported inhibitors preferentially interact with M. genitalium proteins over their human homologs. Thymidylate kinase was the top predicted target and piperidinylthymines were the top compounds. Further work is needed to experimentally validate piperidinylthymines. In summary, combined systems and structural modeling is a powerful tool for drug repositioning.

Asunto(s)

RESUMEN

The relative stability of biologically relevant, hydrogen bonded complexes with shortened distances can be assessed at low cost by the electrostatic multipole term alone more successfully than by ab initio methods. These results imply that atomic multipole moments may help improve ligand-receptor ranking predictions, particularly in cases where accurate structural data are not available.

Asunto(s)

RESUMEN

Fatty acid amide hydrolase (FAAH) is a member of the amidase signature family and is responsible for the hydrolytic deactivation of fatty acid amide neuromodulators, such as anandamide. FAAH carries an unusual catalytic triad consisting of Lys-Ser-Ser, which uniquely enables the enzyme to cleave amides and esters at similar rates. The acylation of 9Z-octadecenamide (oleamide, a FAAH reference substrate) has been widely investigated by computational methods, and those have shown that conformational fluctuations of the active site affect the reaction barrier. Empirical descriptors have been devised to provide a possible mechanistic explanation for such conformational effects, but a first-principles understanding is still missing. A comparison of FAAH acylation with a reference reaction in water suggests that transition-state stabilization is crucial for catalysis because the activation energy barrier falls by 6 kcal/mol in the presence of the active site. With this in mind, we have analyzed the enzymatic reaction using the differential transition-state stabilization (DTSS) approach to determine key active-site residues for lowering the barrier. We examined several QM/MM structures at the MP2 level of theory and analyzed catalytic effects with a variation-perturbation partitioning of the interaction energy into electrostatic multipole and penetration, exchange, delocalization, and correlation terms. Three residues - Thr236, Ser218, and one water molecule - appear to be essential for the stabilization of the transition state, a conclusion that is also reflected by catalytic fields and agrees with site-directed mutagenesis data. An analogous analysis for URB524, URB618, and URB694 (three potent representatives of covalent, carbamate-based FAAH inhibitors) confirms the importance of the residues involved in oleamide acylation, providing insight for future inhibitor design.

Asunto(s)

RESUMEN

External electric fields align nanostructured block copolymers by either rotation of grains or nucleation and growth depending on how strongly the chemically distinct block copolymer components are segregated. In close vicinity to the order-disorder transition, theory and simulations suggest a third mechanism: selective disordering. We present a time-resolved small-angle X-ray scattering study that demonstrates how an electric field can indeed selectively disintegrate ill-aligned lamellae in a lyotropic block copolymer solution, while lamellae with interfaces oriented parallel to the applied field prevail. The present study adds an additional mechanism to the experimentally corroborated suite of mechanistic pathways, by which nanostructured block copolymers can align with an electric field. Our results further unveil the benefit of electric field assisted annealing for mitigating orientational disorder and topological defects in block copolymer mesophases, both in close vicinity to the order-disorder transition and well below it.

Asunto(s)

RESUMEN

BACKGROUND: The Blue Obelisk movement was established in 2005 as a response to the lack of Open Data, Open Standards and Open Source (ODOSOS) in chemistry. It aims to make it easier to carry out chemistry research by promoting interoperability between chemistry software, encouraging cooperation between Open Source developers, and developing community resources and Open Standards. RESULTS: This contribution looks back on the work carried out by the Blue Obelisk in the past 5 years and surveys progress and remaining challenges in the areas of Open Data, Open Standards, and Open Source in chemistry. CONCLUSIONS: We show that the Blue Obelisk has been very successful in bringing together researchers and developers with common interests in ODOSOS, leading to development of many useful resources freely available to the chemistry community.

RESUMEN

We report a systematic analysis of the intermolecular interactions of cationic ethidium intercalated into a UA/AU step of RNA for a single conformation based on crystallographic coordinates. Interaction energies at the MP2/6-31G** level were partitioned into electrostatic, exchange, delocalization, and correlation components. Various pairwise interaction models built from chemically intuitive fragments reproduce within a few percent values obtained when treating the intercalation site as a whole. Gas phase results are very sensitive to the charge state of the two phosphate groups, with the electrostatic term nearly tripling when the counterions are removed. But this is largely compensated by solvation, an effect represented here within the polarizable continuum model. In a few cases, more diffuse and larger basis sets as well as QCISD(T) corrections were applied in an effort to estimate plausible ethidium-nucleobase electron correlation effects.

RESUMEN

There are now a wide variety of packages for electronic structure calculations, each of which differs in the algorithms implemented and the output format. Many computational chemistry algorithms are only available to users of a particular package despite being generally applicable to the results of calculations by any package. Here we present cclib, a platform for the development of package-independent computational chemistry algorithms. Files from several versions of multiple electronic structure packages are automatically detected, parsed, and the extracted information converted to a standard internal representation. A number of population analysis algorithms have been implemented as a proof of principle. In addition, cclib is currently used as an input filter for two GUI applications that analyze output files: PyMOlyze and GaussSum.

Asunto(s)

RESUMEN

Major components of the interaction energy that define several approximate levels starting from second order Möller-Plesset theory were studied for 58 stacked nucleic acid dimers. They included typical B-DNA and A-DNA structures, and selected published geometries. A survey of the various terms yields an unexpected correlation between the Pauli exchange and dispersion or correlation terms, which holds for each class of similar planar geometries and for various basis sets. The geometries that exhibit these correlations span a specific range of molecular overlaps when compared to a model benzene-pyridine stacked dimer. Also, the relationship between electrostatic interactions and MP2 stabilization energies reported earlier is confirmed and a prediction interval of practical relevance is estimated.

Asunto(s)

RESUMEN

On the basis of the crystallographic structures of three nucleic acid intercalation complexes involving ethidium and proflavine, we have analyzed the interaction energies between intercalator chromophores and their four nearest bases, using a hybrid variation-perturbation method at the second-order Møller-Plesset theory level (MP2) with a 6-31G(d,p) basis set. A total MP2 interaction energy minimum precisely reproduces the crystallographic position of the ethidium chromophore in the intercalation plane between UA/AU bases. The electrostatic component constitutes the same fraction of the total energy for all three studied structures. The multipole electrostatic interaction energy, calculated from cumulative atomic multipole moments (CAMMs), was found to converge only after including components above the fifth order. CAMM interaction surfaces, calculated on grids in the intercalation planes of these structures, reasonably reproduce the alignment of intercalators in crystal structures; they exhibit additional minima in the direction of the DNA grooves, however, which also need to be examined at higher theory levels if no crystallographic data are given.

Asunto(s)