Evaluation of quercetin as a potential β-lactamase CTX-M-15 inhibitor via the molecular docking,dynamics simulations,and MMGBSA

Antimicrobial resistance (AMR) threatens millions of people around the world and has been declared a global risk by the World Economic Forum. One of the important AMR mechanisms in Enterobacteriaceae is the production of extended-spectrum β-lactamases. The most common ESBL, CTX-M β-lactamases, is spread to the world by CTX-M-15 and CTX-M-14. Sulbactam, clavulanic acid, and tazobactam are first-generation β-lactamase inhibitors and avibactam is a new non-β-lactam β-lactamase inhibitor. We studied that avibactam, sulbactam, clavulanic acid, tazobactam, and quercetin natural flavonoids were docked to target protein CTXM-15. Subsequently, the complexes were simulated using the molecular dynamics simulations method during 100 ns for determining the final binding positions of ligands. Clavulanic acid left CTX-M-15 and other ligands remained in the binding site after the simulation. The estimated binding energies were calculated during 100 ns simulation by the MMGBSA-MMPBSA method. The estimated free binding energies of avibactam, sulbactam, quercetin, tazobactam, and clavulanic acid were sorted as –33.61 kcal/mol, –16.04 kcal/mol, –14 kcal/mol, –12.68 kcal/mol, and –2.95 kcal/mol. As a result of both final binding positions and free binding energy calculations, Quercetin may be evaluated an alternative candidate and a more potent β-lactamases inhibitor for new antimicrobial combinations to CTX-M-15. The results obtained in silico studies are predicted to be a preliminary study for in vitro studies for quercetin and similar bioactive natural compounds. These studies are notable for the discovery of natural compounds that can be used in the treatment of infections caused by β-lactamase-producing pathogens.     

Mechanistic Insight into Binding of Huperzine A with Human Serum Albumin: Computational and Spectroscopic Approaches

Human serum albumin (HSA) is the most abundant protein in plasma synthesized by the liver and the main modulator of fluid distribution between body compartments. It has an amazing capacity to bind with multiple ligands, offering a store and transporter for various endogenous and exogenous compounds. Huperzine A (HpzA) is a natural sesquiterpene alkaloid found in Huperzia serrata and used in various neurological conditions, including Alzheimer’s disease (AD). This study elucidated the binding of HpzA with HSA using advanced computational approaches such as molecular docking and molecular dynamic (MD) simulation followed by fluorescence-based binding assays. The molecular docking result showed plausible interaction between HpzA and HSA. The MD simulation and principal component analysis (PCA) results supported the stable interactions of the protein–ligand complex. The fluorescence assay further validated the in silico study, revealing significant binding affinity between HpzA and HSA. This study advocated that HpzA acts as a latent HSA binding partner, which may be investigated further in AD therapy in experimental settings.     

Novel and Potential Small Molecule Scaffolds as DYRK1A Inhibitors by Integrated Molecular Docking-Based Virtual Screening and Dynamics Simulation Study

The dual-specificity tyrosine phosphorylation-regulated kinase 1A (DYRK1A) is a novel, promising and emerging biological target for therapeutic intervention in neurodegenerative diseases, especially in Alzheimer’s disease (AD). The molMall database, comprising rare, diverse and unique compounds, was explored for molecular docking-based virtual screening against the DYRK1A protein, in order to find out potential inhibitors. Ligands exhibiting hydrogen bond interactions with key amino acid residues such as Ile165, Lys188 (catalytic), Glu239 (gk+1), Leu241 (gk+3), Ser242, Asn244, and Asp307, of the target protein, were considered potential ligands. Hydrogen bond interactions with Leu241 (gk+3) were considered key determinants for the selection. High scoring structures were also docked by Glide XP docking in the active sites of twelve DYRK1A related protein kinases, viz. DYRK1B, DYRK2, CDK5/p25, CK1, CLK1, CLK3, GSK3β, MAPK2, MAPK10, PIM1, PKA, and PKCα, in order to find selective DYRK1A inhibitors. MM/GBSA binding free energies of selected ligand–protein complexes were also calculated in order to remove false positive hits. Physicochemical and pharmacokinetic properties of the selected six hit ligands were also computed and related with the proposed limits for orally active CNS drugs. The computational toxicity webserver ProTox-II was used to predict the toxicity profile of selected six hits (molmall IDs 9539, 11352, 15938, 19037, 21830 and 21878). The selected six docked ligand–protein systems were exposed to 100 ns molecular dynamics (MD) simulations to validate their mechanism of interactions and stability in the ATP pocket of human DYRK1A kinase. All six ligands were found to be stable in the ATP binding pocket of DYRK1A kinase.     

