In silico modelling of the interaction of flavonoids with human P-glycoprotein nucleotide-binding domain

A three-dimensional model of human ABCB1 nucleotide-binding domain (NBD) was developed by homology modelling using the high-resolution human TAP1 transporter structure as template. Interactions between NBD and flavonoids were investigated using in silico docking studies. Ring-A of unmodified flavonoid was located within the NBD P-loop with the 5-hydroxyl group involved in hydrogen bonding with Lys1076. Ring-B was stabilised by hydrophobic stacking interactions with Tyr1044. The 3-hydroxyl group and carbonyl oxygen were extensively involved in hydrogen bonding interactions with amino acids within the NBD. Addition of prenyl, benzyl or geranyl moieties to ring-A (position-6) and hydrocarbon substituents (O-n-butyl to O-n-decyl) to ring-B (position-4) resulted in a size-dependent decrease in predicted docking energy which reflected the increased binding affinities reported in vitro.

Statistical deconvolution of enthalpic energetic contributions to MHC-peptide binding affinity

Background - MHC Class I molecules present antigenic peptides to cytotoxic T cells, which forms an integral part of the adaptive immune response. Peptides are bound within a groove formed by the MHC heavy chain. Previous approaches to MHC Class I-peptide binding prediction have largely concentrated on the peptide anchor residues located at the P2 and C-terminus positions. Results - A large dataset comprising MHC-peptide structural complexes was created by re-modelling pre-determined x-ray crystallographic structures. Static energetic analysis, following energy minimisation, was performed on the dataset in order to characterise interactions between bound peptides and the MHC Class I molecule, partitioning the interactions within the groove into van der Waals, electrostatic and total non-bonded energy contributions. Conclusion - The QSAR techniques of Genetic Function Approximation (GFA) and Genetic Partial Least Squares (G/PLS) algorithms were used to identify key interactions between the two molecules by comparing the calculated energy values with experimentally-determined BL50 data. Although the peptide termini binding interactions help ensure the stability of the MHC Class I-peptide complex, the central region of the peptide is also important in defining the specificity of the interaction. As thermodynamic studies indicate that peptide association and dissociation may be driven entropically, it may be necessary to incorporate entropic contributions into future calculations.

Mercury and selenium binding biomolecules in terrestrial mammals (Cervus elaphus and Sus scrofa) from a mercury exposed area

Acknowledgements The authors are grateful to Junta de Comunidades de Castilla-La Mancha (PCC-05-004-2, PAI06-0094, PCI-08-0096, PEII09-0032-5329) and the Ministerio de Economía y Competitividad (CTQ2013-48411-P) for financial support. M.J. Patiño Ropero acknowledges the Junta de Comunidades de Castilla-La Mancha for her PhD. fellowship.

Path Optimization in Free Energy Calculations

Free energy calculations are a computational method for determining thermodynamic quantities, such as free energies of binding, via simulation.

Currently, due to computational and algorithmic limitations, free energy calculations are limited in scope.

In this work, we propose two methods for improving the efficiency of free energy calculations.

First, we expand the state space of alchemical intermediates, and show that this expansion enables us to calculate free energies along lower variance paths.

We use Q-learning, a reinforcement learning technique, to discover and optimize paths at low computational cost.

Second, we reduce the cost of sampling along a given path by using sequential Monte Carlo samplers.

We develop a new free energy estimator, pCrooks (pairwise Crooks), a variant on the Crooks fluctuation theorem (CFT), which enables decomposition of the variance of the free energy estimate for discrete paths, while retaining beneficial characteristics of CFT.

Combining these two advancements, we show that for some test models, optimal expanded-space paths have a nearly 80% reduction in variance relative to the standard path.

Additionally, our free energy estimator converges at a more consistent rate and on average 1.8 times faster when we enable path searching, even when the cost of path discovery and refinement is considered.

Adsorption Calorimetry for Energy Conversion Technologies: Applications in Organic Photovoltaics, Catalysis, and Atomic Layer Deposition

Thesis (Ph.D.)--University of Washington, 2016-07

Single-molecule fluorescence resonance energy transfer study of SNARE-mediated membrane fusion

