Search

The presence of toxic dyes in industrial waste dramatically diminishes the beneficial effects of remediation efforts. To overcome the hazardous impacts of dyes on biodiversity and environment, we integrated polymers into nanoparticles to substantially enhance their functionality and performance. 2 and 4 wt% of chitosan (CS) and 3 wt% of polyacrylic acid (PAA) doped cadmium selenide (CdSe) nanostructures (NSs) were prepared by co-precipitation approach. CdSe quantum dots (QDs) exhibit a narrow band gap energy, high solubility, and tunable properties, which are appropriate for redox reactions but show less adsorption and catalytic behavior. In this work, catalytic and antibacterial activities of CdSe QDs enhanced upon the integration of PAA due to increment in surface area confirmed by BET analysis. Furthermore, the addition of CS escalates the dye degradation and microbes evolve to the interaction of CdSe surface with the functional groups of CS. Highly doped CdSe shows significant inhibitory zones (8.65 to 9.30 mm) against gram-positive bacteria Staphylococcus aureus (S. aureus). In addition, the inhibitory activity of CS/PAA-CdSe nanostructures against DNA gyrase and dihydrofolate reductase (DHFR) in S. aureus was elucidated using molecular docking investigations, providing a rationale for their bactericidal action., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier B.V. All rights reserved.)

Integrating machine learning potentials (MLPs) with quantum mechanical/molecular mechanical (QM/MM) free energy simulations has emerged as a powerful approach for studying enzymatic catalysis. However, its practical application has been hindered by the time-consuming process of generating the necessary training, validation, and test data for MLP models through QM/MM simulations. Furthermore, the entire process needs to be repeated for each specific enzyme system and reaction. To overcome this bottleneck, it is required that trained MLPs exhibit transferability across different enzyme environments and reacting species, thereby eliminating the need for retraining with each new enzyme variant. In this study, we explore this potential by evaluating the transferability of a pretrained ΔMLP model across different enzyme mutations within the MM environment using the QM/MM-based ML architecture developed by Pan, X. J. Chem. Theory Comput. 2021, 17(9), 5745-5758. The study includes scenarios such as single point substitutions, a homologous enzyme from different species, and even a transition to an aqueous environment, where the last two systems have MM environment that is substantially different from that used in MLP training. The results show that the ΔMLP model effectively captures and predicts the effects of enzyme mutations on electrostatic interactions, producing reliable free energy profiles of enzyme-catalyzed reactions without the need for retraining. The study also identified notable limitations in transferability, particularly when transitioning from enzyme to water-rich MM environments. Overall, this study demonstrates the robustness of the Pan et al.'s QM/MM-based ML architecture for application to diverse enzyme systems, as well as the need for further research and the development of more sophisticated MLP models and training methods.

Identifying target proteins that interact with bioactive molecules is indispensable for understanding their mechanisms of action. In this study, we developed a uniform ribosome display technology using equal-length DNAs and mRNAs to improve molecular display principle for target identification. The equal-length DNAs were designed to contain various coding sequences for full-length proteins with molecular weights of up to 130 kDa and were used to synthesize equal-length mRNAs, which allowed the formation of full-length protein-ribosome-equal-length mRNA complexes. Uniform ribosome display selections of dihydrofolate reductase and haloalkane dehalogenase mutant were performed against methotrexate and chlorohexane, respectively. Quantitative changes of proteins after each selection indicated that the target protein-displaying ribosomal complexes were specifically selected through non-covalent or covalent interactions with the corresponding bioactive molecules. Furthermore, selection of full-length proteins interacting with methotrexate or anti-DDX46 antibody from protein pools showed that only the target proteins could be precisely identified even though the molar amounts of equal-length mRNAs encoding them were adjusted to 1/20,000 of the total equal-length mRNAs. Thus, the uniform ribosome display technology enabled efficient identification of target proteins that interact with bioactive small and large molecules through simplified operations without deep sequencing., (© 2024 Wiley-VCH GmbH.)

Naga chilli (Capsicum chinense Jacq.) have garnered significant attention due to the plant's possible health benefits and variety of phytochemical components. Utilizing cutting-edge analytical techniques such as gas chromatography-mass spectrometry (GC-MS) and high-performance thin layer chromatography (HPTLC) in conjunction with bioautography, this study conducts a thorough phytochemical profiling and biological activity assessment of the Naga chilli plant. An in silico docking study was performed for all the bioactive compounds identified through GC-MS against dihydrofolate reductase, a critical enzyme for bacterial survival. Many important components were identified and quantified with the help of subsequent GC-MS and HPTLC analysis. Among them, capsaicinoids were found to be the most prevalent. GC-MS results showed nonadecane (21.28%), 1-dimethyl(phenyl)silyloxypentane (14.53%), capsaicin (13.55%) and 2-pentanone, 4-hydroxy-4-methyl- (11.42%) were the most prevalent. HPTLC report showed capsaicin was 0.833 mg/g of fresh weight of Naga chilli. This study showed good docking scores for some of the constituents, particularly capsaicin, indicating that this plant is a good candidate for antimicrobial activity. This activity of the extract confirms the docking results, which needs to be in focus for further antimicrobial drug development., (© 2024 John Wiley & Sons Ltd.)

Recently, P218, a new flexible antifolate targeting Plasmodium falciparum dihydrofolate reductase (PfDHFR), has entered its clinical trial with good safety profile and effective Pf infection prevention. However, it carries a free carboxyl terminal, which is hydrophilic and prone to metabolic glucuronidation. Here, a new series of P218 analogues carrying butyrolactone has been synthesized with the purpose of enhancing lipophilicity and minimizing metabolic instability. The inhibition constants against the mutant PfDHFR enzymes are in sub-nanomolar level and the antimalarial activity against antifolate-resistant parasites are in the low micromolar range. The crystal structure of the most potent analogue LA1 bound enzyme complex indicates interaction with multiple residues, including Arg122 and Phe116 in the active site. In vitro log D 7.4 and kinetic solubility confirmed a higher lipophilicity of this butyrolactone series as compared to P218. These outcomes suggest the possibility to further develop butyrolactone derivatives as non-carboxyl antiplasmodial antifolates., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Sasithorn Decharuangsilp has patent #TH2301005512 pending to National Science and Technology Development Agency (NSTDA). If there are other authors, they declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Tuberculosis (TB) is an obstinate and infectious disease requiring a relatively longer treatment duration than other bacterial infections. The current treatment regime is prolonged and cumbersome, with adverse effects, often leading to nonadherence. The upsurge in TB's multidrug-resistant and extensively drug-resistant strains with evolved resistance to existing drugs has compounded the problems. The last two decades witnessed unprecedented progress in developing TB drugs with better efficacy and reduced toxicity. Of late, inhibitors targeting the dihydrofolate reductase (DHFR) enzyme were being explored and developed as antitubercular drugs. A plethora of diverse molecular cores, such as pteridines, diamino heterocycles, diamino triazoles, and nontraditional cores, were developed recently as Mtb-DHFR targets. Besides the characteristic binding pockets of Mtb-DHFR, an extended hydrophilic binding pocket was also studied for intermolecular interactions with the designed compounds to assess the enzyme specificity. In this study, prominent DHFR inhibitors developed in the last two decades were reported. Key features of the designed compounds, such as the structural similarities with existing pharmacophores, interactions with binding pockets, enzyme selectivity and specificity, and percentage of inhibition, were evaluated. The authors hope the study will help streamline the pharmacological pipeline of Mtb-DHFR inhibitors and bring the investigators one step closer to success., (© 2024 John Wiley & Sons Ltd.)

