Molecular dynamics-based identification of binding pathways and two distinct high-affinity sites for succinate in succinate receptor 1/GPR91.

SUCNR1 is an auto- and paracrine sensor of the metabolic stress signal succinate. Using unsupervised molecular dynamics (MD) simulations (170.400 ns) and mutagenesis across human, mouse, and rat SUCNR1, we characterize how a five-arginine motif around the extracellular pole of TM-VI determines the initial capture of succinate in the extracellular vestibule (ECV) to either stay or move down to the orthosteric site. Metadynamics demonstrate low-energy succinate binding in both sites, with an energy barrier corresponding to an intermediate stage during which succinate, with an associated water cluster, unlocks the hydrogen-bond-stabilized conformationally constrained extracellular loop (ECL)-2b. Importantly, simultaneous binding of two succinate molecules through either a "sequential" or "bypassing" mode is a frequent endpoint. The mono-carboxylate NF-56-EJ40 antagonist enters SUCNR1 between TM-I and -II and does not unlock ECL-2b. It is proposed that occupancy of both high-affinity sites is required for selective activation of SUCNR1 by high local succinate concentrations.

The structure of a C. neoformans polysaccharide motif recognized by protective antibodies: A combined NMR and MD study.

Cryptococcus neoformans is a fungal pathogen responsible for cryptococcosis and cryptococcal meningitis. The C. neoformans' capsular polysaccharide and its shed exopolysaccharide function both as key virulence factors and to protect the fungal cell from phagocytosis. Currently, a glycoconjugate of these polysaccharides is being explored as a vaccine to protect against C. neoformans infection. In this study, NOE and J-coupling values from NMR experiments were consistent with a converged structure of the synthetic decasaccharide, GXM10-Ac3, calculated from MD simulations. GXM10-Ac3 was designed as an extension of glucuronoxylomannan (GXM) polysaccharide motif (M2) which is common in the clinically predominant serotype A strains and is recognized by protective forms of GXM-specific monoclonal antibodies. The M2 motif is a hexasaccharide with a three-residue α-mannan backbone, modified by ß-(1→2)-xyloses (Xyl) on the first two mannoses (Man) and a ß-(1→2)-glucuronic acid (GlcA) on the third Man. Combined NMR and MD analyses reveal that GXM10-Ac3 adopts an extended structure, with Xyl/GlcA branches alternating sides along the α-mannan backbone. O-acetyl esters also alternate sides and are grouped in pairs. MD analysis of a twelve M2-repeating unit polymer supports the notion that the GXM10-Ac3 structure is uniformly represented throughout the polysaccharide. This derived GXM model displays high flexibility while maintaining a structural identity, yielding insights to further explore intermolecular interactions between polysaccharides, interactions with anti-GXM mAbs, and the cryptococcal polysaccharide architecture.

Conformational dynamics underlying atypical chemokine receptor 3 activation.

Atypical Chemokine Receptor 3 (ACKR3) belongs to the G protein-coupled receptor family but it does not signal through G proteins. The structural properties that govern the functional selectivity and the conformational dynamics of ACKR3 activation are poorly understood. Here, we combined hydrogen/deuterium exchange mass spectrometry, site-directed mutagenesis, and molecular dynamics simulations to examine the binding mode and mechanism of action of ACKR3 ligands of different efficacies. Our results show that activation or inhibition of ACKR3 is governed by intracellular conformational changes of helix 6, intracellular loop 2, and helix 7, while the DRY motif becomes protected during both processes. Moreover, we identified the binding sites and the allosteric modulation of ACKR3 upon ß-arrestin 1 binding. In summary, this study highlights the structure-function relationship of small ligands, the binding mode of ß-arrestin 1, the activation dynamics, and the atypical dynamic features in ACKR3 that may contribute to its inability to activate G proteins.

Enhancing protein dynamics analysis with hydrophilic polyethylene glycol cross-linkers.

