Structural modeling and gene expression analysis of phosvitinless vitellogenin (vgc) in the Indian freshwater murrel, Channa punctatus (Bloch, 1793)

Vitellogenin (Vg) is a female-specific egg-yolk precursor protein, synthesized in the liver of fish in response to estrogens. In the present study, complete gene of phosvitinless vitellogenin (vgc) was sequenced, its 3D structure was predicted and validated by web-based softwares. The complete nucleotide sequence of vgc was 4126 bp which encodes for 1272 amino acids and showed the presence of three conserved domains viz. LPD_N, DUF1943 and DUF1944. The retrieved amino acid sequence of VgC protein was subjected to in silico analysis for understanding the structural and functional properties of protein. mRNA levels of multiple vg genes have also been quantified during annual reproductive cycle employing qPCR. A correlation has been observed between seasonal changes in gonadosomatic index with estradiol levels and hepatic expression of three types of vg genes (vga, vgb, vgc) during ovarian cycle of murrel. During preparatory phase, when photoperiod and temperature are low; low titre of E2 in blood induces expression of vgc gene. A rapid increase in the levels of E2 favours induction of vgb and vga genes in liver of murrel during early pre-spawning phase when photoperiod is long and temperature is high in nature. These results suggest that among three vitellogenin proteins, VgC is synthesized earlier than VgA and VgB during oogenesis.

RNA 3D Modeling with FARFAR2, Online

Understanding the three-dimensional structure of an RNA molecule is often essential to understanding its function. Sampling algorithms and energy functions for RNA structure prediction are improving, due to the increasing diversity of structural data available for training statistical potentials and testing structural data, along with a steady supply of blind challenges through the RNA-Puzzles initiative. The recent FARFAR2 algorithm enables near-native structure predictions on fairly complex RNA structures, including automated selection of final candidate models and estimation of model accuracy. Here, we describe the use of a publicly available webserver for RNA modeling for realistic scenarios using FARFAR2, available at https://rosie.rosettacommons.org/farfar2 . We walk through two cases in some detail: a simple model pseudoknot from the frameshifting element of beet western yellows virus modeled using the "basic interface" to the webserver and a replication of RNA-Puzzle 20, a metagenomic twister sister ribozyme, using the "advanced interface." We also describe example runs of FARFAR2 modeling including two kinds of experimental data: a c-di-GMP riboswitch modeled with low-resolution restraints from MOHCA-seq experiments and a tandem GA motif modeled with 1H NMR chemical shifts.

A hybrid approach to study large conformational transitions of biomolecules from single particle XFEL diffraction data

X-ray free-electron laser (XFEL) is the latest generation of the X-ray source that could become an invaluable technique in structural biology. XFEL has ultrashort pulse duration, extreme peak brilliance, and high spatial coherence, which could enable the observation of the biological molecules in near nature state at room temperature without crystallization. However, for biological systems, due to their low diffraction power and complexity of sample delivery, experiments and data analysis are not straightforward, making it extremely challenging to reconstruct three-dimensional (3D) structures from single particle XFEL data. Given the current limitations to the amount and resolution of the data from such XFEL experiments, we propose a new hybrid approach for characterizing biomolecular conformational transitions by using a single 2D low-resolution XFEL diffraction pattern in combination with another known conformation. In our method, we represent the molecular structure with a coarse-grained model, the Gaussian mixture model, to describe large conformational transitions from low-resolution XFEL data. We obtain plausible 3D structural models that are consistent with the XFEL diffraction pattern by deforming an initial structural model to maximize the similarity between the target pattern and the simulated diffraction patterns from the candidate models. We tested the proposed algorithm on two biomolecules of different sizes with different complexities of conformational transitions, adenylate kinase, and elongation factor 2, using synthetic XFEL data. The results show that, with the proposed algorithm, we can successfully describe the conformational transitions by flexibly fitting the coarse-grained model of one conformation to become consistent with an XFEL diffraction pattern simulated from another conformation. In addition, we showed that the incident beam orientation has some effect on the accuracy of the 3D structure modeling and discussed the reasons for the inaccuracies for certain orientations. The proposed method could serve as an alternative approach for retrieving information on 3D conformational transitions from the XFEL diffraction patterns to interpret experimental data. Since the molecules are represented by Gaussian kernels and no atomic structure is needed in principle, such a method could also be used as a tool to seek initial models for 3D reconstruction algorithms.