Protein folding in vivo begins during synthesis on the ribosome and is modulated by molecular chaperones that engage the nascent polypeptide. How these features of protein biogenesis influence the maturation pathway of nascent proteins is incompletely understood. Here, we use hydrogen-deuterium exchange mass spectrometry to define, at peptide resolution, the cotranslational chaperone-assisted folding pathway of Escherichia coli dihydrofolate reductase. The nascent polypeptide folds along an unanticipated pathway through structured intermediates not populated during refolding from denaturant. Association with the ribosome allows these intermediates to form, as otherwise destabilizing carboxy-terminal sequences remain confined in the ribosome exit tunnel. Trigger factor binds partially folded states without disrupting their structure, and the nascent chain is poised to complete folding immediately upon emergence of the C terminus from the exit tunnel. By mapping interactions between the nascent chain and ribosomal proteins, we trace the path of the emerging polypeptide during synthesis. Our work reveals new mechanisms by which cellular factors shape the conformational search for the native state., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Antimicrobial peptides have garnered increasing attention as potential alternatives due to their broad-spectrum antimicrobial activity and low propensity for developing resistance. This is for the first time; proteome sequences of Aegle marmelos were subjected to in-silico digestion and AMP prediction were performed using DBAASP server. After screening the peptides on the basis of different physiochemical property, peptide sequence GKEAATKAIKEWGQPKSKITH (AM1) shows the maximum binding affinity with - 10.2 Kcal/mol in comparison with the standard drug (Trimethoprim) with - 7.4 kcal/mol and - 6.8 Kcal/mol for DHFR and SaTrmK enzyme respectively. Molecular dynamics simulation performed for 300ns, it has been found that peptide was able to stabilize the protein more effectively, analysed by RMSD, RMSF, and other statistical analysis. Free binding energy for DHFR and SaTrmK interaction from MMPBSA analysis with peptide was found to be -47.69 and - 44.32 Kcal/mol and for Trimethoprim to be -13.85 Kcal/mol and - 11.67 Kcal/mol respectively. Further in-vitro study was performed against Methicillin Susceptible Staphylococcus aureus (MSSA), Methicillin Resistant Staphylococcus aureus (MRSA), Multi-Drug Resistant Staphylococcus aureus (MDR-SA) strain, where MIC values found to be 2, 4, and 8.5 µg/ml lesser in comparison to trimethoprim which has higher MIC values 2.5, 5, and 9.5 µg/ml respectively. Thus, our study provides the insight for the further in-vivo study of the peptides against multi-drug resistant S. aureus., (© 2024. The Author(s).)

Dihydrofolate reductase (DHFR), due to its universality and the depth with which it has been studied, is a model system in the study of protein dynamics. Myriad previous works have identified networks of residues in positions near to and remote from the active site that are involved in the dynamics. For example, specific mutations on the Met20 loop in Escherichia coli DHFR (N23PP/S148A) are known to disrupt millisecond-time scale motions as well as reduce catalytic activity. However, how and if networks of dynamically coupled residues influence the evolution of DHFR is still an unanswered question. In this study, we first identify, by statistical coupling analysis and molecular dynamic simulations, a network of coevolving residues that possesses increased correlated motions. We then go on to show that allosteric communication in this network is knocked down in N23PP/S148A mutant E. coli DHFR. We also identify two sites in the human DHFR sector which may accommodate the Met20 loop double proline motif. Finally, we demonstrate a concerted evolutionary change in the human DHFR allosteric networks, which maintains dynamic communication. These findings strongly implicate protein dynamics as a driving force for evolution.

The sensitivity to protein inhibitors is altered by modifications or protein mutations, as represented by drug resistance. The mode of stable drug binding to the protein pocket has been experimentally clarified. However, the nature of the binding of inhibitors with reduced sensitivity remains unclear at the atomic level. In this study, we analyzed the thermodynamics and kinetics of inhibitor binding to the surface of wild-type and mutant dihydrofolate reductase (DHFR) using molecular dynamics simulations combined with Markov state modeling. A strong inhibitor of methotrexate (MTX) showed a preference for the active site of wild-type DHFR with minimal binding to unrelated (secondary) sites. Deletion of a side-chain fragment in MTX largely destabilized the active site-bound state, with clear evidence of binding to secondary sites. Similarly, the F31V mutation in DHFR diminished the specificity of MTX binding to the active site. These results reveal the presence of multiple-bound states whose stabilities are comparable to or higher than those of the unbound state, suggesting that a reduction in the binding affinity for the active site significantly elevates the fractions of these states. This study presents a theoretical model that more accurately interprets the altered drug sensitivity than the traditional two-state model.

