Homology Modeling, Molecular Docking and DNA Binding Studies of Nucleotide Excision Repair UvrC Protein from M. tuberculosis

Molecular modeling, simulation and docking of Rv1250 protein from Mycobacterium tuberculosis

Computational prediction and protein structure modeling have come to the aid of various biological problems in determining the structure of proteins. These technologies have revolutionized the biological world of research, allowing scientists and researchers to gain insights into their biological questions and design experimental research much more efficiently. Pathogenic Mycobacterium spp. is known to stay alive within the macrophages of its host. Mycobacterium tuberculosis is an acid-fast bacterium that is the most common cause of tuberculosis and is considered to be the main cause of resistance of tuberculosis as a leading health issue. The genome of Mycobacterium tuberculosis contains more than 4,000 genes, of which the majority are of unknown function. An attempt has been made to computationally model and dock one of its proteins, Rv1250 (MTV006.22), which is considered as an apparent drug-transporter, integral membrane protein, and member of major facilitator superfamily (MFS). The most widely used techniques, i.e., homology modeling, molecular docking, and molecular dynamics (MD) simulation in the field of structural bioinformatics, have been used in the present work to study the behavior of Rv1250 protein from M. tuberculosis. The structure of unknown TB protein, i.e., Rv1250 was retrived using homology modeling with the help of I-TASSER server. Further, one of the sites responsible for infection was identified and docking was done by using the specific Isoniazid ligand which is an inhibitor of this protein. Finally, the stability of protein model and analysis of stable and static interaction between protein and ligand molecular dynamic simulation was performed at 100 ns The designing of novel Rv1250 enzyme inhibitors is likely achievable with the use of proposed predicted model, which could be helpful in preventing the pathogenesis caused by M. tuberculosis. Finally, the MD simulation was done to evaluate the stability of the ligand for the specific protein.

In vivo developmental studies of Helicoverpa armigera and in silico molecular interactions with trypsin reveal the bio-insecticidal potential of trypsin inhibitor (SSTI) isolated from Solanum surattense

We report the impact of gut protease inhibition on the development of Helicoverpa armigera by trypsin inhibitor and the use of molecular modeling to understand the mechanism of trypsin inhibition. Larvae of H. armigera fed on an artificial diet containing 150 and 300 μg/ml SSTI showed a negative impact on the insects' development in terms of mean larval weight, larval fatality, survival rate, and nutritional indices. Prominent physical abnormalities like curled wings, malformed appendages, and small body size were observed during the development. Gene expression studies revealed down regulation in trypsin (HaTry 1, 2, 3, 4, 6, 8) and chymotrypsin (HaChy 1, 2, 3, 4) genes of the larval gut upon treatment of SSTI. Homology modeling has been used to build the three-dimensional structure of SSTI, which showed β-sheets having a stable canonical inhibitory loop (CIL) with conserved lysine residue. Molecular docking studies showed the strong binding of SSTI at the active site of trypsin. Molecular dynamic (MD) simulation revealed the stable interactions of the rigid CIL of SSTI at the active site of trypsin, leading to its destabilization. Conserved lysine63 of the P1 site in SSTI forms a strong hydrogen bonding network with residues Asp189 and Ser190 of trypsin.

Interrogating the substrate specificity landscape of UvrC reveals novel insights into its non-canonical function

Although it is relatively unexplored, accumulating data highlight the importance of tripartite crosstalk between nucleotide excision repair (NER), DNA replication, and recombination in the maintenance of genome stability; however, elucidating the underlying mechanisms remains challenging. While Escherichia coli uvrA and uvrB can fully complement polAΔ cells in DNA replication, uvrC attenuates this alternative DNA replication pathway, but the exact mechanism by which uvrC suppresses DNA replication is unknown. Furthermore, the identity of bona fide canonical and non-canonical substrates for UvrCs are undefined. Here, we reveal that Mycobacterium tuberculosis UvrC (MtUvrC) strongly binds to, and robustly cleaves, key intermediates of DNA replication/recombination as compared with the model NER substrates. Notably, inactivation of MtUvrC ATPase activity significantly attenuated its endonuclease activity, thus suggesting a causal link between these two functions. We built an in silico model of the interaction of MtUvrC with the Holliday junction (HJ), using a combination of homology modeling, molecular docking, and molecular dynamic simulations. The model predicted residues that were potentially involved in HJ binding. Six of these residues were mutated either singly or in pairs, and the resulting MtUvrC variants were purified and characterized. Among them, residues Glu595 and Arg597 in the helix-hairpin-helix motif were found to be crucial for the interaction between MtUvrC and HJ; consequently, mutations in these residues, or inhibition of ATP hydrolysis, strongly abrogated its DNA-binding and endonuclease activities. Viewed together, these findings expand the substrate specificity landscape of UvrCs and provide crucial mechanistic insights into the interplay between NER and DNA replication/recombination.