Halogen Bonding in Haspin-Halogenated Tubercidin Complexes: Molecular Dynamics and Quantum Chemical Calculations

Haspin, an atypical serine/threonine protein kinase, is a potential target for cancer therapy. 5-iodotubercidin (5-iTU), an adenosine derivative, has been identified as a potent Haspin inhibitor in vitro. In this paper, quantum chemical calculations and molecular dynamics (MD) simulations were employed to identify and quantitatively confirm the presence of halogen bonding (XB), specifically halogen∙∙∙π (aromatic) interaction between halogenated tubercidin ligands with Haspin. Consistent with previous theoretical finding, the site specificity of the XB binding over the ortho-carbon is identified in all cases. A systematic increase of the interaction energy down Group 17, based on both quantum chemical and MD results, supports the important role of halogen bonding in this series of inhibitors. The observed trend is consistent with the experimental observation of the trend of activity within the halogenated tubercidin ligands (F < Cl < Br < I). Furthermore, non-covalent interaction (NCI) plots show that cooperative non-covalent interactions, namely, hydrogen and halogen bonds, contribute to the binding of tubercidin ligands toward Haspin. The understanding of the role of halogen bonding interaction in the ligand–protein complexes may shed light on rational design of potent ligands in the future.     

Pharmacological Characterization of µ-Opioid Receptor Agonists with Biased G Protein or β-Arrestin Signaling,and Computational Study of Conformational Changes during Receptor Activation

In recent years, G protein vs. β-arrestin biased agonism at opioid receptors has been proposed as an opportunity to produce antinociception with reduced adverse effects. However, at present this approach is highly debated, a reason why more information about biased ligands is required. While the practical relevance of bias in the case of µ-opioid receptors (MOP) still needs to be validated, it remains important to understand the basis of this bias of MOP (and other GPCRs). Recently, we reported two cyclopeptides with high affinity for MOP, the G protein biased Dmt-c[d-Lys-Phe-pCF₃-Phe-Asp]NH₂ (F-81), and the β-arrestin 2 biased Dmt-c[d-Lys-Phe-Asp]NH₂ (C-33), as determined by calcium mobilization assay and bioluminescence resonance energy transfer-based assay. The biased character of F-81 and C-33 has been further analyzed in the [³⁵S]GTPγS binding assay in human MOP-expressing cells, and the PathHunter enzyme complementation assay, used to measure β-arrestin 2 recruitment. To investigate the structural features of peptide-MOP complexes, we performed conformational analysis by NMR spectroscopy, molecular docking, and molecular dynamics simulation. These studies predicted that the two ligands form alternative complexes with MOP, engaging specific ligand–receptor contacts. This would induce different displays of the cytosolic side of the seven-helices bundle, in particular by stabilizing different angulations of helix 6, that could favor intracellular coupling to either G protein or β-arrestin.     

A–X⋯σ Interactions—Halogen Bonds with σ-Electrons as the Lewis Base Centre

CCSD(T)/aug-cc-pVTZ//ωB97XD/aug-cc-pVTZ calculations were performed for halogen-bonded complexes. Here, the molecular hydrogen, cyclopropane, cyclobutane and cyclopentane act as Lewis base units that interact through the electrons of the H–H or C–C σ-bond. The FCCH, ClCCH, BrCCH and ICCH species, as well as the F₂, Cl₂, Br₂ and I₂ molecular halogens, act as Lewis acid units in these complexes, interacting through the σ-hole localised at the halogen centre. The Quantum Theory of Atoms in Molecules (QTAIM), the Natural Bond Orbital (NBO) and the Energy Decomposition Analysis (EDA) approaches were applied to analyse these aforementioned complexes. These complexes may be classified as linked by A–X···σ halogen bonds, where A = C, X (halogen). However, distinct properties of these halogen bonds are observed that depend partly on the kind of electron donor: dihydrogen, cyclopropane, or another cycloalkane. Examples of similar interactions that occur in crystals are presented; Cambridge Structural Database (CSD) searches were carried out to find species linked by the A–X···σ halogen bonds.     

