X-ray Crystallography Module in MD Simulation Program Amber 2023. Refining the Models of Protein Crystals.

The MD simulation package Amber offers an attractive platform to refine crystallographic structures of proteins: (i) state-of-the-art force fields help to regularize protein coordinates and reconstruct the poorly diffracting elements of the structure, such as flexible loops; (ii) MD simulations restrained by the experimental diffraction data provide an effective strategy to optimize structural models of protein crystals, including explicitly modeled interstitial solvent as well as crystal contacts. Here, we present the new crystallography module xray, released as a part of the Amber 2023 package. This module contains functions to calculate and scale structure factors (including the contributions from bulk solvent), evaluate the maximum-likelihood-type crystallographic potential, and compute its derivative forces. The X-ray functionality of Amber no longer relies on external dependencies so that the full advantage of GPU acceleration can be taken. This makes it possible to refine in a short time hundreds of crystal models, including supercell models comprised of multiple unit cells. The new automated Amber-based refinement procedure leads to an appreciable improvement in Rfree (in some cases, by as much as 0.067) as well as MolProbity scores.

Enhanced TROSY Effect in [2-19 F, 2-13 C] Adenosine and ATP Analogs Facilitates NMR Spectroscopy of Very Large Biological RNAs in Solution.

Large RNAs are central to cellular functions, but characterizing such RNAs remains challenging by solution NMR. We present two labeling technologies based on [2-19 F, 2-13 C]-adenosine, which allow the incorporation of aromatic 19 F-13 C spin pairs. The labels when coupled with the transverse relaxation optimized spectroscopy (TROSY) enable us to probe RNAs comprising up to 124 nucleotides. With our new [2-19 F, 2-13 C]-adenosine-phosphoramidite, all resonances of the human hepatitis B virus epsilon RNA could be readily assigned. With [2-19 F, 2-13 C]-adenosine triphosphate, the 124 nt pre-miR-17-NPSL1-RNA was produced via in vitro transcription and the TROSY spectrum of this 40 kDa [2-19 F, 2-13 C]-A-labeled RNA featured sharper resonances than the [2-1 H, 2-13 C]-A sample. The mutual cancelation of the chemical-shift-anisotropy and the dipole-dipole-components of TROSY-resonances leads to narrow linewidths over a wide range of molecular weights. With the synthesis of a non-hydrolysable [2-19 F, 2-13 C]-adenosine-triphosphate, we facilitate the probing of co-factor binding in kinase complexes and NMR-based inhibitor binding studies in such systems. Our labels allow a straightforward assignment for larger RNAs via a divide-and-conquer/mutational approach. The new [2-19 F, 2-13 C]-adenosine precursors are a valuable addition to the RNA NMR toolbox and will allow the study of large RNAs/RNA protein complexes in vitro and in cells.

AmberTools.

AmberTools is a free and open-source collection of programs used to set up, run, and analyze molecular simulations. The newer features contained within AmberTools23 are briefly described in this Application note.

Molecular dynamics analysis of a flexible loop at the binding interface of the SARS-CoV-2 spike protein receptor-binding domain.

Since the identification of the SARS-CoV-2 virus as the causative agent of the current COVID-19 pandemic, considerable effort has been spent characterizing the interaction between the Spike protein receptor-binding domain (RBD) and the human angiotensin converting enzyme 2 (ACE2) receptor. This has provided a detailed picture of the end point structure of the RBD-ACE2 binding event, but what remains to be elucidated is the conformation and dynamics of the RBD prior to its interaction with ACE2. In this work, we utilize molecular dynamics simulations to probe the flexibility and conformational ensemble of the unbound state of the receptor-binding domain from SARS-CoV-2 and SARS-CoV. We have found that the unbound RBD has a localized region of dynamic flexibility in Loop 3 and that mutations identified during the COVID-19 pandemic in Loop 3 do not affect this flexibility. We use a loop-modeling protocol to generate and simulate novel conformations of the CoV2-RBD Loop 3 region that sample conformational space beyond the ACE2 bound crystal structure. This has allowed for the identification of interesting substates of the unbound RBD that are lower energy than the ACE2-bound conformation, and that block key residues along the ACE2 binding interface. These novel unbound substates may represent new targets for therapeutic design.

Integral equation models for solvent in macromolecular crystals.

