RESUMEN
INTRODUCTION: Human metabolism is sustained by functional networks that operate at diverse scales. Capturing local and global dynamics in the human body by hierarchically bridging multi-scale functional networks is a major challenge in physiological modeling. OBJECTIVES: To develop an interactive, user-friendly web application that facilitates the simulation and visualization of advection-dispersion transport in three-dimensional (3D) microvascular networks, biochemical exchange, and metabolic reactions in the tissue layer surrounding the vasculature. METHODS: To help modelers combine and simulate biochemical processes occurring at multiple scales, KiPhyNet deploys our discrete graph-based modeling framework that bridges functional networks existing at diverse scales. KiPhyNet is implemented in Python based on Apache web server using MATLAB as the simulator engine. KiPhyNet provides the functionality to assimilate multi-omics data from clinical and experimental studies as well as vascular data from imaging studies to investigate the role of structural changes in vascular topology on the functional response of the tissue. RESULTS: With the network topology, its biophysical attributes, values of initial and boundary conditions, parameterized kinetic constants, biochemical species-specific transport properties such as diffusivity as inputs, a user can use our application to simulate and view the simulation results. The results of steady-state velocity and pressure fields and dynamic concentration fields can be interactively examined. CONCLUSION: KiPhyNet provides barrier-free access to perform time-course simulation experiments by building multi-scale models of microvascular networks in physiology, using a discrete modeling framework. KiPhyNet is freely accessible at http://pallab.cds.iisc.ac.in/kiphynet/ and the documentation is available at https://deepamahm.github.io/kiphynet_docs/ .
Asunto(s)
Simulación por Computador , Programas Informáticos , Humanos , Cinética , Transporte Biológico/fisiología , Modelos Biológicos , InternetRESUMEN
This paper presents an evaluation of predictions submitted for the "HMBS" challenge, a component of the sixth round of the Critical Assessment of Genome Interpretation held in 2021. The challenge required participants to predict the effects of missense variants of the human HMBS gene on yeast growth. The HMBS enzyme, critical for the biosynthesis of heme in eukaryotic cells, is highly conserved among eukaryotes. Despite the application of a variety of algorithms and methods, the performance of predictors was relatively similar, with Kendall's tau correlation coefficients between predictions and experimental scores around 0.3 for a majority of submissions. Notably, the median correlation (≥ 0.34) observed among these predictors, especially the top predictions from different groups, was greater than the correlation observed between their predictions and the actual experimental results. Most predictors were moderately successful in distinguishing between deleterious and benign variants, as evidenced by an area under the receiver operating characteristic (ROC) curve (AUC) of approximately 0.7 respectively. Compared with the recent two rounds of CAGI competitions, we noticed more predictors outperformed the baseline predictor, which is solely based on the amino acid frequencies. Nevertheless, the overall accuracy of predictions is still far short of positive control, which is derived from experimental scores, indicating the necessity for considerable improvements in the field. The most inaccurately predicted variants in this round were associated with the insertion loop, which is absent in many orthologs, suggesting the predictors still heavily rely on the information from multiple sequence alignment.
RESUMEN
The transcription termination process is an important part of the gene expression process in the cell. It has been studied extensively, but many aspects of the mechanism are not well understood. The widespread availability of experimental RNA-seq data from high-throughput experiments provides a unique opportunity to infer the end of the transcription units genome wide. This data is available for both Rho-dependent and Rho-independent termination pathways that drive transcription termination in bacteria. Our book chapter gives an overview of the current knowledge of Rho-independent transcription termination mechanisms and the prediction approaches currently deployed to infer the termination sites. Thereafter, we describe our method that uses cluster hairpins to detect Rho-independent transcription termination sites. These clusters are a group of hairpins that lies at <15 bp from each other and are together capable of enforcing the termination process. The idea of a group of hairpins being extensively used for transcription termination is new, and results show that at least 52% of the total cases are of this type, while in the remaining cases, a single strong hairpin is capable of driving transcription termination. The reads derived from the RNA-seq data for corresponding bacteria have been used to validate the predicted sites. The predictions that match these RNA-seq derived sites have higher confidence, and we find almost 98% of the predicted sites, including alternate termination sites, to match the RNA-seq data. We discuss the features of predicted hairpins in detail for a better understanding of the Rho-independent transcription termination mechanism in bacteria. We also explain how users can use the tools developed by us to do transcription terminator predictions and design their experiments through genome-level visualization of the transcription termination sites from the precomputed INTERPIN database.
Asunto(s)
RNA-Seq , Terminación de la Transcripción Genética , RNA-Seq/métodos , Programas Informáticos , Biología Computacional/métodos , ARN Bacteriano/genética , Bacterias/genética , Análisis de Secuencia de ARN/métodos , Regiones Terminadoras Genéticas/genética , Regulación Bacteriana de la Expresión GénicaRESUMEN
The Endoplasmic Reticulum is a pervasive, dynamic cellular organelle that performs a wide range of functions in the eukaryotic cell, including protein folding and maturation. Upon stress, ER activates an adaptive cellular pathway, namely Unfolded Protein Response, that transduces information from ER to nucleus, restoring homeostasis in the ER milieu. UPR consists of three membrane-tethered sensors; IRE1, PERK and ATF6. Among all the UPR sensors, the IRE1 branch acts as a central pathway that orchestrates several pathways to determine cell fate. However, the detailed knowledge underlying the whole process is not understood yet. Previously, we determined the sMEK1 as one of the interacting partners of IRE1. sMEK1 is a protein phosphatase, which has been indicated in a number of critical cellular functions like apoptosis, cell proliferation, and tumour suppression. In this study, we evaluated the role of sMEK1 on the IRE1 signalling pathway. Our data indicate that sMEK1 can inhibit IRE1 phosphorylation under ER stress. This inhibitory effect of sMEK1 could be reflected in its downstream effectors, Xbp1 and RIDD, which are downregulated in the presence of sMEK1. We also found that the repressing effect of sMEK1 was specific to the IRE1 signalling pathway and could be preserved even under prolonged ER stress. Our findings also indicate that sMEK1 can inhibit IRE1 and its downstream molecules under ER stress irrespective of other UPR sensors. These results help to draw the mechanistic details giving insights into different molecular connections of UPR with other pathways.
RESUMEN
We compare the WebGeSTer and INtrinsic transcription TERmination hairPIN (INTERPIN) databases used for intrinsic transcription termination (ITT) site prediction in bacteria. The former deploys inverted nucleotide repeat detection for identification of RNA hairpin, while the latter a pair-potential function - the hairpin energy score evaluation being identical for both. We find INTERPIN more sensitive than WebGeSTer with about 6% and 51% additional predictions for ITTs in chromosomal and plasmid operons, respectively. INTERPIN hairpins are relatively shorter in length with ungapped stem, and even located in AT-rich segments, compared to GC-rich longer hairpins with a gapped stem in WebGeSTer. The GC%, length, and energy score from INTERPIN transcription units (TUs) are best inter-correlated while the lowest energy single hairpins from WebGeSTer, considered suitable for ITT, being the worst. Around 72% TUs from the two databases overlap, and â¼60% of all alternate ITT sites downstream of TUs overlap, of which 65% are cluster hairpins. This helps highlight hairpin features that can be used to identify termination sites in bacteria across different prediction methods. Overall, the pair-potential-function-based hairpins screened appear to be more consistent with the kinetic and thermodynamics processes of ITT known to date.Communicated by Ramaswamy H. Sarma.
RESUMEN
The value of accurate protein structural models closely conforming to the experimental data is indisputable. DREAMweb deploys an improved DREAM algorithm, DREAMv2, that incorporates a tighter bound in the constraint set of the underlying optimization approach. This reduces the artifacts while modeling the protein structure by solving the distance-geometry problem. DREAMv2 follows a bottom-up strategy of building smaller substructures for regions with a larger concentration of experimental bounds and consolidating them before modeling the rest of the protein structure. This improves secondary structure conformance in the final models consistent with experimental data. The proposed method efficiently models regions with sparse coverage of experimental data by reducing the possibility of artifacts compared to DREAM. To balance performance and accuracy, smaller substructures ( â¼ 200 $\sim 200$ atoms) are solved in this regime, allowing faster builds for the other parts under relaxed conditions. DREAMweb is accessible as an internet resource. The improvements in results are showcased through benchmarks on 10 structures. DREAMv2 can be used in tandem with any NMR-based protein structure determination workflow, including an iterative framework where the NMR assignment for the NOESY spectra is incomplete or ambiguous. DREAMweb is freely available for public use at http://pallab.cds.iisc.ac.in/DREAM/ and downloadable at https://github.com/niladriranjandas/DREAMv2.git.
RESUMEN
Millions of people throughout the world are affected by neurodegenerative disorders like Alzheimer's disease (AD), making them a major public health concern. To create successful medicines, early diagnosis and illness monitoring are required. Emerging as possible diagnostic and treatment tools for neurodegenerative illnesses are biomarkers such as microRNAs (miRNAs). In the realm of neuroscience, miRNAs have been discovered to function as essential regulators of gene expression, with roles spanning development, differentiation, and illness. Several neurodegenerative diseases, including AD, have been linked to miRNA dysregulation. As high-throughput methods have been developed for monitoring miRNA expression and identifying miRNA targets, miRNAs have become a prime candidate for use in diagnostics and therapy. The techniques for isolating miRNAs and the most up-to-date computational methods for finding miRNA target transcripts are both described in this chapter. This chapter will be a helpful reference for anyone investigating the role of miRNAs in AD and related neurodegenerative illnesses.
Asunto(s)
Enfermedad de Alzheimer , MicroARNs , Humanos , MicroARNs/genética , MicroARNs/metabolismo , Enfermedad de Alzheimer/metabolismo , Neuroprotección , Biomarcadores/metabolismo , Perfilación de la Expresión Génica/métodosRESUMEN
The pancreas plays a critical role in maintaining glucose homeostasis through the secretion of hormones from the islets of Langerhans. Glucose-stimulated insulin secretion (GSIS) by the pancreatic ß-cell is the main mechanism for reducing elevated plasma glucose. Here we present a systematic modeling workflow for the development of kinetic pathway models using the Systems Biology Markup Language (SBML). Steps include retrieval of information from databases, curation of experimental and clinical data for model calibration and validation, integration of heterogeneous data including absolute and relative measurements, unit normalization, data normalization, and model annotation. An important factor was the reproducibility and exchangeability of the model, which allowed the use of various existing tools. The workflow was applied to construct a novel data-driven kinetic model of GSIS in the pancreatic ß-cell based on experimental and clinical data from 39 studies spanning 50 years of pancreatic, islet, and ß-cell research in humans, rats, mice, and cell lines. The model consists of detailed glycolysis and phenomenological equations for insulin secretion coupled to cellular energy state, ATP dynamics and (ATP/ADP ratio). Key findings of our work are that in GSIS there is a glucose-dependent increase in almost all intermediates of glycolysis. This increase in glycolytic metabolites is accompanied by an increase in energy metabolites, especially ATP and NADH. One of the few decreasing metabolites is ADP, which, in combination with the increase in ATP, results in a large increase in ATP/ADP ratios in the ß-cell with increasing glucose. Insulin secretion is dependent on ATP/ADP, resulting in glucose-stimulated insulin secretion. The observed glucose-dependent increase in glycolytic intermediates and the resulting change in ATP/ADP ratios and insulin secretion is a robust phenomenon observed across data sets, experimental systems and species. Model predictions of the glucose-dependent response of glycolytic intermediates and biphasic insulin secretion are in good agreement with experimental measurements. Our model predicts that factors affecting ATP consumption, ATP formation, hexokinase, phosphofructokinase, and ATP/ADP-dependent insulin secretion have a major effect on GSIS. In conclusion, we have developed and applied a systematic modeling workflow for pathway models that allowed us to gain insight into key mechanisms in GSIS in the pancreatic ß-cell.
Asunto(s)
Células Secretoras de Insulina , Humanos , Animales , Ratones , Ratas , Secreción de Insulina , Reproducibilidad de los Resultados , Glucosa/farmacología , Adenosina TrifosfatoRESUMEN
The 'Nest' motif plays a functional role in protein owing to its ligand binding potential aided by geometric concavity. The presence of less favored left-handed conformation (L-state) in its structure makes this concavity possible and in shaping the native chemical environment amenable to stable binding interactions. To understand the persistent appearance of L-state torsion in the Nest motif, we analyzed 0.5µs Molecular Dynamics (MD) simulation trajectories of 35 six-residue peptides (out of a total of 50 large Nest sequences of ≥6 residues) identified in our previous study. Analysis of the MD trajectories of the individual peptides reveals initial L-state in 60% of the peptides persists for >40% of the trajectory. Further, Nests with different sequences appear to adopt a specific conformational state driven by the neighboring L-state residues. The sequences also possess short secondary structures and amino acid repeats, suggesting evolutionary conservation and the specific role of amino acids in locally predisposing the torsion angle to the L-state. These findings help us to understand how L-state conformation is an essential prerequisite in stabilizing the Nest motif and shed light on the sequence-structure-function paradigm in the rational design of peptides and peptidomimetics for therapeutics.Communicated by Ramaswamy H. Sarma.
RESUMEN
The large-scale detection of putative intrinsic transcription terminators is limited to only a few bacteria currently. We discovered a group of hairpins, called cluster hairpins, present within 15 nucleotides from each other. These are expected to work in tandem to cause intrinsic transcription termination (ITT), while the single hairpin can do the same alone. Therefore, exploring these ITT sites and the hairpins across bacterial genomes becomes highly desirable. INTERPIN is the largest archived collection of in silico inferred ITT hairpins in bacteria, covering 12745 bacterial genomes and encompassing ten bacterial phyla for â¼25 million hairpins. Users can obtain details on operons, individual cluster, and single ITT hairpins that were screened therein. Integrated Genome Viewer (IGV) software interactively visualizes hairpin secondary and tertiary structures in the genomic context. We also discuss statistics for the occurrence of cluster or single hairpins and other termination alternatives while showing the validation of predicted hairpins against in vivo detected hairpins. The database is freely available at http://pallab.cds.iisc.ac.in/INTERPIN/. INTERPIN (database and software) can make predictions for both AT and GC-rich genomes, which has not been achieved by any other program so far. It can also be used to improve genome annotation as well as to get predictions to improve the understanding of the ITT pathway by further analysis.
Asunto(s)
Bacterias , Transcripción Genética , Bacterias/genética , Genómica , Programas Informáticos , Genoma BacterianoRESUMEN
IRE1 belongs to a type I transmembrane protein family harboring two functional domains, cytoplasmic domain with kinase and RNAse catalytic activity, and the luminal domain, which is involved in the sensing of unfolded proteins. IRE1 molecule undergoes dimerization in the lumenal domain, which functionally activates the catalytic C-terminal domain. IRE1 activation is directly related to transition between monomeric and dimeric forms. We have deduced two quaternary structures from the published crystal structure of IRE1. One structure with a large stable interface that requires large activation and deactivation energy to active IRE1. The other quaternary structure has low dissociation energy and is more suitable for IRE1 oligomeric transition.
RESUMEN
The pancreas plays a critical role in maintaining glucose homeostasis through the secretion of hormones from the islets of Langerhans. Glucose-stimulated insulin secretion (GSIS) by the pancreatic ß -cell is the main mechanism for reducing elevated plasma glucose. Here we present a systematic modeling workflow for the development of kinetic pathway models using the Systems Biology Markup Language (SBML). Steps include retrieval of information from databases, curation of experimental and clinical data for model calibration and validation, integration of heterogeneous data including absolute and relative measurements, unit normalization, data normalization, and model annotation. An important factor was the reproducibility and exchangeability of the model, which allowed the use of various existing tools. The workflow was applied to construct the first consensus model of GSIS in the pancreatic ß -cell based on experimental and clinical data from 39 studies spanning 50 years of pancreatic, islet, and ß -cell research in humans, rats, mice, and cell lines. The model consists of detailed glycolysis and equations for insulin secretion coupled to cellular energy state (ATP/ADP ratio). Key findings of our work are that in GSIS there is a glucose-dependent increase in almost all intermediates of glycolysis. This increase in glycolytic metabolites is accompanied by an increase in energy metabolites, especially ATP and NADH. One of the few decreasing metabolites is ADP, which, in combination with the increase in ATP, results in a large increase in ATP/ADP ratios in the ß -cell with increasing glucose. Insulin secretion is dependent on ATP/ADP, resulting in glucose-stimulated insulin secretion. The observed glucose-dependent increase in glycolytic intermediates and the resulting change in ATP/ADP ratios and insulin secretion is a robust phenomenon observed across data sets, experimental systems and species. Model predictions of the glucose-dependent response of glycolytic intermediates and insulin secretion are in good agreement with experimental measurements. Our model predicts that factors affecting ATP consumption, ATP formation, hexokinase, phosphofructokinase, and ATP/ADP-dependent insulin secretion have a major effect on GSIS. In conclusion, we have developed and applied a systematic modeling workflow for pathway models that allowed us to gain insight into key mechanisms in GSIS in the pancreatic ß -cell.
RESUMEN
IRE1 is a transmembrane signalling protein that activates the unfolded protein response under endoplasmic reticulum stress. IRE1 is endowed with kinase and endoribonuclease activities. The ribonuclease activity of IRE1 can switch substrate specificities to carry out atypical splicing of Xbp1 mRNA or trigger the degradation of specific mRNAs. The mechanisms regulating the distinct ribonuclease activities of IRE1 have yet to be fully understood. Here, we report the Bcl-2 family protein Bid as a novel recruit of the IRE1 complex, which directly interacts with the cytoplasmic domain of IRE1. Bid binding to IRE1 leads to a decrease in IRE1 phosphorylation in a way that it can only perform Xbp1 splicing while mRNA degradation activity is repressed. The RNase outputs of IRE1 have been found to regulate the homeostatic-apoptotic switch. This study, thus, provides insight into IRE1-mediated cell survival.
Asunto(s)
Proteínas Serina-Treonina Quinasas , Respuesta de Proteína Desplegada , Estrés del Retículo Endoplásmico/fisiología , Endorribonucleasas/genética , Endorribonucleasas/metabolismo , Proteínas de la Membrana/metabolismo , Proteínas Serina-Treonina Quinasas/genética , Proteínas Serina-Treonina Quinasas/metabolismo , Ribonucleasas/metabolismo , Proteína 1 de Unión a la X-Box/genética , Proteína 1 de Unión a la X-Box/metabolismo , Proteína Proapoptótica que Interacciona Mediante Dominios BH3/metabolismoRESUMEN
The ß-Coronavirus mouse hepatitis virus (MHV-A59)-RSA59 has a patent stretch of fusion peptide (FP) containing two consecutive central prolines (PP) in the S2 domain of the Spike protein. Our previous studies compared the PP-containing fusogenic-demyelinating strain RSA59(PP) to its one proline-deleted mutant strain RSA59(P) and one proline-containing non-fusogenic non-demyelinating parental strain RSMHV2(P) to its one proline inserted mutant strain RSMHV2(PP). These studies highlighted the crucial role of PP in fusogenicity, hepato-neuropathogenesis, and demyelination. Computational studies combined with biophysical data indicate that PP at the center of the FP provides local rigidity while imparting global fluctuation to the Spike protein that enhances the fusogenic properties of RSA59(PP) and RSMHV2(PP). To elaborate on the understanding of the role of PP in the FP of MHV, the differential neuroglial tropism of the PP and P mutant strains was investigated. Comparative studies demonstrated that PP significantly enhances the viral tropism for neurons, microglia, and oligodendrocytes. PP, however, is not essential for viral tropism for either astroglial or oligodendroglial precursors or the infection of meningeal fibroblasts in the blood-brain and blood-CSF barriers. PP in the fusion domain is critical for promoting gliopathy, making it a potential region for designing antivirals for neuro-COVID therapy.
Asunto(s)
Virus de la Hepatitis Murina , Glicoproteína de la Espiga del Coronavirus , Tropismo Viral , Animales , Ratones , Virus de la Hepatitis Murina/fisiología , Péptidos/metabolismo , Prolina , Glicoproteína de la Espiga del Coronavirus/genética , Glicoproteína de la Espiga del Coronavirus/metabolismo , Proteínas del Envoltorio Viral/metabolismoRESUMEN
Cytochrome b (QcrB) is considered an essential subunit in the electron transport chain that coordinates the action of the entire cytochrome bc1 oxidase. It has been identified as an attractive drug target for a new promising clinical candidate Q203 that depletes the intracellular ATP levels in the bacterium, Mycobacterium tuberculosis. However, single point polymorphism (T313A/I) near the quinol oxidation site of QcrB developed resistance to Q203. In the present study, we analyze the structural changes and drug-resistance mechanism of QcrB due to the point mutation in detail through conformational morphing and molecular docking studies. By morphing, we generated conformers between the open and closed state of the electron transporting cytochrome bc1-aa3 super complex. We clustered them to identify four intermediate structures and relevant intra- and intermolecular motions that may be of functional relevance, especially the binding of Q203 in wild and mutant QcrB intermediate structures and their alteration in developing drug resistance. The difference in the binding score and hydrogen bond interactions between Q203 and the wild-type and mutant intermediate structures of QcrB from molecular docking studies showed that the point mutation T313A severely affected the binding affinity of the candidate drug. Together, the findings provide an in-depth understanding of QcrB inhibition in different conformations, including closed, intermediate, and open states of cytochrome bc1-aa3 super complex in Mycobacterium tuberculosis at the atomic level. We also obtain insights for designing QcrB and cytochrome bc1-aa3 inhibitors as potential therapeutics that may combat drug resistance in tuberculosis.
Asunto(s)
Complejo III de Transporte de Electrones , Mycobacterium tuberculosis , Complejo III de Transporte de Electrones/genética , Complejo III de Transporte de Electrones/metabolismo , Ligandos , Simulación del Acoplamiento Molecular , Mycobacterium tuberculosis/genética , Análisis por ConglomeradosRESUMEN
Nuclear magnetic resonance (NMR) spectroscopy can reveal conformational states of a protein in physiological conditions. However, sparsely available NMR data for a protein with large degrees of freedom can introduce structural artifacts in the built models. Currently used state-of-the-art methods deriving protein structure and conformation from NMR deploy molecular dynamics (MD) coupled with simulated annealing for building models. We provide an alternate graph-based modeling approach, where we first build substructures from NMR-derived distance-geometry constraints combined in one shot to form the core structure. The remaining molecule with inadequate data is modeled using a hybrid approach respecting the observed distance-geometry constraints. One-shot structure building is rarely undertaken for large and sparse data systems, but our data-driven bottom-up approach makes this uniquely feasible by suitable partitioning of the problem. A detailed comparison of select models with state-of-art methods reveals differences in the secondary structure regions wherein the correctness of our models is confirmed by NMR data. Benchmarking of 106 protein-folds covering 38-282 length structures shows minimal experimental-constraint violations while conforming to other structure quality parameters such as the proper folding, steric clash, and torsion angle violation based on Ramachandran plot criteria. Comparative MD studies using select protein models from a state-of-art method and ours under identical experimental parameters reveal distinct conformational dynamics that could be attributed to protein structure-function. Our work is thus useful in building enhanced NMR-evidence-based models that encapsulate the contextual secondary and tertiary structure variations present during the experimentation and expand the scope of functional inference.
Asunto(s)
Proteínas , Conformación Proteica , Modelos Moleculares , Proteínas/química , Espectroscopía de Resonancia Magnética/métodos , Estructura Secundaria de ProteínaRESUMEN
Binding studies of the water-soluble thiadicarbocyanine dye 3,3'-diethylthiadicarbocyanine acetate (DTC) with bovine serum albumin (BSA) were examined under physiological conditions using spectroscopic techniques like fluorescence, UV-Visible, circular dichroism (CD), FT-IR and molecular docking methods. Compiled experimental results envisage that DTC quench the fluorescence intensity of BSA. The increasing binding constants (K) were found to be in the order of 103 Mol-1 as a function of temperature, as calculated from the fluorescence quenching data. The quenching mechanism, thermodynamic parameters (ΔH0, ΔS0 and ΔG0) and the number of binding sites have been explored. CD values showed that the secondary structure of the BSA has been altered upon binding to DTC. Displacement experiments were carried out with different site probes to find out the binding site of DTC on BSA and it was found that binding interaction at site II of sub-domain IIIA. The interference of common metal ions on the interaction of DTC with BSA has also been studied. The experimental data exhibit that DTC interacts with BSA by hydrophobic forces. The experimental findings from BSA binding studies were validated by using in silico molecular docking technique. The results of the investigations were accurately supported by studies on molecular docking. The optimal shape of the molecular probe demonstrated the affinity as a free binding energy release of -7.37 Kcal/mol. The present research report endeavors to the approachable nature of water-soluble DTC dye and paves way for targeted biological interactions.Communicated by Ramaswamy H. Sarma.
Asunto(s)
Colorantes Fluorescentes , Albúmina Sérica Bovina , Simulación del Acoplamiento Molecular , Unión Proteica , Albúmina Sérica Bovina/química , Espectroscopía Infrarroja por Transformada de Fourier , Espectrometría de Fluorescencia , Sitios de Unión , Termodinámica , Dicroismo Circular , Agua , Espectrofotometría UltravioletaRESUMEN
SUMMARY: Molecular dynamics (MD) simulations have become an integral part of biomolecular study. Most MD software suites do not include analysis tools and those which do create very basic visualizations. Molecular Dynamics Data Visualizer (MD DaVis) is a python package developed to facilitate quick comparative analysis of MD trajectories of similar proteins or the same protein under different conditions. MD DaVis can quickly generate interactive visualization from molecular dynamics trajectories with a few simple steps. Interactive plots eliminate the need to make multiple plots for comparison, improving productivity and saving time. AVAILABILITY AND IMPLEMENTATION: MD DaVis is an open-source Python 3 package (https://pypi.org/project/md-davis/) distributed under MIT license. The source code is available at https://github.com/djmaity/md-davis or https://doi.org/10.5281/zenodo.6227047. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Asunto(s)
Visualización de Datos , Simulación de Dinámica Molecular , Programas Informáticos , ProteínasRESUMEN
Combined in silico, in vitro, and in vivo comparative studies between isogenic-recombinant Mouse-Hepatitis-Virus-RSA59 and its proline deletion mutant, revealed a remarkable contribution of centrally located two consecutive prolines (PP) from Spike protein fusion peptide (FP) in enhancing virus fusogenic and hepato-neuropathogenic potential. To deepen our understanding of the underlying factors, we extend our studies to a non-fusogenic parental virus strain RSMHV2 (P) with a single proline in the FP and its proline inserted mutant, RSMHV2 (PP). Comparative in vitro and in vivo studies between virus strains RSA59(PP), RSMHV2 (P), and RSMHV2 (PP) in the FP demonstrate that the insertion of one proline significantly resulted in enhancing the virus fusogenicity, spread, and consecutive neuropathogenesis. Computational studies suggest that the central PP in Spike FP induces a locally ordered, compact, and rigid structure of the Spike protein in RSMHV2 (PP) compared to RSMHV2 (P), but globally the Spike S2-domain is akin to the parental strain RSA59(PP), the latter being the most flexible showing two potential wells in the energy landscape as observed from the molecular dynamics studies. The critical location of two central prolines of the FP is essential for fusogenicity and pathogenesis making it a potential site for designing antiviral.
Asunto(s)
Enfermedades Desmielinizantes , Glicoproteína de la Espiga del Coronavirus , Animales , Ratones , Péptidos/metabolismo , Prolina , Glicoproteína de la Espiga del Coronavirus/genética , Glicoproteína de la Espiga del Coronavirus/metabolismo , Proteínas del Envoltorio Viral/metabolismoRESUMEN
Emerging mutations in the SARS-CoV-2 genome pose a challenge for vaccine development and antiviral therapy. The antiviral efficacy of Azadirachta indica bark extract (NBE) was assessed against SARS-CoV-2 and m-CoV-RSA59 infection. Effects of in vivo intranasal or oral NBE administration on viral load, inflammatory response, and histopathological changes were assessed in m-CoV-RSA59-infection. NBE administered inhibits SARS-CoV-2 and m-CoV-RSA59 infection and replication in vitro, reducing Envelope and Nucleocapsid gene expression. NBE ameliorates neuroinflammation and hepatitis in vivo by restricting viral replication and spread. Isolated fractions of NBE enriched in Nimbin isomers shows potent inhibition of m-CoV-RSA59 infection in vitro. In silico studies revealed that NBE could target Spike and RdRp of m-CoV and SARS-CoV-2 with high affinity. NBE has a triterpenoids origin that may allow them to competitively target panoply of viral proteins to inhibit mouse and different strains of human coronavirus infections, suggesting its potential as an antiviral against pan-ß-Coronaviruses.