Recent work from Nguyen et al. unveils massively parallel measurements of epistatic interactions between two enzymes, dihydrofolate reductase and thymidylate synthase, in their natural cellular context. Almost 3000 mutations of DHFR in three TYMS backgrounds reveal a complex interaction network. The authors capture much of this complexity using a simple model., Competing Interests: Declaration of interests The author declares no competing interests., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Dihydrofolate reductase (DHFR) is ubiquitously present in all living organisms and plays a crucial role in the growth of the fungal pathogen R.solani. Sequence alignment confirmed the evolutionary conservation of the essential lid domain, with the amino acid 'P' within the PEKN lid domain appearing with a frequency of 89.5% in higher organisms and 11.8% in lower organisms. Consequently, a K65P variant was introduced into R.solani DHFR (rDHFR). Subsequent enzymatic kinetics assays were conducted for human DHFR (hDHFR), rDHFR, E. coli DHFR (eDHFR), and the K65P variant. hDHFR exhibited the highest k cat of 0.95 s -1 , followed by rDHFR with 0.14 s -1 , while eDHFR displayed the lowest k cat of 0.09 s -1 . Remarkably, the K65P variant induced a significant reduction in K m , resulting in a 1.8-fold enhancement in catalytic efficiency (k cat /K m ) relative to the wild type. Differential scanning fluorimetry and binding free energy calculations confirmed the enhanced substrate affinity for both folate and NADPH in the K65P variant. These results suggest that the K65P mutation enhances substrate affinity and catalytic efficiency in DHFR, highlighting the evolutionary and functional importance of the K65 residue., (© 2024. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

Dihydrofolate reductase (DHFR), an essential enzyme in folate metabolism, presents a promising target for drug development against various diseases, including cancer and tuberculosis. Herein, we present an integrated approach combining in vitro biochemical assays with in silico molecular docking analysis to evaluate the inhibitory potential of 4-piperidine-based thiosemicarbazones 5(a-s) against DHFR. In our in vitro study, a novel series of 4-piperidine-based thiosemicarbazones 5(a-s) were assessed for their inhibitory activity against DHFR enzyme. The synthesized compounds 5(a-s) exhibited potent inhibition with IC 50 values in the range of 13.70 ± 0.25 µM to 47.30 ± 0.86 µM. Among all the derivatives 5p displayed highest inhibitory activity. Simultaneously, in silico analysis were performed and compared with standard drug (Methotrexate) to predict the binding affinity and interaction pattern of synthesized compounds with DHFR active site. SAR analysis was done to elucidate how structural modifications impact compound's biological activity, guiding the rational design of potent and selective drug candidates for targeted diseases. These findings may provide a comprehensive assessment of 4-piperdine-based thiosemicarbazones as DHFR inhibitors and contribute to the development of novel therapeutics targeting DHFR-associated diseases., (© 2024. The Author(s).)

In this study, the structural and anticancer properties of aminopterin, as well as its antiviral characteristics, were elucidated. The preferred conformations of the title molecule were investigated with semiempirical AM1 method, and the obtained the lowest energy conformer was then optimized by using density functional (DFT/B3LYP) method with 6-311++G(d,p) as basis set. The vibrational frequencies of the optimized structure were calculated by the same level of theory and were compared with the experimental values. The vibrational assignments were performed based on the computed potential energy distribution (PED) of the vibrational modes. The molecular electrostatic potential (MEP) and frontier molecular orbitals (HOMO, LUMO) analyses were carried out for the optimized structure and the chemical reactivity has been scrutinized. To enlighten the biological activity of aminopterin as anticancer and anti-COVID-19 agents, aminopterin was docked into DNA, α IIB β 3 and α 5 β 1 integrins, human dihydrofolate reductase, main protease (M pro ) of SARS-CoV-2 and SARS-CoV-2/ACE2 complex receptor. The binding mechanisms of aminopterin with the receptors were clarified. The molecular docking results revealed the strong interaction of the aminopterin with DNA (-8.2 kcal/mol), α IIB β 3 and α 5 β 1 integrins (-9.0 and -10.8 kcal/mol, respectively), human dihydrofolate reductase (-9.7 kcal/mol), M pro of SARS-CoV-2 (-6.7 kcal/mol), and SARS-CoV-2/ACE2 complex receptor (-8.1 kcal/mol). Moreover, after molecular docking calculations, top-scoring ligand-receptor complexes of the aminopterin with SARS-CoV-2 enzymes (6M03 and 6M0J) were subjected to 50 ns all-atom MD simulations to investigate the ligand-receptor interactions in more detail, and to determine the binding free energies accurately. The predicted results indicate that the aminopterin may significantly inhibit SARS-CoV-2 infection. Thus, in this study, as both anticancer and anti-COVID-19 agents, the versatility of the biological activity of aminopterin was shown.Communicated by Ramaswamy H. Sarma.

Folate enzymes, namely, dihydrofolate reductase (DHFR) and pteridine reductase (PTR1) are acknowledged targets for the development of antiparasitic agents against Trypanosomiasis and Leishmaniasis. Based on the amino dihydrotriazine motif of the drug Cycloguanil (Cyc), a known inhibitor of both folate enzymes, we have identified two novel series of inhibitors, the 2-amino triazino benzimidazoles ( 1 ) and 2-guanidino benzimidazoles ( 2 ), as their open ring analogues. Enzymatic screening was carried out against PTR1, DHFR, and thymidylate synthase (TS). The crystal structures of Tb DHFR and Tb PTR1 in complex with selected compounds experienced in both cases a substrate-like binding mode and allowed the rationalization of the main chemical features supporting the inhibitor ability to target folate enzymes. Biological evaluation of both series was performed against T. brucei and L. infantum and the toxicity against THP-1 human macrophages. Notably, the 5,6-dimethyl-2-guanidinobenzimidazole 2g resulted to be the most potent ( K i = 9 nM) and highly selective Tb DHFR inhibitor, 6000-fold over Tb PTR1 and 394-fold over h DHFR. The 5,6-dimethyl tricyclic analogue 1g , despite showing a lower potency and selectivity profile than 2g , shared a comparable antiparasitic activity against T. brucei in the low micromolar domain. The dichloro-substituted 2-guanidino benzimidazoles 2c and 2d revealed their potent and broad-spectrum antitrypanosomatid activity affecting the growth of T. brucei and L. infantum parasites. Therefore, both chemotypes could represent promising templates that could be valorized for further drug development.

The Escherichia coli GroEL/ES chaperonin system facilitates protein folding in an ATP-driven manner. There are <100 obligate clients of this system in E. coli although GroEL can interact and assist the folding of a multitude of proteins in vitro. It has remained unclear, however, which features distinguish obligate clients from all the other proteins in an E. coli cell. To address this question, we established a system for selecting mutations in mouse dihydrofolate reductase (mDHFR), a GroEL interactor, that diminish its dependence on GroEL for folding. Strikingly, both synonymous and non-synonymous codon substitutions were found to reduce mDHFR's dependence on GroEL. The non-synonymous substitutions increase the rate of spontaneous folding whereas computational analysis indicates that the synonymous substitutions appear to affect translation rates at specific sites., (© 2024 The Author(s). Protein Science published by Wiley Periodicals LLC on behalf of The Protein Society.)

In order to develop novel antimicrobial agents, we prepared quinoline bearing pyrimidine analogues 2-7, 8 a-d and 9 a-d and their structures were elucidated by spectroscopic techniques. Furthermore, our second aim was to predict the interactions between the active compounds and enzymes (DNA gyrase and DHFR). In this work, fourteen pyrimido[4,5-b]quinoline derivatives were prepared and assessed for their antimicrobial potential by estimating zone of inhibition. All the screened candidates displayed antibacterial potential with zone of inhibition range of 9-24 mm compared with ampicillin (20-25 mm) as a reference drug. Moreover, the target derivatives 2 (ZI=16), 9 c (ZI=17 mm) and 9 d (ZI=16 mm) recorded higher antifungal activity against C. albicans to that exhibited by the antifungal drug amphotericin B (ZI=15 mm). Finally, the most potent pyrimidoquinoline compounds (2, 3, 8 c, 8 d, 9 c and 9 d) were docked inside DHFR and DNA gyrase active sites and they recorded excellent fitting within the active regions of DNA gyrase and DHFR. These outcomes revealed us that compounds (2, 3, 8 c, 8 d, 9 c and 9 d) could be lead compounds to discover novel antibacterial candidates., (© 2024 Wiley-VHCA AG, Zurich, Switzerland.)

Secondary amines, due to their reactivity, can transform protein templates into catalytically active entities, accelerating the development of artificial enzymes. However, existing methods, predominantly reliant on modified ligands or N-terminal prolines, impose significant limitations on template selection. In this study, genetic code expansion was used to break this boundary, enabling secondary amines to be incorporated into alternative proteins and positions of choice. Pyrrolysine analogues carrying different secondary amines could be incorporated into superfolder green fluorescent protein (sfGFP), multidrug-binding LmrR and nucleotide-binding dihydrofolate reductase (DHFR). Notably, the analogue containing a D-proline moiety demonstrated both proteolytic stability and catalytic activity, conferring LmrR and DHFR with the desired transfer hydrogenation activity. While the LmrR variants were confined to the biomimetic 1-benzyl-1,4-dihydronicotinamide (BNAH) as the hydride source, the optimal DHFR variant favorably used the pro-R hydride from NADPH for stereoselective reactions (e.r. up to 92 : 8), highlighting that a switch of protein template could broaden the nucleophile option for catalysis. Owing to the cofactor compatibility, the DHFR-based secondary amine catalysis could be integrated into an enzymatic recycling scheme. This established method shows substantial potential in enzyme design, applicable from studies on enzyme evolution to the development of new biocatalysts., (© 2024 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

In this study, new series of N'-(2-(substitutedphenoxy)acetyl)-4-(1H-pyrrol-1-yl)benzohydrazides (3a-j) 4-(2,5-dimethyl-1H-pyrrol-1-yl)-N'-(2-(substitutedphenoxy)acetyl)benzohydrazides (5a-j) were synthesized, characterized and assessed as inhibitors of enoyl ACP reductase and DHFR. Most of the compounds exhibited dual inhibition against the enzymes enoyl ACP reductase and DHFR. Several synthesized substances also demonstrated significant antibacterial and antitubercular properties. A molecular docking analysis was conducted in order to determine the potential mechanism of action of the synthesized compounds. The results indicated that there were binding interactions seen with the active sites of dihydrofolate reductase and enoyl ACP reductase. Additionally, important structural details were identified that play a critical role in sustaining the dual inhibitory activity. These findings were useful for the development of future dual inhibitors. Therefore, this study provided strong evidence that several synthesized molecules could exert their antitubercular properties at the cellular level through multi-target inhibition. By shedding light on the mechanisms through which these compounds exert their inhibitory effects, this research opens up promising avenues for the future development of dual inhibitors with enhanced antibacterial and antitubercular properties. The study's findings underscore the importance of multi-target approaches in drug design, providing a strong foundation for the design and optimization of novel compounds that can effectively target bacterial infections at the cellular level., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Mahnashi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The observation of multiple conformations of a functional loop (termed M20) in the Escherichia coli dihydrofolate reductase ( ec DHFR) enzyme triggered the proposition that large-scale motions of protein structural elements contribute to enzyme catalysis. The transition of the M20 loop from a closed conformation to an occluded conformation was thought to aid the rate-limiting release of the products. However, the influence of charged species in the solution environment on the observed M20 loop conformations, independent of charged ligands bound to the enzyme, had not been considered. Molecular dynamics simulations of ec DHFR in model CaCl 2 solutions of varying molar ionic strengths I M reveal a substantial free energy barrier between occluded and closed M20 loop states at I M exceeding the E. coli threshold (∼0.24 M). This barrier may facilitate crystallization of ec DHFR in the occluded state, consistent with ec DHFR structures obtained at I M exceeding 0.3 M. At lower I M (≤0.15 M), the M20 loop can explore the occluded state, but prefers an open / partially closed conformation, again consistent with ec DHFR structures. Our findings caution against using ec DHFR structures obtained at nonphysiological ionic strengths in interpreting catalytic events or in structure-based drug design.

Pneumocystis jirovecii is a fungal pathogen that causes pneumocystis pneumonia, a disease that mainly affects immunocompromised individuals. This fungus has historically been hard to study because of our inability to grow it in vitro. One of the main drug targets in P. jirovecii is its dihydrofolate reductase (PjDHFR). Here, by using functional complementation of the baker's yeast ortholog, we show that PjDHFR can be inhibited by the antifolate methotrexate in a dose-dependent manner. Using deep mutational scanning of PjDHFR, we identify mutations conferring resistance to methotrexate. Thirty-one sites spanning the protein have at least one mutation that leads to resistance, for a total of 355 high-confidence resistance mutations. Most resistance-inducing mutations are found inside the active site, and many are structurally equivalent to mutations known to lead to resistance to different antifolates in other organisms. Some sites show specific resistance mutations, where only a single substitution confers resistance, whereas others are more permissive, as several substitutions at these sites confer resistance. Surprisingly, one of the permissive sites (F199) is without direct contact to either ligand or cofactor, suggesting that it acts through an allosteric mechanism. Modeling changes in binding energy between F199 mutants and drug shows that most mutations destabilize interactions between the protein and the drug. This evidence points towards a more important role of this position in resistance than previously estimated and highlights potential unknown allosteric mechanisms of resistance to antifolate in DHFRs. Our results offer unprecedented resources for the interpretation of mutation effects in the main drug target of an uncultivable fungal pathogen., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Rouleau et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Dihydrofolate reductase (DHFR) is a ubiquitous enzyme that regulates the biosynthesis of tetrahydrofolate among various species of Plasmodium parasite. It is a validated target of the antifolate drug pyrimethamine (Pyr) in Plasmodium falciparum (Pf), but its clinical efficacy has been hampered due to the emergence of drug resistance. This has made the attempt to screen Food & Drug Administration-approved drugs against wild- and mutant PfDHFR by employing an in-silico pipeline to identify potent candidates. The current study has followed a virtual screening approach for identifying potential DHFR inhibitors from DrugBank database, based on a structure similarity search of candidates, followed by absorption, distribution, metabolism, and excretion estimation. The screened drugs were subjected to various parameters like docking, molecular mechanics with generalized born and surface area solvation calculations, and molecular simulations. We have thus identified two potential drug candidates, duloxetine and guanethidine, which can be repurposed to be tested for their efficacy against wild type and drug resistant falciparum malaria., (© 2024 Wiley Periodicals LLC.)

Many natural products have been shown to possess antiplasmodial activities, but their protein targets are unknown. This work employed molecular docking and molecular dynamics simulations to explore the inhibitory activity of some antiplasmodial natural products against wild-type and mutant strains of Plasmodium falciparum dihydrofolate reductase ( Pf DHFR). From the molecular docking study, 6 ligands preferentially bind at the active site of the DHFR domain with binding energies ranging from -6.4 to -9.5 kcal/mol. Interactions of compounds with MET55 and PHE58 were mostly observed in the molecular docking study. From the molecular dynamics study, the binding of 2 of the ligands-nitidine and oplodiol-was observed to be stable against all tested strains of Pf DHFR. The average binding free energy of oplodiol in complex with the various Pf DHFR strains was -93.701 kJ/mol whereas that of nitidine was -106.206 kJ/mol. The impressive in silico activities of the 2 compounds suggest they could be considered for development as potential antifolate agents.Communicated by Ramaswamy H. Sarma.

Pharmacological inhibition of dihydrofolate reductase (DHFR) is an established approach for treating a variety of human diseases, including foreign infections and cancer. However, treatment with classic DHFR inhibitors, such as methotrexate (MTX), are associated with negative side-effects and resistance mechanisms that have prompted the search for alternatives. The DHFR inhibitor pyrimethamine (Pyr) has compelling anti-cancer activity in in vivo models, but lacks potency compared to MTX, thereby requiring higher concentrations to induce therapeutic responses. The purpose of this work was to investigate structural analogues of Pyr to improve its in vitro and cellular activity. A series of 36 Pyr analogues were synthesized and tested in a sequence of in vitro and cell-based assays to monitor their DHFR inhibitory activity, cellular target engagement, and impact on breast cancer cell viability. Ten top compounds were identified, two of which stood out as potential lead candidates, 32 and 34. These functionalized Pyr analogues potently engaged DHFR in cells, at concentrations as low as 1 nM and represent promising DHFR inhibitors that could be further explored as potential anti-cancer agents., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023. Published by Elsevier Masson SAS.)

Background: Quassinoids are degraded triterpene compounds that can be obtained from various species of the Simaroubaceae plant family, including Eurycoma longifolia. Quassinoids are the major compounds in E. longifolia , and they are known to have various medicinal potentials, such as anticancer and antimalarial properties. Dihydrofolate reductase (DHFR) was reported to be one of the important targets for certain anticancer and antimalarial drugs. Twelve quassinoids from E. longifolia were identified to have anticancer effects based on their IC50 values. This study aimed to evaluate the interactions of these twelve quassinoids with DHFR via Autodock 4.2 software and Biovia Discovery Studio Visualiser., Methods: Twelve quassinoids from E. longifolia and their interactions with DHFR were evaluated via Autodock 4.2 software and Biovia Discovery Studio Visualiser. Their drug-likeness and pharmacokinetic properties were also assessed using the ADMETlab 2.0 program., Results: The molecular docking results showed that eleven quassinoids showed better docking scores than methotrexate, in which the binding energy (BE) of these quassinoids ranged from - 7.87 to -9.58 kcal/mol. Their inhibition constant (Ki) ranged from 0.095 to 1.71 μM. At the same time, the BE and Ki values for methotrexate were -7.80 kcal/mol and 1.64 μM, respectively., Conclusion: From the analysis, 6-dehydrolongilactone and eurycomalide B are among the twelve compounds that showed great potential as hit-to-lead compounds based on the docking score on DHFR, drug-likeness, and ADMET properties. These results suggest a great potential to pursue validation studies via in vitro and in vivo models., (Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.net.)

In-silico techniques offer a fast, accurate, reliable, and economical approach to studying the molecular interactions between compounds and proteins. In this study, our main aim is to use in-silico techniques as a rational approach for the prediction of the molecular drug targets for luteolin against Plasmodium falciparum . Multi-target molecular docking, 100 nanoseconds (ns) molecular dynamics (MD) simulations, and Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) binding free energy calculations were carried out for luteolin against dihydrofolate reductase thymidylate synthase ( Pf DHFR-TS), dihydroorotate dehydrogenase ( Pf DHODH), and falcipain-2. The native ligands of each protein were used as a reference to evaluate the performance of luteolin. Luteolin outperformed the native ligands of all proteins at molecular docking and MD simulations studies. However, in the MM-GBSA calculations, luteolin outperformed the native ligand of only Pf DHFR-TS but not Pf DHODH and falcipain-2. Among the studied proteins, the in-silico approach predicted Pf DHFR-TS as the most favorable drug target for luteolin.Communicated by Ramaswamy H. Sarma.