Cross-linkers play a critical role in capturing protein dynamics in chemical cross-linking mass spectrometry techniques. Various types of cross-linkers with different backbone features are widely used in the study of proteins. However, it is still not clear how the cross-linkers' backbone affect their own structure and their interactions with proteins. In this study, we systematically characterized and compared methylene backbone and polyethylene glycol (PEG) backbone cross-linkers in terms of capturing protein structure and dynamics. The results indicate the cross-linker with PEG backbone have a better ability to capture the inter-domain dynamics of calmodulin, adenylate kinase, maltodextrin binding protein and dual-specificity protein phosphatase. We further conducted quantum chemical calculations and all-atom molecular dynamics simulations to analyze thermodynamic and kinetic properties of PEG backbone and methylene backbone cross-linkers. Solution nuclear magnetic resonance was employed to validate the interaction interface between proteins and cross-linkers. Our findings suggest that the polarity distribution of PEG backbone enhances the accessibility of the cross-linker to the protein surface, facilitating the capture of sites located in dynamic regions. By comprehensively benchmarking with disuccinimidyl suberate (DSS)/bis-sulfosuccinimidyl-suberate(BS3), bis-succinimidyl-(PEG)2 revealed superior advantages in protein dynamic conformation analysis in vitro and in vivo, enabling the capture of a greater number of cross-linking sites and better modeling of protein dynamics. Furthermore, our study provides valuable guidance for the development and application of PEG backbone cross-linkers.

Pressure pushes tRNALys3 into excited conformational states.

Conformational dynamics play essential roles in RNA function. However, detailed structural characterization of excited states of RNA remains challenging. Here, we apply high hydrostatic pressure (HP) to populate excited conformational states of tRNALys3, and structurally characterize them using a combination of HP 2D-NMR, HP-SAXS (HP-small-angle X-ray scattering), and computational modeling. HP-NMR revealed that pressure disrupts the interactions of the imino protons of the uridine and guanosine U-A and G-C base pairs of tRNALys3. HP-SAXS profiles showed a change in shape, but no change in overall extension of the transfer RNA (tRNA) at HP. Configurations extracted from computational ensemble modeling of HP-SAXS profiles were consistent with the NMR results, exhibiting significant disruptions to the acceptor stem, the anticodon stem, and the D-stem regions at HP. We propose that initiation of reverse transcription of HIV RNA could make use of one or more of these excited states.

Road-blocker HSP disease mutation disrupts pre-organization for ATP hydrolysis in kinesin through a second sphere control.

Kinesin motor proteins perform several essential cellular functions powered by the adenosine triphosphate (ATP) hydrolysis reaction. Several single-point mutations in the kinesin motor protein KIF5A have been implicated to hereditary spastic paraplegia disease (HSP), a lethal neurodegenerative disease in humans. In earlier studies, we have shown that a series of HSP-related mutations can impair the kinesin's long-distance displacement or processivity by modulating the order-disorder transition of the linker connecting the heads to the coiled coil. On the other hand, the reduction of kinesin's ATP hydrolysis reaction rate by a distal asparagine-to-serine mutation is also known to cause HSP disease. However, the molecular mechanism of the ATP hydrolysis reaction in kinesin by this distal mutation is still not fully understood. Using classical molecular dynamics simulations combined with quantum mechanics/molecular mechanics calculations, the pre-organization geometry required for optimal hydrolysis in kinesin motor bound to α/ß-tubulin is determined. This optimal geometry has only a single salt-bridge (of the possible two) between Arg203-Glu236, putting a reactive water molecule at a perfect position for hydrolysis. Such geometry is also needed to create the appropriate configuration for proton translocation during ATP hydrolysis. The distal asparagine-to-serine mutation is found to disrupt this optimal geometry. Therefore, the current study along with our previous one demonstrates how two different effects on kinesin dynamics (processivity and ATP hydrolysis), caused by a different set of genotypes, can give rise to the same phenotype leading to HSP disease.

How acidic amino acid residues facilitate DNA target site selection.

Despite the negative charge of the DNA backbone, acidic residues (Asp/Glu) commonly participate in the base readout, with a strong preference for cytosine. In fact, in the solved DNA/protein structures, cytosine is recognized almost exclusively by Asp/Glu through a direct hydrogen bond, while at the same time, adenine, regardless of its amino group, shows no propensity for Asp/Glu. Here, we analyzed the contribution of Asp/Glu to sequence-specific DNA binding using classical and ab initio simulations of selected transcription factors and found that it is governed by a fine balance between the repulsion from backbone phosphates and attractive interactions with cytosine. Specifically, Asp/Glu lower the affinity for noncytosine sites and thus act as negative selectors preventing off-target binding. At cytosine-containing sites, the favorable contribution does not merely rely on the formation of a single H-bond but usually requires the presence of positive potential generated by multiple cytosines, consistently with the observed excess of cytosine in the target sites. Finally, we show that the preference of Asp/Glu for cytosine over adenine is a result of the repulsion from the adenine imidazole ring and a tendency of purine-purine dinucleotides to adopt the BII conformation.

