Including crystallographic symmetry in quantum-based refinement: Q|R#2

Quantum refinement in real and reciprocal space using the Phenix and ORCA software

X-ray and neutron crystallography, as well as cryogenic electron microscopy (cryo-EM), are the most common methods to obtain atomic structures of biological macromolecules. A feature they all have in common is that, at typical resolutions, the experimental data need to be supplemented by empirical restraints, ensuring that the final structure is chemically reasonable. The restraints are accurate for amino acids and nucleic acids, but often less accurate for substrates, inhibitors, small-molecule ligands and metal sites, for which experimental data are scarce or empirical potentials are harder to formulate. This can be solved using quantum mechanical calculations for a small but interesting part of the structure. Such an approach, called quantum refinement, has been shown to improve structures locally, allow the determination of the protonation and oxidation states of ligands and metals, and discriminate between different interpretations of the structure. Here, we present a new implementation of quantum refinement interfacing the widely used structure-refinement software Phenix and the freely available quantum mechanical software ORCA. Through application to manganese superoxide dismutase and V- and Fe-nitrogenase, we show that the approach works effectively for X-ray and neutron crystal structures, that old results can be reproduced and structural discrimination can be performed. We discuss how the weight factor between the experimental data and the empirical restraints should be selected and how quantum mechanical quality measures such as strain energies should be calculated. We also present an application of quantum refinement to cryo-EM data for particulate methane monooxygenase and show that this may be the method of choice for metal sites in such structures because no accurate empirical restraints are currently available for metals.

AQuaRef: Machine learning accelerated quantum refinement of protein structures

Cryo-EM and X-ray crystallography provide crucial experimental data for obtaining atomic-detail models of biomacromolecules. Refining these models relies on library-based stereochemical restraints, which, in addition to being limited to known chemical entities, do not include meaningful noncovalent interactions relying solely on nonbonded repulsions. Quantum mechanical (QM) calculations could alleviate these issues but are too expensive for large molecules. We present a novel AI-enabled Quantum Refinement (AQuaRef) based on AIMNet2 neural network potential mimicking QM at substantially lower computational costs. By refining 41 cryo-EM and 30 X-ray structures, we show that this approach yields atomic models with superior geometric quality compared to standard techniques, while maintaining an equal or better fit to experimental data.

Accelerating reliable multiscale quantum refinement of protein-drug systems enabled by machine learning

Biomacromolecule structures are essential for drug development and biocatalysis. Quantum refinement (QR) methods, which employ reliable quantum mechanics (QM) methods in crystallographic refinement, showed promise in improving the structural quality or even correcting the structure of biomacromolecules. However, vast computational costs and complex quantum mechanics/molecular mechanics (QM/MM) setups limit QR applications. Here we incorporate robust machine learning potentials (MLPs) in multiscale ONIOM(QM:MM) schemes to describe the core parts (e.g., drugs/inhibitors), replacing the expensive QM method. Additionally, two levels of MLPs are combined for the first time to overcome MLP limitations. Our unique MLPs+ONIOM-based QR methods achieve QM-level accuracy with significantly higher efficiency. Furthermore, our refinements provide computational evidence for the existence of bonded and nonbonded forms of the Food and Drug Administration (FDA)-approved drug nirmatrelvir in one SARS-CoV-2 main protease structure. This study highlights that powerful MLPs accelerate QRs for reliable protein-drug complexes, promote broader QR applications and provide more atomistic insights into drug development.

In situ ligand restraints from quantum-mechanical methods

In macromolecular crystallographic structure refinement, ligands present challenges for the generation of geometric restraints due to their large chemical variability, their possible novel nature and their specific interaction with the binding pocket of the protein. Quantum-mechanical approaches are useful for providing accurate ligand geometries, but can be plagued by the number of minima in flexible molecules. In an effort to avoid these issues, the Quantum Mechanical Restraints (QMR) procedure optimizes the ligand geometry in situ, thus accounting for the influence of the macromolecule on the local energy minima of the ligand. The optimized ligand geometry is used to generate target values for geometric restraints during the crystallographic refinement. As demonstrated using a sample of >2330 ligand instances in >1700 protein-ligand models, QMR restraints generally result in lower deviations from the target stereochemistry compared with conventionally generated restraints. In particular, the QMR approach provides accurate torsion restraints for ligands and other entities.

Structural and Vibrational Properties of Amino Acids from Composite Schemes and Double-Hybrid DFT: Hydrogen Bonding in Serine as a Test Case

The structures, relative stabilities, and vibrational wavenumbers of the two most stable conformers of serine, stabilized by the O-H···N, O-H···O═C and N-H···O-H intramolecular hydrogen bonds, have been evaluated by means of state-of-the-art composite schemes based on coupled-cluster (CC) theory. The so-called "cheap" composite approach (CCSD(T)/(CBS+CV)^MP2) allowed determination of accurate equilibrium structures and harmonic vibrational wavenumbers, also pointing out significant corrections beyond the CCSD(T)/cc-pVTZ level. These accurate results stand as a reference for benchmarking selected hybrid and double-hybrid, dispersion-corrected DFT functionals. B2PLYP-D3 and DSDPBEP86 in conjunction with a triple-ζ basis set have been confirmed as effective methodologies for structural and spectroscopic studies of medium-sized flexible biomolecules, also showing intramolecular hydrogen bonding. These best performing double-hybrid functionals have been employed to simulate IR spectra by means of vibrational perturbation theory, also considering hybrid CC/DFT schemes. The best overall agreement with experiment, with mean absolute error of 8 cm^-1, has been obtained by combining CCSD(T)/(CBS+CV)^MP2 harmonic wavenumbers with B2PLYP-D3/maug-cc-pVTZ anharmonic corrections. Finally, a composite scheme entirely based on CCSD(T) calculations (CCSD(T)/CBS+CV) has been employed for energetics, further confirming that serine II is the most stable conformer, also when zero-point vibrational energy corrections are included.

