Romina Oliva

SambVca 2. A Web Tool for Analyzing Catalytic Pockets with Topographic Steric Maps

Developing more efficient catalysts remains one of the primary targets of organometallic chemists. To accelerate reaching this goal, effective molecular descriptors and visualization tools can represent a remarkable aid. Here, we present a Web application for analyzing the catalytic pocket of metal complexes using topographic steric maps as a general and unbiased descriptor that is suitable for every class of catalysts. To show the broad applicability of our approach, we first compared the steric map of a series of transition metal complexes presenting popular mono-, di-, and tetracoordinated ligands and three classic zirconocenes. This comparative analysis highlighted similarities and differences between totally unrelated ligands. Then, we focused on a recently developed Fe(II) catalyst that is active in the asymmetric transfer hydrogenation of ketones and imines. Finally, we expand the scope of these tools to rationalize the inversion of enantioselectivity in enzymatic catalysis, achieved by point mutat...

Structural and kinetic characterization of native laccases from Pleurotus ostreatus, Rigidoporus lignosus, and Trametes trogii.

A comparative study has been performed on five native laccases purified from the three basidiomycete fungi Pleurotus ostreatus, Rigidoporus lignosus, and Trametes trogii to relate their different catalytic capacities to their structural properties. Spectroscopic absorption features and EPR spectra at various pH values of the five enzymes are very similar and typical of the blue oxidases. The analysis of the dependence of kinetic parameters on pH suggested that a histidine residue is involved in the binding of nonphenolic substrates, whereas both a histidine and an acidic residue may be involved in the binding of phenolic compounds. His and an Asp residue are indeed found at the bottom of a cavity which may be regarded as a suitable substrate channel for approaching to type 1 copper in the 3D homology models of the two laccases from Pleuorotus ostreatus (POXC and POXA1b) whose sequences are known.

COCOMAPS: a web application to analyze and visualize contacts at the interface of biomolecular complexes.

SUMMARY Herein we present COCOMAPS, a novel tool for analyzing, visualizing and comparing the interface in protein-protein and protein-nucleic acids complexes. COCOMAPS combines traditional analyses and 3D visualization of the interface with the effectiveness of intermolecular contact maps. AVAILABILITY COCOMAPS is accessible as a public web tool at http://www.molnac.unisa.it/BioTools/cocomaps CONTACT romina.oliva@uniparthenope.it; lcavallo@unisa.it.

Accurate energies of hydrogen bonded nucleic acid base pairs and triplets in tRNA tertiary interactions.

Tertiary interactions are crucial in maintaining the tRNA structure and functionality. We used a combined sequence analysis and quantum mechanics approach to calculate accurate energies of the most frequent tRNA tertiary base pairing interactions. Our analysis indicates that six out of the nine classical tertiary interactions are held in place mainly by H-bonds between the bases. In the remaining three cases other effects have to be considered. Tertiary base pairing interaction energies range from −8 to −38 kcal/mol in yeast tRNAPhe and are estimated to contribute roughly 25% of the overall tRNA base pairing interaction energy. Six analyzed posttranslational chemical modifications were shown to have minor effect on the geometry of the tertiary interactions. Modifications that introduce a positive charge strongly stabilize the corresponding tertiary interactions. Non-additive effects contribute to the stability of base triplets.

Frequency and effect of the binding of Mg2+, Mn2+, and Co2+ ions on the guanine base in Watson-Crick and reverse Watson-Crick base pairs.

We performed MP2 calculations to elucidate the structure and energetics of the Mg(2+), Mn(2+), and Co(2+) hexahydrated aquaions, and the effect of the metal binding to the N7 atom of (i) a single guanine, (ii) a guanine involved in a Watson-Crick pair, and (iii) a guanine involved in a reverse Watson-Crick base pair. Our comparative analysis of the three aquaions indicates a clear inverse correlation between the radius of the cation and the binding energy, that indeed increases in the order Mn(2+) < Co(2+) < Mg(2+). The trend in the binding energies of the pentahydrated cations to the N7 atom of the guanine is instead Mg(2+) < Mn(2+) < Co(2+), suggesting a rather different bonding scheme that, for the two transition metals, involves back-donation from the aromatic ring of the guanine to their empty d orbitals. In the gas phase, the three hydrated metals significantly stabilize both G-C base pair geometries, Watson-Crick and reverse Watson-Crick, we investigated. Inclusion of a continuous solvent model, however, remarkably reduces this additional stabilization, which becomes almost negligible in the case of the Mg(2+) cation coordinated to the guanine in the standard Watson-Crick geometry. Conversely, all three metal ions sensibly stabilize the reverse Watson-Crick geometry, also in water. Our results are supported by a screening of the structures available in the Protein Data Bank, which clearly indicates that the two transition metals we investigated have a tendency greater than Mg(2+) to coordinate to the N7 atom of guanines, and that there is no clear correlation between the number of guanines in experimental structures with a metal bound to N7 atom and their involvement in Watson-Crick base pairs.

Electrostatics of aquaporin and aquaglyceroporin channels correlates with their transport selectivity.