Helical domain of hGBP3 cannot stimulate the second phosphate cleavage of GTP.

Interferon-gamma-inducible large GTPases, hGBPs, possess antipathogenic and antitumor activities in human cells. Like hGBP1, its closest homolog, hGBP3 has two domains; an N-terminal catalytic domain and a C-terminal helical domain, connected by an intermediate region. The biochemical function of this protein and the role of its domains in substrate hydrolysis have not yet been investigated. Here, we report that while hGBP3 can produce both GDP and GMP, GMP is the minor product, 30% (unlike 85% in hGBP1), indicating that hGBP3 is unable to produce enhanced GMP. To understand which domain(s) are responsible for this deficiency, we created hGBP3 truncated variants. Surprisingly, GMP production was similar upon deletion of the helical domain, suggesting that in contrast to hGBP1, the helical domain of hGBP3 cannot stimulate the second phosphate cleavage of GTP. We conducted computational and solution studies to understand the underlying basis. We found that the regulatory residue W79, present in the catalytic domain, forms an H-bond with the backbone carbonyl of K76 (located in the catalytic loop) of the substrate-bound hGBP3. However, after gamma-phosphate cleavage of GTP, the W79-containing region does not undergo a conformational change, failing to redirect the catalytic loop toward the beta-phosphate. This is necessary for efficient GMP formation because hGBP homologs utilize the same catalytic residue for both phosphate cleavages. We suggest that the lack of specific interdomain contacts mediated by the helical domain prevents the catalytic loop movement, resulting in reduced GMP formation. These findings may provide insight into how hGBP3 contributes to immunity.

Bayesian network models identify cooperative GPCR:G protein interactions that contribute to G protein coupling.

Cooperative interactions in protein-protein interfaces demonstrate the interdependency or the linked network-like behavior and their effect on the coupling of proteins. Cooperative interactions also could cause ripple or allosteric effects at a distance in protein-protein interfaces. Although they are critically important in protein-protein interfaces, it is challenging to determine which amino acid pair interactions are cooperative. In this work, we have used Bayesian network modeling, an interpretable machine learning method, combined with molecular dynamics trajectories to identify the residue pairs that show high cooperativity and their allosteric effect in the interface of G protein-coupled receptor (GPCR) complexes with Gα subunits. Our results reveal six GPCR:Gα contacts that are common to the different Gα subtypes and show strong cooperativity in the formation of interface. Both the C terminus helix5 and the core of the G protein are codependent entities and play an important role in GPCR coupling. We show that a promiscuous GPCR coupling to different Gα subtypes, makes all the GPCR:Gα contacts that are specific to each Gα subtype (Gαs, Gαi, and Gαq). This work underscores the potential of data-driven Bayesian network modeling in elucidating the intricate dependencies and selectivity determinants in GPCR:G protein complexes, offering valuable insights into the dynamic nature of these essential cellular signaling components.

Structural and dynamic effects of pseudouridine modifications on noncanonical interactions in RNA.

Pseudouridine is the most frequently naturally occurring RNA modification, found in all classes of biologically functional RNAs. Compared to uridine, pseudouridine contains an additional hydrogen bond donor group and is therefore widely regarded as a structure stabilizing modification. However, the effects of pseudouridine modifications on the structure and dynamics of RNAs have so far only been investigated in a limited number of different structural contexts. Here, we introduced pseudouridine modifications into the U-turn motif and the adjacent U:U closing base pair of the neomycin-sensing riboswitch (NSR)-an extensively characterized model system for RNA structure, ligand binding, and dynamics. We show that the effects of replacing specific uridines with pseudouridines on RNA dynamics crucially depend on the exact location of the replacement site and can range from destabilizing to locally or even globally stabilizing. By using a combination of NMR spectroscopy, MD simulations and QM calculations, we rationalize the observed effects on a structural and dynamical level. Our results will help to better understand and predict the consequences of pseudouridine modifications on the structure and function of biologically important RNAs.

