scholarly journals Electron transfer parameters for Methemoglobin formation in mutant Hemoglobin α-chains

2021 ◽  
Author(s):  
Vaibhav A. Dixit ◽  
Jochen Blumberger ◽  
Shivam Kumar Vyas
Keyword(s):  
Average Energy ◽  
Md Simulations ◽  
Structure Stability ◽  
Redox Potentials ◽  
Free Energies ◽  
Redox States ◽  
Methemoglobin Formation ◽  
Α Chain ◽  
Reorganization Energies ◽  
The Stability

AbstractHemoglobin mediated transport of dioxygen (O2) critically depends on the stability of the reduced (Fe2+) form of the Heme cofactors. Some protein mutations stabilize oxidized (Fe3+) state (Methemoglobin, Hb M) causing methemoglobinemia and can be lethal above 30 %. Majority of the analyses of factors influencing Hb oxidation are retrospective and give insights only for inner sphere mutations of Heme (His58, His87). Herein, we report the first all atom MD simulations on redox states and calculations of the Marcus ET parameters for the α-chain Hb oxidation and reduction rates for Hb M. The Hb (wild type), and most of the studied α-chain variants maintain globin structure except the Hb M Iwate (H87Y). Using linear response approximation we calculated average energy gaps (<ΔE>), total (λ), protein (λprot), solvent (λsolv) reorganization energies, and redox potentials (E°), and oxidation free energies (ΔG°). The total λ ranges from 0.685 – 0.730 eV in agreement with literature on Hb and similar Heme proteins. The mutants forming Hb M tend to lower the E° and thus stabilize the oxidized (Fe3+) state (e.g. the Hb Miyagi variant with K61E mutation). Solvent reorganization (λsolv 73 – 96 %) makes major contributions to λ, while protein reorganization (λprot) accounts for 27 – 30 % except for the Miyagi and J-Buda variants (λprot ∼ 4 %). Analysis of Heme-solvent H-bonding interactions among variants provide insights into the role of Lys61 residue in stabilizing Fe2+ state and ET parameters. The ET parameters provide valuable insights into the Hb oxidation to Hb M in agreement with the experimental data. Thus our methodology explains the effect of mutations on the structure, stability and Hb oxidation, and has potential for the prediction of methemoglobinemia.

Structure, Stability and Water Adsorption on Ultra-Thin TiO2 Supported on TiN

2019 ◽  
Author(s):  
Jose Julio Gutierrez Moreno ◽  
Marco Fronzi ◽  
Pierre Lovera ◽  
alan O'Riordan ◽  
Mike J Ford ◽  
...  
Keyword(s):  
Density Functional ◽  
Water Adsorption ◽  
Chemical Potential ◽  
Energy Gap ◽  
Molecular Water ◽  
Structure Stability ◽  
Ambient Conditions ◽  
Rutile Structure ◽  
The Stability ◽  
The One

<p></p><p>Interfacial metal-oxide systems with ultrathin oxide layers are of high interest for their use in catalysis. In this study, we present a density functional theory (DFT) investigation of the structure of ultrathin rutile layers (one and two TiO<sub>2</sub> layers) supported on TiN and the stability of water on these interfacial structures. The rutile layers are stabilized on the TiN surface through the formation of interfacial Ti–O bonds. Charge transfer from the TiN substrate leads to the formation of reduced Ti<sup>3+</sup> cations in TiO<sub>2.</sub> The structure of the one-layer oxide slab is strongly distorted at the interface, while the thicker TiO<sub>2</sub> layer preserves the rutile structure. The energy cost for the formation of a single O vacancy in the one-layer oxide slab is only 0.5 eV with respect to the ideal interface. For the two-layer oxide slab, the introduction of several vacancies in an already non-stoichiometric system becomes progressively more favourable, which indicates the stability of the highly non-stoichiometric interfaces. Isolated water molecules dissociate when adsorbed at the TiO<sub>2</sub> layers. At higher coverages the preference is for molecular water adsorption. Our ab initio thermodynamics calculations show the fully water covered stoichiometric models as the most stable structure at typical ambient conditions. Interfacial models with multiple vacancies are most stable at low (reducing) oxygen chemical potential values. A water monolayer adsorbs dissociatively on the highly distorted 2-layer TiO<sub>1.75</sub>-TiN interface, where the Ti<sup>3+</sup> states lying above the top of the valence band contribute to a significant reduction of the energy gap compared to the stoichiometric TiO<sub>2</sub>-TiN model. Our results provide a guide for the design of novel interfacial systems containing ultrathin TiO<sub>2</sub> with potential application as photocatalytic water splitting devices.</p><p></p>

