PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “In silico modeling”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7Linked to original sources

Ligand recognition by drug-activated nuclear receptors PXR and CAR: structural, site-directed mutagenesis and molecular modeling studies.

Pregnane X receptor (PXR, NR1I2) and constitutive androstane receptor (CAR, NR1I3) are the principal regulators of drug/xenobiotic disposition and toxicity. These nuclear receptors display considerable cross-regulation of their target genes, and species-specific, yet promiscuous activation by a large number of structurally dissimilar ligands. Activation of PXR and/or CAR will frequently result in enhanced drug metabolism, disturbances in homeostasis of endogenous substances, and increased toxicity. Thus, understanding, measurement and prediction of ligand-elicited activation of PXR and CAR receptors is of utmost importance for the drug development process. In this mini-review, we will review the recent elucidation of structural properties of PXR and CAR, the molecular determinants of their ligand and species specificities and progress made in in silico models for identification of PXR and CAR activators.

Animals↗

Investigations into microsporidian methionine aminopeptidase type 2: a therapeutic target for microsporidiosis.

The Microsporidia have been reported to cause a wide range of clinical diseases particularly in patients that are immunosuppressed. They can infect virtually any organ system and cases of gastrointestinal infection, encephalitis, ocular infection, sinusitis, myositis and disseminated infection are well described in the literature. While benzimidazoles such as albendazole are active against many species of Microsporidia, these drugs do not have significant activity against Enterocytozoon bieneusi. Fumagillin, ovalicin and their analogues have been demonstrated to have antimicrosporidial activity in vitro and in animal models of microsporidiosis. Fumagillin has also been demonstrated to have efficacy in human infections due to E. bieneusi. Fumagillin is an irreversible inhibitor of methionine aminopeptidase type 2 (MetAP2). Homology cloning employing the polymerase chain reaction was used to identify the MetAP2 gene from the human pathogenic microsporidia Encephalitozoon cuniculi, Encephalitozoon hellem, Encephalitozoon intestinalis, Brachiola algerae and E. bieneusi. The full-length MetAP2 coding sequence was obtained for all of the Encephalitozoonidae. Recombinant E. cuniculi MetAP2 was produced in baculovirus and purified using chromatographic techniques. The in vitro activity and effect of the inhibitors bestatin and TNP-470 on this recombinant microsporidian MetAP2 was characterized. An in silico model of E. cuniculi MetAP2 was developed based on crystallographic data on human MetAP2. These reagents provide new tools for the development of in vitro assay systems to screen candidate compounds for use as new therapeutic agents for the treatment of microsporidiosis.

Amino Acid Sequence↗

Amino acid variations in hepatitis C virus p7 and sensitivity to antiviral combination therapy with amantadine in chronic hepatitis C.

BACKGROUND: Formation of transmembrane ion channels by hepatitis C virus (HCV) p7 and abrogation of channel function by amantadine was demonstrated in vitro. The relevance of HCV p7 amino acid (aa) variations for response to antiviral therapy with amantadine is unknown. METHODS: HCV p7 was sequenced in 86 individuals who were infected with HCV genotype 1. Thirty-six of 86 patients received amantadine within an interferon-alpha (IFN-alpha)-based antiviral therapy. Helical wheel modelling for HCV p7 was performed. RESULTS: No significant correlation of overall aa variations within HCV p7 was observed with response to IFN-alpha-based therapy with amantadine in HCV genotype 1alpha/b infected patients. When analysis was restricted to non-conservative aa variations, a higher number of aa substitutions within complete HCV p7 and transmembrane helix 2 was associated with non-response in HCV-1b-infected patients receiving therapy with amantadine (P=0.015 and P=0.037, respectively), without amantadine (P=0.106 and P=0.118, respectively), and in the total cohort of HCV-1b-infected patients (P=0.00007 and P=0.011, respectively). Furthermore, substitution L20F was observed more often in non-responders than responders with HCV-1b infection and therapy with amantadine (P=0.099). By in silico modelling, aa 20 was located toward the p7 channel lumen. Substitution L20F may impair amantadine action by altering the shape of the ion channel pore. CONCLUSION: Substitution L20F within HCV p7 may be associated with non-response to combination therapy specifically with amantadine in HCV-1b-infected patients. Non-responders with HCV-1b infection showed higher numbers of non-conservative aa variations within HCV p7 than responders, irrespective of the application of amantadine.

