Luis del Pozo-Yauner

Combining metagenomics, metatranscriptomics and viromics to explore novel microbial interactions: towards a systems-level understanding of human microbiome

The advances in experimental methods and the development of high performance bioinformatic tools have substantially improved our understanding of microbial communities associated with human niches. Many studies have documented that changes in microbial abundance and composition of the human microbiome is associated with human health and diseased state. The majority of research on human microbiome is typically focused in the analysis of one level of biological information, i.e., metagenomics or metatranscriptomics. In this review, we describe some of the different experimental and bioinformatic strategies applied to analyze the 16S rRNA gene profiling and shotgun sequencing data of the human microbiome. We also discuss how some of the recent insights in the combination of metagenomics, metatranscriptomics and viromics can provide more detailed description on the interactions between microorganisms and viruses in oral and gut microbiomes. Recent studies on viromics have begun to gain importance due to the potential involvement of viruses in microbial dysbiosis. In addition, metatranscriptomic combined with metagenomic analysis have shown that a substantial fraction of microbial transcripts can be differentially regulated relative to their microbial genomic abundances. Thus, understanding the molecular interactions in the microbiome using the combination of metagenomics, metatranscriptomics and viromics is one of the main challenges towards a system level understanding of human microbiome.

Thermodynamic and Kinetic Characterization of a Germ Line Human λ6 Light-Chain Protein: The Relation between Unfolding and Fibrillogenesis

Proteins encoded by the gene segment 6a of the lambda variable light-chain repertoire are strongly associated with amyloid deposition. 6aJL2 is a model protein constructed with the predicted sequences encoded by the 6a and JL2 germ line genes. In this work, we characterized the urea- and temperature-induced unfolding of 6aJL2. In the short time scale, spectroscopic, hydrodynamic and calorimetric experiments were compatible with a two-state transition. Furthermore, DeltaG, m and the midpoint urea concentration obtained from equilibrium experiments were compatible with those obtained from kinetic experiments. Since fibril formation is a slow process, samples were also incubated for longer times. After incubation for several hours at 37 degrees C, spectroscopic, hydrodynamic and calorimetric experiments revealed the presence of a partially unfolded off-pathway intermediate around the midpoint urea concentration (1.5-3.0 M urea). In vitro fibrillogenesis assays show that the maximum growth rate for fibril formation and the minimum lag time were obtained at urea concentrations where the partially unfolded state was populated (2.5 M urea at 37 degrees C). This indicates that this partially unfolded state is critical for in vitro fibril formation. Concentration-dependent kinetics and hydrodynamic properties of the intermediate were consistent with a soluble oligomeric state. The intermediate is formed around the midpoint urea concentration, where the native and unfolded states are equally populated and their rate of interconversion is the slowest. This situation may promote the slow accumulation of an intermediate state that is prone to aggregate.

High-Quality Draft Genomes of Two Vibrio parahaemolyticus Strains Aid in Understanding Acute Hepatopancreatic Necrosis Disease of Cultured Shrimps in Mexico

ABSTRACT The high-quality draft genomes of two Vibrio parahaemolyticus strains, one that causes the acute hepatopancreatic necrosis disease (AHPND) in cultured shrimps (FIM-S1708+), and another that does not (FIM-S1392−) are reported. A chromosome-scale assembly for the FIM-S1392− genome is reported here. The analysis of the two genomes gives some clues regarding the genomic differences between the strains.

The N-terminal strand modulates immunoglobulin light chain fibrillogenesis

It has been suggested that the N-terminal strand of the light chain variable domain (V(L)) protects the molecule from aggregation by hindering spurious intermolecular contacts. We evaluated the impact of mutations in the N-terminal strand on the thermodynamic stability and kinetic of fibrillogenesis of the V(L) protein 6aJL2. Mutations in this strand destabilized the protein in a position-dependent manner, accelerating the fibrillogenesis by shortening the lag time; an effect that correlated with the extent of destabilization. In contrast, the effect on the kinetics of fibril elongation, as assessed in seeding experiments was of different nature, as it was not directly dependant on the degree of destabilization. This finding suggests different factors drive the nucleation-dependent and elongation phases of light chain fibrillogenesis. Finally, taking advantage of the dependence of the Trp fluorescence upon environment, four single Trp substitutions were made in the N-terminal strand, and changes in solvent exposure during aggregation were evaluated by acrylamide-quenching. The results suggest that the N-terminal strand is buried in the fibrillar state of 6aJL2 protein. This finding suggest a possible explanation for the modulating effect exerted by the mutations in this strand on the aggregation behavior of 6aJL2 protein.

