Search

Retinal Hemangioblastoma (RH) is the most frequent manifestation of the von Hippel-Lindau syndrome (VHL), a rare disease associated with the germline mutation of the von Hippel-Lindau protein (pVHL). An emblematic feature of RH is the high vascularity, which is easily explained by the overexpression of angiogenic factors (AFs) arising from the pVHL impairment. The introduction of Optical Coherence Tomography Angiography (OCTA) allowed observing this feature with exceptional detail. However, our understanding of RH is limited by the absence of an animal model fully recapitulating the tumor. Here, we exploit a cancer mathematical model as an alternative way to explore RH development and angiogenesis. We derived our model from the agreed pathology for this tumor and compared our results with patient-specific OCTA images. Our simulations closely resemble the medical images, proving the capability of our model to recapitulate RH pathology. Our results also suggest that angiogenesis in RH occurs suddenly when the tumor reaches a critical mass, with full capillary invasion in the order of days. These findings open a new perspective on the critical role of time in antiangiogenic therapy in RH, which has resulted ineffective. Indeed, it might be that when RH is diagnosed, angiogenesis is already too advanced to be effectively targeted with this mean.

The SARS-CoV-2 pandemic is maintained by the emergence of successive variants, highlighting the flexibility of the protein sequences of the virus. We show that experimentally determined intrinsically disordered regions (IDRs) are abundant in the SARS-CoV-2 viral proteins, making up to 28% of disorder content for the S1 subunit of spike and up to 51% for the nucleoprotein, with the vast majority of mutations occurring in the 13 major variants mapped to these IDRs. Strikingly, antigenic sites are enriched in IDRs, in the receptor-binding domain (RBD) and in the N-terminal domain (NTD), suggesting a key role of structural flexibility in the antigenicity of the SARS-CoV-2 protein surface. Mutations occurring in the S1 subunit and nucleoprotein (N) IDRs are critical for immune evasion and antibody escape, suggesting potential additional implications for vaccines and monoclonal therapeutic strategies. Overall, this suggests the presence of variable regions on S1 and N protein surfaces, which confer sequence and antigenic flexibility to the virus without altering its protein functions.

RING-PyMOL is a plugin for PyMOL providing a set of analysis tools for structural ensembles and molecular dynamic simulations. RING-PyMOL combines residue interaction networks, as provided by the RING software, with structural clustering to enhance the analysis and visualization of the conformational complexity. It combines precise calculation of non-covalent interactions with the power of PyMOL to manipulate and visualize protein structures. The plugin identifies and highlights correlating contacts and interaction patterns that can explain structural allostery, active sites, and structural heterogeneity connected with molecular function. It is easy to use and extremely fast, processing and rendering hundreds of models and long trajectories in seconds. RING-PyMOL generates a number of interactive plots and output files for use with external tools. The underlying RING software has been improved extensively. It is 10 times faster, can process mmCIF files and it identifies typed interactions also for nucleic acids. Availability and implementation https://github.com/BioComputingUP/ring-pymol

Older age is a major risk factor for adverse outcomes of COVID-19, potentially due to immunosenescence and chronic low-grade inflammation, both characteristics of older adults which synergistically contribute to their vulnerability. Furthermore, older age is also associated with decreased kidney function and is consequently associated with an increased risk of cardiovascular disease. All of this in the course of COVID-19 infection can worsen and promote the progression of chronic kidney damage and all its sequelae. Frailty is a condition characterized by the decline in function of several homeostatic systems, leading to increased vulnerability to stressors and risk of adverse health outcomes. Thus, it is very likely that frailty, together with comorbidities, may have contributed to the high vulnerability to severe clinical manifestations and deaths from COVID-19 among older people. The combination of viral infection and chronic inflammation in the elderly could cause multiple unforeseen harmful consequences, affecting overall disability and mortality rates. In post-COVID-19 patients, inflammation has been implicated in sarcopenia progression, functional activity decline, and dementia. After the pandemic, it is imperative to shine a spotlight on these sequelae so that we can be prepared for the future outcomes of the ongoing pandemic. Here, we discuss the potential long-term consequences of SARS-CoV-2 infection and its possibility of causing permanent damage to the precarious balance existing in the frail elderly with multiple pathologies.

Context Mutations in the androgen receptor (AR) gene might be associated with infertility mainly because they cause various degrees of androgen insensitivity. Objective The aim of the study was to evaluate the frequency and type of AR variants in a large cohort of infertile males. Methods A total of 8224 males of Italian idiopathic infertile couples were referred to the University Hospital of Padova. The main outcome measures were mutational screening of AR, computational, and functional analyses. Results We found 131 patients (1.6%) harboring 45 variants in AR gene, of which 18 were novel missense AR variants. Patients with AR gene variants had lower sperm count (P = .048), higher testosterone (T) concentration (P < .0001), and higher androgen sensitivity index (ASI) (luteinizing hormone × T, P < .001) than patients without variants. Statistical analyses found T ≥ 15.38 nmol/L and ASI ≥ 180 IU × nmol/L2 as the threshold values to discriminate with good accuracy patients with AR variants. Patients with oligozoospermia and T ≥ 15.38 nmol/L had a 9-fold increased risk of harboring mutations compared with patients with normal sperm count and T < 15.38 nmol/L (odds ratio 9.29, 95% CI 5.07-17.02). Using computational and functional approaches, we identified 2 novel variants, L595P and L791I, as potentially pathogenic. Conclusion This is the largest study screening AR gene variants in men of idiopathic infertile couples. We found that the prevalence of variants increased to 3.4% in oligozoospermic subjects with T ≥ 15.38 nmol/L. Conversely, more than 80% of men with AR gene variants had low sperm count and high T levels. Based on our findings, we suggest AR sequencing as a routine genetic test in cases of idiopathic oligozoospermia with T ≥ 15.38 nmol/L.

Cerebral sinus venous thrombosis (CSVT) is a relatively rare acute disorder of cerebral circulation, but it can potentially be associated with serious sequelae and a poor prognosis. The neurological manifestations associated with it are often not adequately taken into consideration given the extreme variability and nuances of its clinical presentation and given the need for radiological methods suitable for this type of diagnosis. CSVT is usually more common in women, but so far there are little data available in the literature on sex-specific characteristics regarding this pathology. CSVT is the result of multiple conditions and is therefore to be considered a multifactorial disease where at least one risk factor is present in over 80% of cases. From the literature, we learn that congenital or acquired prothrombotic states are to be considered extremely associated with the occurrence of an acute episode of CSVT and its recurrences. It is, therefore, necessary to fully know the origins and natural history of CSVT, in order to implement the diagnostic and therapeutic pathways of these neurological manifestations. In this report, we summarize the main causes of CSVT considering the possible influence of gender, bearing in mind that most of the causes listed above are pathological conditions closely linked to the female sex.

