Machine Learning-Enhanced Quantum Chemistry-Assisted Refinement of the Active Site Structure of Metalloproteins.

Understanding the fine structural details of inhibitor binding at the active site of metalloenzymes can have a profound impact on the rational drug design targeted to this broad class of biomolecules. Structural techniques such as NMR, cryo-EM, and X-ray crystallography can provide bond lengths and angles, but the uncertainties in these measurements can be as large as the range of values that have been observed for these quantities in all the published structures. This uncertainty is far too large to allow for reliable calculations at the quantum chemical (QC) levels for developing precise structure-activity relationships or for improving the energetic considerations in protein-inhibitor studies. Therefore, the need arises to rely upon computational methods to refine the active site structures well beyond the resolution obtained with routine application of structural methods. In a recent paper, we have shown that it is possible to refine the active site of cobalt(II)-substituted MMP12, a metalloprotein that is a relevant drug target, by matching to the experimental pseudocontact shifts (PCS) those calculated using multireference ab initio QC methods. The computational cost of this methodology becomes a significant bottleneck when the starting structure is not sufficiently close to the final one, which is often the case with biomolecular structures. To tackle this problem, we have developed an approach based on a neural network (NN) and a support vector regression (SVR) and applied it to the refinement of the active site structure of oxalate-inhibited human carbonic anhydrase 2 (hCAII), another prototypical metalloprotein target. The refined structure gives a remarkably good agreement between the QC-calculated and the experimental PCS. This study not only contributes to the knowledge of CAII but also demonstrates the utility of combining machine learning (ML) algorithms with QC calculations, offering a promising avenue for investigating other drug targets and complex biological systems in general.

Molecular Engineering of Self-Immolative Bioresponsive MR Probes.

Real-time detection of bio-event in whole animals provides essential information for understanding biological and therapeutic processes. Magnetic resonance (MR) imaging represents a non-invasive approach to generating three-dimensional anatomic images with high spatial-temporal resolution and unlimited depth penetration. We have developed several self-immolative enzyme-activatable agents that provide excellent in vivo contrast and function as gene expression reporters. Here, we describe a vast improvement in image contrast over our previous generations of these bioresponsive agents based on a new pyridyl-carbamate Gd(III) complex. The pyridyl-carbamate-based agent has a very low MR relaxivity in the "off-state" (r1 = 1.8 mM-1 s-1 at 1.41 T). However, upon enzymatic processing, it generates a significantly higher relaxivity with a Δr1 = 106% versus Δr1 ∼ 20% reported previously. Single X-ray crystal and nuclear magnetic relaxation dispersion analyses offer mechanistic insights regarding MR signal enhancement at the molecular scale. This work demonstrates a pyridyl-carbamate-based self-immolative molecular platform for the construction of enzymatic bio-responsive MR agents, which can be adapted to a wide range of other targets for exploring stimuli-responsive materials and biomedical applications.

Novel Nanogels Loaded with Mn(II) Chelates as Effective and Biologically Stable MRI Probes.

Here it is described nanogels (NG) based on a chitosan matrix, which are covalently stabilized by a bisamide derivative of Mn-t-CDTA (t-CDTA = trans-1,2-diaminocyclohexane-N,N,N',N'-tetraacetic acid). the Mn(II) complex acts both as a contrast medium and as a cross-linking agent. These nanogels are proposed as an alternative to the less stable paramagnetic nanogels obtained by electrostatic interactions between the polymeric matrix and paramagnetic Gd(III) chelates. The present novel nanogels show: i) relaxivity values seven times higher than that of typical monohydrated Mn(II) chelates at the clinical fields, thanks to the combination of a restricted mobility of the complex with a fast exchange of the metal-bound water molecule; ii) high stability of the formulation over time at pH 5 and under physiological conditions, thus excluding metal leaking or particles aggregation; iii) good extravasation and accumulation, with a maximum contrast achieved at 24 h post-injection in mice bearing subcutaneous breast cancer tumor; iv) high T1 contrast (1 T) in the tumor 24 h post-injection. These improved properties pave the way for the use of these paramagnetic nanogels as promising magnetic resonance imaging (MRI) probes for in vitro and in vivo preclinical applications.