Aquaporins are homotetrameric channel proteins, which allow the diffusion of water and small solutes across biological membranes. According to their transport function, aquaporins can be divided into “orthodox aquaporins”, which allow the flux of water molecules only, and “aquaglyceroporins”, which facilitate the diffusion of glycerol and other small solutes in addition to water. The contribution of individual residues in the pore to the selectivity of orthodox aquaporins and aquaglyceroporins is not yet fully understood. To gain insights into aquaporin selectivity, we focused on the sequence variation and electrostatics of their channels. The continuum Poisson-Boltzmann electrostatic potential along the channel was calculated and compared for ten three-dimensional-structures which are representatives of different aquaporin subfamilies, and a panel of functionally characterized mutants, for which high-accuracy three-dimensional-models could be derived. Interestingly, specific electrostatic profiles associated with the main selectivity to water or glycerol could be identified. In particular: (i) orthodox aquaporins showed a distinctive electrostatic potential maximum at the periplasmic side of the channel around the aromatic/Arg (ar/R) constriction site; (ii) aquaporin-0 (AQP0), a mammalian aquaporin with considerably low water permeability, had an additional deep minimum at the cytoplasmic side; (iii) aquaglyceroporins showed a rather flat potential all along the channel; and (iv) the bifunctional protozoan PfAQP had an unusual all negative profile. Evaluation of electrostatics of the mutants, along with a thorough sequence analysis of the aquaporin pore-lining residues, illuminated the contribution of specific residues to the electrostatics of the channels and possibly to their selectivity.

Can yeast be used to study mitochondrial diseases? Biolistic tRNA mutants for the analysis of mechanisms and suppressors.

Base substitutions equivalent to those causing human pathologies have been introduced in yeast mitochondrial tRNA genes. These mutants can be utilized as flexible tools to investigate the molecular aspects of mitochondrial diseases and identify correcting genes. We show that for all studied tRNA mutations (including an homoplasmic one in tRNA(Val)) the severity of phenotypes follows the same trend in four different nuclear backgrounds. Correcting genes include TUF1 and genes encoding aminoacyl-tRNA synthetase. The effect of suppressors was analyzed by Northern blot. Mutated leucyl-tRNA synthetase with highly reduced catalytic activity maintains full suppressing effect, thus suggesting a chaperone-like and/or stabilizing function.

Higher order structural effects stabilizing the reverse Watson–Crick Guanine-Cytosine base pair in functional RNAs

The G:C reverse Watson–Crick (W:W trans) base pair, also known as Levitt base pair in the context of tRNAs, is a structurally and functionally important base pair that contributes to tertiary interactions joining distant domains in functional RNA molecules and also participates in metabolite binding in riboswitches. We previously indicated that the isolated G:C W:W trans base pair is a rather unstable geometry, and that dicationic metal binding to the Guanine base or posttranscriptional modification of the Guanine can increase its stability. Herein, we extend our survey and report on other H-bonding interactions that can increase the stability of this base pair. To this aim, we performed a bioinformatics search of the PDB to locate all the occurencies of G:C trans base pairs. Interestingly, 66% of the G:C trans base pairs in the PDB are engaged in additional H-bonding interactions with other bases, the RNA backbone or structured water molecules. High level quantum mechanical calculations on a data set of representative crystal structures were performed to shed light on the structural stability and energetics of the various crystallographic motifs. This analysis was extended to the binding of the preQ1 metabolite to a preQ1-II riboswitch.

Ranking multiple docking solutions based on the conservation of inter‐residue contacts

Molecular docking is the method of choice for investigating the molecular basis of recognition in a large number of functional protein complexes. However, correctly scoring the obtained docking solutions (decoys) to rank native‐like (NL) conformations in the top positions is still an open problem. Herein we present CONSRANK, a simple and effective tool to rank multiple docking solutions, which relies on the conservation of inter‐residue contacts in the analyzed decoys ensemble. First it calculates a conservation rate for each inter‐residue contact, then it ranks decoys according to their ability to match the more frequently observed contacts. We applied CONSRANK to 102 targets from three different benchmarks, RosettaDock, DOCKGROUND, and Critical Assessment of PRedicted Interactions (CAPRI). The method performs consistently well, both in terms of NL solutions ranked in the top positions and of values of the area under the receiver operating characteristic curve. Its ideal application is to solutions coming from different docking programs and procedures, as in the case of CAPRI targets. For all the analyzed CAPRI targets where a comparison is feasible, CONSRANK outperforms the CAPRI scorers. The fraction of NL solutions in the top ten positions in the RosettaDock, DOCKGROUND, and CAPRI benchmarks is enriched on average by a factor of 3.0, 1.9, and 9.9, respectively. Interestingly, CONSRANK is also able to specifically single out the high/medium quality (HMQ) solutions from the docking decoys ensemble: it ranks 46.2 and 70.8% of the total HMQ solutions available for the RosettaDock and CAPRI targets, respectively, within the top 20 positions. Proteins 2013.

CONS-COCOMAPS: a novel tool to measure and visualize the conservation of inter-residue contacts in multiple docking solutions

BackgroundThe development of accurate protein-protein docking programs is making this kind of simulations an effective tool to predict the 3D structure and the surface of interaction between the molecular partners in macromolecular complexes. However, correctly scoring multiple docking solutions is still an open problem. As a consequence, the accurate and tedious screening of many docking models is usually required in the analysis step.MethodsAll the programs under CONS-COCOMAPS have been written in python, taking advantage of python libraries such as SciPy and Matplotlib. CONS-COCOMAPS is freely available as a web tool at the URL:http://www.molnac.unisa.it/BioTools/conscocomaps/.ResultsHere we presented CONS-COCOMAPS, a novel tool to easily measure and visualize the consensus in multiple docking solutions. CONS-COCOMAPS uses the conservation of inter-residue contacts as an estimate of the similarity between different docking solutions. To visualize the conservation, CONS-COCOMAPS uses intermolecular contact maps.ConclusionsThe application of CONS-COCOMAPS to test-cases taken from recent CAPRI rounds has shown that it is very efficient in highlighting even a very weak consensus that often is biologically meaningful.