This is a comprehensive study of protein-mediated membrane fusion through single-molecule fluorescence resonance energy transfer (smFRET). Membrane fusion is one of the important cellular processes by which two initially distinct lipid bilayers merge their hydrophobic cores, resulting in one interconnected structure. For example, exocytosis, fertilization of an egg by a sperm and communication between neurons are a few among many processes that rely on some form of fusion. Proteins called soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) play a central role in fusion processes which is also regulated by many accessory proteins, such as synaptotagmin, complexin and Munc18. By a new lipid mixing method at the single-vesicle level, we are able to accurately detect different stages of SNARE-mediated membrane fusion including docking, hemi and full fusion via FRET value of single donor/acceptor vesicle pair. Through this single-vesicle lipid mixing assay, we discovered the vesicle aggregation induced by C2AB/Ca2+, the dual function of complexin, and the fusion promotion role of Munc18/SNARE-core binding mode. While this new method provides the information regarding the extent of the ensemble lipid mixing, the fusion pore opening between two vesicular cavities and the interaction between proteins cannot be detected. In order to overcome these limitations, we then developed a single-vesicle content mixing method to reveal the key factor of pore expansion by detecting the FRET change of dual-labeled DNA probes encapsulated in vesicles. Through our single-vesicle content mixing assay, we found the fusion pore expansion role of yeast SNAREs as well as neuronal SNAREs plus synaptotagmin 1.

SOVLENT REACTIVITY AND INTERFACE EVOLUTION AT MODEL ELECTRODES FOR ENERGY APPLICATIONS

The Li-ion rechargeable battery (LIB) is widely used as an energy storage device, but has significant limitations in battery cycle life and safety. During initial charging, decomposition of the ethylene carbonate (EC)-based electrolytes of the LIB leads to the formation of a passivating layer on the anode known as the solid electrolyte interphase (SEI). The formation of an SEI has great impact on the cycle life and safety of LIB, yet mechanistic aspects of SEI formation are not fully understood. In this dissertation, two surface science model systems have been created under ultra-high vacuum (UHV) to probe the very initial stage of SEI formation at the model carbon anode surfaces of LIB. The first model system, Model System I, is an lithium-carbonate electrolyte/graphite C(0001) system. I have developed a temperature programmed desorption/temperature programmed reaction spectroscopy (TPD/TPRS) instrument as part of my dissertation to study Model System I in quantitative detail. The binding strengths and film growth mechanisms of key electrolyte molecules on model carbon anode surfaces with varying extents of lithiation were measured by TPD. TPRS was further used to track the gases evolved from different reduction products in the early-stage SEI formation. The branching ratio of multiple reaction pathways was quantified for the first time and determined to be 70.% organolithium products vs. 30% inorganic lithium product. The obtained branching ratio provides important information on the distribution of lithium salts that form at the very onset of SEI formation. One of the key reduction products formed from EC in early-stage SEI formation is lithium ethylene dicarbonate (LEDC). Despite intensive studies, the LEDC structure in either the bulk or thin-film (SEI) form is unknown. To enable structural study, pure LEDC was synthesized and subject to synchrotron X-ray diffraction measurements (bulk material) and STM measurements (deposited films). To enable studies of LEDC thin films, Model System II, a lithium ethylene dicarbonate (LEDC)-dimethylformamide (DMF)/Ag(111) system was created by a solution microaerosol deposition technique. Produced films were then imaged by ultra-high vacuum scanning tunneling microscopy (UHV-STM). As a control, the dimethylformamide (DMF)-Ag(111) system was first prepared and its complex 2D phase behavior was mapped out as a function of coverage. The evolution of three distinct monolayer phases of DMF was observed with increasing surface pressure — a 2D gas phase, an ordered DMF phase, and an ordered Ag(DMF)2 complex phase. The addition of LEDC to this mixture, seeded the nucleation of the ordered DMF islands at lower surface pressures (DMF coverages), and was interpreted through nucleation theory. A structural model of the nucleation seed was proposed, and the implication of ionic SEI products, such as LEDC, in early-stage SEI formation was discussed.

Novel Mutation in the ATP-Binding Cassette Transporter A3 (ABCA3) Encoding Gene Causes Respiratory Distress Syndrome in A Term Newborn in Southwest Iran