Combining Solid-State NMR with Structural and Biophysical Techniques to Design Challenging Protein-Drug Conjugates.

Several protein-drug conjugates are currently being used in cancer therapy. These conjugates rely on cytotoxic organic compounds that are covalently attached to the carrier proteins or that interact with them via non-covalent interactions. Human transthyretin (TTR), a physiological protein, has already been identified as a possible carrier protein for the delivery of cytotoxic drugs. Here we show the structure-guided development of a new stable cytotoxic molecule based on a known strong binder of TTR and a well-established anticancer drug. This example is used to demonstrate the importance of the integration of multiple biophysical and structural techniques, encompassing microscale thermophoresis, X-ray crystallography and NMR. In particular, we show that solid-state NMR has the ability to reveal effects caused by ligand binding which are more easily relatable to structural and dynamical alterations that impact the stability of macromolecular complexes.

Large Protein Assemblies for High-Relaxivity Contrast Agents: The Case of Gadolinium-Labeled Asparaginase.

Biologics are emerging as the most important class of drugs and are used to treat a large variety of pathologies. Most of biologics are proteins administered in large amounts, either by intramuscular injection or by intravenous infusion. Asparaginase is a large tetrameric protein assembly, currently used against acute lymphoblastic leukemia. Here, a gadolinium(III)-DOTA derivative has been conjugated to asparaginase, and its relaxation properties have been investigated to assess its efficiency as a possible theranostic agent. The field-dependent 1H longitudinal relaxation measurements of water solutions of gadolinium(III)-labeled asparaginase indicate a very large increase in the relaxivity of this paramagnetic protein complex with respect to small gadolinium chelates, opening up the possibility of its use as an MRI contrast agent.

Not only manganese, but fruit component effects dictate the efficiency of fruit juice as an oral magnetic resonance imaging contrast agent.

Several fruit juices are used as oral contrast agents to improve the quality of images in magnetic resonance cholangiopancreatography. They are often preferred to conventional synthetic contrast agents because of their very low cost, natural origin, intrinsic safety, and comparable image qualities. Pineapple and blueberry juices are the most employed in clinical practice due to their higher content of manganese(II) ions. The interest of pharmaceutical companies in these products is testified by the appearance in the market of fruit juice derivatives with improved contrast efficacy. Here, we investigate the origin of the contrast of blueberry juice, analyze the parameters that can effect it, and elucidate the differences with pineapple juice and manganese(II) solutions. It appears that, although manganese(II) is the paramagnetic ion responsible for the contrast, it is the interaction of manganese(II) with other juice components that modulates the efficiency of the juice as a magnetic resonance contrast agent. On these grounds, we conclude that blueberry juice concentrated to the same manganese concentration of pineapple juice would prove a more efficient contrast agent than pineapple juice.

High Relaxivity with No Coordinated Waters: A Seemingly Paradoxical Behavior of [Gd(DOTP)]5- Embedded in Nanogels.

Nanogels (NGs) obtained by electrostatic interactions between chitosan and hyaluronic acid and comprising paramagnetic Gd chelates are gaining increasing attention for their potential application in magnetic resonance bioimaging. Herein, the macrocyclic complexes [Gd(DOTP)]5-, lacking metal-bound water molecules (q = 0), were confined or used as a cross-linker in this type of NG. Unlike the typical behavior of Gd complexes with q = 0, a remarkable relaxivity value of 78.0 mM-1 s-1 was measured at 20 MHz and 298 K, nearly 20 times greater than that found for the free complex. A careful analysis of the relaxation data emphasizes the fundamental role of second sphere water molecules with strong and long-lived hydrogen bonding interactions with the complex. Finally, PEGylated derivatives of nanoparticles were used for the first in vivo magnetic resonance imaging study of this type of NG, revealing a fast renal excretion of paramagnetic complexes after their release from the NGs.

The evolution of paramagnetic NMR as a tool in structural biology.

