Changes in genome architecture and transcriptional dynamics progress independently of sensory experience during post-natal brain development.

Both transcription and three-dimensional (3D) architecture of the mammalian genome play critical roles in neurodevelopment and its disorders. However, 3D genome structures of single brain cells have not been solved; little is known about the dynamics of single-cell transcriptome and 3D genome after birth. Here, we generated a transcriptome (3,517 cells) and 3D genome (3,646 cells) atlas of the developing mouse cortex and hippocampus by using our high-resolution multiple annealing and looping-based amplification cycles for digital transcriptomics (MALBAC-DT) and diploid chromatin conformation capture (Dip-C) methods and developing multi-omic analysis pipelines. In adults, 3D genome "structure types" delineate all major cell types, with high correlation between chromatin A/B compartments and gene expression. During development, both transcriptome and 3D genome are extensively transformed in the first post-natal month. In neurons, 3D genome is rewired across scales, correlated with gene expression modules, and independent of sensory experience. Finally, we examine allele-specific structure of imprinted genes, revealing local and chromosome (chr)-wide differences. These findings uncover an unknown dimension of neurodevelopment.

Atomic resolution structure of full-length human insulin fibrils.

Patients with type 1 diabetes mellitus who are dependent on an external supply of insulin develop insulin-derived amyloidosis at the sites of insulin injection. A major component of these plaques is identified as full-length insulin consisting of the two chains A and B. While there have been several reports that characterize insulin misfolding and the biophysical properties of the fibrils, atomic-level information on the insulin fibril architecture remains elusive. We present here an atomic resolution structure of a monomorphic insulin amyloid fibril that has been determined using magic angle spinning solid-state NMR spectroscopy. The structure of the insulin monomer yields a U-shaped fold in which the two chains A and B are arranged in parallel to each other and are oriented perpendicular to the fibril axis. Each chain contains two ß-strands. We identify two hydrophobic clusters that together with the three preserved disulfide bridges define the amyloid core structure. The surface of the monomeric amyloid unit cell is hydrophobic implicating a potential dimerization and oligomerization interface for the assembly of several protofilaments in the mature fibril. The structure provides a starting point for the development of drugs that bind to the fibril surface and disrupt secondary nucleation as well as for other therapeutic approaches to attenuate insulin aggregation.

What is allosteric regulation? Exploring the exceptions that prove the rule!

"Allosteric" was first introduced to mean the other site (i.e., a site distinct from the active or orthosteric site), an adjective for "regulation" to imply a regulatory outcome resulting from ligand binding at another site. That original idea outlines a system with two ligand-binding events at two distinct locations on a macromolecule (originally a protein system), which defines a four-state energy cycle. An allosteric energy cycle provides a quantifiable allosteric coupling constant and focuses our attention on the unique properties of the four equilibrated protein complexes that constitute the energy cycle. Because many observed phenomena have been referenced as "allosteric regulation" in the literature, the goal of this work is to use literature examples to explore which systems are and are not consistent with the two-ligand thermodynamic energy cycle-based definition of allosteric regulation. We emphasize the need for consistent language so comparisons can be made among the ever-increasing number of allosteric systems. Building on the mutually exclusive natures of an energy cycle definition of allosteric regulation versus classic two-state models, we conclude our discussion by outlining how the often-proposed Rube-Goldberg-like mechanisms are likely inconsistent with an energy cycle definition of allosteric regulation.

Characterization of erythroferrone oligomerization and its impact on BMP antagonism.

Hepcidin, a peptide hormone that negatively regulates iron metabolism, is expressed by bone morphogenetic protein (BMP) signaling. Erythroferrone (ERFE) is an extracellular protein that binds and inhibits BMP ligands, thus positively regulating iron import by indirectly suppressing hepcidin. This allows for rapid erythrocyte regeneration after blood loss. ERFE belongs to the C1Q/TNF-related protein family and is suggested to adopt multiple oligomeric forms: a trimer, a hexamer, and a high molecular weight species. The molecular basis for how ERFE binds BMP ligands and how the different oligomeric states impact BMP inhibition are poorly understood. In this study, we demonstrated that ERFE activity is dependent on the presence of stable dimeric or trimeric ERFE and that larger species are dispensable for BMP inhibition. Additionally, we used an in silico approach to identify a helix, termed the ligand-binding domain, that was predicted to bind BMPs and occlude the type I receptor pocket. We provide evidence that the ligand-binding domain is crucial for activity through luciferase assays and surface plasmon resonance analysis. Our findings provide new insight into how ERFE oligomerization impacts BMP inhibition, while identifying critical molecular features of ERFE essential for binding BMP ligands.

