pyKVFinder: an efficient and integrable Python package for biomolecular cavity detection and characterization in data science.

BACKGROUND: Biomolecular interactions that modulate biological processes occur mainly in cavities throughout the surface of biomolecular structures. In the data science era, structural biology has benefited from the increasing availability of biostructural data due to advances in structural determination and computational methods. In this scenario, data-intensive cavity analysis demands efficient scripting routines built on easily manipulated data structures. To fulfill this need, we developed pyKVFinder, a Python package to detect and characterize cavities in biomolecular structures for data science and automated pipelines. RESULTS: pyKVFinder efficiently detects cavities in biomolecular structures and computes their volume, area, depth and hydropathy, storing these cavity properties in NumPy arrays. Benefited from Python ecosystem interoperability and data structures, pyKVFinder can be integrated with third-party scientific packages and libraries for mathematical calculations, machine learning and 3D visualization in automated workflows. As proof of pyKVFinder's capabilities, we successfully identified and compared ADRP substrate-binding site of SARS-CoV-2 and a set of homologous proteins with pyKVFinder, showing its integrability with data science packages such as matplotlib, NGL Viewer, SciPy and Jupyter notebook. CONCLUSIONS: We introduce an efficient, highly versatile and easily integrable software for detecting and characterizing biomolecular cavities in data science applications and automated protocols. pyKVFinder facilitates biostructural data analysis with scripting routines in the Python ecosystem and can be building blocks for data science and drug design applications.

Novel quaternary structures of the human prion protein globular domain.

Prion disease is caused by the misfolding of the cellular prion protein, PrPC, into a self-templating conformer, PrPSc. Nuclear magnetic resonance (NMR) and X-ray crystallography revealed the 3D structure of the globular domain of PrPC and the possibility of its dimerization via an interchain disulfide bridge that forms due to domain swap or by non-covalent association of two monomers. On the contrary, PrPSc is composed by a complex and heterogeneous ensemble of poorly defined conformations and quaternary arrangements that are related to different patterns of neurotoxicity. Targeting PrPC with molecules that stabilize the native conformation of its globular domain emerged as a promising approach to develop anti-prion therapies. One of the advantages of this approach is employing structure-based drug discovery methods to PrPC. Thus, it is essential to expand our structural knowledge about PrPC as much as possible to aid such drug discovery efforts. In this work, we report a crystallographic structure of the globular domain of human PrPC that shows a novel dimeric form and a novel oligomeric arrangement. We use molecular dynamics simulations to explore its structural dynamics and stability and discuss potential implications of these new quaternary structures to the conversion process.

An Intestine/Liver Microphysiological System for Drug Pharmacokinetic and Toxicological Assessment.

The recently introduced microphysiological systems (MPS) cultivating human organoids are expected to perform better than animals in the preclinical tests phase of drug developing process because they are genetically human and recapitulate the interplay among tissues. In this study, the human intestinal barrier (emulated by a co-culture of Caco-2 and HT-29 cells) and the liver equivalent (emulated by spheroids made of differentiated HepaRG cells and human hepatic stellate cells) were integrated into a two-organ chip (2-OC) microfluidic device to assess some acetaminophen (APAP) pharmacokinetic (PK) and toxicological properties. The MPS had three assemblies: Intestine only 2-OC, Liver only 2-OC, and Intestine/Liver 2-OC with the same media perfusing both organoids. For PK assessments, we dosed the APAP in the media at preset timepoints after administering it either over the intestinal barrier (emulating the oral route) or in the media (emulating the intravenous route), at 12 µM and 2 µM respectively. The media samples were analyzed by reversed-phase high-pressure liquid chromatography (HPLC). Organoids were analyzed for gene expression, for TEER values, for protein expression and activity, and then collected, fixed, and submitted to a set of morphological evaluations. The MTT technique performed well in assessing the organoid viability, but the high content analyses (HCA) were able to detect very early toxic events in response to APAP treatment. We verified that the media flow does not significantly affect the APAP absorption whereas it significantly improves the liver equivalent functionality. The APAP human intestinal absorption and hepatic metabolism could be emulated in the MPS. The association between MPS data and in silico modeling has great potential to improve the predictability of the in vitro methods and provide better accuracy than animal models in pharmacokinetic and toxicological studies.