Mutational and genetic determinants of λ6 light chain amyloidogenesis

Approximately 25% of the λ6 light chains have glycine rather than arginine at position 25, which is an allelic variant of the IGLV6‐57 (6a) locus. The Gly25 variant has been shown to decrease the folding stability of the germline λ6 VL protein 6aJL2 by 1.7 kcal·mol−1. In this work, we compared the thermodynamic and fibrillogenic properties of the amyloidosis (AL) derived recombinant (r) VL protein AR, which contains the allelic variant Gly25, with those of germline rVL 6aJL2‐R25G and the λ6 disease‐associated VL proteins Wil (AL) and Jto (myeloma). Our experiments show that of the four proteins AR is the least stable; forms amyloid fibrils at physiological temperature, pH and ionic strength; has the shortest lag time; and elongates homologous seeds most efficiently. We conclude that the Gly25 allelic variant, together with the somatic mutations, contributes importantly to the extremely low stability and high amyloidogenicity of the AL‐derived protein AR.

Secretome Prediction of Two M. tuberculosis Clinical Isolates Reveals Their High Antigenic Density and Potential Drug Targets

The Excreted/Secreted (ES) proteins play important roles during Mycobacterium tuberculosis invasion, virulence, and survival inside the host and they are a major source of immunogenic proteins. However, the molecular complexity of the bacillus cell wall has made difficult the experimental isolation of the total bacterial ES proteins. Here, we reported the genomes of two Beijing genotype M. tuberculosis clinical isolates obtained from patients from Vietnam (isolate 46) and South Africa (isolate 48). We developed a bioinformatics pipeline to predict their secretomes and observed that ~12% of the genome-encoded proteins are ES, being PE, PE-PGRS, and PPE the most abundant protein domains. Additionally, the Gene Ontology, KEGG pathways and Enzyme Classes annotations supported the expected functions for the secretomes. The ~70% of an experimental secretome compiled from literature was contained in our predicted secretomes, while only the 34–41% of the experimental secretome was contained in the two previously reported secretomes for H37Rv. These results suggest that our bioinformatics pipeline is better to predict a more complete set of ES proteins in M. tuberculosis genomes. The predicted ES proteins showed a significant higher antigenic density measured by Abundance of Antigenic Regions (AAR) value than the non-ES proteins and also compared to random constructed secretomes. Additionally, we predicted the secretomes for H37Rv, H37Ra, and two M. bovis BCG genomes. The antigenic density for BGG and for isolates 46 and 48 was higher than the observed for H37Rv and H37Ra secretomes. In addition, two sets of immunogenic proteins previously reported in patients with tuberculosis also showed a high antigenic density. Interestingly, mice infected with isolate 46 showed a significant lower survival rate than the ones infected with isolate 48 and both survival rates were lower than the one previously reported for the H37Rv in the same murine model. Finally, after a druggability analysis of the secretomes, we found potential drug targets such as cytochrome P450, thiol peroxidase, the Ag85C, and Ribonucleoside Reductase in the secreted proteins that could be used as drug targets for novel treatments against Tuberculosis.

Importance of the interaction protein-protein of the CaM-PDE1A and CaM-MLCK complexes in the development of new anti-CaM drugs.