Paramagnetic NMR data contain extremely accurate long-range information on metalloprotein structures and, when used in the frame of integrative structural biology approaches, they allow for the retrieval of structural details to a resolution that is not achievable using other techniques. Paramagnetic data thus represent an extremely powerful tool to refine protein models in solution, especially when coupled to X-ray or cryoelectron microscopy data, to monitor the formation of complexes and determine the relative arrangements of their components, and to highlight the presence of conformational heterogeneity. More recently, theoretical and computational advancements in quantum chemical calculations of paramagnetic NMR observables are progressively opening new routes in structural biology, because they allow for the determination of the structure within the coordination sphere of the metal center, thus acting as a loupe on sites that are difficult to observe but very important for protein function.

Theoretical analysis of the long-distance limit of NMR chemical shieldings.

After some years of controversy, it was recently demonstrated how to obtain the correct long-distance limit [point-dipole approximation (PDA)] of pseudo-contact nuclear magnetic resonance chemical shifts from rigorous first-principles quantum mechanics [Lang et al., J. Phys. Chem. Lett. 11, 8735 (2020)]. This result confirmed the classical Kurland-McGarvey theory. In the present contribution, we elaborate on these results. In particular, we provide a detailed derivation of the PDA both from the Van den Heuvel-Soncini equation for the chemical shielding tensor and from a spin Hamiltonian approximation. Furthermore, we discuss in detail the PDA within the approximate density functional theory and Hartree-Fock theories. In our previous work, we assumed a relatively crude effective nuclear charge approximation for the spin-orbit coupling operator. Here, we overcome this assumption by demonstrating that the derivation is also possible within the fully relativistic Dirac equation and even without the assumption of a specific form for the Hamiltonian. Crucial ingredients for the general derivation are a Hamiltonian that respects gauge invariance, the multipolar gauge, and functional derivatives of the Hamiltonian, where it is possible to identify the first functional derivative with the electron number current density operator. The present work forms an important foundation for future extensions of the Kurland-McGarvey theory beyond the PDA, including induced magnetic quadrupole and higher moments to describe the magnetic hyperfine field.

Solid-state NMR methods for the characterization of bioconjugations and protein-material interactions.

Protein solid-state NMR has evolved dramatically over the last two decades, with the development of new hardware and sample preparation methodologies. This technique is now ripe for complex applications, among which one can count bioconjugation, protein chemistry and functional biomaterials. In this review, we provide our account on this aspect of protein solid-state NMR.

Detection of Metabolite-Protein Interactions in Complex Biological Samples by High-Resolution Relaxometry: Toward Interactomics by NMR.

Metabolomics, the systematic investigation of metabolites in biological fluids, cells, or tissues, reveals essential information about metabolism and diseases. Metabolites have functional roles in a myriad of biological processes, as substrates and products of enzymatic reactions but also as cofactors and regulators of large numbers of biochemical mechanisms. These functions involve interactions of metabolites with macromolecules. Yet, methods to systematically investigate these interactions are still scarce to date. In particular, there is a need for techniques suited to identify and characterize weak metabolite-macromolecule interactions directly in complex media such as biological fluids. Here, we introduce a method to investigate weak interactions between metabolites and macromolecules in biological fluids. Our approach is based on high-resolution NMR relaxometry and does not require any invasive procedure or separation step. We show that we can detect interactions between small and large molecules in human blood serum and quantify the size of the complex. Our work opens the way for investigations of metabolite (or other small molecules)-protein interactions in biological fluids for interactomics or pharmaceutical applications.

Comparison of Different Reweighting Approaches for the Calculation of Conformational Variability of Macromolecules from Molecular Simulations.

Conformational variability and heterogeneity are crucial determinants of the function of biological macromolecules. The possibility of accessing this information experimentally suffers from severe under-determination of the problem, since there are a few experimental observables to be accounted for by a (potentially) infinite number of available conformational states. Several computational methods have been proposed over the years in order to circumvent this theoretically insurmountable obstacle. A large share of these strategies is based on reweighting an initial conformational ensemble which arises from, for example, molecular simulations of different qualities and levels of theory. In this work, we compare the outcome of three reweighting approaches based on radically different views of the conformational heterogeneity problem, namely Maximum Entropy, Maximum Parsimony and Maximum Occurrence, and we do so using the same experimental data. In this comparison we find both expected as well as unexpected similarities.

