Computational study of conformational changes in human 3-hydroxy-3-methylglutaryl coenzyme reductase induced by substrate binding

Pharmacophore Modeling and Binding Affinity of Secondary Metabolites from Angelica keiskei to HMG Co-A Reductase

Statins are cholesterol-lowering drugs with a mechanism of inhibiting 3-hydroxy-3-methylglutaryl-CoA reductase, but long-term use can cause side effects. An example of a plant capable of reducing cholesterol levels is Angelica keiskei (ashitaba). Therefore, this study aimed to obtain suitable compounds with inhibitory activity against the HMG-CoA reductase enzyme from ashitaba through in silico tests. The experiment began with screening and pharmacophore modeling, followed by molecular docking on ashitaba's compounds, statins groups, and the native ligand was (3R,5R)-7-[4-(benzyl carbamoyl)-2-(4-fluorophenyl)-5-(1-methylethyl)-1H-imidazole-1-yl]-3,5-dihydroxyheptanoic acid (4HI). Based on the results of the molecular docking simulations, 15 hit compounds had a small binding energy (ΔG). Pitavastatin, as the comparator drug (ΔG = -8.24 kcal/mol; Ki = 2.11 µM), had a lower ΔG and inhibition constant (Ki) than the native ligand 4HI (ΔG = -7.84 kcal/mol; Ki = 7.96µM). From ashitaba's compounds, it was found that 4'-O-geranylnaringenin, luteolin, isobavachalcone, dorsmannin A, and 3'-carboxymethyl-4,2'-dihydroxy-4'-methoxychalcone have low ΔG of below -6 kcal/mol. The lowest ΔG value was found in 3'-carboxymethyl-4,2'-dihydroxy-4'-methoxy chalcone with a ΔG of -6.67 kcal/mol and Ki value of 16.66 µM, which was lower than the ΔG value of the other comparator drugs, atorvastatin (ΔG = -5.49 kcal/mol; Ki = 1148.17 µM) and simvastatin (ΔG = -6.50 kcal/mol; Ki = 22.34 µM). This compound also binds to the important amino acid residues, including ASN755D, ASP690C, GLU559D, LYS735D, LYS691C, and SER684C, through hydrogen bonds. Based on the results, the compound effectively binds to six important amino acids with good binding affinity and only requires a small concentration to reduce half of the enzyme activity.

Evolving Strategies for Use of Phytochemicals in Prevention and Long-Term Management of Cardiovascular Diseases (CVD)

This report describes major pathomechanisms of disease in which the dysregulation of host inflammatory processes is a major factor, with cardiovascular disease (CVD) as a primary model, and reviews strategies for countermeasures based on synergistic interaction between various agents, including drugs and generally regarded as safe (GRAS) natural medical material (NMM), such as Ginkgo biloba, spice phytochemicals, and fruit seed flavonoids. The 15 well-defined CVD classes are explored with particular emphasis on the extent to which oxidative stressors and associated ischemia-reperfusion tissue injury contribute to major symptoms. The four major categories of pharmaceutical agents used for the prevention of and therapy for CVD: statins, beta blockers (β-blockers), blood thinners (anticoagulants), and aspirin, are presented along with their adverse effects. Analyses of major cellular and molecular features of drug- and NMM-mediated cardioprotective processes are provided in the context of their development for human clinical application. Future directions of the evolving research described here will be particularly focused on the characterization and manipulation of calcium- and calcineurin-mediated cascades of signaling from cell surface receptors on cardiovascular and immune cells to the nucleus, with the emergence of both protective and pathological epigenetic features that may be modulated by synergistically-acting combinations of drugs and phytochemicals in which phytochemicals interact with cells to promote signaling that reduces the effective dosage and thus (often) toxicity of drugs.

Enhancement of Squalene Production by Constitutive Expression of the 3-Hydroxy-3-Methylglutaryl-CoA Reductase in Aurantiochytrium sp. 18W-13a