The solvent can occupy up to ∼70% of macromolecular crystals, and hence, having models that predict solvent distributions in periodic systems could improve the interpretation of crystallographic data. Yet, there are few implicit solvent models applicable to periodic solutes, and crystallographic structures are commonly solved assuming a flat solvent model. Here, we present a newly developed periodic version of the 3D-reference interaction site model (RISM) integral equation method that is able to solve efficiently and describe accurately water and ion distributions in periodic systems; the code can compute accurate gradients that can be used in minimizations or molecular dynamics simulations. The new method includes an extension of the Ornstein-Zernike equation needed to yield charge neutrality for charged solutes, which requires an additional contribution to the excess chemical potential that has not been previously identified; this is an important consideration for nucleic acids or any other charged system where most or all the counter- and co-ions are part of the "disordered" solvent. We present several calculations of proteins, RNAs, and small molecule crystals to show that x-ray scattering intensities and the solvent structure predicted by the periodic 3D-RISM solvent model are in closer agreement with the experiment than are intensities computed using the default flat solvent model in the refmac5 or phenix refinement programs, with the greatest improvement in the 2 to 4 Šrange. Prospects for incorporating integral equation models into crystallographic refinement are discussed.

Simulations of Kindlin-2 PIP binding domains reveal protonation-dependent membrane binding modes.

Kindlin-2, a member of the Kindlin family of peripheral membrane proteins, is important for integrin activation and stabilization of epidermal growth factor receptor. It associates with the cytoplasmic face of the plasma membrane via dedicated phosphatidylinositol phosphate binding domains located in the N-terminal F0 and Pleckstrin Homology domains. These domains have binding affinity for phosphatidylinositol 4,5-bisphosphate and, to a greater degree, phosphatidylinositol 3,4,5-trisphosphate. The biological significance of the differential binding of these phosphatidylinositol phosphates to Kindlin-2 and the mechanism by which they activate Kindlin-2 are not well understood. Recently, ssNMR identified the predominant protonation states of phosphatidylinositol 4,5-bisphosphate and phosphatidylinositol 3,4,5-trisphosphate near physiological pH in the presence of anionic lipids. Here, we perform atomistic simulation of the bound state of the Pleckstrin Homology and F0 domains of Kindlin-2 at membranes containing phosphatidylinositol 4,5-bisphosphate/phosphatidylinositol 3,4,5-trisphosphate with differing protonation states. This computational approach demonstrates that these two phosphatidylinositol phosphates differently modulate Kindlin-2 subdomain binding in a protonation-state-dependent manner. We speculate these variations in binding mode provide a mechanism for intracellular pH and Ca2+ influx to control the membrane binding behavior and activity of Kindlin-2.

Coupled intra- and interdomain dynamics support domain cross-talk in Pin1.

The functional mechanisms of multidomain proteins often exploit interdomain interactions, or "cross-talk." An example is human Pin1, an essential mitotic regulator consisting of a Trp-Trp (WW) domain flexibly tethered to a peptidyl-prolyl isomerase (PPIase) domain, resulting in interdomain interactions important for Pin1 function. Substrate binding to the WW domain alters its transient contacts with the PPIase domain via means that are only partially understood. Accordingly, we have investigated Pin1 interdomain interactions using NMR paramagnetic relaxation enhancement (PRE) and molecular dynamics (MD) simulations. The PREs show that apo-Pin1 samples interdomain contacts beyond the range suggested by previous structural studies. They further show that substrate binding to the WW domain simultaneously alters interdomain separation and the internal conformation of the WW domain. A 4.5-µs all-atom MD simulation of apo-Pin1 suggests that the fluctuations of interdomain distances are correlated with fluctuations of WW domain interresidue contacts involved in substrate binding. Thus, the interdomain/WW domain conformations sampled by apo-Pin1 may already include a range of conformations appropriate for binding Pin1's numerous substrates. The proposed coupling between intra-/interdomain conformational fluctuations is a consequence of the dynamic modular architecture of Pin1. Such modular architecture is common among cell-cycle proteins; thus, the WW-PPIase domain cross-talk mechanisms of Pin1 may be relevant for their mechanisms as well.

The Interplay between Molten Globules and Heme Disassociation Defines Human Hemoglobin Disassembly.