The Prolyl Hydroxylases (PHDs) are an enzymatic family that regulates cell oxygen-sensing. PHDs hydroxylate hypoxia-inducible transcription factors α (HIFs-α) driving their proteasomal degradation. Hypoxia inhibits PHDs activity, inducing HIFs-α stabilization and cell adaptation to hypoxia. As a hallmark of cancer, hypoxia promotes neo-angiogenesis and cell proliferation. PHD isoforms are thought to have a variable impact on tumor progression. All isoforms hydroxylate HIF-α (HIF-1,2,3α) with different affinities. However, what determines these differences and how they pair with tumor growth is poorly understood. Here, molecular dynamics simulations were used to characterize the PHD2 binding properties in complexes with HIF-1α and HIF-2α. In parallel, conservation analysis and binding free energy calculations were performed to better understand PHD2 substrate affinity. Our data suggest a direct association between the PHD2 C-terminus and HIF-2α that is not observed in the PHD2/HIF-1α complex. Furthermore, our results indicate that phosphorylation of a PHD2 residue, Thr405, causes a variation in binding energy, despite the fact that this PTM has only a limited structural impact on PHD2/HIFs-α complexes. Collectively, our findings suggest that the PHD2 C-terminus may act as a molecular regulator of PHD’s activity.

The mitochondrial ATP synthase is a multi-subunit enzyme complex located in the inner mitochondrial membrane which is essential for oxidative phosphorylation under physiological conditions. In this review, we analyse the enzyme functions involved in cancer progression by dissecting specific conditions in which ATP synthase contributes to cancer development or metastasis. Moreover, we propose the role of ATP synthase in the formation of the permeability transition pore (PTP) as an additional mechanism which controls tumour cell death. We further describe transcriptional and translational modifications of the enzyme subunits and of the inhibitor protein IF1 that may promote adaptations leading to cancer metabolism. Finally, we outline ATP synthase gene mutations and epigenetic modifications associated with cancer development or drug resistance, with the aim of highlighting this enzyme complex as a potential novel target for future anti-cancer therapy.

Mutations in the RYR1 gene, encoding ryanodine receptor 1 (RyR1), are a well-known cause of Central Core Disease (CCD) and Multi-minicore Disease (MmD). We screened a cohort of 153 patients carrying an histopathological diagnosis of core myopathy (cores and minicores) for RYR1 mutation. At least one RYR1 mutation was identified in 69 of them and these patients were further studied. Clinical and histopathological features were collected. Clinical phenotype was highly heterogeneous ranging from asymptomatic or paucisymptomatic hyperCKemia to severe muscle weakness and skeletal deformity with loss of ambulation. Sixty-eight RYR1 mutations, generally missense, were identified, of which 16 were novel. The combined analysis of the clinical presentation, disease progression and the structural bioinformatic analyses of RYR1 allowed to associate some phenotypes to mutations in specific domains. In addition, this study highlighted the structural bioinformatics potential in the prediction of the pathogenicity of RYR1 mutations. Further improvement in the comprehension of genotype–phenotype relationship of core myopathies can be expected in the next future: the actual lack of the human RyR1 crystal structure paired with the presence of large intrinsically disordered regions in RyR1, and the frequent presence of more than one RYR1 mutation in core myopathy patients, require designing novel investigation strategies to completely address RyR1 mutation effect.

Functional impairment of the von Hippel–Lindau tumor suppressor (pVHL) is causative of a familiar increased risk of developing cancer. As an E3 substrate recognition particle, pVHL marks the hypoxia inducible factor 1α (HIF-1α) for degradation in normoxic conditions, thus acting as a key regulator of both acute and chronic cell adaptation to hypoxia. The male mice model carrying VHL gene conditional knockout presents significant abnormalities in testis development paired with defects in spermatogenesis and infertility, indicating that pVHL exerts testis-specific roles. Here we aimed to explore whether pVHL could have a similar role in humans by performing a testis-tissue library screening complemented with in-depth bioinformatics analysis. We identified 55 novel pVHL binding proteins directly involved in spermatogenesis, cell differentiation and reproductive metabolism. In addition, computational investigation of these new interactors identified multiple pVHL-specific binding motifs and demonstrated that somatic mutations described in human cancers reside in these binding regions. Collectively, these findings suggest that, in addition to its role in cancer formation, pVHL may also be pivotal in normal gonadal development in humans.

Summary Mathematical models are effective in studying cancer development at different scales from metabolism to tissue. Phase Field Models (PFMs) have been shown to reproduce accurately cancer growth and other related phenomena, including expression of relevant molecules, extracellular matrix remodeling and angiogenesis. However, implementations of such models are rarely published, reducing access to these techniques. To reduce this gap, we developed Mocafe, a modular open-source Python package that implements some of the most important PFMs reported in the literature. Mocafe is designed to handle both PFMs purely based on differential equations and hybrid agent-based PFMs. Moreover, Mocafe is meant to be extensible, allowing the inclusion of new models in future releases. Availability and implementation Mocafe is a Python package based on FEniCS, a popular computing platform for solving partial differential equations. The source code, extensive documentation and demos are provided on GitHub at URL: https://github.com/BioComputingUP/mocafe. Moreover, we uploaded on Zenodo an archive of the package, which is available at https://doi.org/10.5281/zenodo.6366052.