Introduction: ABCA3 glycoprotein belongs to the ATP-binding cassette (ABC) superfamily of transporters, which utilize the energy derived from hydrolysis of ATP for the translocation of a wide variety of substrates across the plasma membrane. Mutations in the ABCA3 gene are knowingly causative for fatal surfactant deficiency, particularly respiratory distress syndrome (RDS) in term babies. Case Presentation: In this study, Sanger sequencing of the whole ABCA3 gene (NCBI NM_001089) was performed in a neonatal boy with severe RDS. A homozygous mutation has been identified in the patient. Parents were heterozygous for the same missense mutation GGA > AGA at position 202 in exon 6 of the ABCA3 gene (c.604G > A; p.G202R). Furthermore, 70 normal individuals have been analyzed for the mentioned change with negative results. Conclusions: Regarding Human Genome Mutation Database (HGMD) and other literature recherche, the detected change is a novel mutation and has not been reported before. Bioinformatics mutation predicting tools prefer it as pathogenic.

Understanding Co-translational Protein Targeting and Lithium Dendrite Formation through Free Energy Simulations and Coarse-Grained Models

We describe the application of alchemical free energy methods and coarse-grained models to study two key problems: (i) co-translational protein targeting and insertion to direct membrane proteins to the endoplasmic reticulum for proper localization and folding, (ii) lithium dendrite formation during recharging of lithium metal batteries. We show that conformational changes in the signal recognition particle, a central component of the protein targeting machinery, confer additional specificity during the the recognition of signal sequences. We then develop a three-dimensional coarse-grained model to study the long-timescale dynamics of membrane protein integration at the translocon and a framework for the calculation of binding free energies between the ribosome and translocon. Finally, we develop a coarse-grained model to capture the dynamics of lithium deposition and dissolution at the electrode interface with time-dependent voltages to show that pulse plating and reverse pulse plating methods can mitigate dendrite growth.

Beam-energy Dependence Of The Directed Flow Of Protons, Antiprotons, And Pions In Au+au Collisions.

Rapidity-odd directed flow (v1) measurements for charged pions, protons, and antiprotons near midrapidity (y=0) are reported in sNN=7.7, 11.5, 19.6, 27, 39, 62.4, and 200 GeV Au+Au collisions as recorded by the STAR detector at the Relativistic Heavy Ion Collider. At intermediate impact parameters, the proton and net-proton slope parameter dv1/dy|y=0 shows a minimum between 11.5 and 19.6 GeV. In addition, the net-proton dv1/dy|y=0 changes sign twice between 7.7 and 39 GeV. The proton and net-proton results qualitatively resemble predictions of a hydrodynamic model with a first-order phase transition from hadronic matter to deconfined matter, and differ from hadronic transport calculations.

Direct Visualization Of The Action Of Triton X-100 On Giant Vesicles Of Erythrocyte Membrane Lipids.

The raft hypothesis proposes that microdomains enriched in sphingolipids, cholesterol, and specific proteins are transiently formed to accomplish important cellular tasks. Equivocally, detergent-resistant membranes were initially assumed to be identical to membrane rafts, because of similarities between their compositions. In fact, the impact of detergents in membrane organization is still controversial. Here, we use phase contrast and fluorescence microscopy to observe giant unilamellar vesicles (GUVs) made of erythrocyte membrane lipids (erythro-GUVs) when exposed to the detergent Triton X-100 (TX-100). We clearly show that TX-100 has a restructuring action on biomembranes. Contact with TX-100 readily induces domain formation on the previously homogeneous membrane of erythro-GUVs at physiological and room temperatures. The shape and dynamics of the formed domains point to liquid-ordered/liquid-disordered (Lo/Ld) phase separation, typically found in raft-like ternary lipid mixtures. The Ld domains are then separated from the original vesicle and completely solubilized by TX-100. The insoluble vesicle left, in the Lo phase, represents around 2/3 of the original vesicle surface at room temperature and decreases to almost 1/2 at physiological temperature. This chain of events could be entirely reproduced with biomimetic GUVs of a simple ternary lipid mixture, 2:1:2 POPC/SM/chol (phosphatidylcholine/sphyngomyelin/cholesterol), showing that this behavior will arise because of fundamental physicochemical properties of simple lipid mixtures. This work provides direct visualization of TX-100-induced domain formation followed by selective (Ld phase) solubilization in a model system with a complex biological lipid composition.

Hypothalamic Inflammation And The Central Nervous System Control Of Energy Homeostasis.