Amantadine↗

In silico classification of HERG channel blockers: a knowledge-based strategy.

The blockage of the hERG potassium channel by a wide number of diverse compounds has become a major pharmacological safety concern as it can lead to sudden cardiac death. In silico models can be potent tools to screen out potential hERG blockers as early as possible during the drug-discovery process. In this study, predictive models developed using the recursive partitioning method and created using diverse datasets from 203 molecules tested on the hERG channel are described. The first model was built with hERG compounds grouped into two classes, with a separation limit set at an IC50 value of 1 microm, and reaches an overall accuracy of 81%. The misclassification of molecules having a range of activity between 1 and 10 microM led to the generation of a tri-class model able to correctly classify high, moderate, and weak hERG blockers with an overall accuracy of 90%. Another model, constructed with the high and weak hERG-blocker categories, successfully increases the accuracy to 96%. The results reported herein indicate that a combination of precise, knowledge management resources and powerful modeling tools are invaluable to assessing potential cardiotoxic side effects related to hERG blockage.

Drug Design↗

Accurate computer-based design of a new backbone conformation in the second turn of protein L.

The rational design of loops and turns is a key step towards creating proteins with new functions. We used a computational design procedure to create new backbone conformations in the second turn of protein L. The Protein Data Bank was searched for alternative turn conformations, and sequences optimal for these turns in the context of protein L were identified using a Monte Carlo search procedure and an energy function that favors close packing. Two variants containing 12 and 14 mutations were found to be as stable as wild-type protein L. The crystal structure of one of the variants has been solved at a resolution of 1.9 A, and the backbone conformation in the second turn is remarkably close to that of the in silico model (1.1 A RMSD) while it differs significantly from that of wild-type protein L (the turn residues are displaced by an average of 7.2 A). The folding rates of the redesigned proteins are greater than that of the wild-type protein and in contrast to wild-type protein L the second beta-turn appears to be formed at the rate limiting step in folding.

Amino Acid Sequence↗

Efficacy of 2-halogen substituted D-glucose analogs in blocking glycolysis and killing "hypoxic tumor cells".

PURPOSE: Since 2-deoxy-D-glucose (2-DG) is currently in phase I clinical trials to selectively target slow-growing hypoxic tumor cells, 2-halogenated D-glucose analogs were synthesized for improved activity. Given the fact that 2-DG competes with D-glucose for binding to hexokinase, in silico modeling of molecular interactions between hexokinase I and these new analogs was used to determine whether binding energies correlate with biological effects, i.e. inhibition of glycolysis and subsequent toxicity in hypoxic tumor cells. METHODS AND RESULTS: Using a QSAR-like approach along with a flexible docking strategy, it was determined that the binding affinities of the analogs to hexokinase I decrease as a function of increasing halogen size as follows: 2-fluoro-2-deoxy-D-glucose (2-FG) > 2-chloro-2-deoxy-D-glucose (2-CG) > 2-bromo-2-deoxy-D-glucose (2-BG). Furthermore, D-glucose was found to have the highest affinity followed by 2-FG and 2-DG, respectively. Similarly, flow cytometry and trypan blue exclusion assays showed that the efficacy of the halogenated analogs in preferentially inhibiting growth and killing hypoxic vs. aerobic cells increases as a function of their relative binding affinities. These results correlate with the inhibition of glycolysis as measured by lactate inhibition, i.e. ID50 1 mM for 2-FG, 6 mM for 2-CG and > 6 mM for 2-BG. Moreover, 2-FG was found to be more potent than 2-DG for both glycolytic inhibition and cytotoxicity. CONCLUSIONS: Overall, our in vitro results suggest that 2-FG is more potent than 2-DG in killing hypoxic tumor cells, and therefore may be more clinically effective when combined with standard chemotherapeutic protocols.