Hemoglobin functions as a tetrameric oxygen transport protein, with each subunit containing a heme cofactor. Its denaturation, either in vivo or in vitro, involves autoxidation to methemoglobin, followed by cofactor loss and globin unfolding. We have proposed a global disassembly scheme for human methemoglobin, linking hemin (ferric protoporphyrin IX) disassociation and apoprotein unfolding pathways. The model is based on the evaluation of circular dichroism and visible absorbance measurements of guanidine-hydrochloride-induced disassembly of methemoglobin and previous measurements of apohemoglobin unfolding. The populations of holointermediates and equilibrium disassembly parameters were estimated quantitatively for adult and fetal hemoglobins. The key stages are characterized by hexacoordinated hemichrome intermediates, which are important for preventing hemin disassociation from partially unfolded, molten globular species during early disassembly and late-stage assembly events. Both unfolding experiments and independent small angle x-ray scattering measurements demonstrate that heme disassociation leads to the loss of tetrameric structural integrity. Our model predicts that after autoxidation, dimeric and monomeric hemichrome intermediates occur along the disassembly pathway inside red cells, where the hemoglobin concentration is very high. This prediction suggests why misassembled hemoglobins often get trapped as hemichromes that accumulate into insoluble Heinz bodies in the red cells of patients with unstable hemoglobinopathies. These Heinz bodies become deposited on the cell membranes and can lead to hemolysis. Alternatively, when acellular hemoglobin is diluted into blood plasma after red cell lysis, the disassembly pathway appears to be dominated by early hemin disassociation events, which leads to the generation of higher fractions of unfolded apo subunits and free hemin, which are known to damage the integrity of blood vessel walls. Thus, our model provides explanations of the pathophysiology of hemoglobinopathies and other disease states associated with unstable globins and red cell lysis and also insights into the factors governing hemoglobin assembly during erythropoiesis.

Atomistic Simulations of Heme Dissociation Pathways in Human Methemoglobins Reveal Hidden Intermediates.

Heme dissociations disrupt function and structural integrity of human hemoglobin and trigger various cardiovascular complications. These events become significant in methemoglobins that have undergone autoxidation of ferrous into ferric heme. We have structurally characterized the heme disassociation pathways for adult tetrameric methemoglobins using all-atom molecular dynamics simulations. These reveal that bis-histidine hemichromes, characterized here by the coordination of heme iron to both the F8 (proximal) and E7 (distal) histidines, are seen as intermediates following dissociation of the water molecule distally bound to each heme iron. Later, the breaking of coordination between heme iron and proximal histidine disrupts the F helix and pushes it away from the heme cavity, enabling both bulk solvent penetration and disruption of tetramer interface interactions. The interactions inhibiting heme dissociation were then seen to be (i) either a direct or a water-molecule-mediated interaction between distal histidine and heme iron and (ii) stacking between heme and the αCE1/ßCD1 phenylalanine residue. These interactions are less important in the ß than in α subunits due to a more flexible ß subunit CE loop region. The absence of a distal histidine interaction in the H(E7)L mutant and increased heme cavity volume in the V(E11)A mutant both promoted heme escape from the protein interior. Adult and fetal hemoglobins were seen to share a general heme disassociation pathway and intermediates due to the conservation of key heme pocket residues. The intermediates seen here are analyzed in light of experimental studies of heme dissociation and pathways of certain hemoglobinopathies.

Molecular underpinnings of integrin binding to collagen-mimetic peptides containing vascular Ehlers-Danlos syndrome-associated substitutions.

Collagens carry out critical extracellular matrix (ECM) functions by interacting with numerous cell receptors and ECM components. Single glycine substitutions in collagen III, which predominates in vascular walls, result in vascular Ehlers-Danlos syndrome (vEDS), leading to arterial, uterine, and intestinal rupture and an average life expectancy of <50 years. Collagen interactions with integrin α2ß1 are vital for platelet adhesion and activation; however, how these interactions are impacted by vEDS-associated mutations and by specific amino acid substitutions is unclear. Here, we designed collagen-mimetic peptides (CMPs) with previously reported Gly → Xaa (Xaa = Ala, Arg, or Val) vEDS substitutions within a high-affinity integrin α2ß1-binding motif, GROGER. We used these peptides to investigate, at atomic-level resolution, how these amino acid substitutions affect the collagen III-integrin α2ß1 interaction. Using a multitiered approach combining biological adhesion assays, CD, NMR, and molecular dynamics (MD) simulations, we found that these substitutions differentially impede human mesenchymal stem cell spreading and integrin α2-inserted (α2I) domain binding to the CMPs and were associated with triple-helix destabilization. Although an Ala substitution locally destabilized hydrogen bonding and enhanced mobility, it did not significantly reduce the CMP-integrin interactions. MD simulations suggested that bulkier Gly → Xaa substitutions differentially disrupt the CMP-α2I interaction. The Gly → Arg substitution destabilized CMP-α2I side-chain interactions, and the Gly → Val change broke the essential Mg2+ coordination. The relationship between the loss of functional binding and the type of vEDS substitution provides a foundation for developing potential therapies for managing collagen disorders.