Anti-Inflammatory and Anti-Rheumatic Potential of Selective Plant Compounds by Targeting TLR-4/AP-1 Signaling: A Comprehensive Molecular Docking and Simulation Approaches

Plants are an important source of drug development and numerous plant derived molecules have been used in clinical practice for the ailment of various diseases. The Toll-like receptor-4 (TLR-4) signaling pathway plays a crucial role in inflammation including rheumatoid arthritis. The TLR-4 binds with pro-inflammatory ligands such as lipopolysaccharide (LPS) to induce the downstream signaling mechanism such as nuclear factor κappa B (NF-κB) and mitogen activated protein kinases (MAPKs). This signaling activation leads to the onset of various diseases including inflammation. In the present study, 22 natural compounds were studied against TLR-4/AP-1 signaling, which is implicated in the inflammatory process using a computational approach. These compounds belong to various classes such as methylxanthine, sesquiterpene lactone, alkaloid, flavone glycosides, lignan, phenolic acid, etc. The compounds exhibited different binding affinities with the TLR-4, JNK, NF-κB, and AP-1 protein due to the formation of multiple hydrophilic and hydrophobic interactions. With TLR-4, rutin had the highest binding energy (−10.4 kcal/mol), poncirin had the highest binding energy (−9.4 kcal/mol) with NF-κB and JNK (−9.5 kcal/mol), respectively, and icariin had the highest binding affinity (−9.1 kcal/mol) with the AP-1 protein. The root means square deviation (RMSD), root mean square fraction (RMSF), and radius of gyration (RoG) for 150 ns were calculated using molecular dynamic simulation (MD simulation) based on rutin’s greatest binding energy with TLR-4. The RMSD, RMSF, and RoG were all within acceptable limits in the MD simulation, and the complex remained stable for 150 ns. Furthermore, these compounds were assessed for the potential toxic effect on various organs such as the liver, heart, genotoxicity, and oral maximum toxic dose. Moreover, the blood–brain barrier permeability and intestinal absorption were also predicted using SwissADME software (Lausanne, Switzerland). These compounds exhibited promising physico-chemical as well as drug-likeness properties. Consequently, these selected compounds portray promising anti-inflammatory and drug-likeness properties.     

Reactive Acrylamide-Modified DNA Traps for Accurate Cross-Linking with Cysteine Residues in DNA–Protein Complexes Using Mismatch Repair Protein MutS as a Model

Covalent protein capture (cross-linking) by reactive DNA derivatives makes it possible to investigate structural features by fixing complexes at different stages of DNA–protein recognition. The most common cross-linking methods are based on reactive groups that interact with native or engineered cysteine residues. Nonetheless, high reactivity of most of such groups leads to preferential fixation of early-stage complexes or even non-selective cross-linking. We synthesised a set of DNA reagents carrying an acrylamide group attached to the C5 atom of a 2′-deoxyuridine moiety via various linkers and studied cross-linking with MutS as a model protein. MutS scans DNA for mismatches and damaged nucleobases and can form multiple non-specific complexes with DNA that may cause non-selective cross-linking. By varying the length of the linker between DNA and the acrylamide group and by changing the distance between the reactive nucleotide and a mismatch in the duplex, we showed that cross-linking occurs only if the distance between the acrylamide group and cysteine is optimal within the DNA–protein complex. Thus, acrylamide-modified DNA duplexes are excellent tools for studying DNA–protein interactions because of high selectivity of cysteine trapping.     

Anti-Metatype Antibody Screening,Sandwich Immunoassay Development,and Structural Insights for β-Lactams Based on Penicillin Binding Protein

Theoretically, sandwich immunoassay is more sensitive and has a wider working range than that of competitive format. However, it has been thought that small molecules cannot be detected by the sandwich format due to their limited size. In the present study, we proposed a novel strategy for achieving sandwich immunoassay of β-lactams with low molecular weights. Firstly, five β-lactam antibiotics were selected to bind with penicillin binding protein (PBP)2x* to form complexes. Then, monoclonal and polyclonal antibodies against PBP2x*-β-lactams complexes were produced by animal immunization. Subsequently, the optimal pairing antibodies were utilized to establish sandwich immunoassay for detection of 18 PBP2x*-β-lactam complexes. Among them, ceftriaxone could be detected at as low as 1.65 ng/mL with working range of 1–1000 ng/mL in milk. To reveal the detection mechanism, computational chemistry and molecular recognition study were carried out. The results showed that β-lactams with a large size and complex structures maybe conducive to induce conformational changes of PBP2x*, and then exhibit greater possibility of being detected by sandwich immunoassay after combination with PBP2x*. This study provides insights for subsequent investigations of anti-metatype antibody screening and sandwich immunoassay establishment for small-molecule detection.     