Fitness landscape theory predicts that rugged landscapes with multiple peaks impair Darwinian evolution, but experimental evidence is limited. In this study, we used genome editing to map the fitness of >260,000 genotypes of the key metabolic enzyme dihydrofolate reductase in the presence of the antibiotic trimethoprim, which targets this enzyme. The resulting landscape is highly rugged and harbors 514 fitness peaks. However, its highest peaks are accessible to evolving populations via abundant fitness-increasing paths. Different peaks share large basins of attraction that render the outcome of adaptive evolution highly contingent on chance events. Our work shows that ruggedness need not be an obstacle to Darwinian evolution but can reduce its predictability. If true in general, the complexity of optimization problems on realistic landscapes may require reappraisal.

Protein space is a rich analogy for genotype-phenotype maps, where amino acid sequence is organized into a high-dimensional space that highlights the connectivity between protein variants. It is a useful abstraction for understanding the process of evolution, and for efforts to engineer proteins towards desirable phenotypes. Few mentions of protein space consider how protein phenotypes can be described in terms of their biophysical components, nor do they rigorously interrogate how forces like epistasis-describing the nonlinear interaction between mutations and their phenotypic consequences-manifest across these components. In this study, we deconstruct a low-dimensional protein space of a bacterial enzyme (dihydrofolate reductase; DHFR) into "subspaces" corresponding to a set of kinetic and thermodynamic traits [k&#95;{cat}, K&#95;{M}, K&#95;{i}, and T&#95;{m} (melting temperature)]. We then examine how combinations of three mutations (eight alleles in total) display pleiotropy, or unique effects on individual subspace traits. We examine protein spaces across three orthologous DHFR enzymes (Escherichia coli, Listeria grayi, and Chlamydia muridarum), adding a genotypic context dimension through which epistasis occurs across subspaces. In doing so, we reveal that protein space is a deceptively complex notion, and that future applications to bioengineering should consider how interactions between amino acid substitutions manifest across different phenotypic subspaces.

Enzymes are complex macromolecules capable of catalyzing a wide variety of chemical reactions with high efficiency. Nonetheless, biological catalysis can be rudimentary. Here, we describe an enzyme that is built from a simple protein fold. This short protein sequence - almost a peptide - belongs to the ancient SH3 family of binding modules. Surprisingly, this binding module catalyzes the specific reduction of dihydrofolate using NADPH as a reducing cofactor, making this a dihydrofolate reductase. Too small to provide all the required binding and catalytic machinery on its own, it homotetramerizes, thus creating a large, central active site environment. Remarkably, none of the active site residues is essential to the catalytic function. Instead, backbone interactions juxtapose the reducing cofactor proximal to the target imine of the folate substrate, and a specific motion of the substrate promotes formation of the transition state. In this feature article, we describe the features that make this small protein a functional enzyme capable of catalyzing a metabolically essential reaction, highlighting the characteristics that make it a model for the evolution of primitive enzymes from binding modules.

The rapid emergence of various drug resistance and unfavourable aliphatic medication side effects endangers people's health. Phytocompounds with antibacterial activity and less harmful effects are known to be present in medicinal plants. Alkaloids from Eclipta alba were tested for their in vitro antibacterial capabilities and in silico docking studies against pathogenic bacteria and their target proteins in the current investigation. The alkaloid compounds verazine, ecliptine, 4-hydroxyverazine, 20-Epi-4beta-hydroxyverazine and hydroxyverazine were subjected to molecular docking studies to determine the method of binding as well as potential interactions and the docking score. The in vitro antibacterial activity of verazine alkaloid was assessed against two gram-positive and two gram-negative bacteria. Verazine alkaloid has the best inhibitory ability against DNA gyrase of E. coli (ΔG= -8.44 kcal/mol) and dihydrofolate reductase (DHFR) of S. aureus ( ΔG= -10.04 kcal/mol), according to docking studies. Verazine shown substantial in vitro antibacterial activity in this investigation against all test bacteria, with MIC and MBC values of 31.25 and 62.50 µg/mL for S. aureus and 15.63 and 31.25 µg/mL for B. cereus , respectively. The results of this work highlighted the value of unique alkaloid compounds from E. alba , which may offer effective antibacterial agents and DNA gyrase, DHFR inhibitors due to their novel structural properties capable of combating antimicrobial resistance. These findings call for more investigation into the compounds' function as antibacterial agents, as well as their unique-binding locations and mechanisms.

Dihydrofolate reductase (DHFR) is an important drug target and a highly studied model protein for understanding enzyme dynamics. DHFR's crucial role in folate synthesis renders it an ideal candidate to understand protein function and protein evolution mechanisms. In this study, to understand how a newly proposed DHFR inhibitor, 4'-deoxy methyl trimethoprim (4'-DTMP), alters evolutionary trajectories, we studied interactions that lead to its superior performance over that of trimethoprim (TMP). To elucidate the inhibition mechanism of 4'-DTMP, we first confirmed, both computationally and experimentally, that the relative binding free energy cost for the mutation of TMP and 4'-DTMP is the same, pointing the origin of the characteristic differences to be kinetic rather than thermodynamic. We then employed an interaction-based analysis by focusing first on the active site and then on the whole enzyme. We confirmed that the polar modification in 4'-DTMP induces additional local interactions with the enzyme, particularly, the M20 loop. These changes are propagated to the whole enzyme as shifts in the hydrogen bond networks. To shed light on the allosteric interactions, we support our analysis with network-based community analysis and show that segmentation of the loop domain of inhibitor-bound DHFR must be avoided by a successful inhibitor.

We investigated the relationship between mutations and dynamics in Escherichia coli dihydrofolate reductase (DHFR) using computational methods. Our study focused on the M20 and FG loops, which are known to be functionally important and affected by mutations distal to the loops. We used molecular dynamics simulations and developed position-specific metrics, including the dynamic flexibility index (DFI) and dynamic coupling index (DCI), to analyze the dynamics of wild-type DHFR and compared our results with existing deep mutational scanning data. Our analysis showed a statistically significant association between DFI and mutational tolerance of the DHFR positions, indicating that DFI can predict functionally beneficial or detrimental substitutions. We also applied an asymmetric version of our DCI metric (DCI asym ) to DHFR and found that certain distal residues control the dynamics of the M20 and FG loops, whereas others are controlled by them. Residues that are suggested to control the M20 and FG loops by our DCI asym metric are evolutionarily nonconserved; mutations at these sites can enhance enzyme activity. On the other hand, residues controlled by the loops are mostly deleterious to function when mutated and are also evolutionary conserved. Our results suggest that dynamics-based metrics can identify residues that explain the relationship between mutation and protein function or can be targeted to rationally engineer enzymes with enhanced activity., (© 2023 The Protein Society.)

Protein flexible loop regions were once thought to be simple linkers between other more functional secondary structural elements. However, as it becomes clearer that these loop domains are critical players in a plethora of biological processes, accurate conformational sampling of 3D loop structures is vital to the advancement of drug design techniques and the overall growth of knowledge surrounding molecular systems. While experimental techniques provide a wealth of structural information, the resolution of flexible loop domains is sometimes low or entirely absent due to their complex and dynamic nature. This highlights an opportunity for de novo structure prediction using in silico methods with molecular dynamics (MDs). This study evaluates some of the AMBER protein force field's (ffs) ability to accurately model dihydrofolate reductase (DHFR) conformations, a protein complex characterized by specific arrangements and interactions of multiple flexible loops whose conformations are determined by the presence or absence of bound ligands and cofactors. Although the AMBER ffs, including ff19SB, studied well model most protein structures with rich secondary structure, results obtained here suggest the inability to significantly sample the expected DHFR loop-loop conformations - of the six distinct protein-ligand systems simulated, a majority lacked consistent stabilization of experimentally derived metrics definitive the three enzyme conformations. Although under-sampling and the chosen ff parameter combinations could be the cause, given past successes with these MD approaches for many protein systems, this suggests a potential misbalance in available ff parameters required to accurately predict the structure of multiple flexible loop regions present in proteins.Communicated by Ramaswamy H. Sarma.

The present study is conducted to find the solution of rising antimicrobial resistance (AMR) in Escherichia coli which is a pathogen responsible for fatal systemic infections in human and animals. The enzyme dihydrofolate reductase (DHFR) is found in all organisms. In this study DHFR of E. coli (ec-DHFR) and human DHFR (h-DHFR) is targeted by novel chemical entities (NCE) from the Pathogen box of Medicines for Malaria Venture, Switzerland (MMV) using molecular modelling. The in-silico studies were further validated by in-vitro assays. The virtual screening of 400 MMV compounds was conducted using PyRx standard tool followed by manual docking of selected compounds by Autodock vina and Ligplot program. The in-silico studies showed good binding energy and strong hydrogen bond in docking of MMV675968 with ec-DHFR and no hydrogen bond with h-DHFR. This was further validated by the Molecular dynamic studies that revealed high binding free energy in ec-DHFR and in-vitro assays that produced good synergy in combination study of MMV675968 with last line (meropenem) and last resort (colistin) antibiotics. The extensive MD simulation and energetic analysis thus concludes that MMV675968 targets ec-DHFR. The combination studies were conducted with MMV675968 and FDA approved drugs against a panel of multidrug resistant Escherichia coli isolates. The synergistic results obtained in combination studies concluded that in-vitro data is consistent with in-silico data and that MMV675968 is a potential lead for future process of antimicrobial drug development against the multidrug resistance E. coli .Communicated by Ramaswamy H. Sarma.

Plasmodium falciparum dihydrofolate reductase enzyme ( P. falciparum DHFR) is one of the vital drug targets for malaria treatment, as this protein is indispensable for nucleotide metabolic pathways. This research aimed to discover promising phthalide derivatives against both wild and mutant P. falciparum DHFR enzymes through various computational techniques. The binding affinities were investigated using molecular docking, which showed five compounds having the highest affinity scores against both enzymes compared to the reference compounds. MM-GBSA calculations displayed favourable free binding energy. Moreover, the ADMET properties of the compounds are within acceptable ranges. The stability of the ligand-protein complexes was studied by Molecular Dynamics (MD) simulations. Depending on the results obtained from this research, we propose three compounds to be hit against P. falciparum DHFR activity which could be examined experimentally.Communicated by Ramaswamy Sarma.

Plasmodium falciparum dihydrofolate reductase-thymidylate synthase ( Pf DHFR-TS) is an important target enzyme in malarial chemotherapy. An understanding of how novel inhibitors interact with wild-type (wt Pf DHFR), quadruple-mutant (qm Pf DHFR), and human ( h DHFR) enzymes is required for the development of these compounds as antimalarials. This study is focused on a series of des -Cl and m -Cl phenyl analogs of pyrimethamine with various flexible 6-substituents. The interactions of these compounds with DHFR enzymes were investigated by 3 D-QSAR, MD simulations, MM-PBSA, and DFT calculations. CoMFA and CoMSIA models were developed with good predictive abilities for wt Pf DHFR and qm Pf DHFR. For h DHFR, CoMSIA models combined with clogP descriptor were successfully derived. Binding free energy using MM-PBSA and comparison of per residue decomposition energy analyses with the DFT method at M06-2X/6-31G ++(d,p) level of theory indicated that Asp54 and Phe58 play important roles in the binding of the most potent compound in the series (compound 27) with both wt Pf DHFR and qm Pf DHFR, whereas Arg59 and Arg122 were additionally found to interact with this inhibitor in qm Pf DHFR. For h DHFR, the residues Glu30 and Phe34 but not Arg70, equivalent to Asp54, Phe58, and Arg122 in Pf DHFR, also play role in compound 27 binding through strong hydrophobic interactions (Phe34) and hydrogen bond network with Glu30, Ile7, and Val115. From the key interactions identified in the DHFR-inhibitor complexes, a general scheme is proposed for designing new inhibitors selective for Pf DHFR that is important for the development of novel antifolate antimalarials.Communicated by Ramaswamy H. Sarma.

Dihydrofolate reductase (DHFR) is a ubiquitous cellular enzyme involved in the biosynthesis of nucleotide and protein precursors, thus, the inhibition of human DHFR can be a promising strategy in cancer treatment. The design of effective anticancer drugs is an urgent need today according to the high spread of cancer. The indole molecule with diverse mechanisms of action and anticancer properties is one of the efficient pharmacophores in drug design. Hence, a virtual library of indole derivatives as a scaffold was selected for designing safer and more effective anticancer drugs against DHFR in this work. All indole derivatives utilized in the library design were selected regarding appreciable tumor growth inhibition. Structure-activity relationship (SAR), docking energy, ADMET (absorption, distribution, metabolism, excretion, and toxicity) parameters, and effective non-covalent interactions were used to identify potential anticancer with indole scaffold. Results showed a higher number of indole moieties provide a strong attachment to the DHFR binding pocket and therefore more effective anticancer activity. The indole scaffold in combination with dichlorobenzene improves DHFR inhibition whereas barbituric acid weakens inhibition activity. In the following to validate the docking results, Molecular dynamics (MD) simulation and molecular mechanics generalized-Born surface area (MM-GBSA) indicated the permanent stability of the selected ligands into the DHFR binding pocket and the key amino acids. Therefore, promising pharmacophores based on indole-DHFR interactions were discovered, and the outcome could be useful in guiding future in vitro and in vivo drug discovery in cancer medicine.Communicated by Ramaswamy H. Sarma.

We present a potential mechanism for emergence of catalytic activity that is essential for survival, from a non-catalytic protein fold. The type B dihydrofolate reductase (DfrB) family of enzymes were first identified in pathogenic bacteria because their dihydrofolate reductase activity is sufficient to provide trimethoprim (TMP) resistance. DfrB enzymes are described as poorly evolved as a result of their unusual structural and kinetic features. No characterized protein shares sequence homology with DfrB enzymes; how they evolved to emerge in the modern resistome is unknown. In this work, we identify DfrB homologues from a database of putative and uncharacterized proteins. These proteins include an SH3-like fold homologous to the DfrB enzymes, embedded in a variety of additional structural domains. By means of functional, structural and biophysical characterization, we demonstrate that these distant homologues and their extracted SH3-like fold can display dihydrofolate reductase activity and confer TMP resistance. We provide evidence of tetrameric assembly and catalytic mechanism analogous to that of DfrB enzymes. These results contribute, to our knowledge, the first insights into a potential evolutionary path taken by this SH3-like fold to emerge in the modern resistome following introduction of TMP. This article is part of the theme issue 'Reactivity and mechanism in chemical and synthetic biology'.

Selective inhibitors of Escherichia coli dihydrofolate reductase (eDHFR) are crucial chemical biology tools that have widespread clinical applications. We developed a set of eDHFR-selective photoswitchable inhibitors by derivatizing the structure of our previously reported methotrexate (MTX) azolog, azoMTX . Substitution of the skeletal p -phenylene group of azoMTX with bulky bis-alkylated arylazopyrazole moieties significantly increased its selectivity toward eDHFR over human DHFR. Owing to the physical properties of arylazopyrazoles, the new ligands exhibited nearly complete Z-to-E photoconversion and high thermostability of Z-isomers. In addition, real-time photoreversible control of eDHFR activity was achieved by alternatively switching the illumination light wavelengths.

Our previous research suggests an important regulatory role of CK2-mediated phosphorylation of enzymes involved in the thymidylate biosynthesis cycle, i.e., thymidylate synthase (TS), dihydrofolate reductase (DHFR), and serine hydroxymethyltransferase (SHMT). The aim of this study was to show whether silencing of the CK2α gene affects TS and DHFR expression in A-549 cells. Additionally, we attempted to identify the endogenous kinases that phosphorylate TS and DHFR in CCRF-CEM and A-549 cells. We used immunodetection, immunofluorescence/confocal analyses, reverse transcription-quantitative polymerase chain reaction (RT-qPCR), in-gel kinase assay, and mass spectrometry analysis. Our results demonstrate that silencing of the CK2α gene in lung adenocarcinoma cells significantly increases both TS and DHFR expression and affects their cellular distribution. Additionally, we show for the first time that both TS and DHFR are very likely phosphorylated by endogenous CK2 in two types of cancer cells, i.e., acute lymphoblastic leukaemia and lung adenocarcinoma. Moreover, our studies indicate that DHFR is phosphorylated intracellularly by CK2 to a greater extent in leukaemia cells than in lung adenocarcinoma cells. Interestingly, in-gel kinase assay results indicate that the CK2α' isoform was more active than the CK2α subunit. Our results confirm the previous studies concerning the physiological relevance of CK2-mediated phosphorylation of TS and DHFR.

Dihydrofolate reductase (DHFR) is a crucial enzyme that maintains the levels of 5,6,7,8-tetrahydrofolate (THF) required for the biological synthesis of the building blocks of DNA, RNA, and proteins. Over-activation of DHFR results in the progression of multiple pathological conditions such as cancer, bacterial infection, and inflammation. Therefore, DHFR inhibition plays a major role in treating these illnesses. Sesquiterpenes of various types are prime metabolites derived from the marine sponge Dactylospongia elegans and have demonstrated antitumor, anti-inflammation, and antibacterial capacities. Here, we investigated the in silico potential inhibitory effects of 87 D. elegans metabolites on DHFR and predicted their ADMET properties. Compounds were prepared computationally for molecular docking into the selected crystal structure of DHFR (PDB: 1KMV). The docking scores of metabolites 34 , 28 , and 44 were the highest among this series (gscore values of -12.431, -11.502, and -10.62 kcal/mol, respectively), even above the co-crystallized inhibitor SRI-9662 score (-10.432 kcal/mol). The binding affinity and protein stability of these top three scored compounds were further estimated using molecular dynamic simulation. Compounds 34 , 28 , and 44 revealed high binding affinity to the enzyme and could be possible leads for DHFR inhibitors; however, further in vitro and in vivo investigations are required to validate their potential.

There is an urgent need to discover and develop novel drugs to combat Mycobacterium tuberculosis , the causative agent of tuberculosis (TB) in humans. Alkaloids have been shown to have wide-ranging therapeutic application and could be ideal candidates for drug development, and research is underway to develop new anti-tubercular drugs from natural sources. In this regard, the current research deals with finding novel lead compounds from the Withania somnifera (WS) plant. Broad health benefits of WS are due to the presence of diverse chemical constituents which include anaferine and anahygrine and which belong to the alkaloid family. In the present study, these two compounds have been theoretically studied to understand their electronic properties using the density functional theory (DFT) at the B3LYP/6-311 + G (d,p) level. HOMO and LUMO properties and molecular electrostatic potential (MEP) surface were calculated. Further, to understand the mechanism of action of these compounds and to identify their putative drug target, molecular docking and dynamics studies were employed against Mycobacterium tuberculosis dihydrofolate reductase (DHFR). It was determined that NADP + affects stability of the complexes by reducing fluctuations of residues 14-23 and 117-126. It was also found that Ile5 and Gln28 play an important role in complexation. Electron density analysis (using the AlteQ method) of the intermolecular region, analyzing both the anaferin-NADP + and anahygrin-NADP + complexes showed that anaferin and anahygrin complexes are more stable in the presence of NADP + . It has been established that in most intermolecular contacts the contribution of the ligand to the electron density is greater than that of NADP + . The present study thus provides an excellent way to analyze the effect of anaferine and anahygrine in essential processes of M. tuberculosis .Communicated by Ramaswamy H. Sarma.

In all species, dihydrofolate reductase (DHFR) is an essential enzyme that regulates the cellular amount of tetrahydrofolate. Human DHFR (hDHFR) activity inhibition results in tetrahydrofolate depletion and cell death. This property has made hDHFR a therapeutic target for cancer. Methotrexate is a well-known hDHFR inhibitor, but its administration has shown some light to severe adverse effects. Therefore, we aimed to find new potential hDHFR inhibitors using structure-based virtual screening, ADMET prediction, molecular docking, and molecular dynamics simulations. Here, we used the PubChem database to find all compounds with at least 90% structural similarity to known natural DHFR inhibitors. To explore their interaction pattern and estimate their binding affinities, the screened compounds (2023) were subjected to structure-based molecular docking against hDHFR. The fifteen compounds that showed higher binding affinity to the hDHFR than the reference compound (methotrexate) displayed important molecular orientation and interactions with key residues in the enzyme's active site. These compounds were subjected to Lipinski and ADMET prediction. PubChem CIDs: 46886812 and 638190 were identified as putative inhibitors. In addition, molecular dynamics simulations revealed that the binding of compounds (CIDs: 46886812 and 63819) stabilized the hDHFR structure and caused minor conformational changes. Our findings suggest that two compounds (CIDs: 46886812 and 63819) could be promising potential inhibitors of hDHFR in cancer therapy.Communicated by Ramaswamy H. Sarma.

Antibiotic resistance is a global health problem in which mutations occurring in functional proteins render drugs ineffective. The working mechanisms of the arising mutants are seldom apparent; a methodology to decipher these mechanisms systematically would render devising therapies to control the arising mutational pathways possible. Here we utilize C α -C β bond vector relaxations obtained from moderate length MD trajectories to determine conduits for functionality of the resistance conferring mutants of Escherichia coli dihydrofolate reductase. We find that the whole enzyme is synchronized to the motions of the substrate, irrespective of the mutation introducing gain-of-function or loss-of function. The total coordination of the motions suggests changes in the hydrogen bond dynamics with respect to the wild type as a possible route to determine and classify the mode-of-action of individual mutants. As a result, nine trimethoprim-resistant point mutations arising frequently in evolution experiments are categorized. One group of mutants that display the largest occurrence (L28R, W30G) work directly by modifying the dihydrofolate binding region. Conversely, W30R works indirectly by the formation of the E139-R30 salt bridge which releases energy resulting from tight binding by distorting the binding cavity. A third group (D27E, F153S, I94L) arising as single, resistance invoking mutants in evolution experiment trajectories allosterically and dynamically affects a hydrogen bonding motif formed at residues 59-69-71 which in turn modifies the binding site dynamics. The final group (I5F, A26T, R98P) consists of those mutants that have properties most similar to the wild type; these only appear after one of the other mutants is fixed on the protein structure and therefore display clear epistasis. Thus, we show that the binding event is governed by the entire enzyme dynamics while the binding site residues play gating roles. The adjustments made in the total enzyme in response to point mutations are what make quantifying and pinpointing their effect a hard problem. Here, we show that hydrogen bond dynamics recorded on sub-μs time scales provide the necessary fingerprints to decipher the various mechanisms at play.

Conventional quantum mechanical-molecular mechanics (QM/MM) simulation approaches for modeling enzyme reactions often assume that there is one dominant reaction pathway and that this pathway can be sampled starting from an X-ray structure of the enzyme. These assumptions reduce computational cost; however, their validity has not been extensively tested. This is due in part to the lack of a rigorous formalism for integrating disparate pathway information from dynamical QM/MM calculations. Here, we present a way to model ensembles of reaction pathways efficiently using a divide-and-conquer strategy through Hierarchical Markov State Modeling (Hi-MSM). This approach allows information on multiple, distinct pathways to be incorporated into a chemical kinetic model, and it allows us to test these two assumptions. Applying Hi-MSM to the reaction carried out by dihydrofolate reductase (DHFR) we find (i) there are multiple, distinct pathways significantly contributing to the overall flux of the reaction that the conventional approach does not identify and (ii) that the conventional approach does not identify the dominant reaction pathway. Thus, both assumptions underpinning the conventional approach are violated. Since DHFR is a relatively small enzyme, and configuration space scales exponentially with protein size, accounting for multiple reaction pathways is likely to be necessary for most enzymes.

A series of 5-[(phenethylamino)methyl]pyrimidine-2,4-diamines were assessed in silico as potential inhibitors of Plasmodium falciparum dihydrofolate reductase (PfDHFR), synthesised and tested for inhibitory activity against PfDHFR in vitro. The compounds displayed promising inhibitory activity against both wild-type (K i 1.3-243 nM) and quadruple mutant (K i 13-208 nM) PfDHFR in the biochemical enzyme assay, but were less potent in the whole-cell P. falciparum assay (IC 50 (TM4/8.2) 0.4-28 μM; IC 50 (V1S) 3.7-54 μM). Further investigation into the pharmacokinetic properties of these compounds may guide the development of more potent analogues., (© 2022 The Authors. ChemMedChem published by Wiley-VCH GmbH.)

Introduction: Dihydrofolate reductase (DHFR) plays an important role in the biosynthesis of amino acid and folic acid. It participates by reducing dihydrofolate to tetrahydrofolate, in the presence of nicotinamide dinucleotide phosphate cofactor, and has been verified by various clinical studies to use DHFR as a target for the treatment of cancer and various bacterial infections., Area Covered: In this review, we have disclosed patents of synthetics and natural DHFR inhibitors with diaminopyrimidine and quinazoline nucleus from 2001. Additionally, this review highlights the clinical progression of numerous DHFR inhibitors received from the last five years., Expert Opinion: From 2001 to 2021, numerous active chemical scaffolds have been introduced and are exposed as lead candidates that have entered clinical trials as potent DHFR inhibitors. Moreover, researchers have paid considerable attention to the development of a new class of DHFR inhibitors with higher selectivity and potency. This development includes synthesis of synthetic as well as natural compounds that are potent DHFR inhibitors. On the basis of literature review, we can anticipate that there are a huge number of novel active molecules available for the future that could possess superior abilities to target this enzyme with a profound pharmacological profile.

Dihydrofolate reductase (DHFR) is a well-studied, clinically relevant enzyme known for being highly dynamic over the course of its catalytic cycle. However, the role dynamic motions play in the explicit hydride transfer from the nicotinamide cofactor to the dihydrofolate substrate remains unclear because reaction initiation and direct spectroscopic examination on the appropriate time scale for such femtosecond to picosecond motions is challenging. Here, we employ pre-steady-state kinetics to observe the hydride transfer as directly as possible in two different species of DHFR: Escherichia coli and Homo sapiens . While the hydride transfer has been well-characterized in DHFR from E. coli , improvements in time resolution now allow for sub-millisecond dead times for stopped-flow spectroscopy, which reveals that the maximum rate is indeed faster than previously recorded. The rate in the human enzyme, previously only estimated, is also able to be directly observed using cutting-edge stopped-flow instrumentation. In addition to the pH dependence of the hydride transfer rates for both enzymes, we examine the primary H/D kinetic isotope effect to reveal a temperature dependence in the human enzyme that is absent from the E. coli counterpart. This dependence, which appears above a temperature of 15 °C is a shared feature among other hydride transfer enzymes and is also consistent with computational work suggesting the presence of a fast promoting-vibration that provides donor-acceptor compression on the time scale of catalysis to facilitate the chemistry step.

In this issue of Structure, Kampmeyer et al. provide detailed mechanistic insights into how structural changes in disease-associated dihydrofolate reductase (DHFR) missense variants affect their cellular protein abundance and discuss implications for hereditary megaloblastic anemia disease., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022. Published by Elsevier Ltd.)

New antibiotics are urgently needed to counter the emergence of antimicrobial-resistant pathogenic bacteria. A major challenge in antibiotic drug discovery is to turn potent biochemical inhibitors of essential bacterial components into effective antimicrobials. This difficulty is underpinned by a lack of methods to investigate the physicochemical properties needed for candidate antibiotics to permeate the bacterial cell envelope and avoid clearance by the action of bacterial efflux pumps. To address these issues, here we used a target engagement assay to measure the equilibrium and kinetic binding parameters of antibiotics targeting dihydrofolate reductase (DHFR) in live bacteria. We also used this assay to identify novel DHFR ligands having antimicrobial activity. We validated this approach using the Gram-negative bacteria Escherichia coli and the emerging human pathogen Mycobacterium abscessus . We expect the use of target engagement assays in bacteria to expedite the discovery and progression of novel, cell-permeable antibiotics with on-target activity.