Deleterious effects of carbon-carbon dipolar coupling on RNA NMR dynamics.

Many regulatory RNAs undergo dynamic exchanges that are crucial for their biological functions and NMR spectroscopy is a versatile tool for monitoring dynamic motions of biomolecules. Meaningful information on biomolecular dynamics requires an accurate measurement of relaxation parameters such as longitudinal (R1) rates, transverse (R2) rates and heteronuclear Overhauser effect (hNOE). However, earlier studies have shown that the large 13C-13C interactions complicate analysis of the carbon relaxation parameters. To investigate the effect of 13C-13C interactions on RNA dynamic studies, we performed relaxation measurements on various RNA samples with different labeling patterns and compared these measurements with the computational simulations. For uniformly labeled samples, contributions of the neighboring carbon to R1 measurements were observed. These contributions increased with increasing magnetic field and overall correlation time ([Formula: see text]) for R1 rates, necessitating more careful analysis for uniformly labeled large RNAs. In addition, the hNOE measurements were also affected by the adjacent carbon nuclei. Unlike R1 rates, R1ρ rates showed relatively good agreement between uniformly- and site-selectively labeled samples, suggesting no dramatic effect from their attached carbon, in agreement with previous observations. Overall, having more accurate rate measurements avoids complex analysis and will be a key for interpreting 13C relaxation rates for molecular motion that can provide valuable insights into cellular molecular recognition events.

A twist in the road less traveled: The AMBER ff15ipq-m force field for protein mimetics.

We present a new force field, AMBER ff15ipq-m, for simulations of protein mimetics in applications from therapeutics to biomaterials. This force field is an expansion of the AMBER ff15ipq force field that was developed for canonical proteins and enables the modeling of four classes of artificial backbone units that are commonly used alongside natural α residues in blended or "heterogeneous" backbones: chirality-reversed D-α-residues, the Cα-methylated α-residue Aib, homologated ß-residues (ß3) bearing proteinogenic side chains, and two cyclic ß residues (ßcyc; APC and ACPC). The ff15ipq-m force field includes 472 unique atomic charges and 148 unique torsion terms. Consistent with the AMBER IPolQ lineage of force fields, the charges were derived using the Implicitly Polarized Charge (IPolQ) scheme in the presence of explicit solvent. To our knowledge, no general force field reported to date models the combination of artificial building blocks examined here. In addition, we have derived Karplus coefficients for the calculation of backbone amide J-coupling constants for ß3Ala and ACPC ß residues. The AMBER ff15ipq-m force field reproduces experimentally observed J-coupling constants in simple tetrapeptides and maintains the expected conformational propensities in reported structures of proteins/peptides containing the artificial building blocks of interest-all on the µs timescale. These encouraging results demonstrate the power and robustness of the IPolQ lineage of force fields in modeling the structure and dynamics of natural proteins as well as mimetics with protein-inspired artificial backbones in atomic detail.

Characterizing Watson-Crick versus Hoogsteen Base Pairing in a DNA-Protein Complex Using Nuclear Magnetic Resonance and Site-Specifically 13C- and 15N-Labeled DNA.

A( syn)-T and G( syn)-C+ Hoogsteen base pairs in protein-bound DNA duplexes can be difficult to resolve by X-ray crystallography due to ambiguous electron density and by nuclear magnetic resonance (NMR) spectroscopy due to poor chemical shift dispersion and size limitations with solution-state NMR spectroscopy. Here we describe an NMR strategy for characterizing Hoogsteen base pairs in protein-DNA complexes, which relies on site-specifically incorporating 13C- and 15N-labeled nucleotides into DNA duplexes for unambiguous resonance assignment and to improve spectral resolution. The approach was used to resolve the conformation of an A-T base pair in a crystal structure of an ∼43 kDa complex between a 34 bp duplex DNA and the integration host factor (IHF) protein. In the crystal structure (Protein Data Bank entry 1IHF ), this base pair adopts an unusual Hoogsteen conformation with a distorted sugar backbone that is accommodated by a nearby nick used to aid in crystallization. The NMR chemical shifts and interproton nuclear Overhauser effects indicate that this base pair predominantly adopts a Watson-Crick conformation in the intact DNA-IHF complex under solution conditions. Consistent with these NMR findings, substitution of 7-deazaadenine at this base pair resulted in only a small (∼2-fold) decrease in the IHF-DNA binding affinity. The NMR strategy provides a new approach for resolving crystallographic ambiguity and more generally for studying the structure and dynamics of protein-DNA complexes in solution.