Strategy for modeling higher-order G-quadruplex structures recalcitrant to NMR determination.

Guanine-rich nucleic acids can form intramolecularly folded four-stranded structures known as G-quadruplexes (G4s). Traditionally, G4 research has focused on short, highly modified DNA or RNA sequences that form well-defined homogeneous compact structures. However, the existence of longer sequences with multiple G4 repeats, from proto-oncogene promoters to telomeres, suggests the potential for more complex higher-order structures with multiple G4 units that might offer selective drug-targeting sites for therapeutic development. These larger structures present significant challenges for structural characterization by traditional high-resolution methods like multi-dimensional NMR and X-ray crystallography due to their molecular complexity. To address this current challenge, we have developed an integrated structural biology (ISB) platform, combining experimental and computational methods to determine self-consistent molecular models of higher-order G4s (xG4s). Here we outline our ISB method using two recent examples from our lab, an extended c-Myc promoter and long human telomere G4 repeats, that highlights the utility and generality of our approach to characterizing biologically relevant xG4s.

Membrane-bound chemoreception of bitter bile acids and peptides is mediated by the same subset of bitter taste receptors.

The vertebrate sense of taste allows rapid assessment of the nutritional quality and potential presence of harmful substances prior to ingestion. Among the five basic taste qualities, salty, sour, sweet, umami, and bitter, bitterness is associated with the presence of putative toxic substances and elicits rejection behaviors in a wide range of animals including humans. However, not all bitter substances are harmful, some are thought to be health-beneficial and nutritious. Among those compound classes that elicit a bitter taste although being non-toxic and partly even essential for humans are bitter peptides and L-amino acids. Using functional heterologous expression assays, we observed that the 5 dominant human bitter taste receptors responsive to bitter peptides and amino acids are activated by bile acids, which are notorious for their extreme bitterness. We further demonstrate that the cross-reactivity of bitter taste receptors for these two different compound classes is evolutionary conserved and can be traced back to the amphibian lineage. Moreover, we show that the cross-detection by some receptors relies on "structural mimicry" between the very bitter peptide L-Trp-Trp-Trp and bile acids, whereas other receptors exhibit a phylogenetic conservation of this trait. As some bile acid-sensitive bitter taste receptor genes fulfill dual-roles in gustatory and non-gustatory systems, we suggest that the phylogenetic conservation of the rather surprising cross-detection of the two substance classes could rely on a gene-sharing-like mechanism in which the non-gustatory function accounts for the bitter taste response to amino acids and peptides.

Ultrafast photooxidation of protein-bound anionic flavin radicals.

The photophysical properties of anionic semireduced flavin radicals are largely unknown despite their importance in numerous biochemical reactions. Here, we studied the photoproducts of these intrinsically unstable species in five different flavoprotein oxidases where they can be stabilized, including the well-characterized glucose oxidase. Using ultrafast absorption and fluorescence spectroscopy, we unexpectedly found that photoexcitation systematically results in the oxidation of protein-bound anionic flavin radicals on a time scale of less than ∼100 fs. The thus generated photoproducts decay back in the remarkably narrow 10- to 20-ps time range. Based on molecular dynamics and quantum mechanics computations, positively charged active-site histidine and arginine residues are proposed to be the electron acceptor candidates. Altogether, we established that, in addition to the commonly known and extensively studied photoreduction of oxidized flavins in flavoproteins, the reverse process (i.e., the photooxidation of anionic flavin radicals) can also occur. We propose that this process may constitute an excited-state deactivation pathway for protein-bound anionic flavin radicals in general. This hitherto undocumented photochemical reaction in flavoproteins further extends the family of flavin photocycles.

The emergence of machine learning force fields in drug design.