A Quantum Chemistry View on Two Archetypical Paramagnetic Pentacoordinate Nickel(II) Complexes Offers a Fresh Look on Their NMR Spectra.

Quantum chemical methods for calculating paramagnetic NMR observables are becoming increasingly accessible and are being included in the inorganic chemistry practice. Here, we test the performance of these methods in the prediction of proton hyperfine shifts of two archetypical high-spin pentacoordinate nickel(II) complexes (NiSAL-MeDPT and NiSAL-HDPT), which, for a variety of reasons, turned out to be perfectly suited to challenge the predictions to the finest level of detail. For NiSAL-MeDPT, new NMR experiments yield an assignment that perfectly matches the calculations. The slightly different hyperfine shifts from the two "halves" of the molecules related by a pseudo-C2 axis, which are experimentally divided into two well-defined spin systems, are also straightforwardly distinguished by the calculations. In the case of NiSAL-HDPT, for which no X-ray structure is available, the quality of the calculations allowed us to refine its structure using as a starting template the structure of NiSAL-MeDPT.

A High-Resolution View of the Coordination Environment in a Paramagnetic Metalloprotein from its Magnetic Properties.

Metalloproteins constitute a significant fraction of the proteome of all organisms and their characterization is critical for both basic sciences and biomedical applications. A large portion of metalloproteins bind paramagnetic metal ions, and paramagnetic NMR spectroscopy has been widely used in their structural characterization. However, the signals of nuclei in the immediate vicinity of the metal center are often broadened beyond detection. In this work, we show that it is possible to determine the coordination environment of the paramagnetic metal in the protein at a resolution inaccessible to other techniques. Taking the structure of a diamagnetic analogue as a starting point, a geometry optimization is carried out by fitting the pseudocontact shifts obtained from first principles quantum chemical calculations to the experimental ones.

Single Peptide Backbone Surrogate Mutations to Regulate Angiotensin GPCR Subtype Selectivity.

Mutating the side-chains of amino acids in a peptide ligand, with unnatural amino acids, aiming to mitigate its short half-life is an established approach. However, it is hypothesized that mutating specific backbone peptide bonds with bioisosters can be exploited not only to enhance the proteolytic stability of parent peptides, but also to tune its receptor subtype selectivity. Towards this end, four [Y]6 -Angiotensin II analogues are synthesized where amide bonds have been replaced by 1,4-disubstituted 1,2,3-triazole isosteres in four different backbone locations. All the analogues possessed enhanced stability in human plasma in comparison with the parent peptide, whereas only two of them achieved enhanced AT2 R/AT1 R subtype selectivity. This diversification has been studied through 2D NMR spectroscopy and unveiled a putative more structured microenvironment for the two selective ligands accompanied with increased number of NOE cross-peaks. The most potent analogue, compound 2, has been explored regarding its neurotrophic potential and resulted in an enhanced neurite growth with respect to the established agent C21.

The Photocatalyzed Thiol-ene reaction: A New Tag to Yield Fast, Selective and reversible Paramagnetic Tagging of Proteins.

Paramagnetic restraints have been used in biomolecular NMR for the last three decades to elucidate and refine biomolecular structures, but also to characterize protein-ligand interactions. A common technique to generate such restraints in proteins, which do not naturally contain a (paramagnetic) metal, consists in the attachment to the protein of a lanthanide-binding-tag (LBT). In order to design such LBTs, it is important to consider the efficiency and stability of the conjugation, the geometry of the complex (conformational exchanges and coordination) and the chemical inertness of the ligand. Here we describe a photo-catalyzed thiol-ene reaction for the cysteine-selective paramagnetic tagging of proteins. As a model, we designed an LBT with a vinyl-pyridine moiety which was used to attach our tag to the protein GB1 in fast and irreversible fashion. Our tag T1 yields magnetic susceptibility tensors of significant size with different lanthanides and has been characterized using NMR and relaxometry measurements.

