Sequence and Structure-Based Analyses of Human Ankyrin Repeats

Ankyrin is one of the most abundant protein repeat families found across all forms of life. It is found in a variety of multi-domain and single domain proteins in humans with diverse number of repeating units. They are observed to occur in several functionally diverse proteins, such as transcriptional initiators, cell cycle regulators, cytoskeletal organizers, ion transporters, signal transducers, developmental regulators, and toxins, and, consequently, defects in ankyrin repeat proteins have been associated with a number of human diseases. In this study, we have classified the human ankyrin proteins into clusters based on the sequence similarity in their ankyrin repeat domains. We analyzed the amino acid compositional bias and consensus ankyrin motif sequence of the clusters to understand the diversity of the human ankyrin proteins. We carried out network-based structural analysis of human ankyrin proteins across different clusters and showed the association of conserved residues with topologically important residues identified by network centrality measures. The analysis of conserved and structurally important residues helps in understanding their role in structural stability and function of these proteins. In this paper, we also discuss the significance of these conserved residues in disease association across the human ankyrin protein clusters.

Mapping the Substrate-Binding Sites in the Phosphatidylserine Synthase in Candida albicans

The fungal phosphatidylserine (PS) synthase, a membrane protein encoded by the CHO1 gene, is a potential drug target for pathogenic fungi, such as Candida albicans. However, both substrate-binding sites of C. albicans Cho1 have not been characterized. Cho1 has two substrates: cytidyldiphosphate-diacylglycerol (CDP-DAG) and serine. Previous studies identified a conserved CDP-alcohol phosphotransferase (CAPT) binding motif, which is present within Cho1. We tested the CAPT motif for its role in PS synthesis by mutating conserved residues using alanine substitution mutagenesis. PS synthase assays revealed that mutations in all but one conserved amino acid within the CAPT motif resulted in decreased Cho1 function. In contrast, there were no clear motifs in Cho1 for binding serine. Therefore, to identify the serine binding site, PS synthase sequences from three fungi were aligned with sequences of a similar enzyme, phosphatidylinositol (PI) synthase, from the same fungi. This revealed a motif that was unique to PS synthases. Using alanine substitution mutagenesis, we found that some of the residues in this motif are required for Cho1 function. Two alanine substitution mutants, L184A and R189A, exhibited contrasting impacts on PS synthase activity, and were characterized for their Michaelis-Menten kinetics. The L184A mutant displayed enhanced PS synthase activity and showed an increased Vmax. In contrast, R189A showed decreased PS synthase activity and increased Km for serine, suggesting that residue R189 is involved in serine binding. These results help to characterize PS synthase substrate binding, and should direct rational approaches for finding Cho1 inhibitors that may lead to better antifungals.

Finding Aquaporins in Annelids: An Evolutionary Analysis and a Case Study

Aquaporins (AQPs) are a family of membrane channels facilitating diffusion of water and small solutes into and out of cells. Despite their biological relevance in osmoregulation and ubiquitous distribution throughout metazoans, the presence of AQPs in annelids has been poorly investigated. Here, we searched and annotated Aqp sequences in public genomes and transcriptomes of annelids, inferred their evolutionary relationships through phylogenetic analyses and discussed their putative physiological relevance. We identified a total of 401 Aqp sequences in 27 annelid species, including 367 sequences previously unrecognized as Aqps. Similar to vertebrates, phylogenetic tree reconstructions clustered these annelid Aqps in four clades: AQP1-like, AQP3-like, AQP8-like and AQP11-like. We found no clear indication of the existence of paralogs exclusive to annelids; however, several gene duplications seem to have occurred in the ancestors of some Sedentaria annelid families, mainly in the AQP1-like clade. Three of the six Aqps annotated in Alitta succinea, an estuarine annelid showing high salinity tolerance, were validated by RT-PCR sequencing, and their similarity to human AQPs was investigated at the level of “key” conserved residues and predicted three-dimensional structure. Our results suggest a diversification of the structures and functions of AQPs in Annelida comparable to that observed in other taxa.

