Search

Alterations to intestinal microbiota are assumed to occur in the pathogenesis of inflammatory bowel disease (IBD). This study aims to analyze the association of fecal microbiota composition, body composition, and lipid characteristics in patients with Crohn's disease (CD). In our cross-sectional study, patients with CD were enrolled and blood and fecal samples were collected. Clinical and endoscopic disease activity and body composition were assessed and laboratory tests were made. Fecal bacterial composition was analyzed using the shotgun method. Microbiota alterations based on obesity, lipid parameters, and disease characteristics were analyzed. In this study, 27 patients with CD were analyzed, of which 37.0% were obese based on visceral fat area (VFA). Beta diversities were higher in non-obese patients ( p < 0.001), but relative abundances did not differ. C. innocuum had a higher abundance at a high cholesterol level than Bacillota ( p = 0.001, p = 0.0034). Adlercreutzia , B. longum , and Blautia alterations were correlated with triglyceride levels. Higher Clostridia ( p = 0.009) and B. schinkii ( p = 0.032) and lower Lactobacillus ( p = 0.035) were connected to high VFA. Disease activity was coupled with dysbiotic elements. Microbiota alterations in obesity highlight the importance of gut microbiota in diseases with a similar inflammatory background and project therapeutic options.

Background and Aims: Crohn's disease [CD] and ulcerative colitis [UC] require lifelong treatment and patient monitoring. Current biomarkers have several limitations; therefore, there is an unmet need to identify novel biomarkers in inflammatory bowel disease [IBD]. Previously, the role of plasminogen activator inhibitor 1 [PAI-1] was established in the pathogenesis of IBD and suggested as a potential biomarker. Therefore, we aimed to comprehensively analyse the selectivity of PAI-1 in IBD, its correlation with disease activity, and its potential to predict therapeutic response., Methods: Blood, colon biopsy, organoid cultures [OC], and faecal samples were used from active and inactive IBD patients and control subjects. Serpin E1 gene expressions and PAI-1 protein levels and localisation in serum, biopsy, and faecal samples were evaluated by qRT-PCR, ELISA, and immunostaining, respectively., Results: The study population comprised 132 IBD patients [56 CD and 76 UC] and 40 non-IBD patients. We demonstrated that the serum, mucosal, and faecal PAI-1 concentrations are elevated in IBD patients, showing clinical and endoscopic activity. In responders [decrease of eMayo ≥3 in UC; or SES-CD  50% in CD], the initial PAI-1 level decreased significantly upon successful therapy. OCs derived from active IBD patients produced higher concentrations of PAI-1 than the controls, suggesting that epithelial cells could be a source of PAI-1. Moreover, faecal PAI-1 selectively increases in active IBD but not in other organic gastrointestinal diseases., Conclusions: The serum, mucosal, and faecal PAI-1 concentration correlates with disease activity and therapeutic response in IBD, suggesting that PAI-1 could be used as a novel, non-invasive, disease-specific, faecal biomarker in patient follow-up., (© The Author(s) 2023. Published by Oxford University Press on behalf of European Crohn’s and Colitis Organisation.)

Patients with recurrent acute pancreatitis (RAP) are at significant risk of developing early chronic pancreatitis (CP), which progresses into irreversible, end-stage CP with severe symptoms. There is no specific therapy in RAP or in early CP that may hinder disease progression. The pathogenesis of CP is complex and involves interactions among multiple cell types, including pancreatic acinar, ductal, and stellate cells (PSC). Therefore, it is pivotal to identify common pathogenic pathways in these cells that could be targeted pharmacologically. The Orai1-mediated store-operated Ca2+ entry (SOCE) is a ubiquitous signaling mechanism that may become overactivated in pathological states resulting in intracellular Ca2+ overload. In this study, we used ex vivo and in vivo preclinical disease models to demonstrate that Orai1 inhibition prevents progression of RAP and early CP. The selective Orai1 inhibitor CM5480 restored the expression of SOCE-associated regulatory factor in acinar cells, prevented uncontrolled Ca2+ elevation, protected acinar and ductal functions, mitigated immune cell infiltration, and diminished PSC activation, proliferation, and migration. We suggest that the overactivation of Orai1 is a crucial pathogenetic event in the progression of early CP and that inhibition of Orai1 could prevent the development of end-stage CP.

Epithelial ion and fluid secretion determine the physiological functions of a broad range of organs, such as the lung, liver, or pancreas. The molecular mechanism of pancreatic ion secretion is challenging to investigate due to the limited access to functional human ductal epithelia. Patient-derived organoids may overcome these limitations, however direct accessibility of the apical membrane is not solved. In addition, due to the vectorial transport of ions and fluid the intraluminal pressure in the organoids is elevated, which may hinder the study of physiological processes. To overcome these, we developed an advanced culturing method for human pancreatic organoids based on the removal of the extracellular matrix that induced an apical-to-basal polarity switch also leading to reversed localization of proteins with polarized expression. The cells in the apical-out organoids had a cuboidal shape, whereas their resting intracellular Ca 2+ concentration was more consistent compared to the cells in the apical-in organoids. Using this advanced model, we demonstrated the expression and function of two novel ion channels, the Ca 2+ activated Cl - channel Anoctamin 1 (ANO1) and the epithelial Na + channel (ENaC), which were not considered in ductal cells yet. Finally, we showed that the available functional assays, such as forskolin-induced swelling, or intracellular Cl - measurement have improved dynamic range when performed with apical-out organoids. Taken together our data suggest that polarity-switched human pancreatic ductal organoids are suitable models to expand our toolset in basic and translational research., (© 2023. The Author(s).)