Predicting Site-Binding Modes of Ions and Water to Nucleic Acids Using Molecular Solvation Theory.

Site binding of ions and water shapes nucleic acids folding, dynamics, and biological function, complementing the more diffuse, nonspecific "territorial" ion binding. Unlike territorial binding, prediction of site-specific binding to nucleic acids remains an unsolved challenge in computational biophysics. This work presents a new toolset based on the 3D-RISM molecular solvation theory and topological analysis that predicts cation and water site binding to nucleic acids. 3D-RISM is shown to accurately capture alkali cations and water binding to the central channel, transversal loops, and grooves of the Oxytricha nova's telomeres' G-quadruplex ( Oxy-GQ), in agreement with high-resolution crystallographic data. To improve the computed cation occupancy along the Oxy-GQ central channel, it was necessary to refine and validate new cation-oxygen parameters using structural and thermodynamic data available for crown ethers and ion channels. This single set of parameters that describes both localized and delocalized binding to various biological systems is used to gain insight into cation occupancy along the Oxy-GQ channel under various salt conditions. The paper concludes with prospects for extending the method to predict divalent cation binding to nucleic acids. This work advances the forefront of theoretical methods able to provide predictive insight into ion atmosphere effects on nucleic acids function.

High-Resolution ENDOR Spectroscopy Combined with Quantum Chemical Calculations Reveals the Structure of Nitrogenase Janus Intermediate E4(4H).

We have shown that the key state in N2 reduction to two NH3 molecules by the enzyme nitrogenase is E4(4H), the "Janus" intermediate, which has accumulated four [e-/H+] and is poised to undergo reductive elimination of H2 coupled to N2 binding and activation. Initial 1H and 95Mo ENDOR studies of freeze-trapped E4(4H) revealed that the catalytic multimetallic cluster (FeMo-co) binds two Fe-bridging hydrides, [Fe-H-Fe]. However, the analysis failed to provide a satisfactory picture of the relative spatial relationships of the two [Fe-H-Fe]. Our recent density functional theory (DFT) study yielded a lowest-energy form, denoted as E4(4H)(a), with two parallel Fe-H-Fe planes bridging pairs of "anchor" Fe on the Fe2,3,6,7 face of FeMo-co. However, the relative energies of structures E4(4H)(b), with one bridging and one terminal hydride, and E4(4H)(c), with one pair of anchor Fe supporting two bridging hydrides, were not beyond the uncertainties in the calculation. Moreover, a structure of V-dependent nitrogenase resulted in a proposed structure analogous to E4(4H)(c), and additional structures have been proposed in the DFT studies of others. To resolve the nature of hydride binding to the Janus intermediate, we performed exhaustive, high-resolution CW-stochastic 1H-ENDOR experiments using improved instrumentation, Mims 2H ENDOR, and a recently developed pulsed-ENDOR protocol ("PESTRE") to obtain absolute hyperfine interaction signs. These measurements are coupled to DFT structural models through an analytical point-dipole Hamiltonian for the hydride electron-nuclear dipolar coupling to its "anchoring" Fe ions, an approach that overcomes limitations inherent in both experimental interpretation and computational accuracy. The result is the freeze-trapped, lowest-energy Janus intermediate structure, E4(4H)(a).

Parmbsc1: a refined force field for DNA simulations.

We present parmbsc1, a force field for DNA atomistic simulation, which has been parameterized from high-level quantum mechanical data and tested for nearly 100 systems (representing a total simulation time of ∼ 140 µs) covering most of DNA structural space. Parmbsc1 provides high-quality results in diverse systems. Parameters and trajectories are available at http://mmb.irbbarcelona.org/ParmBSC1/.

NGLview-interactive molecular graphics for Jupyter notebooks.