Ion Transport, Selectivity, and Electronic Polarization in Fluoride Channels

Fluoride channels (Fluc) export toxic F- from the cytoplasm. Crystallography and mutagenesis have identified several conserved residues crucial for fluoride transport, but the transport mechanism at the molecular level has remained elusive. Herein we have applied constant-pH molecular dynamics and free energy sampling methods to investigate fluoride transfer through a Fluc protein from Escherichia coli. We find that fluoride is facile to transfer in its charged form, i.e., F-, by traversing through a non-bonded network. The extraordinary F- selectivity is gained by the hydrogen-bonding capability of the central binding site and the Coulombic filter at the channel entrance. The F- transfer rate calculated using an electronically polarizable force field is significantly more accurate compared to the experimental value than that calculated using a more standard additive force field, suggesting an essential role for electronic polarization in the F- - Fluc interactions.

Characterizing disease-associated human proteins without available protein structures or homologues

Mutations in human proteins lead to diseases. The structure of these proteins can help understand the mechanism of such diseases and develop therapeutics against them. With improved deep learning techniques such as RoseTTAFold and AlphaFold, we can predict the structure of these proteins even in the absence of structural homologues. We modeled and extracted the domains from 553 disease-associated human proteins. We noticed that the model quality was higher and the RMSD lower between AlphaFold and RoseTTAFold models for domains that could be assigned to CATH families as compared to those which could be assigned to Pfam families of unknown structure or could not be assigned to either. We predicted ligand-binding sites, protein-protein interfaces, conserved residues and destabilising effects caused by residue mutations in these predicted structures. We then explored whether the disease-associated mutations were in the proximity of these predicted functional sites or if they destabilized the protein structure based on ddG calculations. We could explain 80% of these disease-associated mutations based on proximity to functional sites or structural destabilization. Usage of models from the two state-of-the-art techniques provide better confidence in our predictions, and we explain 93 additional mutations based on RoseTTAFold models which could not be explained based solely on AlphaFold models.

A Study on the Genetics of Primary Ciliary Dyskinesia

Primary ciliary dyskinesia (PCD) is a poorly understood disorder. It is primarily autosomal recessive and is prevalent in tribal communities of the United Arab Emirates due to consanguineous marriages. This retrospective study aimed to assess the pathogenicity of the genetic variants of PCD in indigenous patients with significant clinical respiratory problems. Pathogenicity scores of variants obtained from the chart review were consolidated using the Ensembl Variant Effect Predictor. The multidimensional dataset of scores was clustered into three groups based on their pathogenicity. Sequence alignment and the Jensen–Shannon Divergence (JSD) were generated to evaluate the amino acid conservation at the site of the variation. One-hundred and twelve variants of 28 genes linked to PCD were identified in 66 patients. Twenty-two variants were double heterozygous, two triple heterozygous, and seven homozygous. Of the thirteen novel variants, two, c.11839 + 1G > A in dynein, axonemal, heavy chain 11 (DNAH11) and p.Lys92Trpfs in dynein, axonemal, intermediate chain 1 (DNAI1) were associated with dextrocardia with situs inversus, and one, p.Gly21Val in coiled-coil domain-containing protein 40 (CCDC40), with absent inner dynein arms. Homozygous C1orf127:p.Arg113Ter (rs558323413) was also associated with laterality defects in two related patients. The majority of variants were missense involving conserved residues with a median JSD score of 0.747. Homology models of two deleterious variants in the stalk of DNAH11, p.Gly3102Asp and p.Leu3127Arg, revealed structural importance of the conserved glycine and leucine. These results define potentially damaging PCD variants in the region. Future studies, however, are needed to fully comprehend the genetic underpinnings of PCD.

TwinCons: Conservation score for uncovering deep sequence similarity and divergence