Classification of PTEN missense VUS through exascale simulations.

Phosphatase and tensin homolog (PTEN), a tumor suppressor with dual phosphatase properties, is a key factor in PI3K/AKT signaling pathway. Pathogenic germline variation in PTEN can abrogate its ability to dephosphorylate, causing high cancer risk. Lack of functional evidence lets numerous PTEN variants be classified as variants of uncertain significance (VUS). Utilizing Molecular Dynamics (MD) simulations, we performed a thorough evaluation for 147 PTEN missense VUS, sorting them into 66 deleterious and 81 tolerated variants. Utilizing replica exchange molecular dynamic (REMD) simulations, we further assessed the variants situated in the catalytic core of PTEN's phosphatase domain and uncovered conformational alterations influencing the structural stability of the phosphatase domain. There was a high degree of agreement between our results and the variants classified by Variant Abundance by Massively Parallel Sequencing, saturation mutagenesis, multiplexed functional data and experimental assays. Our extensive analysis of PTEN missense VUS should benefit their clinical applications in PTEN-related cancer. SIGNIFICANCE STATEMENT: Classification of PTEN variants affecting its lipid phosphatase activity is important for understanding the roles of PTEN variation in the pathogenesis of hereditary and sporadic malignancies. Of the 3000 variants identified in PTEN, 1296 (43%) were assigned as VUS. Here, we applied MD and REMD simulations to investigate the effects of PTEN missense VUS on the structural integrity of the PTEN phosphatase domain consisting the WPD, P and TI active sites. We classified a total of 147 missense VUS into 66 deleterious and 81 tolerated variants by referring to the control group comprising 54 pathogenic and 12 benign variants. The classification was largely in concordance with these classified by experimental approaches.

Long-range electron proton coupling in respiratory complex I - insights from molecular simulations of the quinone chamber and antiporter-like subunits.

Respiratory complex I is a redox-driven proton pump. Several high-resolution structures of complex I have been determined providing important information about the putative proton transfer paths and conformational transitions that may occur during catalysis. However, how redox energy is coupled to the pumping of protons remains unclear. In this article, we review biochemical, structural and molecular simulation data on complex I and discuss several coupling models, including the key unresolved mechanistic questions. Focusing both on the quinone-reductase domain as well as the proton-pumping membrane-bound domain of complex I, we discuss a molecular mechanism of proton pumping that satisfies most experimental and theoretical constraints. We suggest that protonation reactions play an important role not only in catalysis, but also in the physiologically-relevant active/deactive transition of complex I.

Leveraging intrinsic flexibility to engineer enhanced enzyme catalytic activity.

Dynamic motions of enzymes occurring on a broad range of timescales play a pivotal role in all steps of the reaction pathway, including substrate binding, catalysis, and product release. However, it is unknown whether structural information related to conformational flexibility can be exploited for the directed evolution of enzymes with higher catalytic activity. Here, we show that mutagenesis of residues exclusively located at flexible regions distal to the active site of Homo sapiens kynureninase (HsKYNase) resulted in the isolation of a variant (BF-HsKYNase) in which the rate of the chemical step toward kynurenine was increased by 45-fold. Mechanistic pre­steady-state kinetic analysis of the wild type and the evolved enzyme shed light on the underlying effects of distal mutations (>10 Å from the active site) on the rate-limiting step of the catalytic cycle. Hydrogen-deuterium exchange coupled to mass spectrometry and molecular dynamics simulations revealed that the amino acid substitutions in BF-HsKYNase allosterically affect the flexibility of the pyridoxal-5'-phosphate (PLP) binding pocket, thereby impacting the rate of chemistry, presumably by altering the conformational ensemble and sampling states more favorable to the catalyzed reaction.

Topological gelation of reconnecting polymers.