Aromatic Rings as Molecular Determinants for the Molecular Recognition of Protein Kinase Inhibitors

Protein kinases are key enzymes in many signal transduction pathways, and play a crucial role in cellular proliferation, differentiation, and various cell regulatory processes. However, aberrant function of kinases has been associated with cancers and many other diseases. Consequently, competitive inhibition of the ATP binding site of protein kinases has emerged as an effective means of curing these diseases. Over the past three decades, thousands of protein kinase inhibitors (PKIs) with varying molecular frames have been developed. Large-scale data mining of the Protein Data Bank resulted in a database of 2139 non-redundant high-resolution X-ray crystal structures of PKIs bound to protein kinases. This provided us with a unique opportunity to study molecular determinants for the molecular recognition of PKIs. A chemoinformatic analysis of 2139 PKIs resulted in findings that PKIs are “flat” molecules with high aromatic ring counts and low fractions of sp³ carbon. All but one PKI possessed one or more aromatic rings. More importantly, it was found that the average weighted hydrogen bond count is inversely proportional to the number of aromatic rings. Based on this linear relationship, we put forward the exchange rule of hydrogen bonding interactions and non-bonded π-interactions. Specifically, a loss of binding affinity caused by a decrease in hydrogen bonding interactions is compensated by a gain in binding affinity acquired by an increase in aromatic ring-originated non-bonded interactions (i.e., π–π stacking interactions, CH–π interactions, cation–π interactions, etc.), and vice versa. The very existence of this inverse relationship strongly suggests that both hydrogen bonding and aromatic ring-originated non-bonded interactions are responsible for the molecular recognition of PKIs. As an illustration, two representative PKI–kinase complexes were employed to examine the relative importance of different modes of non-bonded interactions for the molecular recognition of PKIs. For this purpose, two FDA-approved PKI drugs, ibrutinib and lenvatinib, were chosen. The binding pockets of both PKIs were thoroughly examined to identify all non-bonded intermolecular interactions. Subsequently, the strengths of interaction energies between ibrutinib and its interacting residues in tyrosine kinase BTK were quantified by means of the double hybrid DFT method B2PLYP. The resulting energetics for the binding of ibrutinib in tyrosine kinase BTK showed that CH–π interactions and π–π stacking interactions between aromatic rings of the drug and hydrophobic residues in its binding pocket dominate the binding interactions. Thus, this work establishes that, in addition to hydrogen bonding, aromatic rings function as important molecular determinants for the molecular recognition of PKIs. In conclusion, our findings support the following pharmacophore model for ATP-competitive kinase inhibitors: a small molecule features a scaffold of one or more aromatic rings which is linked with one or more hydrophilic functional groups. The former has the structural role of acting as a scaffold and the functional role of participating in aromatic ring-originated non-bonded interactions with multiple hydrophobic regions in the ATP binding pocket of kinases. The latter ensure water solubility and form hydrogen bonds with the hinge region and other hydrophilic residues of the ATP binding pocket.     

Probing conformational hotspots for the recognition and intervention of protein complexes by lysine reactivity profiling

Probing the conformational and functional hotspot sites within aqueous native protein complexes is still a challenging task. Herein, a mass spectrometry (MS)-based two-step isotope labeling-lysine reactivity profiling (TILLRP) strategy is developed to quantify the reactivities of lysine residues and probe the molecular details of protein–protein interactions as well as evaluate the conformational interventions by small-molecule active compounds. The hotspot lysine sites that are crucial to the SARS-CoV-2 S1–ACE2 combination could be successfully probed, such as S1 Lys⁴¹⁷ and Lys⁴⁴⁴. Significant alteration of the reactivities of lysine residues at the interaction interface of S1-RBD Lys³⁸⁶–Lys⁴⁶² was observed during the formation of complexes, which might be utilized as indicators for investigating the S1-ACE2 dynamic recognition and intervention at the molecular level in high throughput.

A mass spectrometry-based two-step isotope labeling-lysine reactivity profiling strategy is developed to probe the molecular details of protein–protein interactions and evaluate the conformational interventions by small-molecule active compounds.     