Background: In patients with inflammatory bowel disease (IBD), Crohn's disease (CD), and ulcerative colitis (UC), numerous cases of exacerbations could be observed after colonoscopy, raising the possible pathogenetic effect of colonic microbiota alterations in IBD flare., Objectives: We aimed to investigate the changes in the fecal microbiota composition in IBD patients influenced by the bowel preparation with sodium picosulfate., Design: We enrolled patients with IBD undergoing bowel preparation for colonoscopy in the prospective cohort study. The control group (Con) comprised non-IBD patients who underwent colonoscopy. Clinical data, blood, and stool samples were collected before colonoscopy (timepoint A), 3 days later (timepoint B), and 4 weeks later (timepoint C)., Methods: Disease activity and gut microbiota changes were assessed at each timepoint. Fecal microbiota structure - at family level - was determined by sequencing the V4 region of the 16S rRNA gene. Statistical analysis included differential abundance analysis and Mann-Whitney tests., Results: Forty-one patients (9 CD, 13 UC, and 19 Con) were included. After bowel preparation, alpha diversity was lower in the CD group than in the UC ( p  = 0.01) and Con ( p  = 0.02) groups at timepoint B. Alpha diversity was significantly higher in the UC group than in the CD and Con ( p  = 0.03) groups at timepoint C. Beta diversity difference differed between the IBD and Con ( p  = 0.001) groups. Based on the differential abundance analysis, the Clostridiales family was increased, whereas the Bifidobacteriaceae family was decreased in CD patients compared to the Con at timepoint B., Conclusions: Bowel preparation may change the fecal microbial composition in IBD patients, which may have a potential role in disease exacerbation after bowel cleansing., Competing Interests: TM has received speaker’s honoraria from MSD, AbbVie, Egis, Goodwill Pharma, Takeda, Pfizer, Janssen, Sandoz, MundiPharma, Phytotec, Roche, Fresenius, and Teva. KF has received speaker’s honoraria from AbbVie, Janssen, Ferring, Takeda, and Goodwill Pharma. PB, MR, AA, TR, KSZ, BJ, AB, EA, WA, BK, TF, DP, RB, AF, JM, and ZSZ have no conflict of interest to disclose., (© The Author(s), 2023.)

Cationic antimicrobial peptide PGLa gets into close contact with the anionic bacterial cell membrane, facilitating cross-membrane transport phenomena and membrane disruption depending on the concentration. The mechanisms of action are closely associated with the tilted insertion geometry of PGLa. Therefore, we aimed to understand the interaction between the transmembrane potential (TMP) and the orientation of the membrane-bound PGLa helix. Molecular dynamics simulations were performed with TMP, and we found that the PGLa tilt angle relative to the membrane is coupled with the TMP. Elevated TMP increases the population of the tilted state. We observed positive feedback between the tilt angle and the TMP, which occurs due to the electrostatic interaction between the peptidic helix and the Na + cations at the membrane-water interface. These TMP coupled phenomena can contribute to understanding the direct antimicrobial and adjuvant effects of PGLa in combination with regular antibiotics.

Introduction: The relationship between clinical outcomes and serum anti-TNF levels is controversial. The aim of this study was to perform simultaneous analyses of serum, mucosal, and fecal anti-TNF-α levels., Methods: Consecutive IBD patients who received maintenance anti-TNF-α therapy were enrolled. The number of TNF-α positive cells in the mucosa was detected using immunofluorescent labeling on biopsy samples. Serum, mucosal and fecal anti-TNF-α, serum anti-drug antibody, and fecal calprotectin levels were determined using ELISA. Each patient underwent body composition analysis as well., Results: Data of 50 patients were analyzed. The number TNF-α positive cells was significantly higher in the inflamed part of the colon than in the un-inflamed part of the colon. Tissue and fecal drug levels did not show any association with serum drug levels; moreover, serum anti-TNF concentration did not correlate with endoscopic activity. Mucosal anti-TNF levels were higher only in IFX-treated patients in remission and IFX-treated patients with detectable fecal anti-TNF had lower tissue drug levels. Presence of the drug in the feces was significantly different according to disease activity., Conclusion: Fecal drug concentration is suggested to be a better predictor of endoscopic activity and loss of response, and fecal drug monitoring may improve the estimation accuracy of tissue drug levels.