Protein–protein interactions play central roles in physiological and pathological processes. The bases of the mechanisms of drug action are relevant to the discovery of new therapeutic targets. This work focuses on understanding the interactions in protein–protein–ligands complexes, using proteins calmodulin (CaM), human calcium/calmodulin‐dependent 3′,5′‐cyclic nucleotide phosphodiesterase 1A active human (PDE1A), and myosin light chain kinase (MLCK) and ligands αII–spectrin peptide (αII–spec), and two inhibitors of CaM (chlorpromazine (CPZ) and malbrancheamide (MBC)). The interaction was monitored with a fluorescent biosensor of CaM (hCaM M124C–mBBr). The results showed changes in the affinity of CPZ and MBC depending on the CaM–protein complex under analysis. For the Ca2+–CaM, Ca2+–CaM–PDE1A, and Ca2+–CaM–MLCK complexes, CPZ apparent dissociation constants (Kds) were 1.11, 0.28, and 0.55 μM, respectively; and for MBC Kds were 1.43, 1.10, and 0.61 μM, respectively. In competition experiments the addition of calmodulin binding peptide 1 (αII–spec) to Ca2+–hCaM M124C–mBBr quenched the fluorescence (Kd = 2.55 ± 1.75 pM) and the later addition of MBC (up to 16 μM) did not affect the fluorescent signal. Instead, the additions of αII–spec to a preformed Ca2+–hCaM M124C–mBBr–MBC complex modified the fluorescent signal. However, MBC was able to displace the PDE1A and MLCK from its complex with Ca2+–CaM. In addition, docking studies were performed for all complexes with both ligands showing an excellent correlation with experimental data. These experiments may help to explain why in vivo many CaM drugs target prefer only a subset of the Ca2+–CaM regulated proteins and adds to the understanding of molecular interactions between protein complexes and small ligands. Copyright

Evaluation of three different formats of a neutralizing single chain human antibody against toxin Cn2: Neutralization capacity versus thermodynamic stability

The single-chain antibody fragment (scFv) 6009F, obtained by directed evolution, neutralizes the effects of the Cn2 toxin, which is the major toxic component of Centruroides noxius scorpion venom. In this work we compared the neutralization capacity and the thermodynamic stability of scFv 6009F with those of two other derived formats: Fab 6009F and diabody 6009F. Additionally, the affinity constants to Cn2 toxin of the three recombinant antibody fragments were determined by means of BIAcore. We found a correlation between the thermodynamic stability of these antibody fragments with their neutralization capacity. The order of thermodynamic stability determined was Fab≫scFv>diabody. The Fab and scFv were capable of neutralizing the toxic effects of Cn2 and whole venom but the diabody was unable to fully neutralize intoxication. In silico analysis of the diabody format indicates that the reduction of stability and neutralization capacity could be explained by a less cooperative interface between the heavy and the light variable domains.

Insights into molecular interactions between CaM and its inhibitors from molecular dynamics simulations and experimental data

In order to contribute to the structural basis for rational design of calmodulin (CaM) inhibitors, we analyzed the interaction of CaM with 14 classic antagonists and two compounds that do not affect CaM, using docking and molecular dynamics (MD) simulations, and the data were compared to available experimental data. The Ca2+-CaM-Ligands complexes were simulated 20 ns, with CaM starting in the “open” and “closed” conformations. The analysis of the MD simulations provided insight into the conformational changes undergone by CaM during its interaction with these ligands. These simulations were used to predict the binding free energies (ΔG) from contributions ΔH and ΔS, giving useful information about CaM ligand binding thermodynamics. The ΔG predicted for the CaM’s inhibitors correlated well with available experimental data as the r2 obtained was 0.76 and 0.82 for the group of xanthones. Additionally, valuable information is presented here: I) CaM has two preferred ligand binding sites in the open conformation known as site 1 and 4, II) CaM can bind ligands of diverse structural nature, III) the flexibility of CaM is reduced by the union of its ligands, leading to a reduction in the Ca2+-CaM entropy, IV) enthalpy dominates the molecular recognition process in the system Ca2+-CaM-Ligand, and V) the ligands making more extensive contact with the protein have higher affinity for Ca2+-CaM. Despite their limitations, docking and MD simulations in combination with experimental data continue to be excellent tools for research in pharmacology, toward a rational design of new drugs.

1H, 13C and 15N resonance assignment of 6aJL2(R25G), a highly fibrillogenic λVI light chain variable domain

An allotypic variation at position 25 influences the fibrillogenicity of λVI light chains, which are related to humoral immune response and have been associated with AL amyloidosis. The full resonance assignment and a preliminary structural characterization of 6aJL2(R25G) are reported.