scholarly journals Protonation Equilibrium in the Active Site of the Photoactive Yellow Protein

Molecules ◽  
2021 ◽  
Vol 26 (7) ◽  
pp. 2025
Author(s):  
Pablo Campomanes ◽  
Stefano Vanni

The role and existence of low-barrier hydrogen bonds (LBHBs) in enzymatic and protein activity has been largely debated. An interesting case is that of the photoactive yellow protein (PYP). In this protein, two short HBs adjacent to the chromophore, p-coumaric acid (pCA), have been identified by X-ray and neutron diffraction experiments. However, there is a lack of agreement on the chemical nature of these H-bond interactions. Additionally, no consensus has been reached on the presence of LBHBs in the active site of the protein, despite various experimental and theoretical studies having been carried out to investigate this issue. In this work, we perform a computational study that combines classical and density functional theory (DFT)-based quantum mechanical/molecular mechanical (QM/MM) simulations to shed light onto this controversy. Furthermore, we aim to deepen our understanding of the chemical nature and dynamics of the protons involved in the two short hydrogen bonds that, in the dark state of PYP, connect pCA with the two binding pocket residues (E46 and Y42). Our results support the existence of a strong LBHB between pCA and E46, with the H fully delocalized and shared between both the carboxylic oxygen of E46 and the phenolic oxygen of pCA. Additionally, our findings suggest that the pCA interaction with Y42 can be suitably described as a typical short ionic H-bond of moderate strength that is fully localized on the phenolic oxygen of Y42.

Molecules ◽  
2020 ◽  
Vol 25 (12) ◽  
pp. 2794
Author(s):  
Carly C. Carter ◽  
Thomas R. Cundari

In the present density functional theory (DFT) research, nine different molecules, each with different combinations of A (triel) and E (divalent metal) elements, were reacted to effect methane C–H activation. The compounds modeled herein incorporated the triels A = B, Al, or Ga and the divalent metals E = Be, Mg, or Zn. The results show that changes in the divalent metal have a much bigger impact on the thermodynamics and methane activation barriers than changes in the triels. The activating molecules that contained beryllium were most likely to have the potential for activating methane, as their free energies of reaction and free energy barriers were close to reasonable experimental values (i.e., ΔG close to thermoneutral, ΔG‡ ~30 kcal/mol). In contrast, the molecules that contained larger elements such as Zn and Ga had much higher ΔG‡. The addition of various substituents to the A–E complexes did not seem to affect thermodynamics but had some effect on the kinetics when substituted closer to the active site.


Molecules ◽  
2021 ◽  
Vol 26 (20) ◽  
pp. 6129
Author(s):  
Iori Era ◽  
Yasutaka Kitagawa ◽  
Natsumi Yasuda ◽  
Taigo Kamimura ◽  
Naoka Amamizu ◽  
...  

The effect of hydrogen bonds around the active site of Anabaena [2Fe-2S] ferredoxin (Fd) on a vertical ionization potential of the reduced state (IP(red)) is examined based on the density functional theory (DFT) calculations. The results indicate that a single hydrogen bond increases the relative stability of the reduced state, and shifts IP(red) to a reductive side by 0.31–0.33 eV, regardless of the attached sulfur atoms. In addition, the IP(red) value can be changed by the number of hydrogen bonds around the active site. The results also suggest that the redox potential of [2Fe-2S] Fd is controlled by the number of hydrogen bonds because IP(red) is considered to be a major factor in the redox potential. Furthermore, there is a possibility that the redox potentials of artificial iron-sulfur clusters can be finely controlled by the number of the hydrogen bonds attached to the sulfur atoms of the cluster.


2021 ◽  
Vol 13 (2) ◽  
pp. 679-694
Author(s):  
J. Hossen ◽  
T. K. Pal

Efinaconazole (ECZ) is an antifungal drug. Various non-covalent interactions between ECZ and a holoenzyme (protein id: 3idb) has been investigated through computational study. The structure of ECZ was optimized using density functional theory (DFT) applying B3LYP/6-311G+(d,p) method. HOMO, LUMO, chemical hardness and softness, several thermochemical parameters, electrostatic potential surface, vibrational spectrum, total energy, and maximum internal force and maximum internal displacement with respect to optimization step number have been determined. The optimized ECZ ligand was subjected to molecular docking against the protein 3idb in Autodock Vina program. The different non-covalent interactions in the ligand-protein complex were visualized in BIOVIA Discovery Studio Visualizer. Various surface plots such as hydrogen bonds, ionizability, SAS, hydrophobicity, aromatic and charge surfaces were excerpted. ECZ molecule forms three strong hydrogen bonds with amino acid residues of the holoenzyme. In addition to this, it is significantly capable to form several other bonds which strengthen the ligand-protein interaction. The result showed that the ECZ molecule posed considerable binding affinity against the macromolecule.