Polysaccharide Structures in the Outer Mucilage of Arabidopsis Seeds Visualized by AFM

2020 ◽  
Author(s):  
MAK Williams ◽  
V Cornuault ◽  
AH Irani ◽  
VV Symonds ◽  
J Malmström ◽  
...  
Keyword(s):  
Average Length ◽  
American Chemical Society ◽  
Md Simulations ◽  
Chemical Society ◽  
Side Chains ◽  
Rhamnogalacturonan I ◽  
Afm Imaging ◽  
Minor Population ◽  
The Stability ◽  
A Minor

© 2020 American Chemical Society. Evidence is presented that the polysaccharide rhamnogalacturonan I (RGI) can be biosynthesized in remarkably organized branched configurations and surprisingly long versions and can self-assemble into a plethora of structures. AFM imaging has been applied to study the outer mucilage obtained from wild-type (WT) and mutant (bxl1-3 and cesa5-1) Arabidopsis thaliana seeds. For WT mucilage, ordered, multichain structures of the polysaccharide RGI were observed, with a helical twist visible in favorable circumstances. Molecular dynamics (MD) simulations demonstrated the stability of several possible multichain complexes and the possibility of twisted fibril formation. For bxl1-3 seeds, the imaged polymers clearly showed the presence of side chains. These were surprisingly regular and well organized with an average length of ∼100 nm and a spacing of ∼50 nm. The heights of the side chains imaged were suggestive of single polysaccharide chains, while the backbone was on average 4 times this height and showed regular height variations along its length consistent with models of multichain fibrils examined in MD. Finally, in mucilage extracts from cesa5-1 seeds, a minor population of chains in excess of 30 μm long was observed.

In Silico and in Vitro Evaluation of Deamidation Effects on the Stability of the Fusion Toxin DAB389IL-2

2019 ◽  
Vol 16 (4) ◽  
pp. 307-313 ◽  
Author(s):  
Nasrin Zarkar ◽  
Mohammad Ali Nasiri Khalili ◽  
Fathollah Ahmadpour ◽  
Sirus Khodadadi ◽  
Mehdi Zeinoddini
Keyword(s):  
In Silico ◽  
Catalytic Domain ◽  
Md Simulations ◽  
Special Importance ◽  
Native Form ◽  
Vitro Experiment ◽  
In Vitro Experiment ◽  
Pharmaceutical Protein ◽  
The Stability

Background: DAB389IL-2 (Denileukin diftitox) as an immunotoxin is a targeted pharmaceutical protein and is the first immunotoxin approved by FDA. It is used for the treatment of various kinds of cancer such as CTCL lymphoma, melanoma, and Leukemia but among all of these, treatment of CTCL has special importance. DAB389IL-2 consists of two distinct parts; the catalytic domain of Diphtheria Toxin (DT) that genetically fused to the whole IL-2. Deamidation is the most important reaction for chemical instability of proteins occurs during manufacture and storage. Deamidation of asparagine residues occurs at a higher rate than glutamine residues. The structure of proteins, temperature and pH are the most important factors that influence the rate of deamidation. Methods: Since there is not any information about deamidation of DAB389IL-2, we studied in silico deamidation by Molecular Dynamic (MD) simulations using GROMACS software. The 3D model of fusion protein DAB389IL-2 was used as a template for deamidation. Then, the stability of deamidated and native form of the drug was calculated. Results: The results of MD simulations were showed that the deamidated form of DAB389IL-2 is more unstable than the normal form. Also, deamidation was carried by incubating DAB389IL-2, 0.3 mg/ml in ammonium hydrogen carbonate for 24 h at 37o C in order to in vitro experiment. Conclusion: The results of in vitro experiment were confirmed outcomes of in silico study. In silico and in vitro experiments were demonstrated that DAB389IL-2 is unstable in deamidated form.

scholarly journals A Computational Protocol Combining DFT and Cheminformatics for Prediction of pH-Dependent Redox Potentials

