Mauricio Carrillo-Tripp

VIPERdb2: an enhanced and web API enabled relational database for structural virology.

VIPERdb (http://viperdb.scripps.edu) is a relational database and a web portal for icosahedral virus capsid structures. Our aim is to provide a comprehensive resource specific to the needs of the virology community, with an emphasis on the description and comparison of derived data from structural and computational analyses of the virus capsids. In the current release, VIPERdb2, we implemented a useful and novel method to represent capsid protein residues in the icosahedral asymmetric unit (IAU) using azimuthal polar orthographic projections, otherwise known as Φ–Ψ (Phi–Psi) diagrams. In conjunction with a new Application Programming Interface (API), these diagrams can be used as a dynamic interface to the database to map residues (categorized as surface, interface and core residues) and identify family wide conserved residues including hotspots at the interfaces. Additionally, we enhanced the interactivity with the database by interfacing with web-based tools. In particular, the applications Jmol and STRAP were implemented to visualize and interact with the virus molecular structures and provide sequence–structure alignment capabilities. Together with extended curation practices that maintain data uniformity, a relational database implementation based on a schema for macromolecular structures and the APIs provided will greatly enhance the ability to do structural bioinformatics analysis of virus capsids.

A comparative study of the hydration of Na+ and K+ with refined polarizable model potentials

The hydrations of Na+ and K+ were investigated by means of Monte Carlo simulations with refined ab initio based potentials. These interaction potentials include intramolecular relaxation, polarizability and many-body nonadditive effects. Care was taken to ensure proper convergence of the MC runs and that the statistical samples were large enough. As a result, agreement was attained with all experimental data available for the hydration of the ions. The water molecules in the first hydration shell were found to have the same intramolecular geometries and dipole moments as those of the bulk. Furthermore, their dipoles were not aligned to the electric field produced by the ion, but quite tilted. The hydration number for the sodium was found to be 5 or 6 water molecules, whereas the potassium’s hydration number had a probability distribution ranging from 5 to 10. From an analysis of the energetic contributions of each hydration shell to the total enthalpy of hydration we propose that the hydrated ions have a ...

Evolution of Substrate Specificity in a Recipient’s Enzyme Following Horizontal Gene Transfer

Despite the prominent role of horizontal gene transfer (HGT) in shaping bacterial metabolism, little is known about the impact of HGT on the evolution of enzyme function. Specifically, what is the influence of a recently acquired gene on the function of an existing gene? For example, certain members of the genus Corynebacterium have horizontally acquired a whole l-tryptophan biosynthetic operon, whereas in certain closely related actinobacteria, for example, Mycobacterium, the trpF gene is missing. In Mycobacterium, the function of the trpF gene is performed by a dual-substrate (βα)8 phosphoribosyl isomerase (priA gene) also involved in l-histidine (hisA gene) biosynthesis. We investigated the effect of a HGT-acquired TrpF enzyme upon PriAs substrate specificity in Corynebacterium through comparative genomics and phylogenetic reconstructions. After comprehensive in vivo and enzyme kinetic analyses of selected PriA homologs, a novel (βα)8 isomerase subfamily with a specialized function in l-histidine biosynthesis, termed subHisA, was confirmed. X-ray crystallography was used to reveal active-site mutations in subHisA important for narrowing of substrate specificity, which when mutated to the naturally occurring amino acid in PriA led to gain of function. Moreover, in silico molecular dynamic analyses demonstrated that the narrowing of substrate specificity of subHisA is concomitant with loss of ancestral protein conformational states. Our results show the importance of HGT in shaping enzyme evolution and metabolism.

An Amphotericin B Derivative Equally Potent to Amphotericin B and with Increased Safety

Amphotericin B is the most potent antimycotic known to date. However due to its large collateral toxicity, its use, although long standing, had been limited. Many attempts have been made to produce derivatives with reduced collateral damage. The molecular mechanism of polyene has also been closely studied for this purpose and understanding it would contribute to the development of safe derivatives. Our study examined polyene action, including chemical synthesis, electrophysiology, pharmacology, toxicology and molecular dynamics. The results were used to support a novel Amphotericin B derivative with increased selectivity: L-histidine methyl ester of Amphotericin B. We found that this derivative has the same form of action as Amphotericin B, i.e. pore formation in the cell membrane. Its reduced dimerization in solution, when compared to Amphotericin B, is at least partially responsible for its increased selectivity. Here we also present the results of preclinical tests, which show that the derivative is just as potent as Amphotericin B and has increased safety.

Insights into the evolution of enzyme substrate promiscuity after the discovery of (βα)8 isomerase evolutionary intermediates from a diverse metagenome

BackgroundCurrent sequence-based approaches to identify enzyme functional shifts, such as enzyme promiscuity, have proven to be highly dependent on a priori functional knowledge, hampering our ability to reconstruct evolutionary history behind these mechanisms. Hidden Markov Model (HMM) profiles, broadly used to classify enzyme families, can be useful to distinguish between closely related enzyme families with different specificities. The (βα)8-isomerase HisA/PriA enzyme family, involved in L-histidine (HisA, mono-substrate) biosynthesis in most bacteria and plants, but also in L-tryptophan (HisA/TrpF or PriA, dual-substrate) biosynthesis in most Actinobacteria, has been used as model system to explore evolutionary hypotheses and therefore has a considerable amount of evolutionary, functional and structural knowledge available. We searched for functional evolutionary intermediates between the HisA and PriA enzyme families in order to understand the functional divergence between these families.ResultsWe constructed a HMM profile that correctly classifies sequences of unknown function into the HisA and PriA enzyme sub-families. Using this HMM profile, we mined a large metagenome to identify plausible evolutionary intermediate sequences between HisA and PriA. These sequences were used to perform phylogenetic reconstructions and to identify functionally conserved amino acids. Biochemical characterization of one selected enzyme (CAM1) with a mutation within the functionally essential N-terminus phosphate-binding site, namely, an alanine instead of a glycine in HisA or a serine in PriA, showed that this evolutionary intermediate has dual-substrate specificity. Moreover, site-directed mutagenesis of this alanine residue, either backwards into a glycine or forward into a serine, revealed the robustness of this enzyme. None of these mutations, presumably upon functionally essential amino acids, significantly abolished its enzyme activities. A truncated version of this enzyme (CAM2) predicted to adopt a (βα)6-fold, and thus entirely lacking a C-terminus phosphate-binding site, was identified and shown to have HisA activity.ConclusionAs expected, reconstruction of the evolution of PriA from HisA with HMM profiles suggest that functional shifts involve mutations in evolutionarily intermediate enzymes of otherwise functionally essential residues or motifs. These results are in agreement with a link between promiscuous enzymes and intragenic epistasis. HMM provides a convenient approach for gaining insights into these evolutionary processes.

A novel method to map and compare protein–protein interactions in spherical viral capsids

Viral capsids are composed of multiple copies of one or a few chemically distinct capsid proteins and are mostly stabilized by inter subunit protein–protein interactions. There have been efforts to identify and analyze these protein–protein interactions, in terms of their extent and similarity, between the subunit interfaces related by quasi‐ and icosahedral symmetry. Here, we describe a new method to map quaternary interactions in spherical virus capsids onto polar angle space with respect to the icosahedral symmetry axes using azimuthal orthographic diagrams. This approach enables one to map the nonredundant interactions in a spherical virus capsid, irrespective of its size or triangulation number (T), onto the reference icosahedral asymmetric unit space. The resultant diagrams represent characteristic fingerprints of quaternary interactions of the respective capsids. Hence, they can be used as road maps of the protein–protein interactions to visualize the distribution and the density of the interactions. In addition, unlike the previous studies, the fingerprints of different capsids, when represented in a matrix form, can be compared with one another to quantitatively evaluate the similarity (S‐score) in the subunit environments and the associated protein–protein interactions. The S‐score selectively distinguishes the similarity, or lack of it, in the locations of the quaternary interactions as opposed to other well‐known structural similarity metrics (e.g., RMSD, TM‐score). Application of this method on a subset of T = 1 and T = 3 capsids suggests that S‐score values range between 1 and 0.6 for capsids that belong to the same virus family/genus; 0.6–0.3 for capsids from different families with the same T‐number and similar subunit fold; and <0.3 for comparisons of the dissimilar capsids that display different quaternary architectures (T‐numbers). Finally, the sequence conserved interface residues within a virus family, whose spatial locations were also conserved have been hypothesized as the essential residues for self‐assembly of the member virus capsids. Proteins 2008.

Antagonism or synergism between papaya ringspot virus and papaya mosaic virus in Carica papaya is determined by their order of infection.

Antagonism between unrelated plant viruses has not been thoroughly described. Our studies show that two unrelated viruses, papaya ringspot virus (PRSV) and papaya mosaic virus (PapMV) produce different symptomatic outcomes during mixed infection depending on the inoculation order. Synergism occurs in plants infected first with PRSV or in plants infected simultaneously with PRSV and PapMV, and antagonism occurs in plants infected first with PapMV and later inoculated with PRSV. During antagonism, elevated pathogenesis-related (PR-1) gene expression and increased reactive oxygen species production indicated the establishment of a host defense resulting in the reduction in PRSV titers. Polyribosomal fractioning showed that PRSV affects translation of cellular eEF1α, PR-1, β-tubulin, and PapMV RNAs in planta, suggesting that its infection could be related to an imbalance in the translation machinery. Our data suggest that primary PapMV infection activates a defense response against PRSV and establishes a protective relationship with the papaya host.

Immune response to a potyvirus with exposed amino groups available for chemical conjugation

BackgroundThe amino terminus of the tobacco etch virus (TEV) capsid protein is located on the external surface of infectious TEV particles, as proposed by previous studies and an in silico model. The epsilon amino groups on the exposed lysine residues are available for chemical conjugation to any given protein, and can thus act as antigen carriers. The availability of amino groups on the surfaces of TEV particles was determined and the immune response to TEV evaluated.ResultsUsing a biotin-tagged molecule that reacts specifically with amino groups, we found that the TEV capsid protein has amino groups on its surface available for coupling to other molecules via crosslinkers. Intraperitoneal TEV was administered to female BALB/c mice, and both their humoral and cellular responses measured. Different IgG isotypes, particularly IgG2a, directed against TEV were induced. In a cell proliferation assay, only spleen cells from vaccinated mice that were stimulated in vitro with TEV showed significant proliferation of CD3+/CD4+ and CD3+/CD8+ subpopulations and secreted significant amounts of interferon γ.ConclusionsTEV has surface amino groups that are available for chemical coupling. TEV induces both humoral and cellular responses when administered alone intraperitoneally to mice. Therefore, TEV should be evaluated as a vaccine adjuvant when chemically coupled to antigens of choice.

Evolution of substrate specificity in a retained enzyme driven by gene loss

The connection between gene loss and the functional adaptation of retained proteins is still poorly understood. We apply phylogenomics and metabolic modeling to detect bacterial species that are evolving by gene loss, with the finding that Actinomycetaceae genomes from human cavities are undergoing sizable reductions, including loss of L-histidine and L-tryptophan biosynthesis. We observe that the dual-substrate phosphoribosyl isomerase A or priA gene, at which these pathways converge, appears to coevolve with the occurrence of trp and his genes. Characterization of a dozen PriA homologs shows that these enzymes adapt from bifunctionality in the largest genomes, to a monofunctional, yet not necessarily specialized, inefficient form in genomes undergoing reduction. These functional changes are accomplished via mutations, which result from relaxation of purifying selection, in residues structurally mapped after sequence and X-ray structural analyses. Our results show how gene loss can drive the evolution of substrate specificity from retained enzymes. DOI: http://dx.doi.org/10.7554/eLife.22679.001

Multi-Phase Solvation Model for Biological Membranes: Molecular Action Mechanism of Amphotericin B

Amphotericin B (AmB) is still the most effective drug for the treatment of systemic fungal infections in humans. Despite significant theoretical and experimental efforts trying to understand its molecular mechanism of action, the answer has remained elusive. In this work, we present a computational methodology to test the current membrane related hypotheses, namely, transmembrane ion channel, adsorption, and sterol sponge. We use a thermodynamic approach in which we represent the membrane by a multiphase solvation model with atomic detail (MMPSM) and calculate the free energy of transferring the drug between phases with different dielectric properties. Furthermore, we compare AmB to a chemical analogue with increased safety, an l-histidine methyl ester of AmB. Our findings reveal that both drugs dimerize in all solvents studied here. Also, it is energetically unfavorable for the drugs to penetrate into the hydrophobic core of the membrane, unless their concentration is high. Finally, it is thermodynamically possible that the sterols migrate from the membrane into a drug droplet adsorbed at the surface of the bilayer. In light of our results, several effects could take place in the complex antibiotic process. We suggest a molecular mechanism that connects all three hypotheses through a drug concentration dependence and propose that the drug promotes the formation of membrane toroidal pores. Because MMPSM is of general interest, we made it available at http://tripplab.com/tools/mmpsm .