Molecules ◽  
2020 ◽  
Vol 25 (11) ◽  
pp. 2675 ◽  
Author(s):  
Fabián G. Cantú Reinhard ◽  
Yen-Ting Lin ◽  
Agnieszka Stańczak ◽  
Sam P. de Visser

The cytochromes P450 are versatile enzymes found in all forms of life. Most P450s use dioxygen on a heme center to activate substrates, but one class of P450s utilizes hydrogen peroxide instead. Within the class of P450 peroxygenases, the P450 OleTJE isozyme binds fatty acid substrates and converts them into a range of products through the α-hydroxylation, β-hydroxylation and decarboxylation of the substrate. The latter produces hydrocarbon products and hence can be used as biofuels. The origin of these product distributions is unclear, and, as such, we decided to investigate substrate positioning in the active site and find out what the effect is on the chemoselectivity of the reaction. In this work we present a detailed computational study on the wild-type and engineered structures of P450 OleTJE using a combination of density functional theory and quantum mechanics/molecular mechanics methods. We initially explore the wild-type structure with a variety of methods and models and show that various substrate activation transition states are close in energy and hence small perturbations as through the protein may affect product distributions. We then engineered the protein by generating an in silico model of the double mutant Asn242Arg/Arg245Asn that moves the position of an active site Arg residue in the substrate-binding pocket that is known to form a salt-bridge with the substrate. The substrate activation by the iron(IV)-oxo heme cation radical species (Compound I) was again studied using quantum mechanics/molecular mechanics (QM/MM) methods. Dramatic differences in reactivity patterns, barrier heights and structure are seen, which shows the importance of correct substrate positioning in the protein and the effect of the second-coordination sphere on the selectivity and activity of enzymes.


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Ming-Yuan Su ◽  
Simon A. Fromm ◽  
Jonathan Remis ◽  
Daniel B. Toso ◽  
James H. Hurley

AbstractMutation of C9ORF72 is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontal temporal degeneration (FTD), which is attributed to both a gain and loss of function. C9orf72 forms a complex with SMCR8 and WDR41, which was reported to have GTPase activating protein activity toward ARF proteins, RAB8A, and RAB11A. We determined the cryo-EM structure of ARF1-GDP-BeF3- bound to C9orf72:SMCR8:WDR41. The SMCR8longin and C9orf72longin domains form the binding pocket for ARF1. One face of the C9orf72longin domain holds ARF1 in place, while the SMCR8longin positions the catalytic finger Arg147 in the ARF1 active site. Mutations in interfacial residues of ARF1 and C9orf72 reduced or eliminated GAP activity. RAB8A GAP required ~10-fold higher concentrations of the C9orf72 complex than for ARF1. These data support a specific function for the C9orf72 complex as an ARF GAP. The structure also provides a model for the active forms of the longin domain GAPs of FLCN and NPRL2 that regulate the Rag GTPases of the mTORC1 pathway.


2021 ◽  
Author(s):  
Ming-Yuan Su ◽  
Simon H Fromm ◽  
Jonathan Remis ◽  
Daniel Toso ◽  
James H Hurley

Mutation of C9ORF72 is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontal temporal degeneration (FTD), which is attributed to both a gain and loss of function. C9orf72 forms a complex with SMCR8 and WDR41, which was reported to have GTPase activating protein activity toward ARF proteins, RAB8A, and RAB11A. We determined the cryo-EM structure of ARF1-GDP-BeF3- bound to C9orf72:SMCR8:WDR41. The SMCR8longin and C9orf72longin domains form the binding pocket for ARF1. One face of the C9orf72longin domain holds ARF1 in place, while the SMCR8longin positions the catalytic finger Arg147 in the ARF1 active site. Mutations in interfacial residues of ARF1 and C9orf72 reduced or eliminated GAP activity. RAB8A GAP required ~10-fold higher concentrations of the C9orf72 complex than for ARF1. These data support a specific function for the C9orf72 complex as an ARF GAP.


2011 ◽  
Vol 64 (7) ◽  
pp. 916 ◽  
Author(s):  
Tom B. Dupree ◽  
Paul A. Keller ◽  
Renate Griffith