Molecules ◽  
2021 ◽  
Vol 26 (13) ◽  
pp. 3978
Author(s):  
Rocco Peter Fornari ◽  
Piotr de Silva
Keyword(s):  
Organic Molecules ◽  
Ph Value ◽  
Computational Cost ◽  
Legendre Transform ◽  
Good Prediction ◽  
Redox Potentials ◽  
Free Energies ◽  
Hybrid Functional ◽  
Data Set ◽  
State Changes

Discovering new materials for energy storage requires reliable and efficient protocols for predicting key properties of unknown compounds. In the context of the search for new organic electrolytes for redox flow batteries, we present and validate a robust procedure to calculate the redox potentials of organic molecules at any pH value, using widely available quantum chemistry and cheminformatics methods. Using a consistent experimental data set for validation, we explore and compare a few different methods for calculating reaction free energies, the treatment of solvation, and the effect of pH on redox potentials. We find that the B3LYP hybrid functional with the COSMO solvation method, in conjunction with thermal contributions evaluated from BLYP gas-phase harmonic frequencies, yields a good prediction of pH = 0 redox potentials at a moderate computational cost. To predict how the potentials are affected by pH, we propose an improved version of the Alberty-Legendre transform that allows the construction of a more realistic Pourbaix diagram by taking into account how the protonation state changes with pH.

scholarly journals Structural genomics approach to investigate deleterious impact of nsSNPs in conserved telomere maintenance component 1

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Arunabh Choudhury ◽  
Taj Mohammad ◽  
Nikhil Samarth ◽  
Afzal Hussain ◽  
Md. Tabish Rehman ◽  
...  
Keyword(s):  
Noncovalent Interactions ◽  
Md Simulations ◽  
Telomere Maintenance ◽  
Nucleotide Polymorphisms ◽  
Telomeric Dna ◽  
Congenital Diseases ◽  
Ob Fold ◽  
Significant Conformational Change ◽  
The Stability ◽  
The Impact

AbstractConserved telomere maintenance component 1 (CTC1) is an important component of the CST (CTC1-STN1-TEN1) complex, involved in maintaining the stability of telomeric DNA. Several non-synonymous single-nucleotide polymorphisms (nsSNPs) in CTC1 have been reported to cause Coats plus syndrome and Dyskeratosis congenital diseases. Here, we have performed sequence and structure analyses of nsSNPs of CTC1 using state-of-the-art computational methods. The structure-based study focuses on the C-terminal OB-fold region of CTC1. There are 11 pathogenic mutations identified, and detailed structural analyses were performed. These mutations cause a significant disruption of noncovalent interactions, which may be a possible reason for CTC1 instability and consequent diseases. To see the impact of such mutations on the protein conformation, all-atom molecular dynamics (MD) simulations of CTC1-wild-type (WT) and two of the selected mutations, R806C and R806L for 200 ns, were carried out. A significant conformational change in the structure of the R806C mutant was observed. This study provides a valuable direction to understand the molecular basis of CTC1 dysfunction in disease progression, including Coats plus syndrome.

scholarly journals Structural Properties of G,T-Parallel Duplexes

2010 ◽  
Vol 2010 ◽  
pp. 1-11 ◽  
Author(s):  
Anna Aviñó ◽  
Elena Cubero ◽  
Raimundo Gargallo ◽  
Carlos González ◽  
Modesto Orozco ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Structural Properties ◽  
Theoretical Calculations ◽  
Md Simulations ◽  
Theoretical Studies ◽  
Duplex Structure ◽  
Strong Stabilization ◽  
The Stability ◽  
The Impact

The structure of G,T-parallel-stranded duplexes of DNA carrying similar amounts of adenine and guanine residues is studied by means of molecular dynamics (MD) simulations and UV- and CD spectroscopies. In addition the impact of the substitution of adenine by 8-aminoadenine and guanine by 8-aminoguanine is analyzed. The presence of 8-aminoadenine and 8-aminoguanine stabilizes the parallel duplex structure. Binding of these oligonucleotides to their target polypyrimidine sequences to form the corresponding G,T-parallel triplex was not observed. Instead, when unmodified parallel-stranded duplexes were mixed with their polypyrimidine target, an interstrand Watson-Crick duplex was formed. As predicted by theoretical calculations parallel-stranded duplexes carrying 8-aminopurines did not bind to their target. The preference for the parallel-duplex over the Watson-Crick antiparallel duplex is attributed to the strong stabilization of the parallel duplex produced by the 8-aminopurines. Theoretical studies show that the isomorphism of the triads is crucial for the stability of the parallel triplex.

scholarly journals Effects of electron transfer on the stability of hydrogen bonds