Photocontrol of CRISPR/Cas9 function by site-specific chemical modification of guide RNA

The function of CRISPR/Cas9 can be conditionally controlled by the rational engineering of guide RNA (gRNA) to target the gene of choice for precise manipulation of the genome. Particularly, chemically modified gRNA that can be activated by using specific stimuli provides a unique tool to expand the versatility of conditional control. Herein, unlike previous engineering of gRNA that generally focused on the RNA part only but neglected RNA–protein interactions, we aimed at the interactive sites between 2′-OH of ribose in the seed region of gRNA and the Cas9 protein and identified that chemical modifications at specific sites could be utilized to regulate the Cas9 activity. By introducing a photolabile group at these specific sites, we achieved optical control of Cas9 activity without disrupting the Watson–Crick base pairing. We further examined our design through CRISPR-mediated gene activation and nuclease cleavage in living cells and successfully manipulated the gene expression by using light irradiation. Our site-specific modification strategy exhibited a highly efficient and dynamic optical response and presented a new perspective for manipulating gRNA based on the RNA–protein interaction rather than the structure of RNA itself. In addition, these specific sites could also be potentially utilized for modification of other stimuli-responsive groups, which would further enrich the toolbox for conditional control of CRISPR/Cas9 function.

The CRISPR/Cas9 function is optically controlled in living cells by the site-specifically caged guide RNA based on the RNA–protein interaction.     

The Theoretical and Experimental Investigation of the Fluorinated Palladium β-Diketonate Derivatives: Structure and Physicochemical Properties

To search for new suitable Pd precursors for MOCVD/ALD processes, the extended series of fluorinated palladium complexes [Pd(CH₃CXCHCO(R))₂] with β-diketone [tfa−1,1,1-trifluoro-2,4-pentanedionato (1); pfpa−5,5,6,6,6-pentafluoro-2,4-hexanedionato (3); hfba−5,5,6,6,7,7,7-heptafluoro-2,4-heptanedionato (5)] and β-iminoketone [i-tfa−1,1,1-trifluoro-2-imino-4-pentanonato (2); i-pfpa−5,5,6,6,6-pentafluoro-2-imino-4-hexanonato (4); i-hfba-5,5,6,6,7,7,7-heptafluoro-2-imino-4-heptanonato (6)] ligands were synthesized with 70–80% yields and characterized by a set of experimental (SXRD, XRD, IR, NMR spectroscopy, TG) and theoretical (DFT, Hirshfeld surface analysis) methods. Solutions of Pd β-diketonates contained both cis and trans isomers, while only trans isomers were detected in the solutions of Pd β-iminoketonates. The molecules 2–6 and new polymorphs of complexes 3 and 5 were arranged preferentially in stacks, and the distance between molecules in the stack generally increased with elongation of the fluorine chain in ligands. The H…F contacts were the main ones involved in the formation of packages of molecules 1–2, and C…F, F…F, NH…F contacts appeared in the structures of complexes 4–6. The stability of complexes and their polymorphs in the crystal phases were estimated from DFT calculations. The TG data showed that the volatility differences between Pd β-iminoketonates and Pd β-diketonates were minimized with the elongation of the fluorine chain in the ligands.     

The Exchange of Cyclometalated Ligands

Reactions of cyclometalated compounds are numerous. This account is focused on one of such reactions, the exchange of cyclometalated ligands, a reaction between a cyclometalated compound and an incoming ligand that replaces a previously cyclometalated ligand to form a new metalacycle: + H-C*~Z ⇄ + H-C~Y. Originally discovered for Pd^II complexes with Y/Z = N, P, S, the exchange appeared to be a mechanistically challenging, simple, and convenient routine for the synthesis of cyclopalladated complexes. Over four decades it was expanded to cyclometalated derivatives of platinum, ruthenium, manganese, rhodium, and iridium. The exchange, which is also questionably referred to as transcyclometalation, offers attractive synthetic possibilities and assists in disclosing key mechanistic pathways associated with the C–H bond activation by transition metal complexes and C–M bond cleavage. Both synthetic and mechanistic aspects of the exchange are reviewed and discussed.     

Magnetic Resonance Access to Transiently Formed Protein Complexes