Cell Hypoxia↗

Molecular imaging of EGFR kinase activity in tumors with 124I-labeled small molecular tracer and positron emission tomography.

Positron emission tomography (PET) with epidermal growth factor receptor (EGFR) kinase-specific radiolabeled tracers could provide the means for noninvasive and repetitive imaging of heterogeneity of EGFR expression and signaling activity in tumors in individual patients before and during therapy with EGFR signaling inhibitors. We developed the synthesis and (124)I-radiolabeling of the (E)-But-2-enedioic acid [4-(3-[(124)I]iodoanilino)-quinazolin-6-yl]-amide-(3-morpholin-4-yl-propyl)-amide (morpholino-[(124)I]-IPQA), which selectively, irreversibly, and covalently binds the adenosine-triphosphate-binding site to the activated (phosphorylated) EGFR kinase, but not to the inactive EGFR kinase. The latter was demonstrated using in silico modeling with crystal structures of the wild type and different gain-of-function mutants of EGFR kinases. Also, this was demonstrated by selective radiolabeling of the EGFR kinase domain with morpholino-[(131)I]-IPQA in A431 human epidermoid carcinoma cells and Western blot autoradiography. In vitro radiotracer accumulation and washout studies demonstrated a rapid accumulation and progressive retention postwashout of morpholino-[(131)I]-IPQA in A431 epidermoid carcinoma and in U87 human glioma cells genetically modified to express the EGFRvIII mutant receptor, but not in the wild-type U87MG glioma cells under serum-starved conditions. Using morpholino-[(124)I]-IPQA, we obtained noninvasive PET images of EGFR activity in A431 subcutaneous tumor xenografts, but not in subcutaneous tumor xenografts grown from K562 human chronic myeloid leukemia cells in immunocompromised rats and mice. Based on these observations, we suggest that PET imaging with morpholino-[(124)I]-IPQA should allow for identification of tumors with high EGFR kinase signaling activity, including brain tumors expressing EGFRvIII mutants and nonsmall-cell lung cancer expressing gain-of-function EGFR kinase mutants. Because of significant hepatobiliary clearance and intestinal reuptake of the morpholino-[(124)I]-IPQA, additional [(124)I]-IPQA derivatives with improved water solubility may be required to optimize the pharmacokinetics of this class of molecular imaging agents.

Animals↗

Hydrophilically enhanced 3-carboranyl thymidine analogues (3CTAs) for boron neutron capture therapy (BNCT) of cancer.

Five novel 3-carboranyl thymidine analogues (3CTAs) were designed and synthesized for boron neutron capture therapy (BNCT) of cancer. Phosphorylation of all five 3CTAs was catalyzed by recombinant human thymidine kinase (hTK1) using adenosine triphosphate (ATP) as the phosphate donor. The obtained phosphorylation rates ranged from 4% to 64.5% relative to that of thymidine. The compound with the most favorable hTK1 binding properties had a k(cat)/K(M) value of 57.4% relative to that of thymidine and an IC(50) of inhibition of thymidine phosphorylation by hTK1 of 92 microM. Among the five synthesized 3CTAs, this agent had also the overall most favorable physicochemical properties. Therefore, it may have the potential to replace N5-2OH, the current lead 3CTA, in preclinical studies. An in silico model for the binding of this compound to hTK1 was developed.

Binding Sites↗

Chemical modification of the naphthoyl 3-position of JWH-015: in search of a fluorescent probe to the cannabinoid CB2 receptor.

