Theoretical Biology and Medical Modelling最新文献

Background: blaVEB-1 is an integron-located extended-spectrum β-lactamase gene initially detected in Escherichia coli and Pseudomonas aeruginosa strains from south-east Asia. Several recent studies have reported that VEB-1-positive strains are highly resistant to ceftazidime, cefotaxime and aztreonam antibiotics. One strategy to overcome resistance involves administering antibiotics together with β-lactamase inhibitors during the treatment of infectious diseases. During this study, four VEB-1 β-lactamase inhibitors were identified using computer-aided drug design.

Methods: The SWISS-MODEL tool was utilized to generate three dimensional structures of VEB-1 β-lactamase, and the 3D model VEB-1 was verified using PROCHECK, ERRAT and VERIFY 3D programs. Virtual screening was performed by docking inhibitors obtained from the ZINC Database to the active site of the VEB-1 protein using AutoDock Vina software.

Results and conclusion: Homology modeling studies were performed to obtain a three-dimensional structure of VEB-1 β-lactamase. The generated model was validated, and virtual screening of a large chemical ligand library with docking simulations was performed using AutoDock software with the ZINC database. On the basis of the dock-score, four molecules were subjected to ADME/TOX analysis, with ZINC4085364 emerging as the most potent inhibitor of the VEB-1 β-lactamase.

Background: Schizophrenia is a neurodegenerative disorder that occurs worldwide and can be difficult to diagnose. It is the foremost neurological disorder leading to suicide among patients in both developed and underdeveloped countries. D-amino acid oxidase activator (DAOA), also known as G72, is directly implicated in the glutamateric hypothesis of schizophrenia. It activates D-amino acid oxidase, which oxidizes D-serine, leading to modulation of the N-methyl-D-aspartate receptor.

Methods: MODELLER (9v10) was utilized to generate three dimensional structures of the DAOA candidate gene. The HOPE server was used for mutational analysis. The Molecular Evolutionary Genetics Analysis (MEGA5) tool was utilized to reconstruct the evolutionary history of the candidate gene DAOA. AutoDock was used for protein-ligand docking and Gramm-X and PatchDock for protein-protein docking.

Results: A suitable template (1ZCA) was selected by employing BLASTp on the basis of 33% query coverage, 27% identity and E-value 4.9. The Rampage evaluation tool showed 91.1% favored region, 4.9% allowed region and 4.1% outlier region in DAOA. ERRAT demonstrated that the predicted model had a 50.909% quality factor. Mutational analysis of DAOA revealed significant effects on hydrogen bonding and correct folding of the DAOA protein, which in turn affect protein conformation. Ciona was inferred as the outgroup. Tetrapods were in their appropriate clusters with bifurcations. Human amino acid sequences are conserved, with chimpanzee and gorilla showing more than 80% homology and bootstrap value based on 1000 replications. Molecular docking analysis was employed to elucidate the binding mode of the reported ligand complex for DAOA. The docking experiment demonstrated that DAOA is involved in major amino acid interactions: the residues that interact most strongly with the ligand C28H28N3O5PS2 are polar but uncharged (Gln36, Asn38, Thr 122) and non-polar hydrophobic (Ile119, Ser171, Ser21, Ala31). Protein-protein docking simulation demonstrated two ionic bonds and one hydrogen bond involving DAOA. Lys-7 of the receptor protein interacted with Lys-163 and Asp-2037. Tyr-03 interacted with Arg-286 of the ligand protein and formed a hydrogen bond.

Conclusion: The predicted interactions might serve to inhibit the disease-related allele. It is assumed that current bioinformatics methods will contribute significantly to identifying, analyzing and curing schizophrenia. There is an urgent need to develop effective drugs for schizophrenia, and tools for examining candidate genes more accurately and efficiently are required.

Background: Lung cancer is the major cause of mortality worldwide. Major signalling pathways that could play significant role in lung cancer therapy include (1) Growth promoting pathways (Epidermal Growth Factor Receptor/Ras/ PhosphatidylInositol 3-Kinase) (2) Growth inhibitory pathways (p53/Rb/P14ARF, STK11) (3) Apoptotic pathways (Bcl-2/Bax/Fas/FasL). Insilico strategy was implemented to solve the mystery behind selected lung cancer pathway by applying comparative modeling and molecular docking studies.

Results: YASARA [v 12.4.1] was utilized to predict structural models of P16-INK4 and RB1 genes using template 4ELJ-A and 1MX6-B respectively. WHAT CHECK evaluation tool demonstrated overall quality of predicted P16-INK4 and RB1 with Z-score of -0.132 and -0.007 respectively which showed a strong indication of reliable structure prediction. Protein-protein interactions were explored by utilizing STRING server, illustrated that CDK4 and E2F1 showed strong interaction with P16-INK4 and RB1 based on confidence score of 0.999 and 0.999 respectively. In order to facilitate a comprehensive understanding of the complex interactions between candidate genes with their functional interactors, GRAMM-X server was used. Protein-protein docking investigation of P16-INK4 revealed four ionic bonds illustrating Arg47, Arg80,Cys72 and Met1 residues as actively participating in interactions with CDK4 while docking results of RB1 showed four hydrogen bonds involving Glu864, Ser567, Asp36 and Arg861 residues which interact strongly with its respective functional interactor E2F1.

Conclusion: This research may provide a basis for understanding biological insights of P16-INK4 and RB1 proteins which will be helpful in future to design a suitable drug to inhibit the disease pathogenesis as we have determined the interacting amino acids which can be targeted in order to design a ligand in-vitro to propose a drug for clinical trials. Protein -protein docking of candidate genes and their important interacting residues likely to be provide a gateway for developing computer aided drug designing.