The citrus flavonoid naringenin impairs the in vitro infection of human cells by Zika virus.

The Zika virus (ZIKV) is an arthropod-borne virus that belongs to the Flaviviridae family. The ZIKV infection is usually asymptomatic or is associated with mild clinical manifestations; however, increased numbers of cases of microcephaly and birth defects have been recently reported. To date, neither a vaccine nor an antiviral treatment has become available to control ZIKV replication. Among the natural compounds recognized for their medical properties, flavonoids, which can be found in fruits and vegetables, have been found to possess biological activity against a variety of viruses. Here, we demonstrate that the citrus flavanone naringenin (NAR) prevented ZIKV infection in human A549 cells in a concentration-dependent and ZIKV-lineage independent manner. NAR antiviral activity was also observed when primary human monocyte-derived dendritic cells were infected by ZIKV. NAR displayed its antiviral activity when the cells were treated after infection, suggesting that NAR acts on the viral replication or assembly of viral particles. Moreover, a molecular docking analysis suggests a potential interaction between NAR and the protease domain of the NS2B-NS3 protein of ZIKV which could explain the anti-ZIKV activity of NAR. Finally, the results support the potential of NAR as a suitable candidate molecule for developing anti-ZIKV treatments.

Synthesis of novel triazole-derived glycopeptides as analogs of α-dystroglycan mucins.

α-Dystroglycan (α-DG) mucins are essential for maintenance of the structural and functional stability of the muscle fiber and, when hypoglycosylated, they are directly involved in pathological processes such as dystroglycanopathies. Thus, this work reports the synthesis of the novel 1,2,3-triazole-derived glycosyl amino acids αGlcNAc-1-O-triazol-2Manα-ThrOH (1) and Gal-ß1,4-αGlcNAc-1-O-triazol-2Manα-ThrOH (2), followed by solid-phase assembly to get the corresponding glycopeptides NHAcThrVal[αGlcNAc-1-triazol-2Manα]ThrIleArgGlyOH (3) and NHAcThrVal[Gal-ß1,4-αGlcNAc-1-triazol-2Manα]ThrIleArgGlyOH (4) as analogs of α-DG mucins. The glycosyl amino acids 1 (72%) and 2 (35%) were synthesized by Cu(I)-assisted 1,3-dipolar azide-alkyne cycloaddition reactions (CuAAC) between the azide-glycosyl amino acid αManN3-FmocThrOBn (5) and the corresponding alkyne-functionalyzed sugars 2'-propynyl-αGlcNAc (6) and 2'-propynyl-Gal-ß1,4-αGlcNAc (7), followed by hydrogenation reactions. Subsequently, glycopeptides 3 (23%) and 4 (12%) were obtained by solid phase synthesis, involving sequential couplings of Fmoc-protected amino acids or the glycosyl amino acids 1 and 2, followed by cleavage from resin, N-acetylation and O-deacetylation (NaOMe) reactions. Lastly, enzymatic galactosylation of glycopeptide 3 with bovine ß-1,4-GalT showed that it was not a substrate for this enzyme, which could be better elucidated by docking simulations with ß-1,4-GalT.

An analysis and evaluation of the WeFold collaborative for protein structure prediction and its pipelines in CASP11 and CASP12.

Every two years groups worldwide participate in the Critical Assessment of Protein Structure Prediction (CASP) experiment to blindly test the strengths and weaknesses of their computational methods. CASP has significantly advanced the field but many hurdles still remain, which may require new ideas and collaborations. In 2012 a web-based effort called WeFold, was initiated to promote collaboration within the CASP community and attract researchers from other fields to contribute new ideas to CASP. Members of the WeFold coopetition (cooperation and competition) participated in CASP as individual teams, but also shared components of their methods to create hybrid pipelines and actively contributed to this effort. We assert that the scale and diversity of integrative prediction pipelines could not have been achieved by any individual lab or even by any collaboration among a few partners. The models contributed by the participating groups and generated by the pipelines are publicly available at the WeFold website providing a wealth of data that remains to be tapped. Here, we analyze the results of the 2014 and 2016 pipelines showing improvements according to the CASP assessment as well as areas that require further adjustments and research.