Squalene has a wide range of applications in the industry sectors of dietary supplements, cosmetics, immunization, and pharmaceuticals. Yet, suitable organisms as the source of squalene are limited. It is reported that the thraustochytrid Aurantiochytrium sp. strain 18W-13a can accumulate high content of squalene. However, squalene production in this organism is fluctuated under various conditions and is not yet optimized for commercialization. In this organism, the mevalonate pathway supplies isopentenyl pyrophosphate, the initial substrate for squalene production. In this pathway, 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR) is the rate-limiting enzyme. We found that the HMGR activity had a strong positive correlation with the squalene contents in the strain. We constitutively expressed the HMGR in this organism and found that the transformant showed increased and stable production of squalene as well as carotenoids and biomass. These results clearly indicated that the HMGR expression is the bottleneck of squalene synthesis in Aurantiochytrium sp.

Assessment of mutations on RBD in the Spike protein of SARS-CoV-2 Alpha, Delta and Omicron variants

The severe acute respiratory syndrome (SARS) coronavirus 2 (CoV-2) variant Omicron spread more rapid than the other variants of SARS-CoV-2 virus. Mutations on the Spike (S) protein receptor-binding domain (RBD) are critical for the antibody resistance and infectivity of the SARS-CoV-2 variants. In this study, we have used accelerated molecular dynamics (aMD) simulations and free energy calculations to present a systematic analysis of the affinity and conformational dynamics along with the interactions that drive the binding between Spike protein RBD and human angiotensin-converting enzyme 2 (ACE2) receptor. We evaluate the impacts of the key mutation that occur in the RBDs Omicron and other variants in the binding with the human ACE2 receptor. The results show that S protein Omicron has stronger binding to the ACE2 than other variants. The evaluation of the decomposition energy per residue shows the mutations N440K, T478K, Q493R and Q498R observed in Spike protein of SARS-CoV-2 provided a stabilization effect for the interaction between the SARS-CoV-2 RBD and ACE2. Overall, the results demonstrate that faster spreading of SARS-CoV-2 Omicron may be correlated with binding affinity of S protein RBD to ACE2 and mutations of uncharged residues to positively charged residues such as Lys and Arg in key positions in the RBD.

Biosynthesis and regulation of terpenoids from basidiomycetes: exploration of new research

Basidiomycetes, also known as club fungi, consist of a specific group of fungi. Basidiomycetes produce a large number of secondary metabolites, of which sesquiterpenoids, diterpenoids and triterpenoids are the primary components. However, these terpenoids tend to be present in low amounts, which makes it difficult to meet application requirements. Terpenoid biosynthesis improves the quantity of these secondary metabolites. However, current understanding of the biosynthetic mechanism of terpenoids in basidiomycetes is insufficient. Therefore, this article reviews the latest research on the biosynthesis of terpenoids in basidiomycetes and summarizes the CYP450 involved in the biosynthesis of terpenoids in basidiomycetes. We also propose opportunities and challenges for chassis microbial heterologous production of terpenoids in basidiomycetes and provide a reference basis for the better development of basidiomycete engineering.

Assessment of the PETase conformational changes induced by poly(ethylene terephthalate) binding

Recently, a bacterium strain of Ideonella sakaiensis was identified with the uncommon ability to degrade the poly(ethylene terephthalate) (PET). The PETase from I. sakaiensis strain 201-F6 (IsPETase) catalyzes the hydrolysis of PET converting it to mono(2-hydroxyethyl) terephthalic acid (MHET), bis(2-hydroxyethyl)-TPA (BHET), and terephthalic acid (TPA). Despite the potential of this enzyme for mitigation or elimination of environmental contaminants, one of the limitations of the use of IsPETase for PET degradation is the fact that it acts only at moderate temperature due to its low thermal stability. Besides, molecular details of the main interactions of PET in the active site of IsPETase remain unclear. Herein, molecular docking and molecular dynamics (MD) simulations were applied to analyze structural changes of IsPETase induced by PET binding. Results from the essential dynamics revealed that the β1-β2 connecting loop is very flexible. This loop is located far from the active site of IsPETase and we suggest that it can be considered for mutagenesis to increase the thermal stability of IsPETase. The free energy landscape (FEL) demonstrates that the main change in the transition between the unbound to the bound state is associated with the β7-α5 connecting loop, where the catalytic residue Asp206 is located. Overall, the present study provides insights into the molecular binding mechanism of PET into the IsPETase structure and a computational strategy for mapping flexible regions of this enzyme, which can be useful for the engineering of more efficient enzymes for recycling plastic polymers using biological systems.