In the field of molecular simulation for drug design, traditional molecular mechanic force fields and quantum chemical theories have been instrumental but limited in terms of scalability and computational efficiency. To overcome these limitations, machine learning force fields (MLFFs) have emerged as a powerful tool capable of balancing accuracy with efficiency. MLFFs rely on the relationship between molecular structures and potential energy, bypassing the need for a preconceived notion of interaction representations. Their accuracy depends on the machine learning models used, and the quality and volume of training data sets. With recent advances in equivariant neural networks and high-quality datasets, MLFFs have significantly improved their performance. This review explores MLFFs, emphasizing their potential in drug design. It elucidates MLFF principles, provides development and validation guidelines, and highlights successful MLFF implementations. It also addresses potential challenges in developing and applying MLFFs. The review concludes by illuminating the path ahead for MLFFs, outlining the challenges to be overcome and the opportunities to be harnessed. This inspires researchers to embrace MLFFs in their investigations as a new tool to perform molecular simulations in drug design.

Characterization of erythroferrone structural domains relevant to its iron-regulatory function.

Iron delivery to the plasma is closely coupled to erythropoiesis, the production of red blood cells, as this process consumes most of the circulating plasma iron. In response to hemorrhage and other erythropoietic stresses, increased erythropoietin stimulates the production of the hormone erythroferrone (ERFE) by erythrocyte precursors (erythroblasts) developing in erythropoietic tissues. ERFE acts on the liver to inhibit bone morphogenetic protein (BMP) signaling and thereby decrease hepcidin production. Decreased circulating hepcidin concentrations then allow the release of iron from stores and increase iron absorption from the diet. Guided by evolutionary analysis and Alphafold2 protein complex modeling, we used targeted ERFE mutations, deletions, and synthetic ERFE segments together with cell-based bioassays and surface plasmon resonance to probe the structural features required for bioactivity and BMP binding. We define the ERFE active domain and multiple structural features that act together to entrap BMP ligands. In particular, the hydrophobic helical segment 81 to 86 and specifically the highly conserved tryptophan W82 in the N-terminal region are essential for ERFE bioactivity and Alphafold2 modeling places W82 between two tryptophans in its ligands BMP2, BMP6, and the BMP2/6 heterodimer, an interaction similar to those that bind BMPs to their cognate receptors. Finally, we identify the cationic region 96-107 and the globular TNFα-like domain 186-354 as structural determinants of ERFE multimerization that increase the avidity of ERFE for BMP ligands. Collectively, our results provide further insight into the ERFE-mediated inhibition of BMP signaling in response to erythropoietic stress.

Deciphering specificity and cross-reactivity in tachykinin NK1 and NK2 receptors.

The tachykinin receptors neurokinin 1 (NK1R) and neurokinin 2 (NK2R) are G protein-coupled receptors that bind preferentially to the natural peptide ligands substance P and neurokinin A, respectively, and have been targets for drug development. Despite sharing a common C-terminal sequence of Phe-X-Gly-Leu-Met-NH2 that helps direct biological function, the peptide ligands exhibit some degree of cross-reactivity toward each other's non-natural receptor. Here, we investigate the detailed structure-activity relationships of the ligand-bound receptor complexes that underlie both potent activation by the natural ligand and cross-reactivity. We find that the specificity and cross-reactivity of the peptide ligands can be explained by the interactions between the amino acids preceding the FxGLM consensus motif of the bound peptide ligand and two regions of the receptor: the ß-hairpin of the extracellular loop 2 (ECL2) and a N-terminal segment leading into transmembrane helix 1. Positively charged sidechains of the ECL2 (R177 of NK1R and K180 of NK2R) are seen to play a vital role in the interaction. The N-terminal positions 1 to 3 of the peptide ligand are entirely dispensable. Mutated and chimeric receptor and ligand constructs neatly swap around ligand specificity as expected, validating the structure-activity hypotheses presented. These findings will help in developing improved agonists or antagonists for NK1R and NK2R.

Barley Nodulin 26-like intrinsic protein permeates water, metalloids, saccharides, and ion pairs due to structural plasticity and diversification.