We have developed the program TwinCons, to detect noisy signals of deep ancestry of proteins or nucleic acids. As input, the program uses a composite alignment containing pre-defined groups, and mathematically determines a ‘cost’ of transforming one group to the other at each position of the alignment. The output distinguishes conserved, variable and signature positions. A signature is conserved within groups but differs between groups. The method automatically detects continuous characteristic stretches (segments) within alignments. TwinCons provides a convenient representation of conserved, variable and signature positions as a single score, enabling the structural mapping and visualization of these characteristics. Structure is more conserved than sequence. TwinCons highlights alternative sequences of conserved structures. Using TwinCons, we detected highly similar segments between proteins from the translation and transcription systems. TwinCons detects conserved residues within regions of high functional importance for the ribosomal RNA (rRNA) and demonstrates that signatures are not confined to specific regions but are distributed across the rRNA structure. The ability to evaluate both nucleic acid and protein alignments allows TwinCons to be used in combined sequence and structural analysis of signatures and conservation in rRNA and in ribosomal proteins (rProteins). TwinCons detects a strong sequence conservation signal between bacterial and archaeal rProteins related by circular permutation. This conserved sequence is structurally colocalized with conserved rRNA, indicated by TwinCons scores of rRNA alignments of bacterial and archaeal groups. This combined analysis revealed deep co-evolution of rRNA and rProtein buried within the deepest branching points in the tree of life.

An Engineered Thermostable Laccase with Great Ability to Decolorize and Detoxify Malachite Green

Laccases can catalyze the remediation of hazardous synthetic dyes in an eco-friendly manner, and thermostable laccases are advantageous to treat high-temperature dyeing wastewater. A novel laccase from Geothermobacter hydrogeniphilus (Ghlac)was cloned and expressed in Escherichia coli. Ghlac containing 263 residues was characterized as a functional laccase of the DUF152 family. By structural and biochemical analyses, the conserved residues H78, C119, and H136 were identified to bind with one copper atom to fulfill the laccase activity. In order to make it more suitable for industrial use, Ghlac variant Mut2 with enhanced thermostability was designed. The half-lives of Mut2 at 50 °C and 60 °C were 80.6 h and 9.8 h, respectively. Mut2 was stable at pH values ranging from 4.0 to 8.0 and showed a high tolerance for organic solvents such as ethanol, acetone, and dimethyl sulfoxide. In addition, Mut2 decolorized approximately 100% of 100 mg/L of malachite green dye in 3 h at 70 °C. Furthermore, Mut2 eliminated the toxicity of malachite green to bacteria and Zea mays. In summary, the thermostable laccase Ghlac Mut2 could effectively decolorize and detoxify malachite green at high temperatures, showing great potential to remediate the dyeing wastewater.

Effect of metal ions and chemical agents on the activity of endoglucanase GH5 exploited from goats-rumen bacterial metagenomic DNA data

A gene coding for GH5 endoglucanase exploited from metagnomic DNA data of bacteria in Vietnamese goats’ rumen was modularity structure including a catalytic module, a fibronectin-3 like module and an X module. The recombinant enzyme was sucessfully expressed in E. coli and purified. To study the effect of some metal ions and chemicals on enzyme activity, in this study, we used some tools including Swiss-Prot, ProFunc, COFACTOR for prediction of enzyme structure and ligands interaction. The obtained results indicated that the most similar structure with enzyme had two conserved residues (Asp-190 và Asp-192) linked with Mn2+ within a radius of ~ 3.5 Å from the center of ion Mn2+ and enzyme molecule contained a disulphide bond. Experimental results for essessment of the effect of some metal ions (Ca2 +, Mn2 +, Mg2+, Ni2+, K+, Co2+, Cu2+, Zn2+, Fe3+) at the final concentration of 10 mM and of six common chemicals including SDS (1%), urea (1 µM), 2-mercaptoethanol (1 µM), EDTA (1 µM), tween 80 (1mM), triton X-100 (1 µM) showed that only Mn2+ increased enzyme activity slightly at concentration of 10 mM and two times at the concentration of 40 mM Mn2+. The Mn2+ has been identified as a specific binding agent may increase the stability and activity of endoglucanase GH5.