Gonçalves Vianez JL, Santana da Costa K, de Molfetta FA, Nahum Alves C. Analysis of Kojic Acid Derivatives as Competitive Inhibitors of Tyrosinase: A Molecular Modeling Approach

Tyrosinases belong to the functional copper-containing proteins family, and their structure contains two copper atoms, in the active site, which are coordinated by three histidine residues. The biosynthesis of melanin in melanocytes has two stages depending on the actions of the natural substrates L-DOPA and L-tyrosine. The dysregulation of tyrosinase is involved in skin cancer initiation. In the present study, using molecular modeling tools, we analyzed the inhibition activity of tyrosinase activity using kojic acid (KA) derivatives designed from aromatic aldehydes and malononitrile. All derivatives showed conformational affinity to the enzyme active site, and a favorable distance to chelate the copper ion, which is essential for enzyme function. Molecular dynamics simulations revealed that the derivatives formed promising complexes, presenting stable conformations with deviations between 0.2 and 0.35 Å. In addition, the investigated KA derivatives showed favorable binding free energies. The most stable KA derivatives showed the following binding free energies: -17.65 kcal mol-1 (D6), -18.07 kcal mol-1 (D2), -18.13 (D5) kcal mol-1, and -10.31 kcal mol-1 (D4). Our results suggest that these derivatives could be potent competitive inhibitors of the natural substrates of L-DOPA (-12.84 kcal mol-1) and L-tyrosine (-9.04 kcal mol-1) in melanogenesis.

A-family anti-inflammatory cyclopentenone prostaglandins: A novel class of non-statin inhibitors of HMG-CoA reductase

Disruption of the intracellular lipid balance leading to cholesterol accumulation is one of the features of cells that participate in the development of atherosclerotic lesions. Evidence form our laboratory indicates that anti-inflammatory cyclopentenone prostaglandins (cyPGs) of A- and J-family deviate lipid metabolism from the synthesis of cholesterol and cholesteryl esters to the synthesis of phospholipids in foam-cell macrophages. cyPGs possessing an α,β-unsaturated cyclopentane ring are highly electrophilic substances able to promptly react with reactive cysteines of intracellular molecules through Michael addition. On the other hand, HMG-CoA reductase (HMGCR), the enzyme responsible for the rate-limiting step in cholesterol biosynthesis, presents critically reactive cysteines at the entry of catalytic domain, particularly Cys561, that could be target of cyPG inhibition. In the present study, we showed that cyPGs (but not other non-α,β-unsaturated PGs) physically interact with HMGCR, in a dithiothreitol- and β-mercaptoethanol-sensitive way, and block the activity of the catalytic subunit of the enzyme (IC₅₀ for PGA₂ = 0.17 μM). PGA₂ inhibits HMGCR activity in cultured rat and human macrophages/macrophage-foam cells and leads to enhanced expression of HMGCR protein, as observed with statins. In cell culture models, PGA₂ effectively inhibits the reductase at non-toxic doses (e.g., 1 μM) that block cell proliferation thus suggesting that part of the well-known antiproliferative effect of PGA₂ may be due to its ability of blocking HMGCR activity, as cells cannot proliferate without a robust cholesterogenesis. Therefore, besides the powerfully anti-inflammatory and antiproliferative effects, the anticholesterogenic effects of PGA₂ should be exploited in atherosclerosis therapeutics.

Using Accelerated Molecular Dynamics Simulation to elucidate the effects of the T198F mutation on the molecular flexibility of the West Nile virus envelope protein