Background/Aims: The permeability transition pore (PTP) is an unselective, Ca2+-dependent high conductance channel of the inner mitochondrial membrane whose molecular identity has long remained a mystery. The most recent hypothesis is that pore formation involves the F-ATP synthase, which consistently generates Ca2+-activated channels. Available structures do not display obvious features that can accommodate a channel; thus, how the pore can form and whether its activity can be entirely assigned to F-ATP synthase is the matter of debate. In this study, we investigated the role of F-ATP synthase subunits e, g and b in PTP formation. Methods: Yeast null mutants for e, g and the first transmembrane (TM) α-helix of subunit b were generated and evaluated for mitochondrial morphology (electron microscopy), membrane potential (Rhodamine123 fluorescence) and respiration (Clark electrode). Homoplasmic C23S mutant of subunit a was generated by in vitro mutagenesis followed by biolistic transformation. F-ATP synthase assembly was evaluated by BN-PAGE analysis. Cu2+ treatment was used to induce the formation of F-ATP synthase dimers in the absence of e and g subunits. The electrophysiological properties of F-ATP synthase were assessed in planar lipid bilayers. Results: Null mutants for the subunits e and g display dimer formation upon Cu2+ treatment and show PTP-dependent mitochondrial Ca2+ release but not swelling. Cu2+ treatment causes formation of disulfide bridges between Cys23 of subunits a that stabilize dimers in absence of e and g subunits and favors the open state of wild-type F-ATP synthase channels. Absence of e and g subunits decreases conductance of the F-ATP synthase channel about tenfold. Ablation of the first TM of subunit b, which creates a distinct lateral domain with e and g, further affected channel activity. Conclusion: F-ATP synthase e, g and b subunits create a domain within the membrane that is critical for the generation of the high-conductance channel, thus is a prime candidate for PTP formation. Subunits e and g are only present in eukaryotes and may have evolved to confer this novel function to F-ATP synthase.