DNA recombination is a ubiquitous process that ensures genetic diversity. Contrary to textbook pictures, DNA recombination, as well as generic DNA translocations, occurs in a confined and highly entangled environment. Inspired by this observation, here, we investigate a solution of semiflexible polymer rings undergoing generic cutting and reconnection operations under spherical confinement. Our setup may be realized using engineered DNA in the presence of recombinase proteins or by considering micelle-like components able to form living (or reversibly breakable) polymer rings. We find that in such systems, there is a topological gelation transition, which can be triggered by increasing either the stiffness or the concentration of the rings. Flexible or dilute polymers break into an ensemble of short, unlinked, and segregated rings, whereas sufficiently stiff or dense polymers self-assemble into a network of long, linked, and mixed loops, many of which are knotted. We predict that the two phases should behave qualitatively differently in elution experiments monitoring the escape dynamics from a permeabilized container. Besides shedding some light on the biophysics and topology of genomes undergoing DNA reconnection in vivo, our findings could be leveraged in vitro to design polymeric complex fluids-e.g., DNA-based complex fluids or living polymer networks-with desired topologies.

Uneven Strain Distribution Induces Consecutive Dislocation Slipping, Plane Gliding, and Subsequent Detwinning of Penta-Twinned Nanoparticles.

Twin structures possess distinct physical and chemical properties by virtue of their specific twin configuration. However, twinning and detwinning processes are not fully understood on the atomic scale. Integrating in situ high resolution transmission electron microscopy and molecular dynamic simulations, we find tensile strain in the asymmetrical 5-fold twins of Au nanoparticles leads to twin boundary migration through dislocation sliding (slipping of an atomic layer) along twin boundaries and dislocation reactions at the 5-fold axis under an electron beam. Migration of one or two layers of twin planes is governed by energy barriers, but overall, the total energy, including surface, lattice strain, and twin boundary energy, is relaxed after consecutive twin boundary migration, leading to a detwinning process. In addition, surface rearrangement of 5-fold twinned nanoparticles can aid in the detwinning process.

Insights Into the Molecular Interactions of MIC2 and M2AP: Role of TSR6 and Conservation Across Species.

Microneme protein 2 (MIC2) and its associated protein M2AP are pivotal for the gliding motility and host cell invasion by Toxoplasma gondii. In our prior work, we showed that M2AP binds specifically to the sixth TSR domain of MIC2, with this interaction mediated dominantly by the hotspot residue H620 situated at the center of TSR6. To delve deeper into the functional significance of H620 and explore the dynamic behavior of Y602, we conducted molecular dynamic (MD) simulations of the Toxoplasma TSR6-M2AP complex, encompassing both wild-type and mutant forms. Our findings underscore the critical role of H620 within TSR6, particularly its hydrogen bond interaction with K72 of M2AP. The H620A mutation disrupts the nearby hydrophobic network while minimally affecting other hydrophilic interactions. Furthermore, our data reveal a highly conserved binding pose between M2AP and TSR6 across different species, consistent with previous trans-genera studies, thereby offering insights for future strategies in infection control development.

Identification of potent BRD4-BD1 inhibitors using classical and steered molecular dynamics based free energy analysis.

In the present work a combination of traditional and steered molecular dynamics based techniques were employed to identify potential inhibitors against the human BRD4 protein (BRD4- BD1); an established drug target for multiple illnesses including various malignancies. Quinoline derivatives that were synthesized in-house were tested for their potential as new BRD4-BD1 inhibitors. Initially molecular docking experiments were performed to determine the binding poses of BRD4-BD1 inhibitors. To learn more about the thermodynamics of inhibitor binding to the BRD4-BD1 active site, the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) free energy calculations were conducted afterwards. The findings of the MM-PBSA analysis were further reinforced by performing steered umbrella sampling simulations which revealed crucial details about the binding/unbinding process of the most potent quinoline derivatives at the BRD4-BD1 active site. We report a novel quinoline derivative which can be developed into a fully functional BRD4-BD1 inhibitor after experimental validation. The identified compound (4 g) shows better properties than the standard BRD4-BD1 inhibitors considered in the study. The study also highlights the crucial role of Gln78, Phe79, Trp81, Pro82, Phe83, Gln84, Gln85, Val87, Leu92, Leu94, Tyr97, Met105, Cys136, Asn140, Ile146 and Met149 in inhibitor binding. The study provides a possible lead candidate and key amino acids involved in inhibitor recognition and binding at the active site of BRD4-BD1 protein. The findings might be of significance to medicinal chemists involved in the development of potent BRD4-BD1 inhibitors.

Nanoparticle-mediated cancer cell therapy: basic science to clinical applications.