Aquaporins can facilitate the passive movement of water, small polar molecules, and some ions. Here, we examined solute selectivity for the barley Nodulin 26-like Intrinsic Protein (HvNIP2;1) embedded in liposomes and examined through stopped-flow light scattering spectrophotometry and Xenopus laevis oocyte swelling assays. We found that HvNIP2;1 permeates water, boric and germanic acids, sucrose, and lactose but not d-glucose or d-fructose. Other saccharides, such as neutral (d-mannose, d-galactose, d-xylose, d-mannoheptaose) and charged (N-acetyl d-glucosamine, d-glucosamine, d-glucuronic acid) aldoses, disaccharides (cellobiose, gentiobiose, trehalose), trisaccharide raffinose, and urea, glycerol, and acyclic polyols, were permeated to a much lower extent. We observed apparent permeation of hydrated KCl and MgSO4 ions, while CH3COONa and NaNO3 permeated at significantly lower rates. Our experiments with boric acid and sucrose revealed no apparent interaction between solutes when permeated together, and AgNO3 or H[AuCl4] blocked the permeation of all solutes. Docking of sucrose in HvNIP2;1 and spinach water-selective SoPIP2;1 aquaporins revealed the structural basis for sucrose permeation in HvNIP2;1 but not in SoPIP2;1, and defined key residues interacting with this permeant. In a biological context, sucrose transport could constitute a novel element of plant saccharide-transporting machinery. Phylogenomic analyses of 164 Viridiplantae and 2993 Archaean, bacterial, fungal, and Metazoan aquaporins rationalized solute poly-selectivity in NIP3 sub-clade entries and suggested that they diversified from other sub-clades to acquire a unique specificity of saccharide transporters. Solute specificity definition in NIP aquaporins could inspire developing plants for food production.

Apo-metallothionein-3 cooperatively forms tightly compact structures under physiological conditions.

Metallothioneins (MTs) are essential mammalian metal chaperones. MT isoform 1 (MT1) is expressed in the kidneys and isoform 3 (MT3) is expressed in nervous tissue. For MTs, the solution-based NMR structure was determined for metal-bound MT1 and MT2, and only one X-ray diffraction structure on a crystallized mixed metal-bound MT2 has been reported. The structure of solution-based metalated MT3 is partially known using NMR methods; however, little is known about the fluxional de novo apo-MT3 because the structure cannot be determined by traditional methods. Here, we used cysteine modification coupled with electrospray ionization mass spectrometry, denaturing reactions with guanidinium chloride, stopped-flow methods measuring cysteine modification and metalation, and ion mobility mass spectrometry to reveal that apo-MT3 adopts a compact structure under physiological conditions and an extended structure under denaturing conditions, with no intermediates. Compared with apo-MT1, we found that this compact apo-MT3 binds to a cysteine modifier more cooperatively at equilibrium and 0.5 times the rate, providing quantitative evidence that many of the 20 cysteines of apo-MT3 are less accessible than those of apo-MT1. In addition, this compact apo-MT3 can be identified as a distinct population using ion mobility mass spectrometry. Furthermore, proposed structural models can be calculated using molecular dynamics methods. Collectively, these findings provide support for MT3 acting as a noninducible regulator of the nervous system compared with MT1 as an inducible scavenger of trace metals and toxic metals in the kidneys.

Dissecting the geometric and hydrophobic constraints of stapled peptides.

Stapled peptides are a promising class of molecules with potential as highly specific probes of protein-protein interactions and as therapeutics. Hydrocarbon stapling affects the peptide properties through the interplay of two factors: enhancing the overall hydrophobicity and constraining the conformational flexibility. By constructing a series of virtual peptides, we study the role of each factor in modulating the structural properties of a hydrocarbon-stapled peptide PM2, which has been shown to enter cells, engage its target Mouse Double Minute 2 (MDM2), and activate p53. Hamiltonian replica exchange molecular dynamics (HREMD) simulations suggest that hydrocarbon stapling favors helical populations of PM2 through a combination of the geometric constraints and the enhanced hydrophobicity of the peptide. To further understand the conformational landscape of the stapled peptides along the binding pathway, we performed HREMD simulations by restraining the peptide at different distances from MDM2. When the peptide approaches MDM2, the binding pocket undergoes dehydration which appears to be greater in the presence of the stapled peptide compared with the linear peptide. In the binding pocket, the helicity of the stapled peptide is increased due to the favorable interactions between the peptide residues as well as the staple and the microenvironment of the binding pocket, contributing to enhanced affinity. The dissection of the multifaceted mechanism of hydrocarbon stapling into individual factors not only deepens fundamental understanding of peptide stapling, but also provides guidelines for the design of new stapled peptides.

Validating Fractalkine receptor as a target and identifying candidates for drug discovery against type 2 diabetes.