The envelope (E) protein is an important target for antibodies in flavivirus. Literature reports that the mutation T198F, located at the domain I-II hinge of the E protein, regulates viral breathing and increases the accessibility of a distal cryptic epitope located on the fusion loop, having a direct impact in the neutralization of West Nile virus (WNV). Our study aimed to describe, using accelerated molecular dynamics simulations, the effects of the T198F mutation in the flexibility of the E protein of WNV and to elucidate the mechanism that regulates epitope accessibility. The simulation results revealed that the mutation favors the formation of alternative hydrogen bonds, hampering the bending movement between domains I and II. We hypothesized that this is the mechanism by which the T198F mutation, located at the middle of the protein, locks the distal cryptc epitope near a single preferred conformation, rendering it more prone to recognition by antibodies.

Unraveling the conformational dynamics of glycerol 3-phosphate dehydrogenase, a nicotinamide adenine dinucleotide-dependent enzyme of Leishmania mexicana

Allosteric changes modulate the enzymatic activity, leading to activation or inhibition of the molecular target. Understanding the induced fit accommodation mechanism of a ligand in its lowest-free energy state and the subsequent conformational changes induced in the protein are important questions for drug design. In the present study, molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA) were applied to analyze the glycerol-3-phosphate dehydrogenase of Leishmania mexicana (LmGPDH) conformational changes induced by its cofactor and substrate binding. GPDH is a nicotinamide adenine dinucleotide (NAD)-dependent enzyme, which has been reported as an interesting target for drug discovery and development against leishmaniasis. Despite its relevance for glycolysis and pentose phosphate pathways, the structural flexibility and conformational motions of LmGPDH in complex with NADH and dihydroxyacetone phosphate (DHAP) remain unexplored. Here, we analyzed the conformational dynamics of the enzyme-NADH complex (cofactor), and the enzyme-NADH-DHAP complex (adduct), mapped the hydrogen-bond interactions for the complexes and pointed some structural determinants of the enzyme that emerge from these contacts to NADH and DHAP. Finally, we proposed a consistent mechanism for the conformational changes on the first step of the reversible redox conversion of dihydroxyacetone phosphate to glycerol 3-phosphate, indicating key residues and interactions that could be further explored in drug discovery.

Investigation of the target-site resistance of EPSP synthase mutants P106T and T102I/P106S against glyphosate

The shikimate pathway enzyme 5-enolpyruvyl shikimate-3-phosphate synthase (EPSPS) catalyzes the reaction involved in the production of amino acids essential for plant growth and survival. Thus, EPSPS is the main target of various herbicides, including glyphosate, a broad-spectrum herbicide that acts as a competitive inhibitor of phosphoenolpyruvate (PEP), which is the natural substrate of EPSPS. However, punctual mutations in the EPSPS gene have led to glyphosate resistance in some plants. Here, we investigated the mechanism of EPSPS resistance to glyphosate in mutants of two weed species, Conyza sumatrensis (mutant, P106T) and Eleusine indica (mutant, T102I/P106S), both of which have an economic impact on industrial crops. Molecular dynamics (MD) simulations and binding free energy calculations revealed the influence of the mutations on the affinity of glyphosate in the PEP-binding site. The amino acid residues of the EPSPS protein in both species involved in glyphosate resistance were elucidated as well as other residues that could be useful for protein engineering. In addition, during MD simulations, we identified conformational changes in glyphosate when complexed with resistant EPSPS, related to loss of herbicide activity and binding affinity. Our computational findings are consistent with previous experimental results and clarify the inhibitory activity of glyphosate as well as the structural target-site resistance of EPSPS against glyphosate.

Single or double EPSP synthase mutations lead glyphosate to undergo conformational changes that limit its inhibitory action.

Exploring the Potentiality of Natural Products from Essential Oils as Inhibitors of Odorant-Binding Proteins: A Structure- and Ligand-Based Virtual Screening Approach To Find Novel Mosquito Repellents