Rad5 HIRAN domain: Structural insights into its interaction with ssDNA through molecular modeling approaches

The Rad5 protein is an SWI/SNF family ubiquitin ligase that contains an N-terminal HIRAN domain and a RING C3HC4 motif. The HIRAN domain is critical for recognition of the stalled replication fork during the replication process and acts as a sensor to initiate the damaged DNA checkpoint. It is a conserved domain widely distributed in eukaryotic organisms and is present in several DNA-binding proteins from all kingdoms. Here we showed that distant species have important differences in key residues that affect affinity for ssDNA. Based on these findings, we hypothesized that different HIRAN domains might affect fork reversal and translesion synthesis through different metabolic processes. To address this question, we predicted the tertiary structure of both yeast and human HIRAN domains using molecular modeling. Structural dynamics experiments showed that the yeast HIRAN domain exhibited higher structural denaturation than its human homolog, although both domains became stable in the presence of ssDNA. Analysis of atomic contacts revealed that a greater number of interactions between the ssDNA nucleotides and the Rad5 domain are electrostatic. Taken together, these results provide new insights into the molecular mechanism of the HIRAN domain of Rad5 and may guide us to further elucidate differences in the ancient eukaryotes HIRAN sequences and their DNA affinity. Communicated by Ramaswamy H. Sarma.

Characterization of squalene synthase gene from Gymnema sylvestre R. Br

Background

Squalene synthase (SQS) is a rate-limiting enzyme necessary to produce pentacyclic triterpenes in plants. It is an important enzyme producing squalene molecules required to run steroidal and triterpenoid biosynthesis pathways working in competitive inhibition mode. Reports are available on information pertaining to SQS gene in several plants, but detailed information on SQS gene in Gymnema sylvestre R. Br. is not available. G. sylvestre is a priceless rare vine of central eco-region known for its medicinally important triterpenoids. Our work aims to characterize the GS-SQS gene in this high-value medicinal plant.

Results

Coding DNA sequences (CDS) with 1245 bp length representing GS-SQS gene predicted from transcriptome data in G. sylvestre was used for further characterization. The SWISS protein structure modeled for the GS-SQS amino acid sequence data had MolProbity Score of 1.44 and the Clash Score 3.86. The quality estimates and statistical score of Ramachandran plots analysis indicated that the homology model was reliable. For full-length amplification of the gene, primers designed from flanking regions of CDS encoding GS-SQS were used to get amplification against genomic DNA as template which resulted in approximately 6.2-kb sized single-band product. The sequencing of this product through NGS was carried out generating 2.32 Gb data and 3347 number of scaffolds with N50 value of 457 bp. These scaffolds were compared to identify similarity with other SQS genes as well as the GS-SQSs of the transcriptome. Scaffold_3347 representing the GS-SQS gene harbored two introns of 101 and 164 bp size. Both these intronic regions were validated by primers designed from adjoining outside regions of the introns on the scaffold representing GS-SQS gene. The amplification took place when the template was genomic DNA and failed when the template was cDNA confirmed the presence of two introns in GS-SQS gene in Gymnema sylvestre R. Br.

Conclusion

This study shows GS-SQS gene was very closely related to Coffea arabica and Gardenia jasminoides and this gene harbored two introns of 101 and 164 bp size.