The control of energy homeostasis relies on robust neuronal circuits that regulate food intake and energy expenditure. Although the physiology of these circuits is well understood, the molecular and cellular response of this program to chronic diseases is still largely unclear. Hypothalamic inflammation has emerged as a major driver of energy homeostasis dysfunction in both obesity and anorexia. Importantly, this inflammation disrupts the action of metabolic signals promoting anabolism or supporting catabolism. In this review, we address the evidence that favors hypothalamic inflammation as a factor that resets energy homeostasis in pathological states.

Changes In Liver Cell Dna Methylation Status In Diabetic Mice Affect Its Ft-ir Characteristics.

Lower levels of cytosine methylation have been found in the liver cell DNA from non-obese diabetic (NOD) mice under hyperglycemic conditions. Because the Fourier transform-infrared (FT-IR) profiles of dry DNA samples are differently affected by DNA base composition, single-stranded form and histone binding, it is expected that the methylation status in the DNA could also affect its FT-IR profile. The DNA FT-IR signatures obtained from the liver cell nuclei of hyperglycemic and normoglycemic NOD mice of the same age were compared. Dried DNA samples were examined in an IR microspectroscope equipped with an all-reflecting objective (ARO) and adequate software. Changes in DNA cytosine methylation levels induced by hyperglycemia in mouse liver cells produced changes in the respective DNA FT-IR profiles, revealing modifications to the vibrational intensities and frequencies of several chemical markers, including νas -CH3 stretching vibrations in the 5-methylcytosine methyl group. A smaller band area reflecting lower energy absorbed in the DNA was found in the hyperglycemic mice and assumed to be related to the lower levels of -CH3 groups. Other spectral differences were found at 1700-1500 cm(-1) and in the fingerprint region, and a slight change in the DNA conformation at the lower DNA methylation levels was suggested for the hyperglycemic mice. The changes that affect cytosine methylation levels certainly affect the DNA-protein interactions and, consequently, gene expression in liver cells from the hyperglycemic NOD mice.

Beam-energy Dependence Of Charge Separation Along The Magnetic Field In Au+au Collisions At Rhic.

Local parity-odd domains are theorized to form inside a quark-gluon plasma which has been produced in high-energy heavy-ion collisions. The local parity-odd domains manifest themselves as charge separation along the magnetic field axis via the chiral magnetic effect. The experimental observation of charge separation has previously been reported for heavy-ion collisions at the top RHIC energies. In this Letter, we present the results of the beam-energy dependence of the charge correlations in Au+Au collisions at midrapidity for center-of-mass energies of 7.7, 11.5, 19.6, 27, 39, and 62.4 GeV from the STAR experiment. After background subtraction, the signal gradually reduces with decreased beam energy and tends to vanish by 7.7 GeV. This implies the dominance of hadronic interactions over partonic ones at lower collision energies.

Rapidly Switching Multidirectional Defibrillation: Reversal Of Ventricular Fibrillation With Lower Energy Shocks.

Cardiac arrest after open surgery has an incidence of approximately 3%, of which more than 50% of the cases are due to ventricular fibrillation. Electrical defibrillation is the most effective therapy for terminating cardiac arrhythmias associated with unstable hemodynamics. The excitation threshold of myocardial microstructures is lower when external electrical fields are applied in the longitudinal direction with respect to the major axis of cells. However, in the heart, cell bundles are disposed in several directions. Improved myocardial excitation and defibrillation have been achieved by applying shocks in multiple directions via intracardiac leads, but the results are controversial when the electrodes are not located within the cardiac chambers. This study was designed to test whether rapidly switching shock delivery in 3 directions could increase the efficiency of direct defibrillation. A multidirectional defibrillator and paddles bearing 3 electrodes each were developed and used in vivo for the reversal of electrically induced ventricular fibrillation in an anesthetized open-chest swine model. Direct defibrillation was performed by unidirectional and multidirectional shocks applied in an alternating fashion. Survival analysis was used to estimate the relationship between the probability of defibrillation and the shock energy. Compared with shock delivery in a single direction in the same animal population, the shock energy required for multidirectional defibrillation was 20% to 30% lower (P < .05) within a wide range of success probabilities. Rapidly switching multidirectional shock delivery required lower shock energy for ventricular fibrillation termination and may be a safer alternative for restoring cardiac sinus rhythm.