Type 2 diabetes mellitus (T2DM) is one of the most common chronic diseases employing abnormal levels of insulin. Enhancing the insulin production is greatly aided by the regulatory mechanisms of the Fractalkine receptor (CX3CR1) system in islet ß-cell function. However, elements including a high-fat diet, obesity, and ageing negatively impact the expression of CX3CR1 in islets. CX3CL1/CX3CR1 receptor-ligand complex is now recognized as a novel therapeutic target. It suggests that T2DM-related ß-cell dysfunction may result from lower amount of these proteins. We analyzed the differential expression of CX3CR1 gene samples taken from persons with T2DM using data obtained from the Gene Expression Omnibus database. Homology modeling enabled us to generate the three-dimensional structure of CX3CR1 and a possible binding pocket. The optimized CX3CR1 structure was subjected to rigorous screening against a massive library of 693 million drug-like molecules from the ZINC15 database. This screening process led to the identification of three compounds with strong binding affinity at the identified binding pocket of CX3CR1. To further evaluate the potential of these compounds, molecular dynamics simulations were conducted over a 50 ns time scale to assess the stability of the protein-ligand complexes. These simulations revealed that ZINC000032506419 emerged as the most promising drug-like compound among the three potent molecules. The discovery of ZINC000032506419 holds exciting promise as a potential therapeutic agent for T2D and other related metabolic disorders. These findings pave the way for the development of effective medications to address the complexities of T2DM and its associated metabolic diseases.

Using Multiscale Molecular Modeling to Analyze Possible NS2b-NS3 Protease Inhibitors from Philippine Medicinal Plants.

Within the field of Philippine folkloric medicine, the utilization of indigenous plants like Euphorbia hirta (tawa-tawa), Carica papaya (papaya), and Psidium guajava (guava) as potential dengue remedies has gained attention. Yet, limited research exists on their comprehensive effects, particularly their anti-dengue activity. This study screened 2944 phytochemicals from various Philippine plants for anti-dengue activity. Absorption, distribution, metabolism, excretion, and toxicity (ADMET) profiling provided 1265 compounds demonstrating pharmacokinetic profiles suitable for human use. Molecular docking targeting the dengue virus NS2b-NS3 protease's catalytic triad (Asp 75, Ser 135, and His 51) identified ten ligands with higher docking scores than reference compounds idelalisib and nintedanib. Molecular dynamics simulations confirmed the stability of eight of these ligand-protease complexes. Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) analysis highlighted six ligands, including veramiline (-80.682 kJ/mol), cyclobranol (-70.943 kJ/mol), chlorogenin (-63.279 kJ/mol), 25beta-Hydroxyverazine (-61.951 kJ/mol), etiolin (-59.923 kJ/mol), and ecliptalbine (-56.932 kJ/mol) with favorable binding energies, high oral bioavailability, and drug-like properties. This integration of traditional medical knowledge with advanced computational drug discovery methods paves new pathways for the development of treatments for dengue.

Ground and excited state charge transfer at aqueous nanodiamonds.

Nanodiamonds (NDs) are unique carbonaceous materials with exceptionally high stability, hardness, and notable electronic properties. Their applications in photocatalysis, biomedicine, and energy materials are usually carried out in aqueous environments, where they interact with aqueous adsorbates. Especially, electron density may rearrange from the diamond material toward oxidative adsorbates such as oxygen, which is known as charge transfer doping. In this article, we quantify the charge transfer doping for NDs with inhomogeneous surface coverings (hydroxyl, fluorine, and amorphous carbon), as well as NDs doped with heteroatoms (B, Si, N) using hybrid density functional theory (DFT) calculations. The transfer doping magnitude is largely determined by the NDs' highest occupied molecular orbital energies, which can in turn be modified by the surface covering and doping. However, local modifications of the ND structures do not have any local effects on the magnitude of the charge transfer. We furthermore analyze the impact of aqueous adsorbates on the excited states of an aqueous ND in the context of photocatalysis via time-dependent DFT. Here, we find that the excited electrons are biased to move in the direction of the respective oxidative adsorbate. Surprisingly, we find that also unreactive species such as nitrous oxide may attract the excited electrons, which is probably due to the positive partial charge that is induced by the local N 2 O solvation geometry.

Structure-based improvement of the binding affinity and recognition specificity of peptide competitors to target pediatric IL-5R/IL-5 interaction by gluing halogen bonds at their complex interface.

