Chemical Shifts of the Carbohydrate Binding Domain of Galectin-3 from Magic Angle Spinning NMR and Hybrid Quantum Mechanics/Molecular Mechanics Calculations

Predicted and Experimental NMR Chemical Shifts at Variable Temperatures: The Effect of Protein Conformational Dynamics

NMR chemical shifts provide a sensitive probe of protein structure and dynamics but remain challenging to predict and interpret. We examine the effect of protein conformational distributions on 15N chemical shifts for dihydrofolate reductase (DHFR), comparing QM/MM predicted shifts with experimental shifts in solution as well as frozen distributions. Representative snapshots from MD trajectories exhibit variation in predicted 15N chemical shifts of up to 25 ppm. The average over the fluctuations is in significantly better agreement with room temperature solution experimental values than the prediction for any single optimal conformations. Meanwhile, solid-state NMR (SSNMR) measurements of frozen solutions at 105 K exhibit broad lines whose widths agree well with the widths of distributions of predicted shifts for samples from the trajectory. The backbone torsion angle ψi-1 varies over 60° on the picosecond time scale, compensated by φi. These fluctuations can explain much of the shift variation.

Cholesterol Biases the Conformational Landscape of the Chemokine Receptor CCR3: A MAS SSNMR-Filtered Molecular Dynamics Study

Cholesterol directs the pathway of ligand-induced G protein-coupled receptor (GPCR) signal transduction. The GPCR C-C motif chemokine receptor 3 (CCR3) is the principal chemotactic receptor for eosinophils, with roles in cancer metastasis and autoinflammatory conditions. Recently, we discovered a direct correlation between bilayer cholesterol and increased agonist-triggered CCR3 signal transduction. However, the allosteric molecular mechanism escalating ligand affinity and G protein coupling is unknown. To study cholesterol-guided CCR3 conformational selection, we implement comparative, objective measurement of protein architectures by scoring shifts (COMPASS) to grade model structures from molecular dynamics simulations. In this workflow, we scored predicted chemical shifts against 2-dimensional solid-state NMR 13C-13C correlation spectra of U-15N,13C-CCR3 samples prepared with and without cholesterol. Our analysis of trajectory model structures uncovers that cholesterol induces site-specific conformational restraint of extracellular loop (ECL) 2 and conserved motion in transmembrane helices and ECL3 not observed in simulations of bilayers with only phosphatidylcholine lipids. PyLipID analysis implicates direct cholesterol agency in CCR3 conformational selection and dynamics. Residue-residue contact scoring shows that cholesterol biases the conformational selection of the orthosteric pocket involving Y411.39, Y1133.32, and E2877.39. Lastly, we observe contact remodeling in activation pathway residues centered on the initial transmission switch, Na+ pocket, and R3.50 in the DRY motif. Our observations have unique implications for understanding of CCR3 ligand recognition and specificity and provide mechanistic insight into how cholesterol functions as an allosteric regulator of CCR3 signal transduction.

Predicted and Experimental NMR Chemical Shifts at Variable Temperatures: The Effect of Protein Conformational Dynamics

NMR chemical shifts provide a sensitive probe of protein structure and dynamics. Prediction of shifts, and therefore interpretation of shifts, particularly for the frequently measured amidic 15 N sites, remains a tall challenge. We demonstrate that protein 15 N chemical shift prediction from QM/MM predictions can be improved if conformational variation is included via MD sampling, focusing on the antibiotic target, E. coli Dihydrofolate reductase (DHFR). Variations of up to 25 ppm in predicted 15 N chemical shifts are observed over the trajectory. For solution shifts the average of fluctuations on the low picosecond timescale results in a superior prediction to a single optimal conformation. For low temperature solid state measurements, the histogram of predicted shifts for locally minimized snapshots with specific solvent arrangements sampled from the trajectory explains the heterogeneous linewidths; in other words, the conformations and associated solvent are 'frozen out' at low temperatures and result in inhomogeneously broadened NMR peaks. We identified conformational degrees of freedom that contribute to chemical shift variation. Backbone torsion angles show high amplitude fluctuations during the trajectory on the low picosecond timescale. For a number of residues, including I60, ψ varies by up to 60º within a conformational basin during the MD simulations, despite the fact that I60 (and other sites studied) are in a secondary structure element and remain well folded during the trajectory. Fluctuations in ψ appear to be compensated by other degrees of freedom in the protein, including φ of the succeeding residue, resulting in "rocking" of the amide plane with changes in hydrogen bonding interactions. Good agreement for both room temperature and low temperature NMR spectra provides strong support for the specific approach to conformational averaging of computed chemical shifts.

A Computational Perspective on Molecular Recognition by Galectins

This article presents an overview of recent computational studies dedicated to the analysis of binding between galectins and small-molecule ligands. We first present a summary of the most popular simulation techniques adopted for calculating binding poses and binding energies, and then discuss relevant examples reported in the literature for the three main classes of galectins (dimeric, tandem and chimera). We show that simulation of galectin-ligand interactions is a mature field which has proven invaluable for completing and unraveling experimental observations. Future perspectives to further improve the accuracy and cost-effectiveness of existing computational approaches will involve the development of new schemes to account for solvation and entropy effects, which represent the main current limitations to the accuracy of computational results.