Identification of productive and futile encounters in an electron transfer protein complex.

Well-defined, stereospecific states in protein complexes are often in exchange with an ensemble of more dynamic orientations: the encounter states. The structure of the stereospecific complex between cytochrome P450cam and putidaredoxin was solved recently by X-ray diffraction as well as paramagnetic NMR spectroscopy. Other than the stereospecific complex, the NMR data clearly show the presence of additional states in the complex in solution. In these encounter states, populated for a small percentage of the time, putidaredoxin assumes multiple orientations and samples a large part of the surface of cytochrome P450cam. To characterize the nature of the encounter states, an extensive paramagnetic NMR dataset has been analyzed using the Maximum Occurrence of Regions methodology. The analysis reveals the location and maximal spatial extent of the additional states needed to fully explain the NMR data. Under the assumption of sparsity of the size of the conformational ensemble, several minor states can be located quite precisely. The distribution of these minor states correlates with the electrostatic potential map around cytochrome P450cam. Whereas some minor states are on isolated positively charged patches, others are connected to the stereospecific site via positively charged paths. The existence of electrostatically favorable pathways between the stereospecific interaction site and the different minor states or lack thereof suggests a means to discriminate between productive and futile encounter states.

Assessing Structural Preferences of Unstructured Protein Regions by NMR.

Biomacromolecules, such as proteins, often exhibit significant motions intimately associated with their function. Intrinsically disordered proteins and proteins with intrinsically disordered regions, although extremely important for a plethora of cellular functions, are difficult to structurally characterize at the atomic level because the experimental parameters report on ensemble and time averages. Here, we demonstrate for the C-terminal domain of the human immunodeficiency virus type 1 capsid protein that NMR and, in particular, residual dipolar couplings (RDCs) measured for the folded portion of the protein can inform on the structural preferences of the unstructured portion using RDC-prediction tools and the maximum occurrence approach.

Metal centers in biomolecular solid-state NMR.

Solid state NMR (SSNMR) has earned a substantial success in the characterization of paramagnetic systems over the last decades. Nowadays, the resolution and sensitivity of solid state NMR in biological molecules has improved significantly and these advancements can be translated into the study of paramagnetic biomolecules. However, the electronic properties of different metal centers affect the quality of their SSNMR spectra differently, and not all systems turn out to be equally easy to approach by this technique. In this review we will try to give an overview of the properties of different paramagnetic centers and how they can be used to increase the chances of experimental success.

Bimodal Fluorescence-Magnetic Resonance Contrast Agent for Apoptosis Imaging.

Effective cancer therapy largely depends on inducing apoptosis in cancer cells via chemotherapy and/or radiation. Monitoring apoptosis in real-time provides invaluable information for evaluating cancer therapy response and screening preclinical anticancer drugs. In this work, we describe the design, synthesis, characterization, and in vitro evaluation of caspase probe 1 (CP1), a bimodal fluorescence-magnetic resonance (FL-MR) probe that exhibits simultaneous FL-MR turn-on response to caspase-3/7. Both caspases exist as inactive zymogens in normal cells but are activated during apoptosis and are unique biomarkers for this process. CP1 has three distinct components: a DOTA-Gd(III) chelate that provides the MR signal enhancement, tetraphenylethylene as the aggregation induced emission luminogen (AIEgen), and DEVD peptide which is a substrate for caspase-3/7. In response to caspase-3/7, the water-soluble peptide DEVD is cleaved and the remaining Gd(III)-AIEgen (Gad-AIE) conjugate aggregates leading to increased FL-MR signals. CP1 exhibited sensitive and selective dual FL-MR turn-on response to caspase-3/7 in vitro and was successfully tested by fluorescence imaging of apoptotic cells. Remarkably, we were able to use the FL response of CP1 to quantify the exact concentrations of inactive and active agents and accurately predict the MR signal in vitro. We have demonstrated that the aggregation-driven FL-MR probe design is a unique method for MR signal quantification. This probe design platform can be adapted for a variety of different imaging targets, opening new and exciting avenues for multimodal molecular imaging.