Summary: NGLview is a Jupyter/IPython widget to interactively view molecular structures as well as trajectories from molecular dynamics simulations. Fast and scalable molecular graphics are provided through the NGL Viewer. The widget supports showing data from the file-system, online data bases and from objects of many popular analysis libraries including mdanalysis, mdtraj, pytraj, rdkit and more. Availability and implementation: The source code is freely available under the MIT license at https://github.com/arose/nglview. Python packages are available from PyPI and bioconda. NGLview uses Python on the server-side and JavaScript on the client. The integration with Jupyter is done through the ipywidgets package. The NGL Viewer is embedded client-side to provide WebGL accelerated molecular graphics. Contact: asr.moin@gmail.com.

Atomic structures of excited state A-T Hoogsteen base pairs in duplex DNA by combining NMR relaxation dispersion, mutagenesis, and chemical shift calculations.

NMR relaxation dispersion studies indicate that in canonical duplex DNA, Watson-Crick base pairs (bps) exist in dynamic equilibrium with short-lived low abundance excited state Hoogsteen bps. N1-methylated adenine (m1A) and guanine (m1G) are naturally occurring forms of damage that stabilize Hoogsteen bps in duplex DNA. NMR dynamic ensembles of DNA duplexes with m1A-T Hoogsteen bps reveal significant changes in sugar pucker and backbone angles in and around the Hoogsteen bp, as well as kinking of the duplex towards the major groove. Whether these structural changes also occur upon forming excited state Hoogsteen bps in unmodified duplexes remains to be established because prior relaxation dispersion probes provided limited information regarding the sugar-backbone conformation. Here, we demonstrate measurements of C3' and C4' spin relaxation in the rotating frame (R1ρ) in uniformly 13C/15N labeled DNA as sensitive probes of the sugar-backbone conformation in DNA excited states. The chemical shifts, combined with structure-based predictions using an automated fragmentation quantum mechanics/molecular mechanics method, show that the dynamic ensemble of DNA duplexes containing m1A-T Hoogsteen bps accurately model the excited state Hoogsteen conformation in two different sequence contexts. Formation of excited state A-T Hoogsteen bps is accompanied by changes in sugar-backbone conformation that allow the flipped syn adenine to form hydrogen-bonds with its partner thymine and this in turn results in overall kinking of the DNA toward the major groove. Results support the assignment of Hoogsteen bps as the excited state observed in canonical duplex DNA, provide an atomic view of DNA dynamics linked to formation of Hoogsteen bps, and lay the groundwork for a potentially general strategy for solving structures of nucleic acid excited states.

GPU-Accelerated Molecular Dynamics and Free Energy Methods in Amber18: Performance Enhancements and New Features.

We report progress in graphics processing unit (GPU)-accelerated molecular dynamics and free energy methods in Amber18. Of particular interest is the development of alchemical free energy algorithms, including free energy perturbation and thermodynamic integration methods with support for nonlinear soft-core potential and parameter interpolation transformation pathways. These methods can be used in conjunction with enhanced sampling techniques such as replica exchange, constant-pH molecular dynamics, and new 12-6-4 potentials for metal ions. Additional performance enhancements have been made that enable appreciable speed-up on GPUs relative to the previous software release.

Determination of accurate backbone chemical shift tensors in microcrystalline proteins by integrating MAS NMR and QM/MM.

Chemical shifts are highly sensitive probes of local conformation and overall structure. Both isotropic shifts and chemical shift tensors are readily accessible from NMR experiments but their quantum mechanical calculations remain challenging. In this work, we report and compare accurately measured and calculated 15NH and 13Cα chemical shift tensors in proteins, using the microcrystalline agglutinin from Oscillatoria agardhii (OAA). Experimental 13Cα and 15NH chemical tensors were obtained by solid-state NMR spectroscopy, employing tailored recoupling sequences, and for their quantum mechanics/molecular mechanics (QM/MM) calculations different sets of functionals were evaluated. We show that 13Cα chemical shift tensors are primarily determined by backbone dihedral angles and dynamics, while 15NH tensors mainly depend on local electrostatic contributions from solvation and hydrogen bonding. In addition, the influence of including crystallographic waters, the molecular mechanics geometry optimization protocol, and the level of theory on the accuracy of the calculated chemical shift tensors is discussed. Specifically, the power of QM/MM calculations in accurately predicting the unusually upfield shifted 1HN G26 and G93 resonances is highlighted. Our integrated approach is expected to benefit structure refinement of proteins and protein assemblies.