In silico modelling was used to guide the positioning of the fluorescent dye NBD-F on the cannabinoid CB2 receptor agonist JWH-015. While the ultimate fluorescent conjugate lost extensive binding affinity to the cannabinoid CB2 receptor, affinity and efficacy studies on the naphthoyl 3-position modified precursor molecules have provided new insight into structure-activity relationships associated with this position.

Binding, Competitive↗

Reconstruction of microbial transcriptional regulatory networks.

Although metabolic networks can be readily reconstructed through comparative genomics, the reconstruction of regulatory networks has been hindered by the relatively low level of evolutionary conservation of their molecular components. Recent developments in experimental techniques have allowed the generation of vast amounts of data related to regulatory networks. This data together with literature-derived knowledge has opened the way for genome-scale reconstruction of transcriptional regulatory networks. Large-scale regulatory network reconstructions can be converted to in silico models that allow systematic analysis of network behavior in response to changes in environmental conditions. These models can further be combined with genome-scale metabolic models to build integrated models of cellular function including both metabolism and its regulation.

Bacteria↗

Genome semantics, in silico multicellular systems and the Central Dogma.

Genomes with their complexity and size present what appears to be an impossible challenge. Scientists speak in terms of decades or even centuries before we will understand how genomes and their hosts the cell and the city of cells that make up the multicellular context function. We believe that there will be surprisingly quick progress made in our understanding of genomes. The key is to stop taking the Central Dogma as the only direction in which genome research can scale the semantics of genomes. Instead a top-down approach coupled with a bottom-up approach may snare the unwieldy beast and make sense of genomes. The method we propose is to take in silico biology seriously. By developing in silico models of genomes cells and multicellular systems, we position ourselves to develop a theory of meaning for artificial genomes. Then using that develop a natural semantics of genomes.

Cells↗

Elucidating the clinical and genetic spectrum of inositol polyphosphate phosphatase INPP4A-related neurodevelopmental disorder.

PURPOSE: Biallelic INPP4A variants have recently been associated with severe neurodevelopmental disease in single-case reports. Here, we expand and elucidate the clinical-genetic spectrum and provide a pathomechanistic explanation for genotype-phenotype correlations. METHODS: Clinical and genomic investigations of 30 individuals were undertaken alongside molecular and in silico modelling and translation reinitiation studies. RESULTS: We characterize a clinically variable disorder with cardinal features, including global developmental delay, severe-profound intellectual disability, microcephaly, limb weakness, cerebellar signs, and short stature. A more severe presentation associated with biallelic INPP4A variants downstream of exon 4 has additional features of (ponto)cerebellar hypoplasia, reduced cerebral volume, peripheral spasticity, contractures, intractable seizures, and cortical visual impairment. Our studies identify the likely pathomechanism of this genotype-phenotype correlation entailing translational reinitiation in exon 4 resulting in an N-terminal truncated INPP4A protein retaining partial functionality, associated with less severe disease. We also identified identical reinitiation site conservation in Inpp4a-/- mouse models displaying similar genotype-phenotype correlation. Additionally, we show fibroblasts from a single affected individual exhibit disrupted endocytic trafficking pathways, indicating the potential biological basis of the condition. CONCLUSION: Our studies comprehensively characterize INPP4A-related neurodevelopmental disorder and suggest genotype-specific clinical assessment guidelines. We propose that the potential mechanistic basis of observed genotype-phenotype correlations entails exon 4 translation reinitiation.

Humans↗

Rare variants in MIR184 are a novel genetic cause of Fuchs endothelial corneal dystrophy.