Structural analysis of SARS-Cov-2 nonstructural protein 1 polymorphisms found in the Brazilian Amazon

The coronavirus disease COVID-19 has been the cause of millions of deaths worldwide. Among the SARS-CoV-2 proteins, the non-structural protein 1 (NSP1) has great importance during the virus infection process and is present in both alpha and beta-CoVs. Therefore, monitoring of NSP1 polymorphisms is crucial in order to understand their role during infection and virus-induced pathogenicity. Herein, we analyzed how mutations detected in the circulating SARS-CoV-2 in the population of the city of Manaus, Amazonas state, Brazil could modify the tertiary structure of the NSP1 protein. Three mutations were detected in the SARS-CoV-2 NSP1 gene: deletion of the amino acids KSF from positions 141 to 143 (delKSF), SARS-CoV-2, lineage B.1.195; and two substitutions, R29H and R43C, SARS-CoV-2 lineage B.1.1.28 and B.1.1.33, respectively. The delKSF was found in 47 samples, whereas R29H and R43C were found in two samples, one for each mutation. The NSP1 structures carrying the mutations R43C and R29H on the N-terminal portion (e.g. residues 10 to 127) showed minor backbone divergence compared to the Wuhan model. However, the NSP1 C-terminal region (residues 145 to 180) was severely affected in the delKSF and R29H mutants. The intermediate variable region (residues 144 to 148) leads to changes in the C-terminal region, particularly in the delKSF structure. New investigations must be carried out to analyze how these changes affect NSP1 activity during the infection. Our results reinforce the need for continuous genomic surveillance of SARS-CoV-2 to better understand virus evolution and assess the potential impact of the viral mutations on the approved vaccines and future therapies.

SRD5A3-CDG: 3D structure modeling, clinical spectrum, and computer-based dysmorphic facial recognition