Nonuniformly sampled exclusively-13 C/15 N 4D solid-state NMR experiments: Assignment and characterization of IKe phage capsid

An important step in the process of protein research by NMR is the assignment of chemical shifts. In the coat protein of IKe bacteriophage, there are 53 residues making up a long helix resulting in relatively high spectral ambiguity. Assignment thus requires the collection of a set of three-dimensional (3D) experiments and the preparation of sparsely labeled samples. Increasing the dimensionality can facilitate fast and reliable assignment of IKe and of larger proteins. Recent progress in nonuniform sampling techniques made the application of multidimensional NMR solid-state experiments beyond 3D more practical. 4D ¹ H-detected experiments have been demonstrated in high-fields and at spinning speeds of 60 kHz and higher but are not practical at spinning speeds of 10-20 kHz for fully protonated proteins. Here, we demonstrate the applicability of a nonuniformly sampled 4D ¹³ C/¹⁵ N-only correlation experiment performed at a moderate field of 14.1 T, which can incorporate sufficiently long acquisition periods in all dimensions. We show how a single CANCOCX experiment, supported by several 2D carbon-based correlation experiments, is utilized for the assignment of heteronuclei in the coat protein of the IKe bacteriophage. One sparsely labeled sample was used to validate sidechain assignment of several hydrophobic-residue sidechains. A comparison to solution NMR studies of isolated IKe coat proteins embedded in micelles points to key residues involved in structural rearrangement of the capsid upon assembly of the virus. The benefits of 4D to a quicker assignment are discussed, and the method may prove useful for studying proteins at relatively low fields.

Recent advances in computational methods for biosensor design

Biosensors are analytical tools with a great application in healthcare, food quality control, and environmental monitoring. They are of considerable interest to be designed by using cost-effective and efficient approaches. Designing biosensors with improved functionality or application in new target detection has been converted to a fast-growing field of biomedicine and biotechnology branches. Experimental efforts have led to valuable successes in the field of biosensor design; however, some deficiencies restrict their utilization for this purpose. Computational design of biosensors is introduced as a promising key to eliminate the gap. A set of reliable structure prediction of the biosensor segments, their stability, and accurate descriptors of molecular interactions are required to computationally design biosensors. In this review, we provide a comprehensive insight into the progress of computational methods to guide the design and development of biosensors, including molecular dynamics simulation, quantum mechanics calculations, molecular docking, virtual screening, and a combination of them as the hybrid methodologies. By relying on the recent advances in the computational methods, an opportunity emerged for them to be complementary or an alternative to the experimental methods in the field of biosensor design.

Accurate Backbone 13 C and 15 N Chemical Shift Tensors in Galectin-3 Determined by MAS NMR and QM/MM: Details of Structure and Environment Matter

Chemical shift tensors obtained from solid-state NMR spectroscopy are very sensitive reporters of structure and dynamics in proteins. While accurate 13 C and 15 N chemical shift tensors are accessible by magic angle spinning (MAS) NMR, their quantum mechanical calculations remain challenging, particularly for 15 N atoms. Here we compare experimentally determined backbone 13 Cα and 15 NH chemical shift tensors by MAS NMR with hybrid quantum mechanics/molecular mechanics/molecular dynamics (MD-QM/MM) calculations for the carbohydrate-binding domain of galectin-3. Excellent agreement between experimental and computed 15 NH chemical shift anisotropy values was obtained using the Amber ff15ipq force field when solvent dynamics was taken into account in the calculation. Our results establish important benchmark conditions for improving the accuracy of chemical shift calculations in proteins and may aid in the validation of protein structure models derived by MAS NMR.

Accuracy and precision of protein structures determined by magic angle spinning NMR spectroscopy: for some 'with a little help from a friend'

We present a systematic investigation into the attainable accuracy and precision of protein structures determined by heteronuclear magic angle spinning solid-state NMR for a set of four proteins of varied size and secondary structure content. Structures were calculated using synthetically generated random sets of C-C distances up to 7 Å at different degrees of completeness. For single-domain proteins, 9-15 restraints per residue are sufficient to derive an accurate model structure, while maximum accuracy and precision are reached with over 15 restraints per residue. For multi-domain proteins and protein assemblies, additional information on domain orientations, quaternary structure and/or protein shape is needed. As demonstrated for the HIV-1 capsid protein assembly, this can be accomplished by integrating MAS NMR with cryoEM data. In all cases, inclusion of TALOS-derived backbone torsion angles improves the accuracy for small number of restraints, while no further increases are noted for restraint completeness above 40%. In contrast, inclusion of TALOS-derived torsion angle restraints consistently increases the precision of the structural ensemble at all degrees of distance restraint completeness.