PURPOSE: To identify novel genetic causes of Fuchs endothelial corneal dystrophy (FECD) within a genetically unsolved patient cohort lacking repeat expansions in the TCF4 gene (Exp-). METHODS: A rare variant analysis framework (CoCoRV) was applied to exome data, in combination with in silico modeling, luciferase reporter, and RNA-seq analysis to characterize transcriptome-wide consequences of identified variants. RESULTS: A gene burden analysis identified MIR184, a microRNA encoding gene, to be enriched for rare pathogenic variants within the studied Exp- FECD cohort. In total, 2 noncoding rare variants were identified in 4 unrelated FECD probands: NR_029705.1:n.58G>A and n.73G>T. Both variants altered highly conserved mature sequence residues, were predicted to induce hairpin structural changes, and were experimentally determined to disrupt microRNA-mRNA interactions. RNA-seq of transfected human corneal endothelial cells revealed that the mutants elicited distinct transcriptomic profiles. Enriched KEGG pathways included PI3K-Akt signaling, focal adhesion, and immune response, revealing shared pathogenic mechanisms between MIR184-associated FECD and the more common TCF4 repeat expansion-mediated form of disease. CONCLUSION: MIR184 variants are a novel rare genetic cause of FECD, and common pathways of transcriptomic dysregulation are shared across genetically distinct subtypes of the disease. These pathways may serve as future gene agnostic targets for therapeutic interventions.

Humans↗

Crystal structure of the phosphoenolpyruvate-binding enzyme I-domain from the Thermoanaerobacter tengcongensis PEP: sugar phosphotransferase system (PTS).

Enzyme I (EI), the first component of the phosphoenolpyruvate (PEP):sugar phosphotransferase system (PTS), consists of an N-terminal protein-binding domain (EIN) and a C-terminal PEP-binding domain (EIC). EI transfers phosphate from PEP by double displacement via a histidine residue on EIN to the general phosphoryl carrier protein HPr. Here, we report the 1.82A crystal structure of the homodimeric EIC domain from Thermoanaerobacter tengcongensis, a saccharolytic eubacterium that grows optimally at 75 degrees C. EIC folds into a (betaalpha)(8) barrel with three large helical insertions between beta2/alpha2, beta3/alpha3 and beta6/alpha6. The large amphipathic dimer interface buries 3750A(2) of accessible surface area per monomer. A comparison with pyruvate phosphate dikinase (PPDK) reveals that the active-site residues in the empty PEP-binding site of EIC and in the liganded PEP-binding site of PPDK have almost identical conformations, pointing to a rigid structure of the active site. In silico models of EIC in complex with the Z and E-isomers of chloro-PEP provide a rational explanation for their difference as substrates and inhibitors of EI. The EIC domain exhibits 54% amino acid sequence identity with Escherichia coli and 60% with Bacillus subtilis EIC, has the same amino acid composition but contains additional salt-bridges and a more complex salt-bridge network than the homology model of E.coli EIC. The easy crystallization of EIC suggests that T.tengcongensis can serve as source for stable homologs of mesophilic proteins that are too labile for crystallization.

Binding Sites↗

Disulfide bridges in the mesophilic triosephosphate isomerase from Giardia lamblia are related to oligomerization and activity.

Triosephosphate isomerase from the mesophile Giardia lamblia (GlTIM) is the only known TIM with natural disulfide bridges. We previously found that oxidized and reduced thiol states of GlTIM are involved in the interconversion between native dimers and higher oligomeric species, and in the regulation of enzymatic activity. Here, we found that trophozoites and cysts have different oligomeric species of GlTIM and complexes of GlTIM with other proteins. Our data indicate that the internal milieu of G. lamblia is favorable for the formation of disulfide bonds. Enzyme mutants of the three most solvent exposed Cys of GlTIM (C202A, C222A, and C228A) were prepared to ascertain their contribution to oligomerization and activity. The data show that the establishment of a disulfide bridge between two C202 of two dimeric GlTIMs accounts for multimerization. In addition, we found that the establishment of an intramonomeric disulfide bond between C222 and C228 abolishes catalysis. Multimerization and inactivation are both reversed by reducing conditions. The 3D structure of the C202A GlTIM was solved at 2.1 A resolution, showing that the environment of the C202 is prone to hydrophobic interactions. Molecular dynamics of an in silico model of GlTIM when the intramonomeric disulfide bond is formed, showed that S216 is displaced 4.6 A from its original position, causing loss of hydrogen bonds with residues of the active-site loop. This suggests that this change perturb the conformational state that aligns the catalytic center with the substrate, inducing enzyme inactivation.