2017 ◽  
Vol 8 (11) ◽  
pp. 7324-7329 ◽  
Author(s):  
Tyler M. Porter ◽  
Gavin P. Heim ◽  
Clifford P. Kubiak
Keyword(s):  
Electron Transfer ◽  
Hydrogen Bonds ◽  
Ruthenium Complexes ◽  
Complementary Pair ◽  
Redox States ◽  
The Stability ◽  
Hydrogen Bonded

The measurement of the dimerization constants of hydrogen-bonded ruthenium complexes (12, 22, 32) linked by a self-complementary pair of 4-pyridylcarboxylic acid ligands in different redox states is reported.

scholarly journals Atomistic Analysis of ToxN and ToxI Complex Unbinding Mechanism

2018 ◽  
Vol 19 (11) ◽  
pp. 3524 ◽  
Author(s):  
Guodong Hu ◽  
Xiu Yu ◽  
Yunqiang Bian ◽  
Zanxia Cao ◽  
Shicai Xu ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Noncoding Rna ◽  
Md Simulations ◽  
Normal Function ◽  
Toxin Complex ◽  
Potentials Of Mean Force ◽  
Protein Toxins ◽  
Function Relationship ◽  
The Stability ◽  
Smd Simulations

ToxIN is a triangular structure formed by three protein toxins (ToxNs) and three specific noncoding RNA antitoxins (ToxIs). To respond to stimuli, ToxI is preferentially degraded, releasing the ToxN. Thus, the dynamic character is essential in the normal function interactions between ToxN and ToxI. Here, equilibrated molecular dynamics (MD) simulations were performed to study the stability of ToxN and ToxI. The results indicate that ToxI adjusts the conformation of 3′ and 5′ termini to bind to ToxN. Steered molecular dynamics (SMD) simulations combined with the recently developed thermodynamic integration in 3nD (TI3nD) method were carried out to investigate ToxN unbinding from the ToxIN complex. The potentials of mean force (PMFs) and atomistic pictures suggest the unbinding mechanism as follows: (1) dissociation of the 5′ terminus from ToxN, (2) missing the interactions involved in the 3′ terminus of ToxI without three nucleotides (G31, A32, and A33), (3) starting to unfold for ToxI, (4) leaving the binding package of ToxN for three nucleotides of ToxI, (5) unfolding of ToxI. This work provides information on the structure-function relationship at the atomistic level, which is helpful for designing new potent antibacterial drugs in the future.

scholarly journals Molecular simulations of lipid membrane partitioning and translocation by bacterial quorum sensing modulators

PLoS ONE ◽  
2021 ◽  
Vol 16 (2) ◽  
pp. e0246187
Author(s):  
Tianyi Jin ◽  
Samarthaben J. Patel ◽  
Reid C. Van Lehn
Keyword(s):  
Quorum Sensing ◽  
Lipid Bilayers ◽  
Lipid Membrane ◽  
Computational Cost ◽  
Md Simulations ◽  
Umbrella Sampling ◽  
Communication Process ◽  
Free Energies ◽  
Membrane Partitioning ◽  
Solvent Model

Quorum sensing (QS) is a bacterial communication process mediated by both native and non-native small-molecule quorum sensing modulators (QSMs), many of which have been synthesized to disrupt QS pathways. While structure-activity relationships have been developed to relate QSM structure to the activation or inhibition of QS receptors, less is known about the transport mechanisms that enable QSMs to cross the lipid membrane and access intracellular receptors. In this study, we used atomistic MD simulations and an implicit solvent model, called COSMOmic, to analyze the partitioning and translocation of QSMs across lipid bilayers. We performed umbrella sampling at atomistic resolution to calculate partitioning and translocation free energies for a set of naturally occurring QSMs, then used COSMOmic to screen the water-membrane partition and translocation free energies for 50 native and non-native QSMs that target LasR, one of the LuxR family of quorum-sensing receptors. This screening procedure revealed the influence of systematic changes to head and tail group structures on membrane partitioning and translocation free energies at a significantly reduced computational cost compared to atomistic MD simulations. Comparisons with previously determined QSM activities suggest that QSMs that are least likely to partition into the bilayer are also less active. This work thus demonstrates the ability of the computational protocol to interrogate QSM-bilayer interactions which may help guide the design of new QSMs with engineered membrane interactions.

Sign in / Sign up

Export Citation Format

Share Document