Every sixth person in the world dies due to cancer, making it the second leading severe cause of death after cardiovascular diseases. According to WHO, cancer claimed nearly 10 million deaths in 2020. The most common types of cancers reported have been breast (lung, colon and rectum, prostate cases), skin (non-melanoma) and stomach. In addition to surgery, the most widely used traditional types of anti-cancer treatment are radio- and chemotherapy. However, these do not distinguish between normal and malignant cells. Additional treatment methods have evolved over time for early detection and targeted therapy of cancer. However, each method has its limitations and the associated treatment costs are quite high with adverse effects on the quality of life of patients. Use of individual atoms or a cluster of atoms (nanoparticles) can cause a paradigm shift by virtue of providing point of sight sensing and diagnosis of cancer. Nanoparticles (1-100 nm in size) are 1000 times smaller in size than the human cell and endowed with safer relocation capability to attack mechanically and chemically at a precise location which is one avenue that can be used to destroy cancer cells precisely. This review summarises the extant understanding and the work done in this area to pave the way for physicians to accelerate the use of hybrid mode of treatments by leveraging the use of various nanoparticles.

Structural and Dynamic Features of the Recognition of 8-oxoguanosine Paired with an 8-oxoG-clamp by Human 8-oxoguanine-DNA Glycosylase.

8-oxoguanine (oxoG) is formed in DNA by the action of reactive oxygen species. As a highly mutagenic and the most common oxidative DNA lesion, it is an important marker of oxidative stress. Human 8-oxoguanine-DNA glycosylase (OGG1) is responsible for its prompt removal in human cells. OGG1 is a bifunctional DNA glycosylase with N-glycosylase and AP lyase activities. Aspects of the detailed mechanism underlying the recognition of 8-oxoguanine among numerous intact bases and its subsequent interaction with the enzyme's active site amino acid residues are still debated. The main objective of our work was to determine the effect (structural and thermodynamic) of introducing an oxoG-clamp in model DNA substrates on the process of 8-oxoG excision by OGG1. Towards that end, we used DNA duplexes modeling OGG1-specific lesions: 8-oxoguanine or an apurinic/apyrimidinic site with either cytidine or the oxoG-clamp in the complementary strand opposite to the lesion. It was revealed that there was neither hydrolysis of the N-glycosidic bond at oxoG nor cleavage of the sugar-phosphate backbone during the reaction between OGG1 and oxoG-clamp-containing duplexes. Possible structural reasons for the absence of OGG1 enzymatic activity were studied via the stopped-flow kinetic approach and molecular dynamics simulations. The base opposite the damage was found to have a critical effect on the formation of the enzyme-substrate complex and the initiation of DNA cleavage. The oxoG-clamp residue prevented the eversion of the oxoG base into the OGG1 active site pocket and impeded the correct convergence of the apurinic/apyrimidinic site of DNA and the attacking nucleophilic group of the enzyme. An obtained three-dimensional model of the OGG1 complex with DNA containing the oxoG-clamp, together with kinetic data, allowed us to clarify the role of the contact of amino acid residues with DNA in the formation of (and rearrangements in) the enzyme-substrate complex.

Crystal structure of a thiolase from Archaeal Pyrococcus furiosus and its in silico functional annotation.

In most of the eukaryotes and archaea, isopentenyl pyrophosphate (IPP) and dimethyl allyl pyrophosphate (DMAPP) essential building blocks of all isoprenoids synthesized in the mevalonate pathway. Here, the first enzyme of this pathway, acetoacetyl CoA thiolase (PFC_04095) from an archaea Pyrococcus furiosus is structurally characterized. The crystal structure of PFC_04095 is determined at 2.7 Å resolution, and the crystal structure reveals the absence of catalytic acid/base cysteine in its active site, which is uncommon in thiolases. In place of cysteine, His285 of HDAF motif performs both protonation and abstraction of proton during the reaction. The crystal structure shows that the distance between Cys83 and His335 is 5.4 Å. So, His335 could not abstract a proton from nucleophilic cysteine (Cys83), resulting in the loss of enzymatic activity of PFC_04095. MD simulations of the docked PFC_04095-acetyl CoA complex show substrate binding instability to the active site pocket. Here, we have reported that the stable binding of acetyl CoA to the PFC_04095 pocket requires the involvement of three protein complexes, i.e., thiolase (PFC_04095), DUF35 (PFC_04100), and HMGCS (PFC_04090).