Superimposition of 125 non-nucleoside inhibitors from human immunodeficiency virus reverse transcriptase structures reveals a novel binding space deeper into the enzyme for some of these inhibitors, allowing access to the polymerase active site. This may enable us to design new inhibitors of this enzyme with better mutation resistance profiles. We have analysed this new binding space and have docked our in-house scaffolds into this region, highlighting the possibility of the formation of new hydrogen bonds with residues of the active site.


2020 ◽  
Author(s):  
Jon Uranga ◽  
Lukas Hasecke ◽  
Jonny Proppe ◽  
Jan Fingerhut ◽  
Ricardo A. Mata

The 20S Proteasome is a macromolecule responsible for the chemical step in the ubiquitin-proteasome system of degrading unnecessary and unused proteins of the cell. It plays a central role both in the rapid growth of cancer cells as well as in viral infection cycles. Herein, we present a computational study of the acid-base equilibria in an active site of the human proteasome, an aspect which is often neglected despite the crucial role protons play in the catalysis. As example substrates, we take the inhibition by epoxy and boronic acid containing warheads. We have combined cluster quantum mechanical calculations, replica exchange molecular dynamics and Bayesian optimization of non-bonded potential terms in the inhibitors. In relation to the latter, we propose an easily scalable approach to the reevaluation of non-bonded potentials making use of QM/MM dynamics information. Our results show that coupled acid-base equilibria need to be considered when modeling the inhibition mechanism. The coupling between a neighboring lysine and the reacting threonine is not affected by the presence of the inhibitor.


2019 ◽  
Author(s):  
Michael Fischer

<p>The silicoaluminophosphate zeotype ECR-40, which has the MEI topology, contains linkages of AlO<sub>4</sub> tetrahedra via a common oxygen atom, thereby violating the famous “Löwenstein’s rule”. Due to the proven existence of Al-O-Al linkages in this material, it constitutes an ideal model system to study the acidity and mobility of protons associated with such unusual linkages. In addition, their properties can be directly compared to those of protons associated with more common Si-O-Al linkages, which are also present in ECR-40. In this work, static density functional theory (DFT) calculations including a dispersion correction were employed to study the preferred proton sites as well as the Brønsted acidity of the framework protons, followed by DFT-based ab-initio molecular dynamics (AIMD) to investigate the proton mobility in guest-free and hydrated ECR-40. Initially, two different proton arrangements were compared, one containing both H[O6] protons associated with Al-O-Al linkages and H[O10] protons at Si-O-Al linkages, the other one containing only H[O10] protons. The former model was found to be thermodynamically favoured, as a removal of protons from the Al-O-Al linkages causes a local accumulation of negative charge. Calculations of the deprotonation energy showed a moderately higher Brønsted acidity of the H[O10] protons, at variance with previous empirical explanations, which attributed the exceptional performance of ECR-40 as acid catalyst to the presence of Al‑O‑Al linkages. The AIMD simulations (<i>T</i> = 298 K) delivered no appreciable proton mobility for guest-free ECR-40 and for low levels of hydration (one H<sub>2</sub>O per framework proton). Under saturation conditions, framework deprotonation occurred, leading to the formation of protonated water clusters in the pores. Pronounced differences between the two types of framework protons were observed: While the H[O10] protons were always removed from the Si-O-Al linkages, the Al-O-Al linkages remained mostly protonated, but deprotonation did occur to a minor extent. The observation of a degree of framework deprotonation of Al-O-Al linkages differs from the findings reported in a recent computational study of hydrated aluminosilicate zeolites with such linkages (Heard et al., <i>Chem. Sci.</i> <b>2019</b>, <i>10</i>, 5705), pointing to an influence of the overall framework composition. Further inspection of the AIMD results showed that a coordination of water molecules to framework Al atoms occurred in many cases, especially in the vicinity of the Al-O-Al linkages, sometimes resulting in a pronounced modification of the linkages through additional bridging oxygen atoms. Given the changes in the local structure, it can be expected that such modified linkages are especially prone to break upon dehydration. Thus, in addition to elucidating the deprotonation behaviour of protons associated with different types of linkages, the calculations also provide insights into possible reasons for the instability of Al-O-Al linkages, clarifying why Löwenstein’s rule is mostly obeyed in materials that are formed via a hydrothermal route.</p>


Sign in / Sign up

Export Citation Format

Share Document