scholarly journals A QM/MM Evaluation of the Missing Step in the Reduction Mechanism of HMG-CoA by Human HMG-CoA Reductase

Processes ◽  
2021 ◽  
Vol 9 (7) ◽  
pp. 1085
Author(s):  
Paula Mihaljević-Jurič ◽  
Sérgio F. Sousa

Statins are important drugs in the regulation of cholesterol levels in the human body that have as a primary target the enzyme β-hydroxy-β-methylglutaryl-CoA reductase (HMGR). This enzyme plays a crucial role in the mevalonate pathway, catalyzing the four-electron reduction of HMG-CoA to mevalonate. A second reduction step of this reaction mechanism has been the subject of much speculation in the literature, with different conflicting theories persisting to the present day. In this study, the different mechanistic hypotheses were evaluated with atomic-level detail through a combination of molecular dynamics simulations (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations. The obtained Gibbs free activation and Gibbs free reaction energy (15 kcal mol−1 and −40 kcal mol−1) show that this hydride step takes place with the involvement of a cationic His405 and Lys639, and a neutral Glu98, while Asp715 remains in an anionic state. The results provide an atomic-level portrait of this step, clearly demonstrating the nature and protonation state of the amino acid residues involved, the energetics associated, and the structure and charge of the key participating atoms in the several intermediate states, finally elucidating this missing step.

1998 ◽  
Vol 9 (9) ◽  
pp. 2611-2626 ◽  
Author(s):  
Richard Gardner ◽  
Stephen Cronin ◽  
Benjamin Leader ◽  
Jasper Rine ◽  
Randolph Hampton

The degradation rate of 3-hydroxy-3-methylglutaryl CoA reductase (HMG-R), a key enzyme of the mevalonate pathway, is regulated through a feedback mechanism by the mevalonate pathway. To discover the intrinsic determinants involved in the regulated degradation of the yeast HMG-R isozyme Hmg2p, we replaced small regions of the Hmg2p transmembrane domain with the corresponding regions from the other, stable yeast HMG-R isozyme Hmg1p. When the first 26 amino acids of Hmg2p were replaced with the same region from Hmg1p, Hmg2p was stabilized. The stability of this mutant was not due to mislocalization, but rather to an inability to be recognized for degradation. When amino acid residues 27–54 of Hmg2p were replaced with those from Hmg1p, the mutant was still degraded, but its degradation rate was poorly regulated. The degradation of this mutant was still dependent on the first 26 amino acid residues and on the function of the HRD genes. These mutants showed altered ubiquitination levels that were well correlated with their degradative phenotypes. Neither determinant was sufficient to impart regulated degradation to Hmg1p. These studies provide evidence that there are sequence determinants in Hmg2p necessary for degradation and optimal regulation, and that independent processes may be involved in Hmg2p degradation and its regulation.


2020 ◽  
Vol 16 (4) ◽  
pp. 451-459 ◽  
Author(s):  
Fortunatus C. Ezebuo ◽  
Ikemefuna C. Uzochukwu

Background: Sulfotransferase family comprises key enzymes involved in drug metabolism. Oxamniquine is a pro-drug converted into its active form by schistosomal sulfotransferase. The conformational dynamics of side-chain amino acid residues at the binding site of schistosomal sulfotransferase towards activation of oxamniquine has not received attention. Objective: The study investigated the conformational dynamics of binding site residues in free and oxamniquine bound schistosomal sulfotransferase systems and their contribution to the mechanism of oxamniquine activation by schistosomal sulfotransferase using molecular dynamics simulations and binding energy calculations. Methods: Schistosomal sulfotransferase was obtained from Protein Data Bank and both the free and oxamniquine bound forms were subjected to molecular dynamics simulations using GROMACS-4.5.5 after modeling it’s missing amino acid residues with SWISS-MODEL. Amino acid residues at its binding site for oxamniquine was determined and used for Principal Component Analysis and calculations of side-chain dihedrals. In addition, binding energy of the oxamniquine bound system was calculated using g_MMPBSA. Results: The results showed that binding site amino acid residues in free and oxamniquine bound sulfotransferase sampled different conformational space involving several rotameric states. Importantly, Phe45, Ile145 and Leu241 generated newly induced conformations, whereas Phe41 exhibited shift in equilibrium of its conformational distribution. In addition, the result showed binding energy of -130.091 ± 8.800 KJ/mol and Phe45 contributed -9.8576 KJ/mol. Conclusion: The results showed that schistosomal sulfotransferase binds oxamniquine by relying on hybrid mechanism of induced fit and conformational selection models. The findings offer new insight into sulfotransferase engineering and design of new drugs that target sulfotransferase.


RSC Advances ◽  
2014 ◽  
Vol 4 (89) ◽  
pp. 48621-48631 ◽  
Author(s):  
Eleanor R. Turpin ◽  
Sam Mulholland ◽  
Andrew M. Teale ◽  
Boyan B. Bonev ◽  
Jonathan D. Hirst

2017 ◽  
Vol 114 (42) ◽  
pp. E8830-E8836 ◽  
Author(s):  
Chang Yun Son ◽  
Arun Yethiraj ◽  
Qiang Cui

Cytochrome c oxidase (CcO) is a transmembrane protein that uses the free energy of O2 reduction to generate the proton concentration gradient across the membrane. The regulation of competitive proton transfer pathways has been established to be essential to the vectorial transport efficiency of CcO, yet the underlying mechanism at the molecular level remains lacking. Recent studies have highlighted the potential importance of hydration-level change in an internal cavity that connects the proton entrance channel, the site of O2 reduction, and the putative proton exit route. In this work, we use atomistic molecular dynamics simulations to investigate the energetics and timescales associated with the volume fluctuation and hydration-level change in this central cavity. Extensive unrestrained molecular dynamics simulations (accumulatively ∼4 μs) and free energy computations for different chemical states of CcO support a model in which the volume and hydration level of the cavity are regulated by the protonation state of a propionate group of heme a3 and, to a lesser degree, the redox state of heme a and protonation state of Glu286. Markov-state model analysis of ∼2-μs trajectories suggests that hydration-level change occurs on the timescale of 100–200 ns before the proton-loading site is protonated. The computed energetic and kinetic features for the cavity wetting transition suggest that reversible hydration-level change of the cavity can indeed be a key factor that regulates the branching of proton transfer events and therefore contributes to the vectorial efficiency of proton transport.


Blood ◽  
2003 ◽  
Vol 102 (9) ◽  
pp. 3354-3362 ◽  
Author(s):  
Niels W. C. J. van de Donk ◽  
Marloes M. J. Kamphuis ◽  
Berris van Kessel ◽  
Henk M. Lokhorst ◽  
Andries C. Bloem

AbstractHMG-CoA reductase is the rate-limiting enzyme of the mevalonate pathway leading to the formation of cholesterol and isoprenoids such as farnesylpyrophosphate (FPP) and geranylgeranylpyrophosphate (GGPP). The inhibition of HMG-CoA reductase by lovastatin induced apoptosis in plasma cell lines and tumor cells from patients with multiple myeloma. Here we show that cotreatment with mevalonate or geranylgeranyl moieties, but not farnesyl groups, rescued myeloma cells from lovastatin-induced apoptosis. In addition, the inhibition of geranylgeranylation by specific inhibition of geranylgeranyl transferase I (GGTase I) induced the apoptosis of myeloma cells. Apoptosis triggered by the inhibition of geranylgeranylation was associated with reduction of Mcl-1 protein expression, collapse of the mitochondrial transmembrane potential, expression of the mitochondrial membrane protein 7A6, cytochrome c release from mitochondria into the cytosol, and stimulation of caspase-3 activity. These results imply that protein geranylgeranylation is critical for regulating myeloma tumor cell survival, possibly through regulating Mcl-1 expression. Our results show that pharmacologic agents such as lovastatin or GGTase inhibitors may be useful in the treatment of multiple myeloma.


2019 ◽  
Vol 24 (9) ◽  
pp. 928-938 ◽  
Author(s):  
Luca Palazzolo ◽  
Chiara Paravicini ◽  
Tommaso Laurenzi ◽  
Sara Adobati ◽  
Simona Saporiti ◽  
...  

SLC6A14 (ATB0,+) is a sodium- and chloride-dependent neutral and dibasic amino acid transporter that regulates the distribution of amino acids across cell membranes. The transporter is overexpressed in many human cancers characterized by an increased demand for amino acids; as such, it was recently acknowledged as a novel target for cancer therapy. The knowledge on the molecular mechanism of SLC6A14 transport is still limited, but some elegant studies on related transporters report the involvement of the 12 transmembrane α-helices in the transport mechanism, and describe structural rearrangements mediated by electrostatic interactions with some pivotal gating residues. In the present work, we constructed a SLC6A14 model in outward-facing conformation via homology modeling and used molecular dynamics simulations to predict amino acid residues critical for substrate recognition and translocation. We docked the proteinogenic amino acids and other known substrates in the SLC6A14 binding site to study both gating regions and the exposed residues involved in transport. Interestingly, some of these residues correspond to those previously identified in other LeuT-fold transporters; however, we could also identify a novel relevant residue with such function. For the first time, by combined approaches of molecular docking and molecular dynamics simulations, we highlight the potential role of these residues in neutral amino acid transport. This novel information unravels new aspects of the human SLC6A14 structure–function relationship and may have important outcomes for cancer treatment through the design of novel inhibitors of SLC6A14-mediated transport.


2019 ◽  
Vol 21 (1) ◽  
pp. 1-31 ◽  
Author(s):  
Lindong Weng ◽  
Shannon L. Stott ◽  
Mehmet Toner

Successful stabilization and preservation of biological materials often utilize low temperatures and dehydration to arrest molecular motion. Cryoprotectants are routinely employed to help the biological entities survive the physicochemical and mechanical stresses induced by cold or dryness. Molecular interactions between biomolecules, cryoprotectants, and water fundamentally determine the outcomes of preservation. The optimization of assays using the empirical approach is often limited in structural and temporal resolution, whereas classical molecular dynamics simulations can provide a cost-effective glimpse into the atomic-level structure and interaction of individual molecules that dictate macroscopic behavior. Computational research on biomolecules, cryoprotectants, and water has provided invaluable insights into the development of new cryoprotectants and the optimization of preservation methods. We describe the rapidly evolving state of the art of molecular simulations of these complex systems, summarize the molecular-scale protective and stabilizing mechanisms, and discuss the challenges that motivate continued innovation in this field.


Sign in / Sign up

Export Citation Format

Share Document