The fields of molecular biology and computer science have cooperated over recent years to create a synergy between the cybernetic and biosemiotic relationship found in cellular genomics to that of information and language found in computational systems. Biological information frequently manifests its "meaning" through instruction or actual production of formal bio-function. Such information is called prescriptive information (PI). PI programs organize and execute a prescribed set of choices. Closer examination of this term in cellular systems has led to a dichotomy in its definition suggesting both prescribed data and prescribed algorithms are constituents of PI. This paper looks at this dichotomy as expressed in both the genetic code and in the central dogma of protein synthesis. An example of a genetic algorithm is modeled after the ribosome, and an examination of the protein synthesis process is used to differentiate PI data from PI algorithms.

Background: Based on a literature review, the current study aimed to construct mathematical models of lactate production and removal in both muscles and blood during steady state and at varying intensities during whole-body exercise. In order to experimentally test the models in dynamic situations, a cross-country skier performed laboratory tests while treadmill roller skiing, from where work rate, aerobic power and blood lactate concentration were measured. A two-compartment simulation model for blood lactate production and removal was constructed.

Results: The simulated and experimental data differed less than 0.5 mmol/L both during steady state and varying sub-maximal intensities. However, the simulation model for lactate removal after high exercise intensities seems to require further examination.

Conclusions: Overall, the simulation models of lactate production and removal provide useful insight into the parameters that affect blood lactate response, and specifically how blood lactate concentration during practical training and testing in dynamical situations should be interpreted.

Background: In clinical practice, the common strategy for immunotherapy of nasopharyngeal carcinoma (NPC) is to infuse cytotoxic T-lymphocyte (CTL) lines several times by intravenous injection, but it is difficult by laboratory research to investigate the relationship between treatment time-point, the amount of CTL added and the therapeutic effect. The objective of this study is to establish a mathematical model to study the therapeutic effect of different treatment time-points and amounts of CTL, and to predict the change in therapeutic effect when the percentage of EBV LMP2-specific CTL is increased from 10% to 20%.

Results: The concentration of epidermal growth factor receptor (EGFR) in the tumor cell cytomembranes increases after CTL is added. Concurrently, there is a marked downward trend of the phosphorylated transforming growth factor-α (TGFα)-EGFR complex in the tumor cell cytomembranes, which indicates restriction of tumor growth after CTL immunotherapy. The relationships among the time of addition of CTL, the amount of CTL added, different CTL specificities for LMP2 and the increment rate k of the total number of tumor cells were evaluated.

Conclusions: The simulation results quantify the relationships among treatment time-points, amount of CTL added, and the corresponding therapeutic effect of immunotherapy for NPC.

Background: Solid tumors are heterogeneous in composition. Cancer stem cells (CSCs) are believed to drive tumor progression, but the relative frequencies of CSCs versus non-stem cancer cells span wide ranges even within tumors arising from the same tissue type. Tumor growth kinetics and composition can be studied through an agent-based cellular automaton model using minimal sets of biological assumptions and parameters. Herein we describe a pivotal role for the generational life span of non-stem cancer cells in modulating solid tumor progression in silico.

Results: We demonstrate that although CSCs are necessary for progression, their expansion and consequently tumor growth kinetics are surprisingly modulated by the dynamics of the non-stem cancer cells. Simulations reveal that slight variations in non-stem cancer cell proliferative capacity can result in tumors with distinctly different growth kinetics. Longer generational life spans yield self-inhibited tumors, as the emerging population of non-stem cancer cells spatially impedes expansion of the CSC compartment. Conversely, shorter generational life spans yield persistence-limited tumors, with symmetric division frequency of CSCs determining tumor growth rate. We show that the CSC fraction of a tumor population can vary by multiple orders of magnitude as a function of the generational life span of the non-stem cancer cells.

Conclusions: Our study suggests that variability in the growth rate and CSC content of solid tumors may be, in part, attributable to the proliferative capacity of the non-stem cancer cell population that arises during asymmetric division of CSCs. In our model, intermediate proliferative capacities give rise to the fastest-growing tumors, resulting in self-metastatic expansion driven by a balance between symmetric CSC division and expansion of the non-stem cancer population. Our results highlight the importance of non-stem cancer cell dynamics in the CSC hypothesis, and may offer a novel explanation for the large variations in CSC fractions reported in vivo.

Background: With the development of experimental techniques and bioinformatics, the quantity of data available from protein-protein interactions (PPIs) is increasing exponentially. Functional modules can be identified from protein interaction networks. It follows that the investigation of functional modules will generate a better understanding of cellular organization, processes, and functions. However, experimental PPI data are still limited, and no modularity analysis of PPIs in pathogens has been published to date.

Results: In this study, we predict and analyze the functional modules of E. coli O157:H7 systemically by integrating several bioinformatics methods. After evaluation, most of the predicted modules are found to be biologically significant and functionally homogeneous. Six pathogenicity-related modules were discovered and analyzed, including novel modules. These modules provided new information on the pathogenicity of O157:H7. The modularity of cellular function and cooperativity between modules are also discussed. Moreover, modularity analysis of O157:H7 can provide possible candidates for biological pathway extension and clues for discovering new pathways of cross-talk.

Conclusions: This article provides the first modularity analysis of a pathogen and sheds new light on the study of pathogens and cellular processes. Our study also provides a strategy for applying modularity analysis to any sequenced organism.