Metabotropic glutamate receptor 5 (mGluR5) is a class C G protein-coupled receptor (GPCR) with both an extracellular ligand binding site and an allosteric intrahelical chamber located similarly to the orthosteric ligand binding site of Class A GPCRs. Ligands binding to this ancestral site of mGluR5 can act as positive (PAM), negative (NAM) or silent (SAM) allosteric modulators, and their medicinal chemistry optimization is notoriously difficult, as subtle structural changes may cause significant variation in activity and switch in the functional response. Here we present all atom molecular dynamics simulations of NAM, SAM and PAM complexes formed by closely related ligands and analyse the structural differences of the complexes. Several residues involved in the activation are identified and the formation of a continuous water channel in the active complex but not in the inactive ones is recognized. Our results suggest that the mechanism of mGluR5 activation is similar to that of class A GPCRs.Communicated by Ramaswamy H. Sarma.

EDTA is commonly used as an efficient chelator of metal ion enzyme cofactors. It is highly soluble, optically inactive and does not interfere with most chemicals used in standard buffers making EDTA a common choice to generate metal-free conditions for biochemical and biophysical investigations. However, the controversy in the literature on metal-free enzyme activities achieved using EDTA or by other means called our attention to a putative effect of EDTA beyond chelation. Here, we show that EDTA competes for the nucleotide binding site of the nucleotide hydrolase dUTPase by developing an interaction network within the active site similar to that of the substrate. To achieve these findings, we applied kinetics and molecular docking techniques using two different dUTPases. Furthermore, we directly measured the binding of EDTA to dUTPases and to two other dNTPases, the Taq polymerase and MutT using isothermal titration calorimetry. EDTA binding proved to be exothermic and mainly enthalpy driven with a submicromolar dissociation constant considerably lower than that of the enzyme:substrate or the Mg:EDTA complexes. Control proteins, including an ATPase, did not interact with EDTA. Our findings indicate that EDTA may act as a selective inhibitor against dNTP hydrolyzing enzymes and urge the rethinking of the utilization of EDTA in enzymatic experiments., Competing Interests: The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

1,4-Dioxane is a cytotoxic B2-type human carcinogen, a serious water pollutant produced solely by industrial activity. The effect of 1,4-dioxane on phospholipid membrane models composed of dipalmitoyl-phosphatidylcholine (DPPC) and its branched isomer (isodipalmitoyl-phosphatidylcholine, IPPC) was investigated using MD simulations. Clear and polluted membranes were compared by membrane parameters such as area per lipid (APL), volume per lipid (VPL), compressibility modulus, membrane thickness, and orderliness of lipid tails. While neat systems significantly differ from each other, the presence of the pollutant has the same effect on both types of lipid membranes. High density of dioxane appears in the vicinity of ester groups, which pushes away lipid headgroups from each other, leading to an overall change in lipid structure: APL and VPL grows, while the orderliness of lipid tails, membrane thickness, and compressibility modulus decrease. Orientational preferences of water and dioxane molecules were also investigated and different membrane regions have been specified according to the stance of water molecules. Free-energy profile for 1,4-dioxane penetration mechanism into DPPC membranes was carried out using metadynamics for two different concentrations of the pollutant ( c 1 = 7.51 g/dm 3 , c 2 = 75.10 g/dm 3 ), which showed that the higher the concentration is, the lower the free energy of penetration gets. Only a small free-energy barrier was found in the headgroup region and accumulation of dioxane is thermodynamically unfavored in the middle of the bilayer. The penetration mechanism has been described in detail based on the orientational preference of 1,4-dioxane molecules and the free-energy profiles.

In order to study the effects of peptide exposure to oxidative attack, we chose a model reaction in which the hydroxyl radical discretely abstracts a hydrogen atom from the α-carbon of each residue of a highly amyloidogenic region in the human calcitonin hormone, hCT 15-19 . Based on a combined Molecular Mechanics / Quantum Mechanics approach, the extended and folded L- and D-configuration and radical intermediate hCT 15-19 peptides were optimized to obtain their compactness, secondary structure and relative thermodynamic data. The results suggest that the epimerization of residues is generally an exergonic process that can explain the cumulative nature of molecular aging. Moreover, the configurational inversion induced conformational changes can cause protein dysfunction. The epimerization of the central residue to the D-configuration induced a hairpin structure in hCT 15-19 , concomitant with a possible oligomerization of human calcitonin into Aβ(1-42)-like amyloid fibrils present in patients suffering from Alzheimer's disease., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Computer simulations of four lipid membranes of different compositions, namely neat DPPC and PSM, and equimolar DPPC-cholesterol and PSM-cholesterol mixtures, are performed in the presence and absence of the general anesthetics diethylether and sevoflurane both at 1 and 600 bar. The results are analyzed in order to identify membrane properties that are potentially related to the molecular mechanism of anesthesia, namely that change in the same way in any membrane with any anesthetics, and change oppositely with increasing pressure. We find that the lateral lipid density satisfies both criteria: it is decreased by anesthetics and increased by pressure. This anesthetic-induced swelling is attributed to only those anesthetic molecules that are located close to the boundary of the apolar phase. This lateral expansion is found to lead to increased lateral mobility of the lipids, an effect often thought to be related to general anesthesia; to an increased fraction of the free volume around the outer preferred position of anesthetics; and to the decrease of the lateral pressure in the nearby range of the ester and amide groups, a region into which anesthetic molecules already cannot penetrate. All these changes are reverted by the increase of pressure. Another important finding of this study is that cholesterol has an opposite effect on the membrane properties than anesthetics, and, correspondingly, these changes are less marked in the presence of cholesterol. Therefore, changes in the membrane that can lead to general anesthesia are expected to occur in the membrane domains of low cholesterol content., (Copyright © 2018 Elsevier B.V. All rights reserved.)