Animals↗

A basic peptide within the juxtamembrane region is required for EGF receptor dimerization.

The epidermal growth factor receptor (EGFR) is fundamental for normal cell growth and organ development, but has also been implicated in various pathologies, notably tumors of epithelial origin. We have previously shown that the initial 13 amino acids (P13) within the intracellular juxtamembrane region (R645-R657) are involved in the interaction with calmodulin, thus indicating an important role for this region in EGFR function. Here we show that P13 is required for proper dimerization of the receptor. We expressed either the intracellular domain of EGFR (TKJM) or the intracellular domain lacking P13 (DeltaTKJM) in COS-7 cells that express endogenous EGFR. Only TKJM was immunoprecipitated with an antibody directed against the extracellular part of EGFR, and only TKJM was tyrosine phosphorylated by endogenous EGFR. Using SK-N-MC cells, which do not express endogenous EGFR, that were stably transfected with either wild-type EGFR or recombinant full-length EGFR lacking P13 demonstrated that P13 is required for appropriate receptor dimerization. Furthermore, mutant EGFR lacking P13 failed to be autophosphorylated. P13 is rich in basic amino acids and in silico modeling of the EGFR in conjunction with our results suggests a novel role for the juxtamembrane domain (JM) of EGFR in mediating intracellular dimerization and thus receptor kinase activation and function.

Amino Acid Sequence↗

Somatic mutations and ageing in silico.

Considerable evidence points to an accumulation of somatic mutations in older cells and organisms but the causal role of mutations in the ageing process is still unclear. In addition to demonstrating that mutations accumulate, it is important to address the question of whether they do so at a sufficient rate and with a dynamic profile that is consistent with them playing a causative role. We describe the development of in silico models that can be used to explore the role of somatic mutations in ageing and which form a part of a growing effort to build predictive mathematical and computer models that can help unravel the complexity of the functional genomics of ageing. Our models address, in particular, how mutations affect populations of dividing cells like human fibroblasts, in which the challenge to the somatic mutation theory is greatest, since selection at the cellular level will tend to suppress the accumulation of mutations.

Aging↗

An Acetyltransferase Conferring Self-Resistance of the Producer to Lasso Peptide Antibiotic Lariocidin.

The soil microbiome, a reservoir of antibiotic-producing bacteria, also harbors resistance determinants encoded within antibiotic biosynthetic gene clusters (BGCs). Studying self-resistance mechanisms, which have evolved in producers to protect against their own toxic metabolites, provides critical insights into the evolution of resistance and the potential vulnerabilities of new antibiotics and can facilitate the production of natural products in heterologous hosts. Here, we describe the self-resistance mechanism to lariocidin (LAR), a recently discovered lasso peptide antibiotic that inhibits the ribosomal machinery and exhibits antibacterial activity against key pathogens. We identified and characterized an N-acetyltransferase enzyme (LrcE) encoded within the LAR BGC that mediates self-resistance in LAR-producing Paenibacillus sp. M2. LrcE is a member of the GCN5-related N-acetyltransferase (GNAT) superfamily and performs site-specific acetylation of LAR at a critical lysine residue. This modification disrupts ribosomal binding, thereby reducing LAR's antibacterial activity. Using in silico modeling, we predicted a conserved acetyl-CoA-binding motif and an LAR-binding region on LrcE. Bioinformatic analysis revealed LrcE homologues in environmental but not clinically relevant pathogens, suggesting a limited risk of horizontal gene transfer and, therefore, supporting the further development of LAR as a next-generation antibiotic.

Anti-Bacterial Agents↗