Protein–protein interactions are of utmost importance to an understanding of biological phenomena since non-covalent and therefore reversible couplings between basic proteins leads to the formation of complex regulatory and adaptive molecular systems. Such systems are capable of maintaining their integrity and respond to external stimuli, processes intimately related to living organisms. These interactions, however, span a wide range of dissociation constants, from sub-nanomolar affinities in tight complexes to high-micromolar or even millimolar affinities in weak, transiently formed protein complexes. Herein, we demonstrate how novel NMR and EPR techniques can be used for the characterization of weak protein–protein (ligand) complexes. Applications to intrinsically disordered proteins and transiently formed protein complexes illustrate the potential of these novel techniques to study hitherto unobserved (and unobservable) higher-order structures of proteins.     

Enhancement of Haloperidol Binding Affinity to Dopamine Receptor via Forming a Charge-Transfer Complex with Picric Acid and 7,7,8,8-Tetracyanoquinodimethane for Improvement of the Antipsychotic Efficacy

Haloperidol (HPL) is a typical antipsychotic drug used to treat acute psychotic conditions, delirium, and schizophrenia. Solid charge transfer (CT) products of HPL with 7,7,8,8-tetracyanoquinodimethane (TCNQ) and picric acid (PA) have not been reported till date. Therefore, we conducted this study to investigate the donor–acceptor CT interactions between HPL (donor) and TCNQ and PA (π-acceptors) in liquid and solid states. The complete spectroscopic and analytical analyses deduced that the stoichiometry of these synthesized complexes was 1:1 molar ratio. Molecular docking calculations were performed for HPL as a donor and the resulting CT complexes with TCNQ and PA as acceptors with two protein receptors, serotonin and dopamine, to study the comparative interactions among them, as they are important neurotransmitters that play a large role in mental health. A molecular dynamics simulation was ran for 100 ns with the output from AutoDock Vina to refine docking results and better examine the molecular processes of receptor–ligand interactions. When compared to the reactant donor, the CT complex [(HPL)(TCNQ)] interacted with serotonin and dopamine more efficiently than HPL only. CT complex [(HPL)(TCNQ)] with dopamine (CTtD) showed the greatest binding energy value among all. Additionally, CTtD complex established more a stable interaction with dopamine than HPL–dopamine.     

Nitrogen fixation catalyzed by ferrocene-substituted dinitrogen-bridged dimolybdenum–dinitrogen complexes: unique behavior of ferrocene moiety as redox active site

A series of dinitrogen-bridged dimolybdenum–dinitrogen complexes bearing metallocene-substituted PNP-pincer ligands is synthesized by the reduction of the corresponding monomeric molybdenum–trichloride complexes under 1 atm of molecular dinitrogen. Introduction of ferrocene as a redox-active moiety to the pyridine ring of the PNP-pincer ligand increases the catalytic activity for the formation of ammonia from molecular dinitrogen, up to 45 equiv. of ammonia being formed based on the catalyst (22 equiv. of ammonia based on each molybdenum atom of the catalyst). The time profile for the catalytic reaction reveals that the presence of the ferrocene unit in the catalyst increases the rate of ammonia formation. Electrochemical measurement and theoretical studies indicate that an interaction between the Fe atom of the ferrocene moiety and the Mo atom in the catalyst may play an important role to achieve a high catalytic activity.     

Ab Initio Molecular Dynamics Simulations of the Interaction between Organic Phosphates and Goethite

Today’s fertilizers rely heavily on mining phosphorus (P) rocks. These rocks are known to become exhausted in near future, and therefore effective P use is crucial to avoid food shortage. A substantial amount of P from fertilizers gets adsorbed onto soil minerals to become unavailable to plants. Understanding P interaction with these minerals would help efforts that improve P efficiency. To this end, we performed a molecular level analysis of the interaction of common organic P compounds (glycerolphosphate (GP) and inositol hexaphosphate (IHP)) with the abundant soil mineral (goethite) in presence of water. Molecular dynamics simulations are performed for goethite–IHP/GP–water complexes using the multiscale quantum mechanics/molecular mechanics method. Results show that GP forms monodentate (M) and bidentate mononuclear (B) motifs with B being more stable than M. IHP interacts through multiple phosphate groups with the 3M motif being most stable. The order of goethite–IHP/GP interaction energies is GP M < GP B < IHP M < IHP 3M. Water is important in these interactions as multiple proton transfers occur and hydrogen bonds are formed between goethite–IHP/GP complexes and water. We also present theoretically calculated infrared spectra which match reasonably well with frequencies reported in literature.