Membrane vesiculation induced by proteins of the dengue virus envelope studied by molecular dynamics simulations.

Biological membranes are continuously remodeled in the cell by specific membrane-shaping machineries to form, for example, tubes and vesicles. We examine fundamental mechanisms involved in the vesiculation processes induced by a cluster of envelope (E) and membrane (M) proteins of the dengue virus (DENV) using molecular dynamics simulations and a coarse-grained model. We show that an arrangement of three E-M heterotetramers (EM3) works as a bending unit and an ordered cluster of five such units generates a closed vesicle, reminiscent of the virus budding process. In silico mutagenesis of two charged residues of the anchor helices of the envelope proteins of DENV shows that Arg-471 and Arg-60 are fundamental to produce bending stress on the membrane. The fine-tuning between the size of the EM3 unit and its specific bending action suggests this protein unit is an important factor in determining the viral particle size.

Synthetic 1,2,3-triazole-linked glycoconjugates bind with high affinity to human galectin-3.

This work describes the synthesis of the 1,2,3-triazole amino acid-derived-3-O-galactosides 1-6 and the 1,2,3-triazole di-lactose-derived glycoconjugate 7 as potential galectin-3 inhibitors. The target compounds were synthesized by Cu(I)-catalyzed azide-alkyne cycloaddition reaction ('click chemistry') between the azido-derived amino acids N3-ThrOBn, N3-PheOBn, N3-N-Boc-TrpOBn, N3-N-Boc-LysOBn, N3-O-tBu-AspOBn and N3-l-TyrOH, and the corresponding alkyne-based sugar 3-O-propynyl-GalOMe, as well as by click chemistry reaction between the azido-lactose and 2-propynyl lactose. Surface plasmon resonance (SPR) assays showed that all synthetic glycoconjugates 1-7 bound to galectin-3 with high affinity, but the highest binders were the amino acids-derived glycoconjugates 2 (KD 7.96µM) and 4 (KD 4.56µM), and the divalent lactoside 7 (KD1 0.15µM/KD2 19µM). Molecular modeling results were in agreement with SPR assays, since more stable interactions with galectin-3 were identified for glycoconjugates 2, 4 and 7. Regarding compounds 2 and 4, they established specific cation-π (Arg144) and ionic (Asp148) interactions, whereas glycoconjugate 7 was capable to bridge two independent galectin-3 CRDs, creating a non-covalent cross-link between two monomers and, thus, reaching a submicromolar affinity towards galectin-3.

In silico assessment of S100A12 monomer and dimer structural dynamics: implications for the understanding of its metal-induced conformational changes.

Changes in the concentration of different ions modulate several cellular processes, such as Ca(2+) and Zn(2+) in inflammation. Upon activation of immune system effector cells, the intracellular Ca(2+) concentration rises propagating the activation signal, leading to degranulation and generation of reactive oxygen species, which increases the Zn(2+) intracellular concentration as a consequence of the cellular antioxidant machinery. In this context, S100A12 is of special interest because it is a pro-inflammatory protein expressed in neutrophils whose structure and function are modulated by both Ca(2+) and Zn(2+). The current hypothesis about its mechanism of action was built based on biochemical and crystallographic data. However, there are missing connections between molecular structure and the way in which many events are concatenated at the triggering and along the inflammatory process. In this work we use molecular dynamics simulations to describe how variations in Zn(2+) and Ca(2+) concentrations modulate the structural dynamics of the calcium-free S100A12 dimer and monomer, which was not considered a part of the mechanism of action before. Our results suggest that (i) Zn(2+) have a determinant role in the dimerization step, as well as in the unbinding of the Na(+) complexed to the N-terminal EF-hand; (ii) the N-terminal EF-hand domain is the first to bind Ca(2+), and not the C-terminal, as usually accepted; and that (iii) Ca(2+) modulates the structural dynamics of H-III.