Pathogenic variants in Steroid 5 alpha reductase type 3 (SRD5A3) cause rare inherited congenital disorder of glycosylation known as SRD5A3-CDG (MIM# 612379). To date, 43 affected individuals have been reported. Despite the development of various dysmorphic features in significant number of patients, facial recognition entity has not yet been established for SRD5A3-CDG. Herein, we reported a novel SRD5A3 missense pathogenic variant c.460 T > C p.(Ser154Pro). The 3D structural modeling of the SRD5A3 protein revealed additional transmembrane α-helices and predicted that the p.(Ser154Pro) variant is located in a potential active site and is capable of reducing its catalytic efficiency. Based on phenotypes of our patients and all published SRD5A3-CDG cases, we identified the most common clinical features as well as some recurrent dysmorphic features such as arched eyebrows, wide eyes, shallow nasal bridge, short nose, and large mouth. Based on facial digital 2D images, we successfully designed and validated a SRD5A3-CDG computer based dysmorphic facial analysis, which achieved 92.5% accuracy. The current work integrates genotypic, 3D structural modeling and phenotypic characteristics of CDG-SRD5A3 cases with the successful development of computer tool for accurate facial recognition of CDG-SRD5A3 complex cases to assist in the diagnosis of this particular disorder globally.

Characterization of the Mycobacterial MSMEG-3762/63 Efflux Pump in Mycobacterium smegmatis Drug Efflux

Multi-drug resistant tuberculosis (MDR-TB) represents a major health problem worldwide. Drug efflux and the activity of efflux transporters likely play important roles in the development of drug-tolerant and drug-resistant mycobacterial phenotypes. This study is focused on the action of a mycobacterial efflux pump as a mechanism of drug resistance. Previous studies demonstrated up-regulation of the TetR-like transcriptional regulator MSMEG_3765 in Mycobacterium smegmatis and its ortholog Rv1685c in Mycobacterium tuberculosis (Mtb) in acid-nitrosative stress conditions. MSMEG-3765 regulates the expression of the MSMEG_3762/63/65 operon, and of the orthologous region in Mtb (Rv1687c/86c/85c). MSMEG-3762 and Rv1687c are annotated as ATP-binding proteins, while MSMEG-3763 and Rv1686c are annotated as trans-membrane polypeptides, defining an ABC efflux pump in both M. smegmatis and Mtb. The two putative efflux systems share a high percentage of identity. To examine the role of the putative efflux system MSMEG-3762/63, we constructed and characterized a MSMEG-3763 deletion mutant in M. smegmatis (∆MSMEG_3763). By comparative analysis of wild type, knockout, and complemented strains, together with structural modeling and molecular docking bioinformatics analyses of the MSMEG-3763 trans-membrane protein, we define the protein complex MSMEG-3762/63 as an efflux pump. Moreover, we demonstrate involvement of this pump in biofilm development and in the extrusion of rifampicin and ciprofloxacin (CIP), antimicrobial drugs used in first- and second-line anti-TB therapies.

Defensins of Grasses: A Systematic Review

The grass family (Poaceae) is one of the largest families of flowering plants, growing in all climatic zones of all continents, which includes species of exceptional economic importance. The high adaptability of grasses to adverse environmental factors implies the existence of efficient resistance mechanisms that involve the production of antimicrobial peptides (AMPs). Of plant AMPs, defensins represent one of the largest and best-studied families. Although wheat and barley seed γ-thionins were the first defensins isolated from plants, the functional characterization of grass defensins is still in its infancy. In this review, we summarize the current knowledge of the characterized defensins from cultivated and selected wild-growing grasses. For each species, isolation of defensins or production by heterologous expression, peptide structure, biological activity, and structure–function relationship are described, along with the gene expression data. We also provide our results on in silico mining of defensin-like sequences in the genomes of all described grass species and discuss their potential functions. The data presented will form the basis for elucidation of the mode of action of grass defensins and high adaptability of grasses to environmental stress and will provide novel potent molecules for practical use in medicine and agriculture.

The protein-binding N-terminal domain of human translation elongation factor 1Bβ possesses a dynamic α-helical structural organization

Translation elongation factor 1Bβ (eEF1Bβ) is a metazoan-specific protein involved into the macromolecular eEF1B complex, containing also eEF1Bα and eEF1Bγ subunits. Both eEF1Bα and eEF1Bβ ensure the guanine nucleotide exchange on eEF1A while eEF1Bγ is thought to have a structural role. The structures of the eEF1Bβ catalytic C-terminal domain and neighboring central acidic region are known while the structure of the protein-binding N-terminal domain remains unidentified which prevents clear understanding of architecture of the eEF1B complex. Here we show that the N-terminal domain comprising initial 77 amino acids of eEF1Bβ, eEF1Bβ(1-77), is a monomer in solution with increased hydrodynamic volume. This domain binds eEF1Bγ in equimolar ratio. The CD spectra reveal that the secondary structure of eEF1Bβ(1-77) consists predominantly of α-helices and a portion of disordered region. Very rapid hydrogen/deuterium exchange for all eEF1Bβ(1-77) peptides favors a flexible tertiary organization of eEF1Bβ(1-77). Computational modeling of eEF1Bβ(1-77) suggests several conformation states each composed of three α-helices connected by flexible linkers. Altogether, the data imply that the protein-binding domain of eEF1Bβ shows flexible spatial organization which may be needed for interaction with eEF1Bγ or other protein partners.

Cloning and molecular modeling of a thermostable carboxylesterase from the chicken uropygial glands

Starting from total uropygial glands mRNAs, chicken uropygial carboxylesterase (cuCES) cDNA was synthesized by RT-PCR and cloned into the PGEM-T vector. Amino acid sequence of the cuCES is compared to that of human liver carboxylesterase 1 (hCES1). Given the high amino acid sequence homology between the two enzymes, a 3-D structure model of the chicken carboxylesterase was built using the structure of hCES1 as template. By following this model and utilizing molecular dynamics (MD) simulations, the resistance of the chicken carboxylesterase at high temperatures could be explained. The docking of substrate analogs into the cuCES active site was used to explain the fact that the chicken carboxylesterase cannot hydrolyze efficiently large substrate molecules.

Modeling of mammalian olfactory receptors and docking of odorants

Olfactory receptors (ORs) belong to the superfamily of G protein-coupled receptors (GPCRs), the second largest class of genes after those related to immunity, and account for about 3 % of mammalian genomes. ORs are present in all multicellular organisms and represent more than half the GPCRs in mammalian species (e.g., the mouse OR repertoire contains >1,000 functional genes). ORs are mainly expressed in the olfactory epithelium where they detect odorant molecules, but they are also expressed in a number of other cells, such as sperm cells, although their functions in these cells remain mostly unknown. It has recently been reported that ORs are present in tumoral tissues where they are expressed at different levels than in healthy tissues. A specific OR is over-expressed in prostate cancer cells, and activation of this OR has been shown to inhibit the proliferation of these cells. Odorant stimulation of some of these receptors results in inhibition of cell proliferation. Even though their biological role has not yet been elucidated, these receptors might constitute new targets for diagnosis and therapeutics. It is important to understand the activation mechanism of these receptors at the molecular level, in particular to be able to predict which ligands are likely to activate a particular receptor ('deorphanization') or to design antagonists for a given receptor. In this review, we describe the in silico methodologies used to model the three-dimensional (3D) structure of ORs (in the more general framework of GPCR modeling) and to dock ligands into these 3D structures.

Molecular cloning and 3D structure modeling of APEX1, DNA base excision repair enzyme from the Camel, Camelus dromedarius

The domesticated one-humped camel, Camelus dromedarius, is one of the most important animals in the Arabian Desert. It is exposed most of its life to both intrinsic and extrinsic genotoxic factors that are known to cause gross DNA alterations in many organisms. Ionic radiation and sunlight are known producers of Reactive Oxygen Species (ROS), one of the causes for DNA lesions. The damaged DNA is repaired by many enzymes, among of them Base Excision Repair enzymes, producing the highly mutagenic apurinic/apyrimidinicsites (AP sites). Therefore, recognition of AP sites is fundamental to cell/organism survival. In the present work, the full coding sequence of a putative cAPEX1 gene was amplified for the first time from C. dromedarius by RT-PCR and cloned (NCBI accession number are HM209828 and ADJ96599 for nucleotides and amino acids, respectively). cDNA sequencing was deduced to be 1041 nucleotides, of which 954 nucleotides encode a protein of 318 amino acids, similar to the coding region of the APEX1 gene and the protein from many other species. The calculated molecular weight and isoelectric point of cAPEX1 using Bioinformatics tools was 35.5 kDa and 8.11, respectively. The relative expressions of cAPEX1 in camel kidney, spleen, lung and testis were examined using qPCR and compared with that of the liver using a 18S ribosomal subunit as endogenous control. The highest level of cAPEX1 transcript was found in the testis; 325% higher than the liver, followed by spleen (87%), kidney (20%) and lung (5%), respectively. The cAPEX1 is 94%–97% similar to their mammalian counterparts. Phylogenetic analysis revealed that cAPEX1 is grouped together with that of S. scrofa. The predicted 3D structure of cAPEX1 has similar folds and topology with the human (hAPEX1). The root-mean-square deviation (rmsd) between cAPEX1 and hAPEX1 was 0.582 and the Q-score was 0.939.