Combining crystallography with quantum mechanics

In standard crystallographic refinement of biomacromolecules, the crystallographic raw data are supplemented by empirical restraints that ensure that the structure makes chemical sense. These restraints are typically accurate for amino acids and nucleic acids, but less so for cofactors, substrates, inhibitors, ligands and metal sites. In quantum refinement, this potential is replaced by more accurate quantum mechanical (QM) calculations. Several implementations have been presented, differing in the level of QM and whether it is used for the entire structure or only for a site of particular interest. It has been shown that the method can improve and correct errors in crystal structures and that it can be used to determine protonation and tautomeric states of various ligands and to decide what is really seen in the structure by refining different interpretations and using standard crystallographic and QM quality measures to decide which fits the structure best.

lamaGOET: an interface for quantum crystallography

In quantum crystallography, theoretical calculations and crystallographic refinements are closely intertwined. This means that the employed software must be able to perform both quantum-mechanical calculations and crystallographic least-squares refinements. So far, the program Tonto is the only one able to do that. The lamaGOET interface described herein deals with this issue since it interfaces dedicated quantum-chemical software (the widely used Gaussian package and the specialized ELMOdb program) with the refinement capabilities of Tonto. Three different flavours of quantum-crystallographic refinements of the dipetide glycyl-l-threonine dihydrate are presented to showcase the capabilities of lamaGOET: Hirshfeld atom refinement (HAR), HAR-ELMO, namely HAR coupled with extremely localized molecular orbitals, and X-ray constrained wavefunction fitting.

Macromolecular refinement of X-ray and cryoelectron microscopy structures with Phenix/OPLS3e for improved structure and ligand quality

With the advent of the resolution revolution in cryoelectron microscopy (cryo-EM), low-resolution refinement is common, and likewise increases the need for a reliable force field. Here, we report on the incorporation of the OPLS3e force field with the VSGB2.1 solvation model in the structure determination package Phenix. Our results show significantly improved structure quality and reduced ligand strain at lower resolution for X-ray refinement. For refinement of cryo-EM-based structures, we find comparable quality structures, goodness-of-fit, and reduced ligand strain. We also show how structure quality and ligand strain are related to the map-model cross-correlation as a function of data weight, and how that can detect overfitting. Signs of overfitting are found in over half of our cryo-EM dataset, which can be remedied by a re-refinement at a lower data weight. Finally, a start-to-end script for refining structures with Phenix/OPLS3e is available in the Schrödinger 2020-3 distribution.

Real-space quantum-based refinement for cryo-EM: Q|R#3

Electron cryo-microscopy (cryo-EM) is rapidly becoming a major competitor to X-ray crystallography, especially for large structures that are difficult or impossible to crystallize. While recent spectacular technological improvements have led to significantly higher resolution three-dimensional reconstructions, the average quality of cryo-EM maps is still at the low-resolution end of the range compared with crystallography. A long-standing challenge for atomic model refinement has been the production of stereochemically meaningful models for this resolution regime. Here, it is demonstrated that including accurate model geometry restraints derived from ab initio quantum-chemical calculations (HF-D3/6-31G) can improve the refinement of an example structure (chain A of PDB entry 3j63). The robustness of the procedure is tested for additional structures with up to 7000 atoms (PDB entry 3a5x and chain C of PDB entry 5fn5) using the less expensive semi-empirical (GFN1-xTB) model. The necessary algorithms enabling real-space quantum refinement have been implemented in the latest version of qr.refine and are described here.

Quantum Chemical Calculation of Molecular and Periodic Peptide and Protein Structures

Special-purpose classical force fields (FFs) provide good accuracy at very low computational cost, but their application is limited to systems for which potential energy functions are available. This excludes most metal-containing proteins or those containing cofactors. In contrast, the GFN2-xTB semiempirical quantum chemical method is parametrized for almost the entire periodic table. The accuracy of GFN2-xTB is assessed for protein structures with respect to experimental X-ray data. Furthermore, the results are compared with those of two special-purpose FFs, HF-3c, PM6-D3H4X, and PM7. The test sets include proteins without any prosthetic groups as well as metalloproteins. Crystal packing effects are examined for a set of smaller proteins to validate the molecular approach. For the proteins without prosthetic groups, the special purpose FF OPLS-2005 yields the smallest overall RMSD to the X-ray data but GFN2-xTB provides similarly good structures with even better bond-length distributions. For the metalloproteins with up to 5000 atoms, a good overall structural agreement is obtained with GFN2-xTB. The full geometry optimizations of protein structures with on average 1000 atoms in wall-times below 1 day establishes the GFN2-xTB method as a versatile tool for the computational treatment of various biomolecules with a good accuracy/computational cost ratio.