Human interleukin-5 (IL-5) cytokine mediates the development of eosinophils and is involved in a variety of immune inflammatory responses that play a major role in the pathogenesis of childhood asthma, leukemia, and other pediatric allergic diseases. The immunomodulatory cytokine functions by binding to its cognate cell surface receptor IL-5R in a sheet-by-sheet manner, which can be conformationally mimicked and competitively disrupted by a double-stranded cyclic AF18748 peptide. In this study, we systematically examined the co-crystallized complex structure of human IL-5R with AF18748 peptide and rationally designed a halogen bond to glue at the protein-peptide complex interface by substituting the indole moiety of AF18748 Trp13 residue with a halogen atom (X = F, Cl, Br, or I). High-level theoretical calculations imparted presence of the halogen bond between the oxygen atom (O) of IL-5R Glu58 backbone and the halogen atom (X) of AF18748 Trp13 side chain. Experimental assays confirmed that the halogen bond can promote peptide binding moderately or considerably. More importantly, the halogen bond not only enhances peptide affinity to IL-5R, but also improves peptide selectivity for its cognate IL-5R over other noncognate IL-R proteins. As might be expected, the affinity and selectivity conferred by halogen bond increase consistently in the order: H < F < Cl < Br < I. Structural modeling revealed that the halogen bond plus its vicinal π-cation-π stacking co-define a ringed noncovalent system at the complex interface, which involves a synergistic effect to effectively improve the peptide binding potency and recognition specificity.

PKLR mutations in pyruvate kinase deficient Polish patients: Functional characteristics of c.101-1G > A and c.1058delAAG variants.

Pyruvate kinase (PK) deficiency is a rare autosomal recessive disorder characterized by chronic hemolytic anemia of variable severity. Nine Polish patients with severe hemolytic anemia but normal PK activity were found to carry mutations in the PKLR gene encoding PK, five already known ones and one novel (c.178C > T). We characterized two of the known variants by molecular modeling (c.1058delAAG) and minigene splicing analysis (c.101-1G > A). The former gives a partially destabilized PK tetramer, likely of suboptimal activity, and the c.101-1G > A variant gives alternatively spliced mRNA carrying a premature stop codon, encoding a severely truncated PK and likely undergoing nonsense-mediated decay.

From One-Dimensional Disordered Racemate to Ordered Racemic Conglomerates through Metal-Coordination-Driven Self-Assembly at the Liquid-Solid Interface.

In recent years, there has been significant focus on investigating and controlling chiral self-assembly, specifically in the context of enantiomeric separation. This study explores the self-assembly behavior of 4-dodecyl-3,6-di(2-pyridyl)pyridazine (DPP-C12) at the interface between heptanoic acid (HA) and highly oriented pyrolytic graphite (HOPG) using a combination of scanning tunneling microscopy (STM) and multiscale molecular modeling. The self-assembled monolayer structure formed by DPP-C12 is periodic in one direction, but aperiodic in the direction orthogonal to it. These structures resemble 1D disordered racemic compounds. Upon introducing palladium [Pd(II)] ions, complexing with DPP-C12, these 1D disordered racemic compounds spontaneously transform into 2D racemic conglomerates, which is rationalized with the assistance of force-field simulations. Our findings provide insights into the regulation of two-dimensional chirality.

Rationalization of a Streptavidin Based Enantioselective Artificial Suzukiase: An Integrative Computational Approach.

An Artificial Metalloenzyme (ArM) built employing the streptavidin-biotin technology has been used for the enantioselective synthesis of binaphthyls by means of asymmetric Suzuki-Miyaura cross-coupling reactions. Despite its success, it remains a challenge to understand how the length of the biotin cofactors or the introduction of mutations to streptavidin leads the preferential synthesis of one atropisomer over the other. In this study, we apply an integrated computational modeling approach, including DFT calculations, protein-ligand dockings and molecular dynamics to rationalize the impact of mutations and length of the biotion cofactor on the enantioselectivities of the biaryl product. The results unravel that the enantiomeric differences found experimentally can be rationalized by the disposition of the first intermediate, coming from the oxidative addition step, and the entrance of the second substrate. The work also showcases the difficulties facing to control the enantioselection when engineering ArM to catalyze enantioselective Suzuki-Miyaura couplings and how the combination of DFT calculations, molecular dockings and MD simulations can be used to rationalize artificial metalloenzymes.