Structural insights and inhibition mechanism of TMPRSS2 by experimentally known inhibitors Camostat mesylate, Nafamostat and Bromhexine hydrochloride to control SARS-coronavirus-2: A molecular modeling approach

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) has been responsible for the cause of global pandemic Covid-19 and to date, there is no effective treatment available. The spike ‘S’ protein of SARS-CoV-2 and ACE2 of the host cell are being targeted to design new drugs to control Covid-19. Similarly, a transmembrane serine protease, TMPRSS2 of the host cell plays a significant role in the proteolytic cleavage of viral ‘S’ protein helpful for the priming of ACE2 receptors and viral entry into human cells. However, three-dimensional structural information and the inhibition mechanism of TMPRSS2 is yet to be explored experimentally. Hence, we have used a molecular dynamics (MD) simulated homology model of TMPRSS2 to study the inhibition mechanism of experimentally known inhibitors Camostat mesylate, Nafamostat and Bromhexine hydrochloride (BHH) using molecular modeling techniques. Prior to docking, all three inhibitors were geometry optimized by semi-empirical quantum chemical RM1 method. Molecular docking analysis revealed that Camostat mesylate and its structural analogue Nafamostat interact strongly with residues His296 and Ser441 present in the catalytic triad of TMPRSS2, whereas BHH binds with Ala386 along with other residues. Comparative molecular dynamics simulations revealed the stable behavior of all the docked complexes. MM-PBSA calculations also revealed the stronger binding of Camostat mesylate to TMPRSS2 active site residues as compared to Nafamostat and BHH. Thus, this structural information could be useful to understand the mechanistic approach of TMPRSS2 inhibition, which may be helpful to design new lead compounds to prevent the entry of SARS-Coronavirus 2 in human cells.

The intrinsic ATPase activity of Mycobacterium tuberculosis UvrC is crucial for its damage-specific DNA incision function

To ensure genome stability, bacteria have evolved a network of DNA repair mechanisms; among them, the UvrABC-dependent nucleotide excision repair (NER) pathway is essential for the incision of a variety of bulky adducts generated by exogenous chemicals, UV radiation and by-products of cellular metabolism. However, very little is known about the enzymatic properties of Mycobacterium tuberculosis UvrABC excinuclease complex. Furthermore, the biochemical properties of Escherichia coli UvrC (EcUvrC) are not well understood (compared to UvrA and UvrB), perhaps due to its limited availability and/or activity instability in vitro. In addition, homology modelling of M. tuberculosis UvrC (MtUvrC) revealed the presence of a putative ATP-binding pocket, although its function remains unknown. To elucidate the biochemical properties of UvrC, we constructed and purified wild-type MtUvrC and its eight variants harbouring mutations within the ATP-binding pocket. The data from DNA-binding studies suggest that MtUvrC exhibits high-affinity for duplex DNA containing a bubble or fluorescein-dT moiety, over fluorescein-adducted single-stranded DNA. Most notably, MtUvrC has an intrinsic UvrB-independent ATPase activity, which drives dual incision of the damaged DNA strand. In contrast, EcUvrC is devoid of ATPase activity; however, it retains the ability to bind ATP at levels comparable to that of MtUvrC. The ATPase-deficient variants map to residues lining the MtUvrC ATP-binding pocket. Further analysis of these variants revealed separation of function between ATPase and DNA-binding activities in MtUvrC. Altogether, these findings reveal functional diversity of the bacterial NER machinery and a paradigm for the evolution of a catalytic scaffold in UvrC.

Role of Trace Elements as Cofactor: An Efficient Strategy toward Enhanced Biobutanol Production

Metabolic engineering has the potential to steadily enhance product titers by inducing changes in metabolism. Especially, availability of cofactors plays a crucial role in improving efficacy of product conversion. Hence, the effect of certain trace elements was studied individually or in combinations, to enhance butanol flux during its biological production. Interestingly, nickel chloride (100 mg L^-1) and sodium selenite (1 mg L^-1) showed a nearly 2-fold increase in solvent titer, achieving 16.13 ± 0.24 and 12.88 ± 0.36 g L^-1 total solvents with yields of 0.30 and 0.33 g g^-1, respectively. Subsequently, the addition time (screened entities) was optimized (8 h) to further increase solvent production up to 18.17 ± 0.19 and 15.5 ± 0.13 g L^-1 by using nickel and selenite, respectively. A significant upsurge in butanol dehydrogenase (BDH) levels was observed, which reflected in improved solvent productions. Additionally, a three-dimensional structure of BDH was also constructed using homology modeling and subsequently docked with substrate, cofactor, and metal ion to investigate proper orientation and molecular interactions.