DJ-1 (PARK7) is a multifunctional protein linked to the onset and progression of a number of diseases, most of which are associated with high oxidative stress. The Cys106 of DJ-1 is unusually reactive and thus sensitive to oxidation, and due to high oxidative stress it was observed to be in various oxidized states in disease condition. The oxidation state of Cys106 of DJ-1 is believed to determine the specific functions of the protein in normal and disease conditions. Here we report molecular dynamics simulation and biophysical experimental studies on DJ-1 in reduced (Cys106, S - ), oxidized (Cys106, SO 2 - ), and over-oxidized (Cys106, SO 3 - ) states. To simulate the different oxidation states of Cys106 in DJ-1, AMBER related force field parameters were developed and reported for 3-sulfinoalanine and cysteine sulfonic acid. Our studies found that the overall structure of DJ-1 in different oxidation states was similar globally, while it differed locally significantly, which have implications on its stability, function and its link to disease on-set. Importantly, the results suggest that over-oxidation may trigger loss of functions due to local structural modification in the Cys106 containing pocket of DJ-1 and structurally destabilize the dimeric state of DJ-1, which is believed to be its bioactive conformation. Such loss of functions would result in reduced ability of DJ-1 to protect from oxidative stress insults and may lead to increased progression of disease., (Copyright © 2017 Elsevier B.V. All rights reserved.)

dUTPases catalyze the hydrolysis of dUTP into dUMP and pyrophosphate to maintain the proper nucleotide pool for DNA metabolism. Recent evidence suggests that dUTPases may also represent a selective drug target in mycobacteria because of the crucial role of these enzymes in maintaining DNA integrity. Nucleotide-hydrolyzing enzymes typically harbor a buried ligand-binding pocket at interdomain or intersubunit clefts, facilitating proper solvent shielding for the catalyzed reaction. The mechanism by which substrate binds this hidden pocket and product is released in dUTPases is unresolved because of conflicting crystallographic and spectroscopic data. We sought to resolve this conflict by using a combination of random acceleration molecular dynamics (RAMD) methodology and structural and biochemical methods to study the dUTPase from Mycobacterium tuberculosis In particular, the RAMD approach used in this study provided invaluable insights into the nucleotide dissociation process that reconciles all previous experimental observations. Specifically, our data suggest that nucleotide binding takes place as a small stretch of amino acids transiently slides away and partially uncovers the active site. The in silico data further revealed a new dUTPase conformation on the pathway to a relatively open active site. To probe this model, we developed the Trp 21 reporter and collected crystallographic, spectroscopic, and kinetic data that confirmed the interaction of Trp 21 with the active site shielding C-terminal arm, suggesting that the RAMD method is effective. In summary, our computational simulations and spectroscopic results support the idea that small loop movements in dUTPase allow the shuttlingof the nucleotides between the binding pocket and the solvent., (© 2016 by The American Society for Biochemistry and Molecular Biology, Inc.)

In order to identify molecular models of the human 5-HT6 receptor suitable for virtual screening, homology modeling and membrane-embedded molecular dynamics simulations were performed. Structural requirements for robust enrichment were assessed by an unbiased chemometric analysis of enrichments from retrospective virtual screening studies. The two main structural features affecting enrichment are the outward movement of the second extracellular loop and the formation of a hydrophobic cavity deep in the binding site. These features appear transiently in the trajectories and furthermore the stretches of uniformly high enrichment may only last 4-10 ps. The formation of the inner hydrophobic cavity was also linked to the active-like to inactive-like transition of the receptor, especially the so-called connector region. The best structural models provided significant and robust enrichment over three independent ligand sets.

The energetics of peptide bond formation is an important factor not only in the design of chemical peptide synthesis, but it also has a role in protein biosynthesis. In this work, quantum chemical calculations at 10 different levels of theory including G3MP2B3 were performed on the energetics of glutathione formation. The strength of the peptide bond is found to be closely related to the acid strength of the to-be N-terminal and the basicity of the to-be C-terminal amino acid. It is shown that the formation of the first peptide activates the amino acid for the next condensation step, manifested in bacterial protein synthesis where the first step is the formation of an N-formylmethionine dipeptide. The possible role of glutathione in prebiotic molecular evolution is also analyzed. The implications of the thermodynamics of peptide bond formation in prebiotic peptide formation as well as in the preference of α- instead of β- or γ-amino acids are discussed. An empirical correction is proposed for the compensation of the error due to the incapability of continuum solvation models in describing the change of the first solvation shell when a peptide bond is formed from two zwitterions accompanied by the disappearance of one ion pair.

