Nature Catalysis · 2026
A tailored CoA ligase–N-acyltransferase cascade for on-DNA amide bond formation gives access to broad substrate scope
Author 5 of 10 · cited 1×
10.1038/s41929-026-01576-x
Research
I work on halohydrin dehalogenases, light-controlled enzymes and enzymes that build molecules on DNA. I simulate them with QM, QM/MM and MM, develop the methods those simulations need, and work with experimental collaborators who test the designs.
Research
Quantum mechanics for the electrons, molecular mechanics for the protein, QM/MM to join them, and the methods I develop at each level. The ruler on the left tracks length and time; the panels on the right run live in your browser. Panels marked real data play my own calculations; the others are schematics of the method.
DFT · xTB / g-xTB · ML potentials (AIQM2, MACE, UMA)
Where are the electrons, and where do they go when a bond breaks?
Semi-empirical and machine-learned potentials make it cheap enough to scan, optimise and search for transition states on hundreds of atoms, then check the key points with DFT.
Methods I builtgxtb-pymol-interfaceMLatom-Sella-TS-search
θ = 0°
QM cluster models · relaxed scans · TS search (Sella) · EVB
What does the reaction look like at the atomic level, and what does it cost?
A 303-atom cluster model of a halohydrin dehalogenase active site, computed with AIQM2 through ORCA. Cyanide is pulled towards the epoxide carbon one step at a time. At the last step the epoxide C–O bond breaks and the catalytic tyrosine hands its proton to the leaving oxygen, while the arginine re-protonates the tyrosine.
Methods I builtMLatom-ORCA-interfacehow the scan was run
Loading…
QM/MM electrostatic embedding · environment-aware RESP
How does the enzyme polarise the substrate before anything reacts?
Charges fitted in the gas phase miss the field of the protein. Fitting them inside the enzyme, with the protein as point charges, gives ligand parameters that match where the ligand actually sits.
Methods I builtQMMESPprotonation-optimizer
q(O) = −0.58 e
MD (OpenMM, Amber) · enhanced sampling · MSMs · tunnel analysis
Which conformations does the enzyme visit, and which ones are active?
Long simulations, projected onto a few collective variables, show which states an enzyme samples and how mutations shift the balance, such as the flexible N-terminal loop that controls HheG's activity and enantioselectivity.
Methods I builtMDTunnel
PaperHheG loop
Live Langevin dynamics at 300 K
AlphaFold ensembles · LigandMPNN · BindCraft · Boltz-2
Which positions should we change, and to what?
Simulations point to the residues that matter; design tools and the lab test what to put there. In HheC, the single change F12Y in sequence motif 1 raised specific activity up to 5-fold and the apparent melting temperature by 10 K.
In-silico scan
Nature Catalysis · 2026
Author 5 of 10 · cited 1×
10.1038/s41929-026-01576-x
Angew. Chem. Int. Ed. · 2026
Author 3 of 7 · cited 1×
10.1002/anie.202518896
ACS Catalysis · 2025
Author 2 of 5 · cited 13×
10.1021/acscatal.5c00819
ACS Catalysis · 2024
Author 3 of 6 · cited 10×
10.1021/acscatal.4c04815
Protein Science · 2022
Author 3 of 4 · cited 21×
10.1002/pro.4426
Angew. Chem. Int. Ed. · 2020
Author 3 of 10 · cited 103×
10.1002/anie.202008339