Differentially expressed novel alternatively spliced transcript variant of tumor suppressor Stk11 gene in mouse

Serine/threonine kinase 11 (STK11) is a protein kinase that is encoded by Stk11 gene located on chromosome 19 and 10 in humans and mouse respectively. It acts as a master kinase of adenine monophosphate-activated protein kinase (AMPK) pathway that coordinates the regulation of cellular energy metabolism and cell division. STK11 exerts effect by activating more than 14 kinases including AMPK and AMPK-related kinases. It is also known to regulate cell polarity and acts as tumor suppressor. Alternative splicing of pre-mRNA is a mechanism which results in multiple transcript variants of a single gene. In human, two STK11 isoforms have been reported, an alternatively spliced isoform which has variation at its C-terminal and mostly expressed in testis (LKB1_S). Another isoform exhibiting oncogenic properties lacks few residues at its N-terminal (ΔN-LKB1). In the present study, we report the identification of a new transcript variant Stk11N which is generated through alternative splicing. The new variant was found to have differential and tissue specific expression at Postnatal-7 and adult stages of mouse. As compared to the known variant Stk11C, the conceptually translated amino acid sequences of the new variant differ from exon-E2 onwards. In silico post translational studies of the new and published variant show similarity in some of the properties while differ in properties like nuclear export signals, phosphorylation, glycosylation, etc. Thus, alternative splicing of Stk11 gene generating new variant with heterogeneous properties suggests for complex regulation of these variants in controlling the AMPK pathway and other functions.

Molecular Docking Technique to Understand Enzyme-Ligand Interactions

Molecular docking has advanced to such an extent that one can rapidly and accurately identify pharmaceutically useful lead compounds. It is being used routinely to understand molecular interactions between enzyme and ligand molecules. Several computational approaches are combined with experimental work to investigate molecular mechanisms in detail at the atomic level. Molecular docking method is also useful to investigate proper orientation and interactions between receptor and ligand. In this chapter we have discussed protein-protein approach to understand interactions between enzyme and amyloid beta (Aß) peptide. The Aß peptide is a causative agent of Alzheimer's disease. The Aß peptides can be cleaved specifically by several enzymes. Their interactions with Aß peptide and specific enzyme can be investigated using molecular docking. Thus, the molecular information obtained from docking studies might be useful to design new therapeutic approaches in treatment of Alzheimer's as well as several other diseases.

Molecular modeling approach to explore the role of cathepsin B from Hordeum vulgare in the degradation of Aβ peptides

The pathological hallmark of Alzheimer's disease is the accumulation of Aβ peptides in human brains. These Aβ peptides can be degraded by several enzymes such as hACE, hECE, hIDE and cathepsin B. Out of which cathepsin B also belongs to the papain super family and has been found in human brains, it has a role in Aβ peptide degradation through limited proteolysis. The Aβ concentrations are maintained properly by its production and clearance via receptor-mediated cellular uptake and direct enzymatic degradation. However, the reduced production of Aβ degrading enzymes as well as their Aβ degrading activity in human brains initiate the process of accumulation of Aβ peptides. So it becomes essential to investigate the molecular interactions involved in the process of Aβ degradation in detail at the atomic level. Hence, homology modeling, molecular docking and molecular dynamics simulation techniques have been used to explore the possible role of cathepsin B from Hordeum vulgare in the degradation of amyloid beta (Aβ) peptides. The homology model of cathepsin B from Hordeum vulgare shows good similarity with human cathepsin B. Molecular docking and MD simulation results revealed that the active site residues Cys32, HIS112, HIS113 are involved in the catalytic activity of cathepsin B. The sulfhydryl group of the Cys32 residue of cathepsin B from Hordeum vulgare cleaves the Aβ peptide from the carboxylic end of Glu11. Hence, this structural study might be helpful in designing alternative strategies for the treatment of AD.