Micellisation process for sodium dodecyl sulphate and sodium cholate in 1∶1 molar ratio was investigated in a combined approach, including several experimental methods and coarse grained molecular dynamics simulation. The critical micelle concentration (cmc) of mixed micelle was determined by spectrofluorimetric and surface tension measurements in the temperature range of 0-50°C and the values obtained agreed with each other within the statistical error of the measurements. In range of 0-25°C the cmc values obtained are temperature independent while cmc values were increased at higher temperature, which can be explained by the intensive motion of the monomers due to increased temperature. The evidence of existing synergistic effect among different constituent units of the micelle is indicated clearly by the interaction parameter (β1,2) calculated from cmc values according to Rubingh. As the results of the conductivity measurements showed the negative surface charges of the SDS-NaCA micelle are not neutralized by counterions. Applying a 10 µs long coarse-grained molecular dynamics simulation for system including 30-30 SDS and CA (with appropriate number of Na+ cations and water molecules) we obtained semi-quantitative agreement with the experimental results. Spontaneous aggregation of the surfactant molecules was obtained and the key steps of the micelle formation are identified: First a stable SDS core was formed and thereafter due to the entering CA molecules the size of the micelle increased and the SDS content decreased. In addition the size distribution and composition as well as the shape and structure of micelles are also discussed.

Bursopentin (BP(5), H-Cys(1)-Lys(2)-Arg(3)-Val(4)-Tyr(5)-OH), found in the bursa Fabricius of the chicken, is a pentapeptide that protects the organism from oxidative stress by reducing the intracellular generation of reactive oxygen species. Hydrogen abstraction, a common oxidative reaction occurring in proteins, often results in the formation of d amino acid residues. To study the effect of this phenomenon on the structure of bursopentin, each of its residues were converted from the l configuration to the d configuration, and the structures of these peptide epimers were compared to that of the wild-type bursopentin. The conformations, secondary structures, compactness and hydrogen bonding of bursopentin were compared to its epimers using molecular dynamics simulations and first principles quantum chemical computations. It was discovered that the repulsion between the side chains of Lys(2) and Arg(3) influenced the conformation of the peptide regardless of the configuration of these residues. Epimerisation of the Val(4) and Tyr(5) caused a reduction in the compactness of bursopentin. In all cases, the occurrence of a turn structure was relatively high, especially when Arg(3) was in the d configuration. Thermodynamic analysis of the epimerisation process showed that the formation of d amino acid residues is favourable.

The epimerization of amino acid residues increases with age in living organisms. In the present study, the structural consequences and thermodynamic functions of the epimerization of thymopentin (TP-5), the active site of the thymic hormone thymopoietin, were studied using molecular dynamics and density functional theory methods. The results show that free radical-initiated D-amino acid formation is energetically favoured (-130 kJmol(-1)) for each residue and induces significant changes to the peptide structure. In comparison to the wild-type (each residue in the L-configuration), the radius of gyration of the D-Asp(3) epimer of the peptide decreased by 0.5 Å, and disrupted the intramolecular hydrogen bonding of the native peptide. Beyond establishing important structural, energetic and thermodynamic benchmarks and reference data for the structure of TP-5, these results disseminate the understanding of molecular ageing, the epimerization of amino acid residues.

In spite of recent investigations into the potential pharmaceutical importance of bile acids as drug carriers, the structure of bile acid aggregates is largely unknown. Here, we used global optimization techniques to find the lowest energy configurations for clusters composed between 2 and 10 cholate molecules, and evaluated the relative stabilities of the global minima. We found that the energetically most preferred geometries for small aggregates are in fact reverse micellar arrangements, and the classical micellar behaviour (efficient burial of hydrophobic parts) is achieved only in systems containing more than five cholate units. Hydrogen bonding plays a very important part in keeping together the monomers, and among the size range considered, the most stable structure was found to be the decamer, having 17 hydrogen bonds. Molecular dynamics simulations showed that the decamer has the lowest dissociation propensity among the studied aggregation numbers.