The Protein Ensemble Database (PED) (https://proteinensemble.org), which holds structural ensembles of intrinsically disordered proteins (IDPs), has been significantly updated and upgraded since its last release in 2016. The new version, PED 4.0, has been completely redesigned and reimplemented with cutting-edge technology and now holds about six times more data (162 versus 24 entries and 242 versus 60 structural ensembles) and a broader representation of state of the art ensemble generation methods than the previous version. The database has a completely renewed graphical interface with an interactive feature viewer for region-based annotations, and provides a series of descriptors of the qualitative and quantitative properties of the ensembles. High quality of the data is guaranteed by a new submission process, which combines both automatic and manual evaluation steps. A team of biocurators integrate structured metadata describing the ensemble generation methodology, experimental constraints and conditions. A new search engine allows the user to build advanced queries and search all entry fields including cross-references to IDP-related resources such as DisProt, MobiDB, BMRB and SASBDB. We expect that the renewed PED will be useful for researchers interested in the atomic-level understanding of IDP function, and promote the rational, structure-based design of IDP-targeting drugs.

Alstrom syndrome (OMIM#203800) is an ultra-rare autosomal recessive monogenic disease presenting pathogenic variants in ALMS1 (chromosome 2p13). It is characterized by early onset of blindness, hearing loss and systemic comorbidities, with delayed development without cognitive impairment. We aimed to investigate the cognitive functions and describe new pathogenic variants in Alstrom syndrome patients. Nineteen patients (13 adults, 6 children) underwent a thorough clinical, genetic, laboratory, instrumental, and neurocognitive assessment. Six new pathogenic variants in ALMS1 including the first described in exon 6 were identified. Four patients displayed a "mild phenotype" characterized by slow disease onset or absence of complications, including childhood obesity and association with at least one pathogenic variant in exon 5 or 6. At neurocognitive testing, a significant proportion of patients had deficits in three neurocognitive domains: similarities, phonological memory, and apraxia. In particular, 53% of patients showed difficulties in the auditory working memory test. We found ideomotor and buccofacial apraxia in 74% of patients. "Mild phenotype" patients performed better on auditory working memory and ideomotor apraxia test than "typical phenotype" ones (91.9 + 16.3% vs. 41.7 + 34.5% of correct answers, Z = 64.5, p

The mitochondrial paralog of the Hsp90 chaperone family TRAP1 is often induced in tumors, but the mechanisms controlling its expression, as well as its physiological functions remain poorly understood. Here, we find that TRAP1 is highly expressed in the early stages of Zebrafish development, and its ablation delays embryogenesis while increasing mitochondrial respiration of fish larvae. TRAP1 expression is enhanced by hypoxic conditions both in developing embryos and in cancer models of Zebrafish and mammals. The TRAP1 promoter contains evolutionary conserved hypoxic responsive elements, and HIF1α stabilization increases TRAP1 levels. TRAP1 inhibition by selective compounds or by genetic knock-out maintains a high level of respiration in Zebrafish embryos after exposure to hypoxia.Our data identify TRAP1 as a primary regulator of mitochondrial bioenergetics in highly proliferating cells following reductions in oxygen tension and HIF1α stabilization.

Mutations of the von Hippel-Lindau (pVHL) tumor suppressor are causative of a familiar predisposition to develop different types of cancer. pVHL is mainly known for its role in regulating hypoxia-inducible factor 1-α (HIF-1α) degradation, thus modulating the hypoxia response. There are different pVHL isoforms, including pVHL30 and pVHL19. However, little is known about isoform-specific functions and protein-protein interactions. Integrating in silico predictions with in vitro and in vivo assays, we describe a novel interaction between pVHL and mouse double minute 2 homolog (MDM2). Importantly, we found that pVHL30, and not pVHL19, forms a complex with MDM2, and that the N-terminal acidic tail of pVHL30 is required for its association with MDM2. Further, we demonstrate that an intrinsically disordered region upstream of the tetramerization domain of MDM2 is responsible for its isoform-specific association with pVHL30. This region is highly conserved in higher mammals, including primates, similarly to what has been already proposed for the N-terminal tail of pVHL30. Finally, we show that overexpression of pVHL30 and MDM2 together reduces cell proliferation, suggesting a synergistic effect of these E3 ubiquitin ligases. Collectively, our data support the idea that pVHL30 plays a role in MDM2 regulation, suggesting a wider interplay among hypoxia sensing and cell cycle regulation.

Summary The f subunit is localized at the base of the ATP synthase peripheral stalk. Its function in the human enzyme is poorly characterized. Because full disruption of its ATP5J2 gene with the CRISPR-Cas9 strategy in the HAP1 human model has been shown to cause alterations in the amounts of other ATP synthase subunits, here we investigated the role of the f subunit in HeLa cells by regulating its levels through RNA interference. We confirm the role of the f subunit in ATP synthase dimer stability and observe that its downregulation per se does not alter the amounts of the other enzyme subunits or ATP synthase synthetic/hydrolytic activity. We show that downregulation of the f subunit causes abnormal crista organization and decreases permeability transition pore (PTP) size, whereas its re-expression in f subunit knockdown cells rescues mitochondrial morphology and PTP-dependent swelling.

The von Hippel–Lindau protein (pVHL) is a tumour suppressor mainly known for its role as master regulator of hypoxia-inducible factor (HIF) activity. Functional inactivation of pVHL is causative of the von Hippel–Lindau disease, an inherited predisposition to develop different cancers. Due to its impact on human health, pVHL has been widely studied in the last few decades. However, investigations mostly focus on its role in degrading HIFs, whereas alternative pVHL protein–protein interactions and functions are insistently surfacing in the literature. In this review, we analyse these almost neglected functions by dissecting specific conditions in which pVHL is proposed to have differential roles in promoting cancer. We reviewed its role in regulating phosphorylation as a number of works suggest pVHL to act as an inhibitor by either degrading or promoting downregulation of specific kinases. Further, we summarize hypoxia-dependent and -independent pVHL interactions with multiple protein partners and discuss their implications in tumorigenesis.

ObjectiveTo identify and characterize patients with calsequestrin 1 (CASQ1)–related myopathy.MethodsPatients selected according to histopathologic features underwent CASQ1 genetic screening. CASQ1-mutated patients were clinically evaluated and underwent muscle MRI. Vacuole morphology and vacuolated fiber type were characterized.ResultsTwenty-two CASQ1-mutated patients (12 families) were identified, 21 sharing the previously described founder mutation (p.Asp244Gly) and 1 with the p.Gly103Asp mutation. Patients usually presented in the sixth decade with exercise intolerance and myalgias and later developed mild to moderate, slowly progressive proximal weakness with quadriceps atrophy and scapular winging. Muscle MRI (n = 11) showed a recurrent fibrofatty substitution pattern. Three patients presented subclinical cardiac abnormalities. Muscle histopathology in patients with p.Asp244Gly showed vacuoles in type II fibers appearing empty in hematoxylin-eosin, Gomori, and nicotinamide adenine dinucleotide (NADH) tetrazolium reductase stains but strongly positive for sarcoplasmic reticulum proteins. The muscle histopathology of p.Gly103Asp mutation was different, showing also NADH-positive accumulation consistent with tubular aggregates.ConclusionsWe report the clinical and molecular details of the largest cohort of CASQ1-mutated patients. A possible heart involvement is presented, further expanding the phenotype of the disease. One mutation is common due to a founder effect, but other mutations are possible. Because of a paucity of symptoms, it is likely that CASQ1 mutations may remain undiagnosed if a muscle biopsy is not performed.

Osteoporosis is the most common bone disease characterized by reduced bone mass and increased bone fragility. Genetic contribution is one of the main causes of primary osteoporosis, therefore, both genders are affected by this skeletal disorder. Nonetheless, osteoporosis in men has received little attention, thus being underestimated and undertreated. The aim of this study was to identify novel genetic variants in a cohort of 128 males with idiopathic low bone mass using a next-generation sequencing (NGS) panel including genes whose mutations could result in reduced bone mineral density (BMD). Genetic analysis detected in eleven patients ten rare heterozygous variants within the LRP5 gene, which were categorized as VUS (variant of uncertain significance), likely pathogenic and benign variants according to American College of Medical Genetics and Genomics (ACMG) guidelines. Protein structural and Bayesian analysis performed on identified LRP5 variants pointed out p.R1036Q and p.R1135C as pathogenic, therefore suggesting the likely association of these two variants with the low bone mass phenotype. In conclusion, this study expands our understanding on the importance of a functional LRP5 protein in bone formation and highlights the necessity to sequence this gene in subjects with idiopathic low BMD.

The Gram-positive bacterium Streptococcus mutans is the principal causative agent of human tooth decay, an oral disease that affects the majority of the world’s population. Although the complete S. mutans genome is known, approximately 700 proteins are still annotated as hypothetical proteins, as no three-dimensional structure or homology with known proteins exists for them. Thus, the significant portion of genomic sequences coding for unknown-function proteins makes the knowledge of pathogenicity and survival mechanisms of S. mutans still incomplete. Plasmids are found in virtually every species of Streptococcus, and some of these mediate resistance to antibiotics and pathogenesis. However, there are strains of S. mutans that contain plasmids, such as LM7 and UA140, to which no function has been assigned yet. In this work, we describe an in silico study of the structure and function of all the S. mutans proteins encoded by pLM7 and pUA140 plasmids to gain insight into their biological function. A combination of different structural bioinformatics methodologies led to the identification of plasmidic proteins potentially required for the bacterial survival and pathogenicity. The structural information obtained on these proteins can be used to select novel targets for the design of innovative therapeutic agents towards S. mutans.

Background/Aims: In patients with colon cancer who undergo resection for potential cure, 40–60% have advanced locoregional disease (stage III). Those who are suitable for adjuvant treatment had a definite disease-free-survival benefit. The aim of the present study was to demonstrate whether the presence of desmoplasia influenced the mortality rate of stage III colorectal cancer (CRC) within 5 years from the surgery and adjuvant therapy. Patients and Methods: Sixty-five patients with stage III CRC underwent resection and adjuvant therapy. Qualitative categorization of desmoplasia was obtained using Ueno's stromal CRC classification. Desmoplasia was related to mortality using Spearman correlation and stratified with other histological variables (inflammation, grading) that concurred to the major determinant of malignancy (venous invasion and lymph nodes) using the Chi-square test. Result: The 5-year survival rate was 65% and the relapse rate was 37%. The mortality rate in patients with immature desmoplasia was 86%, 27% in intermediate desmoplasia, and 0% in mature desmoplasia (Spearman correlation coefficient: −0.572, P = 0.05). Conclusion: Immature desmoplasia appears to be associated with disease recurrence and mortality in stage III CRC patients.

We present an in silico characterization of the von Hippel-Lindau-like protein (VLP), the only known human paralog of the von Hippel-Lindau tumor suppressor protein (pVHL). Phylogenetic investigation showed VLP to be mostly conserved in upper mammals and specifically expressed in brain and testis. Structural analysis and molecular dynamics simulations show VLP to be very similar to pVHL three-dimensional organization and binding dynamics. In particular, conservation of elements at the protein interfaces suggests VLP to be a functional pVHL homolog potentially possessing multiple functions beyond HIF-1α-dependent binding activity. Our findings show that VLP may share at least seven interactors with pVHL, suggesting novel functional roles for this understudied human protein. These may occur at precise hypoxia levels where functional overlap with pVHL may permit a finer modulation of pVHL functions.

Helicobacter pylori (HP) is a Gram-negative bacterium that chronically infects the stomach of more than 50% of human population and represents a major cause of gastric cancer, gastric lymphoma, gastric autoimmunity, and peptic ulcer. It still remains to be elucidated, which HP virulence factors are important in the development of gastric disorders. Here, we analysed the role of the HP protein HP1454 in the host-pathogen interaction. We found that a significant proportion of T cells isolated from HP patients with chronic gastritis and gastric adenocarcinoma proliferated in response to HP1454. Moreover, we demonstrated in vivo that HP1454 protein drives Th1/Th17 inflammatory responses. We further analysed the in vitro response of human T cells exposed either to an HP wild-type strain or to a strain with a deletion of the hp1454 gene, and we revealed that HP1454 triggers the T-cell antigen receptor-dependent signalling and lymphocyte proliferation, as well as the CXCL12-dependent cell adhesion and migration. Our study findings prove that HP1454 is a crucial bacterial factor that exerts its proinflammatory activity by directly modulating the T-cell response. The relevance of these results can be appreciated by considering that compelling evidence suggest that chronic gastric inflammation, a condition that paves the way to HP-associated diseases, is dependent on T cells.

The mitochondrial F-ATP synthase is a complex molecular motor arranged in V-shaped dimers that is responsible for most cellular ATP synthesis in aerobic conditions. In the yeast F-ATP synthase, subunits e and g of the F(O) sector constitute a lateral domain, which is required for dimer stability and cristae formation. Here, by using site-directed mutagenesis, we identified Arg-8 of subunit e as a critical residue in mediating interactions between subunits e and g, most likely through an interaction with Glu-83 of subunit g. Consistent with this hypothesis, (i) the substitution of Arg-8 in subunit e (eArg-8) with Ala or Glu or of Glu-83 in subunit g (gGlu-83) with Ala or Lys destabilized the digitonin-extracted F-ATP synthase, resulting in decreased dimer formation as revealed by blue-native electrophoresis; and (ii) simultaneous substitution of eArg-8 with Glu and of gGlu-83 with Lys rescued digitonin-stable F-ATP synthase dimers. When tested in lipid bilayers for generation of Ca(2+)-dependent channels, WT dimers displayed the high-conductance channel activity expected for the mitochondrial megachannel/permeability transition pore, whereas dimers obtained at low digitonin concentrations from the Arg-8 variants displayed currents of strikingly small conductance. Remarkably, double replacement of eArg-8 with Glu and of gGlu-83 with Lys restored high-conductance channels indistinguishable from those seen in WT enzymes. These findings suggest that the interaction of subunit e with subunit g is important for generation of the full-conductance megachannel from F-ATP synthase.

Familiar cancers represent a privileged point of view for studying the complex cellular events inducing tumor transformation. Von Hippel-Lindau syndrome, a familiar predisposition to develop cancer is a clear example. Here, we present our efforts to decipher the role of von Hippel-Lindau tumor suppressor protein (pVHL) in cancer insurgence. We collected high quality information about both pVHL mutations and interactors to investigate the association between patient phenotypes, mutated protein surface and impaired interactions. Our data suggest that different phenotypes correlate with localized perturbations of the pVHL structure, with specific cell functions associated to different protein surfaces. We propose five different pVHL interfaces to be selectively involved in modulating proteins regulating gene expression, protein homeostasis as well as to address extracellular matrix (ECM) and ciliogenesis associated functions. These data were used to drive molecular docking of pVHL with its interactors and guide Petri net simulations of the most promising alterations. We predict that disruption of pVHL association with certain interactors can trigger tumor transformation, inducing metabolism imbalance and ECM remodeling. Collectively taken, our findings provide novel insights into VHL-associated tumorigenesis. This highly integrated in silico approach may help elucidate novel treatment paradigms for VHL disease.Author summaryCancer is generally caused by a series of mutations accumulating over time in a healthy tissue, which becomes re-programmed to proliferate at the expense of the hosting organism. This process is difficult to follow and understand as events in a multitude of different genes can lead to similar outcomes without apparent cause. The von Hippel-Lindau (VHL) tumor suppressor is one of the few genes harboring a familiar cancer syndrome, i.e. VHL mutations are known to cause a predictable series of events leading cancer in the kidneys and a few selected other tissues. This article describes a large-scale analysis to relate known VHL mutations to specific cancer pathways by looking at the molecular interactions. Different cancer types appear to be caused by mutations changing the surface of specific parts of the VHL protein. By looking at the VHL interactors involved, it is therefore possible to identify other candidate genes for mutations leading to very similar cancer types.

Residue interaction networks (RINs) are an alternative way of representing protein structures where nodes are residues and arcs physico–chemical interactions. RINs have been extensively and successfully used for analysing mutation effects, protein folding, domain–domain communication and catalytic activity. Here we present RING 2.0, a new version of the RING software for the identification of covalent and non-covalent bonds in protein structures, including π–π stacking and π–cation interactions. RING 2.0 is extremely fast and generates both intra and inter-chain interactions including solvent and ligand atoms. The generated networks are very accurate and reliable thanks to a complex empirical re-parameterization of distance thresholds performed on the entire Protein Data Bank. By default, RING output is generated with optimal parameters but the web server provides an exhaustive interface to customize the calculation. The network can be visualized directly in the browser or in Cytoscape. Alternatively, the RING-Viz script for Pymol allows visualizing the interactions at atomic level in the structure. The web server and RING-Viz, together with an extensive help and tutorial, are available from URL: http://protein.bio.unipd.it/ring.

Motivation The behavior of a protein is encoded in its sequence, which can be used to predict distinct features such as secondary structure, intrinsic disorder or amphipathicity. Integrating these and other features can help explain the context-dependent behavior of proteins. However, most tools focus on a single aspect, hampering a holistic understanding of protein structure. Here, we present Fast Estimator of Latent Local Structure (FELLS) to visualize structural features from the protein sequence. FELLS provides disorder, aggregation and low complexity predictions as well as estimated local propensities including amphipathicity. A novel fast estimator of secondary structure (FESS) is also trained to provide a fast response. The calculations required for FELLS are extremely fast and suited for large-scale analysis while providing a detailed analysis of difficult cases. Availability and Implementation The FELLS web server is available from URL: http://protein.bio.unipd.it/fells/. The server also exposes RESTful functionality allowing programmatic prediction requests. An executable version of FESS for Linux can be downloaded from URL: protein.bio.unipd.it/download/. Supplementary information Supplementary data are available at Bioinformatics online.

Post-translational modification (PTM) sites have become popular for predictor development. However, with the exception of phosphorylation and a handful of other examples, PTMs suffer from a limited number of available training examples and sparsity in protein sequences. Here, proline hydroxylation is taken as an example to compare different methods and evaluate their performance on new experimentally determined sites. As a guide for effective experimental design, predictors require both high specificity and sensitivity. However, the self-reported performance may often not be indicative of prediction quality and detection of new sites is not guaranteed. We have benchmarked seven published hydroxylation site predictors on two newly constructed independent datasets. The self-reported performance is found to widely overestimate the real accuracy measured on independent datasets. No predictor performs better than random on new examples, indicating the refined models do not sufficiently generalize to detect new sites. The number of false positives is high and precision low, in particular for non-collagen proteins whose motifs are not conserved. As hydroxylation site predictors do not generalize for new data, caution is advised when using PTM predictors in the absence of independent evaluations, in particular for highly specific sites involved in signalling., Author summary Machine learning methods are extensively used by biologists to design and interpret experiments. Predictors which take the only sequence as input are of particular interest due to the large amount of available sequence data and high self-reported performance. In this work, we evaluated post-translational modification (PTM) predictors for hydroxylation sites and found that they perform no better than random, in strong contrast to performances reported in their original publications. PTMs are chemical amino acid alterations providing the cell with conditional mechanisms to fine tune protein function, regulating complex biological processes such as signalling and cell cycle. Hydroxylation sites are a good PTM test case due to the availability of a range of predictors and an abundance of newly experimentally detected modification sites. Poor performances in our results highlight the overlooked problem of predicting PTMs when best practices are not followed and training data are likely incomplete. Experimentalists should be careful when using PTM predictors blindly and more independent assessments are needed to establish their usefulness in practice.

Background: Arrhythmogenic cardiomyopathy (ACM) is an inherited cardiac disease characterized by progressive fibro-fatty myocardial replacement, ventricular arrhythmia, heart failure, and sudden death. Causative mutations can be identified in 60% of patients, and most of them are found in genes encoding mechanical junction proteins of the intercalated disk. Methods: Whole-exome sequencing was performed on the proband of an ACM family. Sanger sequencing was used to screen for mutations the tight junction protein 1 ( TJP1 ) gene in unrelated patients. Predictions of local structure content and molecular dynamics simulations were performed to investigate the structural impact of the variants. Results: A novel c.2006A>G p.(Y669C) variant in TJP1 gene was identified by whole-exome sequencing in a patient with ACM. TJP1 encodes zonula occludens 1, an intercalated disk protein interacting with proteins of gap junctions and area composita. Additional rare TJP1 variants have been identified in 1 of 40 Italian probands (c.793C>T p.(R265W)) with arrhythmogenic right ventricular cardiomyopathy and in 2 of 43 Dutch/German patients (c. 986C>T, p.(S329L) and c.1079A>T, p.(D360V)) with dilated cardiomyopathy and recurrent ventricular tachycardia. The p.(D360V) variant was identified in a proband also carrying the p.(I156N) pathogenic variant in DSP . All 4 TJP1 variants are predicted to be deleterious and affect highly conserved amino acids, either at the GUK (guanylate kinase)–like domain (p.(Y669C)) or at the disordered region of the protein between the PDZ2 and PDZ3 domains (p.(R265W), p.(S329L), and p.(D360V)). The local unfolding induced by the former promotes structural rearrangements of the GUK domain, whereas the others are predicted to impair the function of the disordered region. Furthermore, rare variants in TJP1 are statistically enriched in patients with ACM relative to controls. Conclusions: We provide here the first evidence linking likely pathogenic TJP1 variants to ACM. Prevalence and pathogenic mechanism of TJP1 -mediated ACM remain to be determined.

We present an in silico characterization of the von Hippel-Lindau like protein (VLP), the only known human paralogue of the von Hippel-Lindau tumor suppressor protein (pVHL). Phylogenetic investigation showed VLP to be mostly conserved in upper mammals and specifically expressed in brain and testis. Structural analysis and molecular dynamics simulation show VLP to be very similar to pVHL three-dimensional organization and binding dynamics. In particular, conservation of elements at the protein interfaces suggests VLP to be a functional pVHL homolog potentially possessing multiple functions beyond HIF-1α-dependent binding activity. Our findings show that VLP may share at least seven interactors with pVHL, suggesting novel functional roles for this under-studied human protein. These may occur at precise hypoxia levels where functional overlap with pVHL may permit a finer modulation of pVHL functions.

Light is the main environmental stimulus that synchronizes the endogenous timekeeping systems in most terrestrial organisms. Drosophila cryptochrome (dCRY) is a light-responsive flavoprotein that detects changes in light intensity and wavelength around dawn and dusk. We have previously shown that dCRY acts through Inactivation No Afterpotential D (INAD) in a light-dependent manner on the Signalplex, a multiprotein complex that includes visual-signaling molecules, suggesting a role for dCRY in fly vision. Here, we predict and demonstrate a novel Ca2+-dependent interaction between dCRY and calmodulin (CaM). Through yeast two hybrid, coimmunoprecipitation (Co-IP), nuclear magnetic resonance (NMR) and calorimetric analyses we were able to identify and characterize a CaM binding motif in the dCRY C-terminus. Similarly, we also detailed the CaM binding site of the scaffold protein INAD and demonstrated that CaM bridges dCRY and INAD to form a ternary complex in vivo. Our results suggest a process whereby a rapid dCRY light response stimulates an interaction with INAD, which can be further consolidated by a novel mechanism regulated by CaM.

BACKGROUND Arrhythmogenic cardiomyopathy (ACM) is associated with arrhythmias and risk of sudden death. Mutations in genes encoding proteins of cardiac intercalated discs account for similar to 60% of ACM cases, but the remaining 40% is still genetically elusive.OBJECTIVE The purpose of this study was to identify the underlying genetic cause in probands with ACM.METHODS DNA samples from 40 probands with ACM, negative for mutations in the 3 major ACM genes-DSP, PKP2, and DSG2, were screened by using a targeted gene panel consisting of 15 known ACM genes and 53 candidate genes.RESULTS About half of patients were found to carry rare variant(s) predicted to be damaging; specifically, 9 (22.5%) showed >= 1 variants in genes associated with ACM and/or with other inherited heart diseases and 10 (25%) showed variants in candidate genes. Among the latter, we focused on 2 novel variants in TP63 and PPP1R13L candidate genes (c.796C>T, p.(R266*) and c.1858G>C, p.(A620P), respectively). The encoded proteins p63 and inhibitor of apoptosis stimulating p53 protein are known to be interacting partners. Inhibitor of apoptosis stimulating p53 protein is a shuttling multifunctional protein: in the nucleus it is critical for inhibiting p63 function, whereas in the cytoplasm it regulates desmosome integrity. According to the American College of Medical Genetics and Genomics guidelines, the variant in TP63 has been scored as likely pathogenic and the variant in PPP1R13L as a variant of uncertain significance. Importantly, the mutant TP63 allele leads to nonsense-mediated messenger RNA decay, causing haploinsufficiency.CONCLUSION Our findings identify TP63 as a putative novel disease gene for ACM, while the possible involvement of PPP1R13L remains to be determined.

Modification with arginine-specific glyoxals modulates the permeability transition (PT) of rat liver mitochondria, with inhibitory or inducing effects that depend on the net charge of the adduct(s). Here, we show that phenylglyoxal (PGO) affects the PT in a species-specific manner (inhibition in mouse and yeast, induction in human and Drosophila mitochondria). Following the hypotheses (i) that the effects are mediated by conserved arginine(s) and (ii) that the PT is mediated by the F-ATP synthase, we have narrowed the search to 60 arginines. Most of these residues are located in subunits alpha, beta, gamma, epsilon, a, and c and were excluded because PGO modification did not significantly affect enzyme catalysis. On the other hand, yeast mitochondria lacking subunit g or bearing a subunit g R107A mutation were totally resistant to PT inhibition by PGO. Thus, the effect of PGO on the PT is specifically mediated by Arg-107, the only subunit g arginine that has been conserved across species. These findings are evidence that the PT is mediated by F-ATP synthase.

PHDs (proline hydroxylases) are a small protein family found in all organisms, considered the central regulator of the molecular hypoxia response due to PHDs being completely inactivated under low oxygen concentration. At physiological oxygen concentration, PHDs drive the degradation of the HIF-1α (hypoxia-inducible factor 1-α), which is responsible for upregulating the expression of genes involved in the cellular response to hypoxia. Hypoxia is a common feature of most tumors, in particular during metastasis development. Indeed, cancer reacts by activating pathways promoting new blood vessel formation and activating strategies aimed to improve survival. In this scenario, the PHD family regulates the activation of HIF-1α and cell-cycle regulation. Several PHD mutations were found in cancer patients, underlining their importance for human health. Here, we propose a Bayesian model able to predict the pathological effect of human PHD mutations and their correlation with cancer outcome. The model was developed through an integrative in silico approach, where data collected from the literature has been coupled with sequence evolution and structural analysis. The model was used to assess 135 human PHD variants. Finally, bioinformatics characterization was used to demonstrate how few amino acid changes are able to explain the functional specialization of PHD family members and their physiological role in human health.

Correct phenotypic interpretation of variants of unknown significance for cancer-associated genes is a diagnostic challenge as genetic screenings gain in popularity in the next-generation sequencing era. The Critical Assessment of Genome Interpretation (CAGI) experiment aims to test and define the state of the art of genotype-phenotype interpretation. Here, we present the assessment of the CAGI p16INK4a challenge. Participants were asked to predict the effect on cellular proliferation of ten variants for the p16INK4a tumor suppressor, a cyclin-dependent kinase inhibitor encoded by the CDKN2A gene. Twenty-two pathogenicity predictors were assessed with a variety of accuracy measures for reliability in a medical context. Different assessment measures were combined in an overall ranking to provide more robust results. The R scripts used for assessment are publicly available from a GitHub repository for future use in similar assessment exercises. Despite a limited test-set size, our findings show a variety of results, with some methods performing significantly better. Methods combining different strategies frequently outperform simpler approaches. The best predictor, Yang&Zhou lab, uses a machine learning method combining an empirical energy function measuring protein stability with an evolutionary conservation term. The p16INK4a challenge highlights how subtle structural effects can neutralize otherwise deleterious variants. This article is protected by copyright. All rights reserved.

F‐ATP synthases convert the electrochemical energy of the H+ gradient into the chemical energy of ATP with remarkable efficiency. Mitochondrial F‐ATP synthases can also undergo a Ca2+‐dependent transformation to form channels with properties matching those of the permeability transition pore (PTP), a key player in cell death. The Ca2+ binding site and the mechanism(s) through which Ca2+ can transform the energy‐conserving enzyme into a dissipative structure promoting cell death remain unknown. Through in vitro, in vivo and in silico studies we (i) pinpoint the “Ca2+‐trigger site” of the PTP to the catalytic site of the F‐ATP synthase β subunit and (ii) define a conformational change that propagates from the catalytic site through OSCP and the lateral stalk to the inner membrane. T163S mutants of the β subunit, which show a selective decrease in Ca2+‐ATP hydrolysis, confer resistance to Ca2+‐induced, PTP‐dependent death in cells and developing zebrafish embryos. These findings are a major advance in the molecular definition of the transition of F‐ATP synthase to a channel and of its role in cell death.

Germline inactivation of the von Hippel-Lindau (VHL) tumor suppressor predisposes patients to develop different highly vascularized cancers. pVHL targets the hypoxia-inducible transcription factor (HIF-1α) for degradation, modulating the activation of various genes involved in hypoxia response. Hypoxia plays a relevant role in regulating cell cycle progression, inducing growth arrest in cells exposed to prolonged oxygen deprivation. However, the exact molecular details driving this transition are far from understood. Here, we present novel interactions between pVHL and the cyclin-dependent kinase inhibitor family CDKN1 (p21, p27 and p57). Bioinformatics analysis, yeast two-hybrid screening and co-immunoprecipitation assays were used to predict, dissect and validate the interactions. We found that the CDKN1 proteins share a conserved region mimicking the HIF-1α motif responsible for pVHL binding. Intriguingly, a p27 site-specific mutation associated to cancer is shown to modulate this novel interaction. Our findings suggest a new connection between the pathways regulating hypoxia and cell cycle progression.

Photoreceptors are crucial components for circadian rhythm entrainment in animals, plants, fungi and cyanobacteria. Cryptochromes (CRYs) are flavin adenine dinucleotide (FAD) containing photoreceptors, and FAD is responsible for signal transduction, in contrast to photolyases where it promotes DNA-damage repair. In this work, we investigated an alternative role for FAD in CRY. We analyzed the Drosophila melanogaster CRY crystal structure by means of molecular dynamics, elucidating how this large co-factor within the receptor could be crucial for CRY structural stability. The co-factor appears indeed to improve receptor motility, providing steric hindrance. Moreover, multiple sequence alignments revealed that conserved motifs in the C-terminal tail could be necessary for functional stability.

RepeatsDB (http://repeatsdb.bio.unipd.it/) is a database of annotated tandem repeat protein structures. Tandem repeats pose a difficult problem for the analysis of protein structures, as the underlying sequence can be highly degenerate. Several repeat types haven been studied over the years, but their annotation was done in a case-by-case basis, thus making large-scale analysis difficult. We developed RepeatsDB to fill this gap. Using state-of-the-art repeat detection methods and manual curation, we systematically annotated the Protein Data Bank, predicting 10 745 repeat structures. In all, 2797 structures were classified according to a recently proposed classification schema, which was expanded to accommodate new findings. In addition, detailed annotations were performed in a subset of 321 proteins. These annotations feature information on start and end positions for the repeat regions and units. RepeatsDB is an ongoing effort to systematically classify and annotate structural protein repeats in a consistent way. It provides users with the possibility to access and download high-quality datasets either interactively or programmatically through web services.

Hypoxia inducible factor 1α (HIF-1α) regulates oxygen homeostasis in the cell through a sensing mechanism involving its hydroxylation and binding to the von Hippel–Lindau (VHL) tumor suppressor. This mechanism is mediated through hydroxylation of HIF-1α proline 564, although in vitro tests have previously shown an alternative hydroxylation at proline 567 by PHD-3. Here, molecular dynamics simulations were used to investigate the structural effect of this alternative hydroxylation. A specific hydrogen bond network rearrangement and improved electrostatic energy for hydroxylated P567 are compatible with an increase in HIF-1α binding affinity. Sequence analysis also confirms P567 to be vastly conserved during evolution, indicating a possible role for this alternative, PHD-3 driven, post translational modification in pVHL–HIF-1α complex formation.

Mutations in von Hippel-Lindau tumor suppressor protein (pVHL) predispose to develop tumors affecting specific target organs, such as the retina, epididymis, adrenal glands, pancreas and kidneys. Currently, more than 400 pVHL interacting proteins are either described in the literature or predicted in public databases. This data is scattered among several different sources, slowing down the comprehension of pVHL’s biological role. Here we present VHLdb, a novel database collecting available interaction and mutation data on pVHL to provide novel integrated annotations. In VHLdb, pVHL interactors are organized according to two annotation levels, manual and automatic. Mutation data are easily accessible and a novel visualization tool has been implemented. A user-friendly feedback function to improve database content through community-driven curation is also provided. VHLdb presently contains 478 interactors, of which 117 have been manually curated, and 1,074 mutations. This makes it the largest available database for pVHL-related information. VHLdb is available from URL: http://vhldb.bio.unipd.it/.

Human pendrin (SLC26A4) is an anion transporter mostly expressed in the inner ear, thyroid and kidney. SLC26A4 gene mutations are associated with a broad phenotypic spectrum, including Pendred Syndrome and non-syndromic hearing loss with enlarged vestibular aqueduct (ns-EVA). No experimental structure of pendrin is currently available, making phenotype-genotype correlations difficult as predictions of transmembrane (TM) segments vary in number. Here, we propose a novel three-dimensional (3D) pendrin transmembrane domain model based on the SLC26Dg transporter. The resulting 14 TM topology was found to include two non-canonical transmembrane segments crucial for pendrin activity. Mutation mapping of 147 clinically validated pathological mutations shows that most affect two previously undescribed TM regions.

The Database of Protein Disorder (DisProt, URL: www.disprot.org) has been significantly updated and upgraded since its last major renewal in 2007. The current release holds information on more than 800 entries of IDPs/IDRs, i.e. intrinsically disordered proteins or regions that exist and function without a well-defined three-dimensional structure. We have re-curated previous entries to purge DisProt from conflicting cases, and also upgraded the functional classification scheme to reflect continuous advance in the field in the past 10 years or so. We define IDPs as proteins that are disordered along their entire sequence, i.e. entirely lack structural elements, and IDRs as regions that are at least five consecutive residues without well-defined structure. We base our assessment of disorder strictly on experimental evidence, such as X-ray crystallography and nuclear magnetic resonance (primary techniques) and a broad range of other experimental approaches (secondary techniques). Confident and ambiguous annotations are highlighted separately. DisProt 7.0 presents classified knowledge regarding the experimental characterization and functional annotations of IDPs/IDRs, and is intended to provide an invaluable resource for the research community for a better understanding structural disorder and for developing better computational tools for studying disordered proteins., SCOPUS: ar.j, info:eu-repo/semantics/published

Idiopathic erythrocytosis is a rare disease characterized by an increase in red blood cell mass due to mutations in proteins of the oxygen-sensing pathway, such as prolyl hydroxylase 2 (PHD2). Here, we present a bioinformatics investigation of the pathological effect of twelve PHD2 mutations related to polycythemia insurgence. We show that few mutations impair the PHD2 catalytic site, while most localize to non-enzymatic regions. We also found that most mutations do not overlap the substrate recognition site, suggesting a novel PHD2 binding interface. After a structural analysis of both binding partners, we suggest that this novel interface is responsible for PHD2 interaction with the LIMD1 tumor suppressor.

Graphical abstractDisplay Omitted HighlightsA novel design optimization technique is proposed based on the ant colony optimization (ACO) framework.The proposed approach combines the properties and peculiarities of ACO and Nave Bayes classifier.The proposed approach is applied to a problem concerning enzyme engineering.The simulation and experimental results demonstrate a remarkable scalability of the proposed approach to the overall size of the search space. In a large number of experimental problems, high dimensionality of the search area and economical constraints can severely limit the number of experimental points that can be tested. Within these constraints, classical optimization techniques perform poorly, in particular, when little a priori knowledge is available. In this work we investigate the possibility of combining approaches from statistical modeling and bio-inspired algorithms to effectively explore a huge search space, sampling only a limited number of experimental points. To this purpose, we introduce a novel approach, combining ant colony optimization (ACO) and nave Bayes classifier (NBC) that is, the nave Bayes ant colony optimization (NACO) procedure. We compare NACO with other similar approaches developing a simulation study. We then derive the NACO procedure with the goal to design artificial enzymes with no sequence homology to the extant one. Our final aim is to mimic the natural fold of 200 amino acids 1AGY serine esterase from Fusarium solani.