Odorant-binding proteins (OBPs) are the main olfactory proteins of mosquitoes, and their structures have been widely explored to develop new repellents. In the present study, we combined ligand- and structure-based virtual screening approaches using as a starting point 1633 compounds from 71 botanical families obtained from the Essential Oil Database (EssOilDB). Using as reference the crystallographic structure of N,N-diethyl-meta-toluamide interacting with the OBP1 homodimer of Anopheles gambiae (AgamOBP1), we performed a structural and pharmacophoric similarity search to select potential natural products from the library. Thymol acetate, 4-(4-methyl phenyl)-pentanal, thymyl isovalerate, and p-cymen-8-yl demonstrated a favorable chemical correlation with DEET and also had high-affinity interactions with the OBP binding pocket that molecular dynamics simulations showed to be stable. To the best of our knowledge, this is the first study to evaluate on a large scale the potentiality of NPs from essential oils as inhibitors of the mosquito OBP1 using in silico approaches. Our results could facilitate the design of novel repellents with improved selectivity and affinity to the protein binding pocket and can shed light on the mechanism of action of these compounds against insect olfactory recognition.

Studying effects of different protonation states of His11 and His102 in ribose-5-phosphate isomerase of Trypanosoma cruzi: an example of cooperative behavior

The Trypanosoma cruzi ribose-5-phosphate isomerase B (TcRpiB) is a crucial piece in the pentose phosphate pathway and thus is a potential drug target for treatment of Chagas' disease. TcRpiB residues, such as Cys69, Asp45, Glu149 and Pro47, have confirmed their roles in substrate recognition, catalytic reaction and binding site conformation. However, the joint performance of His11 and His102, in the D-ribose-5-phosphate (R5P) in the catalysis is not well understood. In this work, we probed the influence of different protonation states of His11 and His102 on the behavior of the ligand R5P using molecular dynamics simulations, network analysis and thermodynamic integration. Simulations revealed that a protonated His11 combined with a neutral His102 (His11⁺‒His102) was able to stabilize the ligand R5P in the binding site. Moreover, calculated relative free energy differences showed that when protonated His11 was coupled to a neutral His102 an exergonic process takes place. On the other hand, neutral His11 combined with a protonated His102 (His11‒His102⁺), sampled conformations that resembled the catalyzed product D-ribulose-5-phosphate (Ru5P). Network analysis also demonstrated some peculiarities for these systems with some negatively correlated nodes in the binding site for His11‒His102⁺, and exclusive suboptimal paths for His11⁺‒His102. Therefore, the combined approach presented in this paper proposes two suitable protonation states for the TcRpiB catalytic mechanism, where an extra proton in either histidines might favor R5P binding or influence isomerization reaction to Ru5P. Our results may guide further in silico drug discovery studies. Communicated by Ramaswamy H. Sarma.

Measuring the structural impact of mutations on cytochrome P450 21A2, the major steroid 21-hydroxylase related to congenital adrenal hyperplasia

Congenital adrenal hyperplasia is an inherited autosomal recessive disorder related to deficient cortisol synthesis. The deficiency of steroid 21-hydroxylase (cytochrome P450 21A2), an enzyme involved in cortisol synthesis, is responsible for ∼95% of cases of congenital adrenal hyperplasia. This metabolic disease exhibits three clinical forms: salt-wasting, simple virilizing, and non-classical form, which are divided according to the degree of severity. In the present study, structural and mutational analyses were performed in order to identify the structural impact of mutations on cytochrome P450 21A2 and correlate them with patient clinical severity. The following mutations were selected: arginine-356 to tryptophan (R356W), proline-30 to leucine (P30L), isoleucine-172 to asparagine (I172N), valine-281 to leucine (V281L), and the null mutation glutamine-318 (Q318X). Our computational approach mapped the location of residues on P450 and identified their implications on enzyme electrostatic potential mapping to progesterone and heme binding pockets. Using molecular dynamics simulations, we analyzed the structural stability of ligand binding and protein structure, as well as possible conformational changes at the catalytic pocket that leads to impairment of enzymatic activity. Our study sheds light on the impact structural mutations have over steroid 21-hydroxylase structure-function in the cell.Communicated by Ramaswamy H. Sarma.