We report the identification of novel histamine H4 receptor ligands by ensemble docking on homology model conformers derived from molecular dynamics simulations. Selected receptor models from the trajectories demonstrated superior virtual screening performance compared to the initial models. The ensemble of the best models was able to retrieve a diverse set of known H4 ligands. Prospective virtual screening against these models and subsequent in vitro experimental validation identified novel H4 ligands. Compound 3 showing highest affinity and ligand efficiency represents an interesting scaffold for further medicinal chemistry exploration., (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Penicillin, travels through bodily fluids, targeting and acylatively inactivating enzymes responsible for cell-wall synthesis in gram-positive bacteria. Somehow, it avoids metabolic degradation remaining inactive en route. To resolve this ability to switch from a non-active, to a highly reactive form, we investigated the dynamic structure-activity relationship of penicillin by inelastic neutron spectroscopy, reaction kinetics, NMR and multi-scale theoretical modelling (QM/MM and post-HF ab initio). Results show that by a self-activating physiological pH-dependent two-step proton-mediated process, penicillin changes geometry to activate its irreversibly reactive acylation, facilitated by systemic intramolecular energy management and cooperative vibrations. This dynamic mechanism is confirmed by the first ever reported characterisation of an antibiotic by neutrons, achieved on the TOSCA instrument (ISIS facility, RAL, UK).

The formation of ligand-protein complexes requires simultaneous adaptation of the binding partners. In structure based virtual screening, high throughput docking approaches typically consider the ligand flexibility, but the conformational freedom of the protein is usually taken into account in a limited way. The goal of this study is to elaborate a methodology for incorporating protein flexibility to improve the virtual screening enrichments on GPCRs. Explicit-solvated molecular dynamics simulations (MD) were carried out in lipid bilayers to generate an ensemble of protein conformations for the X-ray structures and homology models of both aminergic and peptidergic GPCRs including the chemokine CXCR4, dopamine D3, histamine H4, and serotonin 5HT6 holo receptor complexes. The quality of the receptor models was assessed by enrichment studies to compare X-ray structures, homology models, and snapshots from the MD trajectory. According to our results, selected frames from the MD trajectory can outperform X-ray structures and homology models in terms of enrichment factor and AUC values. Significant changes were observed considering EF1% values: comparing the original CXCR4, D3, and H4 targets and the additional 5HT6 initial models to that of the best MD frame resulted in 0 to 6.7, 0.32 to 3.5 (10×), 13.3 to 26.7 (2×), and 0 to 14.1 improvements, respectively. It is worth noting that rank-average based ensemble evaluation calculated for different ensemble sizes could not improve the results further. We propose here that MD simulation can capture protein conformations representing the key interacting points of the receptor but less biased toward one specific chemotype. These conformations are useful for the identification of a "consensus" binding site with improved performance in virtual screening.

Non-reactive, comparative (2 × 1.2 μs) molecular dynamics simulations were carried out to characterize the interactions between glutathione (GSH, host molecule) and hydroxyl radical (OH(•), guest molecule). From this analysis, two distinct steps were identified in the recognition process of hydroxyl radical by glutathione: catching and steering, based on the interactions between the host-guest molecules. Over 78% of all interactions are related to the catching mechanism via complex formation between anionic carboxyl groups and the OH radical, hence both terminal residues of GSH serve as recognition sites. The glycine residue has an additional role in the recognition of OH radical, namely the steering. The flexibility of the Gly residue enables the formation of further interactions of other parts of glutathione (e.g. thiol, α- and β-carbons) with the lone electron pair of the hydroxyl radical. Moreover, quantum chemical calculations were carried out on selected GSH/OH(•) complexes and on appropriate GSH conformers to describe the energy profile of the recognition process. The relative enthalpy and the free energy changes of the radical recognition of the strongest complexes varied from -42.4 to -27.8 kJ/mol and from -21.3 to 9.8 kJ/mol, respectively. These complexes, containing two or more intermolecular interactions, would be the starting configurations for the hydrogen atom migration to quench the hydroxyl radical via different reaction channels.

Galactokinase is responsible for the phosphorylation of α-D: -galactose, which is an important step in the metabolism of the latter. Malfunctioning of galactokinase due to a single point mutation causes cataracts and, in serious cases, blindness. This paper reports a study of the Pro(28)Thr point mutation using a variety of theories including molecular dynamics (MD), MM-PBSA/GBSA calculations and AIM analysis. Altered H-bonding networks were detected based on geometric and electron density criteria that resulted in local unfolding of the β-sheet secondary structure. Another consequence was the decrease in stability (5-7 kcal mol(-1)) around this region, as confirmed by ΔG(bind) calculations for the extracted part of the whole system. Local unfolding was verified by several other MD simulations performed with different duration, initial velocities and force field. Based on the results, we propose a possible mechanism for the unfolding caused by the Pro(28)Thr point mutation.

6-Methoxykaempferol-3-O-glucoside (6-MKG) was isolated from a Sudanese herb (El-hazha). The pharmacological effects of 6-MKG were tested on isolated non-pregnant or late-pregnant rat uteri in vitro, whilst docking studies were carried out modelling of the binding of 6-MKG to the rat β(2)-adrenoceptor in silico. In vitro studies revealed that 6-MKG was able to relax both the non-pregnant and the late-pregnant uterine contractility with 50% of the E(max) of terbutaline, whilst the EC(50) for 6-MKG was at least half than that of terbutaline. The β(2)-adrenoceptors antagonist 3-(isopropylamino)-1-[(7-methyl-4-indanyl)oxy]butan-2-ol(ICI118,551) competitively antagonised the relaxing effect of 6-MKG. Radioligand binding and cAMP studies confirmed the β(2)-adrenoceptors agonistic property of the compound. In in silico docking studies, 6-MKG bound to rat β(2)-adrenoceptors with low ∆G(bind) value (-11.53±0.06 kcal/mol) and it interacted with four residues of the active site (Asp(113), Asn(312), Cys(191)and Tyr(316)). It is concluded that 6-MKG exerts weak β(2)-adrenoceptor agonistic activity and can be considered a natural compound with potential therapeutic significance in the field of premature pregnant uterine contractions and asthmatic problems., (Copyright © 2011 Elsevier B.V. All rights reserved.)

All possible X-H (where X can be C, N, O or S) bond dissociation energies (BDEs) of glutathione (γ-L-glutamyl-L-cysteinyl-glycine, GSH) and its fragments have been calculated by first principle methods, and the antioxidant potential of GSH was revealed to be higher than expected in earlier studies. Electron delocalization was found to have an important influence on the antioxidant potential. All structures were optimized and their harmonic vibrational frequencies were calculated in the gas phase at the B3LYP/6-31G(d) level of theory. Solvent effects were taken into account for optimizations at the same level of theory by applying the conductor-like polarizable continuum model (CPCM). Hydrogen cleavage from glutathione proved that the G3MP2B3 composite method provides results consistent with the experimental values for bond dissociation enthalpies (DH(298)) of S-H, O-H, C-H, and N-H bonds. In order to replace the G3MP2B3 energies with accurate single point calculations, six density functionals, namely, MPWKCIS, MPWKCIS1K, M06, TPSS1KCIS, TPSSh, and B3LYP, were tested against G3MP2B3 for obtaining accurate bond dissociation energies. The MPWKCIS1K/6-311++G(3df,2p)//B3LYP/6-31G(d) level of theory provides the best correlation with the G3MP2B3 method for BDEs in both phases, and therefore, it is recommended for similar calculations. Gas phase results showed that the O-H bond was the weakest, while in aqueous phase the N-H bond in the ammonium group proved to have the smallest BDE value in the studied system. In both cases, the cleavage of the X-H bond was followed by decarboxylation which was responsible for the energetic preference of these processes over the S-H dissociation, which was regarded as the most favorable one until now. The calculated BDE values showed that in aqueous phase the most preferred H-abstraction site is at the weakest N-H bond (BDE(aq) = 349.3 kJ mol(-1)) in the glutamine fragment near the α-carbon. In water, the formation of N-centered radicals compared to S-centered ones (BDE(aq) = 351.7 kJ mol(-1)) is more endothermic by 2.54 kJ mol(-1), due to decarboxylation. Hydrogen dissociation energies from the α-carbons are also comparable in energy with those of the thiol hydrogen, within computational error. The higher stability of the radicals--except the S-centered ones--is due to various degrees of electron delocalization. In aqueous phase, four quasi-equivalent stable radical centers (the α-carbons, the N-centered radical of the NH(2) group, and the S-centered radical) were found which provide the antioxidant behavior of glutathione., (© 2011 American Chemical Society)

It has always been a question of considerable scientific interest why amino acids (and other biomolecule building blocks) formed and accumulated in the prebiotic ocean. In this study, we suggest an answer to this question for the simplest amino acid, glycine. We have shown for the first time that classical equilibrium thermodynamics can explain the most likely selection of glycine (and the derivative of its dipeptide) in aqueous media, although glycine is not the lowest free energy structure among all (404) possible constitutional isomers. Species preceding glycine in the free energy order are either supramolecular complexes of small molecules or such molecules likely to dissociate and thus get back to the gas phase. Then, 2-hydroxyacetamide condensates yielding a thermodynamically favored derivative of glycine dipeptide providing an alternative way for peptide formation. It is remarkable that a simple equilibrium thermodynamic model can explain the accumulation of glycine and provide a reason for the importance of water in the formation process.

Molecular dynamics (MD) simulations in a membrane-embedded environment were carried out on the homology model of the human histamine H4 receptor (hH4R) alone and in complex with its endogenous activator histamine and with the first reported selective hH4R antagonist JNJ7777120. During the simulation of the histamine-hH4R complex, considerable changes occurred in the hH4R structure as well as in the interaction pattern of histamine at the binding site. These changes are in agreement with experimental data published on GPCR activation. In particular, the intracellular side of TM helix VI moved significantly away from TM helices III and VII. Moreover, histamine formed an interaction with Asn147 (4.57) that was previously proved to be important in hH4R activation. Results of the MD simulations of the native hH4R and the JNJ7777120-hH4R complex suggest that these models represent an inactive conformation of hH4R. MD simulation in the presence of JNJ7777120 resulted in the movement of the intracellular side of TM helix VI in the direction of TM helix III. Snapshots of the simulations may serve as functionally relevant models in the development of novel hH4R ligands in the future.

The oxytocin receptor belongs to class A receptors within the great family of G protein-coupled receptors. The endogenous ligand oxytocin is a nonapeptide hormone that induces the uterine contractions at parturition and is used to induce the labor. The peptide oxytocin and, even more, its non-peptide antagonist, could be valuable tools in tocolysis. The knowledge of the three-dimensional structure of the oxytocin receptor and the determination of the main interaction points between the receptor and the ligands may help to develop selective oxytocin agonists and antagonist. This review summarizes the knowledge about the mapping of the binding domain of the oxytocin receptor and the efforts in the field of molecular modeling studies related to oxytocin receptor-ligand interactions.

Three-dimensional structure-activity relationship studies of alpha2a-adrenoceptor agonists were carried out by Distance Comparison (DISCOthech) and Comparative Molecular Field Analysis (CoMFA) methods to define the pharmacophore and a quantitative model, respectively, of this class of compounds. The statistical validation of the CoMFA model indicates its high predictive performance for binding affinities of new alpha2a-adrenoceptor agonists.

The aim of this work was to answer the question of whether the general amber force field (GAFF) is good enough to simulate fully hydrated POPC membrane bilayers. The test system contained 128 POPC and 2985 TIP3P water molecules. The equilibration was carried out in a nonarbitrary manner to reach the stable liquid-crystalline phase. The simulations were performed by using particle mesh Ewald electrostatics implemented in molecular dynamics packages Amber8 (NPT ensembles) and NAMD2 (NPgammaT ensembles). The computational results were assessed against the following experimental membrane properties: (i) area per lipid, (ii) area compressibility modulus, (iii) order parameter, (iv) gauche conformations per acyl chain, (v) lateral diffusion coefficients, (vi) electron density profile, and (vii) bound water at the lipid/water interface. The analyses revealed that the tested force field combination approximates the experimental values at an unexpectedly good level when the NPgammaT ensemble is applied with a surface tension of 60 mN m(-1) per bilayer. It is concluded that the GAFF/TIP3P combination can be utilized for aqueous membrane bilayer simulations, as it provides acceptable accuracy for biomolecular modeling., (Copyright 2007 Wiley Periodicals, Inc.)

In the present study, QSAR calculations were performed on the receptor-based alignment of 58 non-peptide human oxytocin receptor antagonists. With the aid of different scoring functions (AutoDock 3.05 built-in and X-Score 1.2) the evolved receptor-ligand complexes were characterized. By means of various datasets it was confirmed that the scoring functions were not capable to predict the biological activity correctly in compounds containing a rigid derivative in the variable region. To improve the pKi prediction 3D-QSAR calculation was performed. The regions related to the biological activity were determined by using cross-validated r2(q2)-guided region selection (q2-GRS) method. The predictive power of the CoMFA model [r(pred)2=0.89, q2(LMO, five groups)=0.695+/-0.034] allowed prediction of the biological activities of newly synthesized compounds and confirmed the receptor-based alignment.

In this study, we performed a molecular docking and dynamics simulation for a benzoxazinone-human oxytocin receptor system to determine the possible hydrophobic and electrostatic interaction points in the dynamic complex. After the homology modeling, the ligand was docked into the putative active using AutoDock 3.05. After the application of energetic and structural filters, the complexes obtained were further refined with a simulated annealing protocol (AMBER8) to remove steric clashes. Three complexes were selected for subjection to the molecular dynamics simulation (5 ns), and the results on the occurrence of average anchor points showed a stable complex between the benzoxazinone derivative and the receptor. The complex could be used as a good starting point for further analysis with site-directed mutagenesis, or further computational research.

In this study, eight oxytocin (OT) antagonists were subjected to conformational analysis. By means of simulated annealing, 2,000 structures were generated to sample the conformational space. The preferred main-chain and side-chain dihedral angles, intramolecular H-bond interactions (types 1<--4 and 1<--3), and secondary-structure elements (beta-turns) were determined in the series of compounds. The influence of the modified amino acid on these structural elements was established; the values determined will be used as filtration criterion in preliminary research.

Molecular docking and 3D-QSAR studies were performed to determine the binding mode for a series of benzoxazine oxytocin antagonists taken from the literature. Structural hypotheses were generated by docking the most active molecule to the rigid receptor by means of AutoDock 3.05. The cluster analysis yielded seven possible binding conformations. These structures were refined by using constrained simulated annealing, and the further ligands were aligned in the refined receptor by molecular docking. A good correlation was found between the estimated deltaG(bind) and the pKi values for complex F. The Connolly-surface analysis, CoMFA and CoMSIA models q2(CoMFA) = 0.653, q2(CoMSA) = 0.630 and r2(pred,CoMFA) = 0.852 , r2(pred,CoMSIA) = 0.815) confirmed the scoring function results. The structural features of the receptor-ligand complex and the CoMFA and CoMSIA fields are in closely connected. These results suggest that receptor-ligand complex F is the most likely binding hypothesis for the studied benzoxazine analogs.