India's Computational Biology Growth and Challenges

India's computational science is growing swiftly due to the outburst of internet and information technology services. The bioinformatics sector of India has been transforming rapidly by creating a competitive position in global bioinformatics market. Bioinformatics is widely used across India to address a wide range of biological issues. Recently, computational researchers and biologists are collaborating in projects such as database development, sequence analysis, genomic prospects and algorithm generations. In this paper, we have presented the Indian computational biology scenario highlighting bioinformatics-related educational activities, manpower development, internet boom, service industry, research activities, conferences and trainings undertaken by the corporate and government sectors. Nonetheless, this new field of science faces lots of challenges.

Evaluation of a possible role of Stigmatella aurantiaca ACE in Aβ peptide degradation: a molecular modeling approach

Amyloid-β (Aβ)-degrading enzymes are known to degrade Aβ peptides, a causative agent of Alzheimer's disease. These enzymes are responsible for maintaining Aβ concentration. However, loss of such enzymes or their Aβ-degrading activity because of certain genetic as well as nongenetic reasons initiates the accumulation of Aβ peptides in the human brain. Considering the limitations of the human enzymes in clearing Aβ peptide, the search for microbial enzymes that could cleave Aβ is necessary. Hence, we built a three-dimensional model of angiotensin-converting enzyme (ACE) from Stigmatella aurantiaca using homology modeling technique. Molecular docking and molecular dynamics simulation techniques were used to outline the possible cleavage mechanism of Aβ peptide. These findings suggest that catalytic residue Glu 434 of the model could play a crucial role to degrade Aβ peptide between Asp 7 and Ser 8. Thus, ACE from S. aurantiaca might cleave Aβ peptides similar to human ACE and could be used to design new therapeutic strategies against Alzheimer's disease.

Combined pharmacophore modeling, 3D-QSAR, homology modeling and docking studies on CYP11B1 inhibitors

The mitochondrial cytochrome P450 enzymes inhibitor steroid 11β-hydroxylase (CYP11B1) can decrease the production of cortisol. Therefore, these inhibitors have an effect in the treatment of Cushing’s syndrome. A pharmacophore model generated by Genetic Algorithm with Linear Assignment for Hypermolecular Alignment of Datasets (GALAHAD) was used to align the compounds and perform comparative molecular field analysis (CoMFA) with Q² = 0.658, R² = 0.959. The pharmacophore model contained six hydrophobic regions and one acceptor atom, and electropositive and bulky substituents would be tolerated at the A and B sites, respectively. A three-dimensional quantitative structure-activity relationship (3D-QSAR) study based on the alignment with the atom root mean square (RMS) was applied using comparative molecular field analysis (CoMFA) with Q² = 0.666, R² = 0.978, and comparative molecular similarity indices analysis (CoMSIA) with Q² = 0.721, R² = 0.972. These results proved that all the models have good predictability of the bioactivities of inhibitors. Furthermore, the QSAR models indicated that a hydrogen bond acceptor substituent would be disfavored at the A and B groups, while hydrophobic groups would be favored at the B site. The three-dimensional (3D) model of the CYP11B1 was generated based on the crystal structure of the CYP11B2 (PDB code 4DVQ). In order to probe the ligand-binding modes, Surflex-dock was employed to dock CYP11B1 inhibitory compounds into the active site of the receptor. The docking result showed that the imidazolidine ring of CYP11B1 inhibitors form H bonds with the amino group of residue Arg155 and Arg519, which suggested that an electronegative substituent at these positions could enhance the activities of compounds. All the models generated by GALAHAD QSAR and Docking methods provide guidance about how to design novel and potential drugs for Cushing’s syndrome treatment.

Insight into molecular interactions of Aβ peptide and gelatinase from Enterococcus faecalis: a molecular modeling approach

Alzheimer's disease is characterized by the presence of extracellular deposition of amyloid beta (Aβ) peptides.