Mirco Zerbetto

General Theoretical/Computational Tool for Interpreting NMR Spin Relaxation in Proteins

We developed in recent years the slowly relaxing local structure (SRLS) approach for analyzing NMR spin relaxation in proteins. SRLS is a two-body coupled rotator model which accounts rigorously for mode-coupling between the global motion of the protein and the local motion of the spin-bearing probe and allows for general properties of the second rank tensors involved. We showed that a general tool of data analysis requires both capabilities. Several important functionalities were missing in our previous implementations of SRLS in data fitting schemes, and in some important cases, the calculations were tedious. Here we present a general implementation which allows for asymmetric local and global diffusion tensors, distinct local ordering and local diffusion frames, and features a rhombic local potential which includes Wigner matrix element terms of ranks 2 and 4. A recently developed hydrodynamics-based approach for calculating global diffusion tensors has been incorporated into the data-fitting scheme. The computational efficiency of the latter has been increased significantly through object-oriented programming within the scope of the C++ programming language, and code parallelization. A convenient graphical user interface is provided. Currently autocorrelated (15)N spin relaxation data can be analyzed effectively. Adaptation to any autocorrelated and cross-correlated relaxation analysis is straightforward. New physical insight is gleaned on largely preserved local structure in solution, even in chain segments which experience slow local motion. Prospects associated with improved dynamic models, and new applications made possible by the current implementation of SRLS, are delineated.

Hydrodynamic modeling of diffusion tensor properties of flexible molecules.

We present a computationally efficient implementation of hydrodynamic modeling for the evaluation of diffusion tensors of molecules with internal degrees of freedom, adapted to take into account information from linear scaling computations of solvent accessible surfaces implemented in the framework of last generation continuum solvent models. Torsional angles are taken also explicitly into account, while retaining correct hydrodynamic interactions. A comparison with literature data is presented to prove the effectiveness of the approach for a wide range of molecular dimensions and solvent environments.

Spectroscopic Insights into Carbon Dot Systems

The controversial nature of the fluorescent properties of carbon dots (CDs), ascribed either to surface states or to small molecules adsorbed onto the carbon nanostructures, is an unresolved issue. To date, an accurate picture of CDs and an exhaustive structure-property correlation are still lacking. Using two unconventional spectroscopic techniques, fluorescence correlation spectroscopy (FCS) and time-resolved electron paramagnetic resonance (TREPR), we contribute to fill this gap. Although electron micrographs indicate the presence of carbon cores, FCS reveals that the emission properties of CDs are based neither on those cores nor on molecular species linked to them, but rather on free molecules. TREPR provides deeper insights into the structure of carbon cores, where C sp2 domains are embedded within C sp3 scaffolds. FCS and TREPR prove to be powerful techniques, characterizing CDs as inherently heterogeneous systems, providing insights into the nature of such systems and paving the way to standardization of these nanomaterials.

Integrated Computational Approach to the Analysis of NMR Relaxation in Proteins: Application to ps−ns Main Chain 15N−1H and Global Dynamics of the Rho GTPase Binding Domain of Plexin-B1

An integrated computational methodology for interpreting NMR spin relaxation in proteins has been developed. It combines a two-body coupled-rotator stochastic model with a hydrodynamics-based approach for protein diffusion, together with molecular dynamics based calculations for the evaluation of the coupling potential of mean force. The method is applied to ¹⁵N relaxation of N-H bonds in the Rho GTPase binding (RBD) domain of plexin-B1, which exhibits intricate internal mobility. Bond vector dynamics are characterized by a rhombic local ordering tensor, S, with principal values S₀² and S₂², and an axial local diffusion tensor, D₂, with principal values D(2,||) and D(2,⊥). For α-helices and β-sheets we find that S₀² ~ -0.5 (strong local ordering), -1.2 < S₂² < -0.8 (large S tensor anisotropy), D(2,⊥) ~ D₁ = 1.93 × 10⁷ s⁻¹ (D₁ is the global diffusion rate), and log(D(2,||)/D₁) ~ 4. For α-helices the z-axis of the local ordering frame is parallel to the C(α)-C(α) axis. For β-sheets the z-axes of the S and D₂ tensors are parallel to the N-H bond. For loops and terminal chain segments the local ordering is generally weaker and more isotropic. On average, D(2,⊥) ~ D₁ also, but log(D(2,||)/D₁) is on the order of 1-2. The tensor orientations are diversified. This study sets forth an integrated computational approach for treating NMR relaxation in proteins by combining stochastic modeling and molecular dynamics. The approach developed provides new insights by its application to a protein that experiences complex dynamics.

Simulation of electron spin resonance spectroscopy in diverse environments: An integrated approach ☆

We discuss in this work a new software tool, named E-SpiReS (Electron Spin Resonance Simulations), aimed at the interpretation of dynamical properties of molecules in fluids from electron spin resonance (ESR) measurements. The code implements an integrated computational approach (ICA) for the calculation of relevant molecular properties that are needed in order to obtain spectral lines. The protocol encompasses information from atomistic level (quantum mechanical) to coarse grained level (hydrodynamical), and evaluates ESR spectra for rigid or flexible single or multi-labeled paramagnetic molecules in isotropic and ordered phases, based on a numerical solution of a stochastic Liouville equation. E-SpiReS automatically interfaces all the computational methodologies scheduled in the ICA in a way completely transparent for the user, who controls the whole calculation flow via a graphical interface. Parallelized algorithms are employed in order to allow running on calculation clusters, and a web applet Java has been developed with which it is possible to work from any operating system, avoiding the problems of recompilation. E-SpiReS has been used in the study of a number of different systems and two relevant cases are reported to underline the promising applicability of the ICA to complex systems and the importance of similar software tools in handling a laborious protocol.

An integrated approach to NMR spin relaxation in flexible biomolecules: Application to β-D-glucopyranosyl-(1 →6)-α-D-mannopyranosyl-OMe

The description of the reorientational dynamics of flexible molecules is a challenging task, in particular when the rates of internal and global motions are comparable. The commonly used simple mode-decoupling models are based on the assumption of statistical independence between these motions. This assumption is not valid when the time scale separation between their rates is small, a situation that was found to arise in oligosaccharides in the context of certain internal motions. To make possible the interpretation of NMR spin relaxation data from such molecules, we developed a comprehensive approach generally applicable to flexible rotators with one internal degree of freedom. This approach integrates a stochastic description of coupled global tumbling and internal torsional motion, quantum chemical calculations of the local potential and the local geometry at the site of the restricted torsion, and hydrodynamics-based calculations of the diffusive properties. The method is applied to the disaccharide beta-D-Glcp-(1-->6)-alpha-D-[6-(13)C]-Manp-OMe dissolved in a DMSO-d(6)/D(2)O cryosolvent. The experimental NMR relaxation parameters, associated with the (13)CH(2) probe residing at the glycosidic linkage, include (13)C T(1) and T(2) and (13)C-{(1)H} nuclear Overhauser enhancement (NOE) as well as longitudinal and transverse dipole-dipole cross-correlated relaxation rates, acquired in the temperature range of 253-293 K. These data are predicted successfully by the new theory with only the H-C-H angle allowed to vary. Previous attempts to fit these data using mode-decoupling models failed.

On the interpretation of continuous wave electron spin resonance spectra of tempo-palmitate in 5-cyanobiphenyl

Electron spin resonance (ESR) measurements are highly informative on the dynamic behavior of molecules, which is of fundamental importance to understand their stability, biological functions and activities, and catalytic action. The wealth of dynamic information which can be extracted from a continuous wave electron spin resonance (cw-ESR) spectrum can be inferred by a basic theoretical approach defined within the stochastic Liouville equation formalism, i.e., the direct inclusion of motional dynamics in the form of stochastic (Fokker-Planck/diffusive) operators in the super Hamiltonian H governing the time evolution of the system. Modeling requires the characterization of magnetic parameters (e.g., hyperfine and Zeeman tensors) and the calculation of ESR observables in terms of spectral densities. The magnetic observables can be pursued by the employment of density functional theory which is apt, provided that hybrid functionals are employed, for the accurate computation of structural properties of molecular systems. Recently, an ab initio integrated computational approach to the in silico interpretation of cw-ESR spectra of multilabeled systems in isotropic fluids has been discussed. In this work we present the extension to the case of nematic liquid crystalline environments by performing simulations of the ESR spectra of the prototypical nitroxide probe 4-(hexadecanoyloxy)-2,2,6,6-tetramethylpiperidine-1-oxy in isotropic and nematic phases of 5-cyanobiphenyl. We first discuss the basic ingredients of the integrated approach, i.e., (1) determination of geometric and local magnetic parameters by quantum-mechanical calculations, taking into account the solvent and, when needed, the vibrational averaging contributions; (2) numerical solution of a stochastic Liouville equation in the presence of diffusive rotational dynamics, based on (3) parameterization of diffusion rotational tensor provided by a hydrodynamic model. Next we present simulated spectra with minimal resorting to fitting procedures, proving that the combination of sensitive ESR spectroscopy and sophisticated modeling can be highly helpful in providing three-dimensional structural and dynamic information on molecular systems in anisotropic environments.

SRLS analysis of 15N spin relaxation from E. coli ribonuclease HI: the tensorial perspective.

15N–H relaxation parameters from ribonuclease HI (RNase H), acquired in previous work at magnetic fields of 14.1 and 18.8 T, and at 300 K, are analyzed with the mode-coupling slowly relaxing local structure (SRLS) approach. In accordance with standard theoretical treatments of restricted motions, SRLS approaches N-H bond dynamics from a tensorial perspective. As shown previously, a physically adequate description of this phenomenon has to account for the asymmetry of the local spatial restrictions. So far, we used rhombic local ordering tensors; this is straightforward but computationally demanding. Here, we propose substantiating the asymmetry of the local spatial restrictions in terms of tilted axial local ordering (S) and local diffusion (D2) tensors. Although less straightforward, this description provides physically sound structural and dynamic information and is efficient computationally. We find that the local order parameter, S(0)2, is on average 0.89 (0.84, and may be as small as 0.6) for the secondary structure elements (loops). The main local ordering axis deviates from the C(i-1)α-C(i)α axis by less than 6°. At 300 K, D(2,perpendicular) is virtually the same as the global diffusion rate, D1 = 1.8 × 10(7) s(-1). The correlation time 1/6D(2,parallel) ranges from 3-125 (208-344) ps for the secondary structure elements (loops) and is on average 125 ps for the C-terminal segment. The main local diffusion axis deviates from the N-H bond by less than 2° (10°) for the secondary structure elements (loops). An effective data-fitting protocol, which leads in most cases to unambiguous results with limited uncertainty, has been devised. A physically sound and computationally effective methodology for analyzing 15N relaxation in proteins, that provides a new picture of N–H bond structural dynamics in proteins, has been set forth.

Analysis of 15N-1H NMR relaxation in proteins by a combined experimental and molecular dynamics simulation approach: picosecond-nanosecond dynamics of the Rho GTPase binding domain of plexin-B1 in the dimeric state indicates allosteric pathways.

We investigate picosecond–nanosecond dynamics of the Rho-GTPase Binding Domain (RBD) of plexin-B1, which plays a key role in plexin-mediated cell signaling. Backbone 15N relaxation data of the dimeric RBD are analyzed with the model-free (MF) method, and with the slowly relaxing local structure/molecular dynamics (SRLS-MD) approach. Independent analysis of the MD trajectories, based on the MF paradigm, is also carried out. MF is a widely popular and simple method, SRLS is a general approach, and SRLS-MD is an integrated approach we developed recently. Corresponding parameters from the RBD dimer, a previously studied RBD monomer mutant, and the previously studied complex of the latter with the GTPase Rac1, are compared. The L2, L3, and L4 loops of the plexin-B1 RBD are involved in interactions with other plexin domains, GTPase binding, and RBD dimerization, respectively. Peptide groups in the loops of both the monomeric and dimeric RBD are found to experience weak and moderately asymmetric local ordering centered approximately at the C(i–1)(α)–C(i)(α) axes, and nanosecond backbone motion. Peptide groups in the α-helices and the β-strands of the dimer (the β-strands of the monomer) experience strong and highly asymmetric local ordering centered approximately at the C(i–1)(α)–C(i)(α) axes (N–H bonds). N–H fluctuations occur on the picosecond time scale. An allosteric pathway for GTPase binding, providing new insights into plexin function, is delineated.

Stochastic Modeling of Flexible Biomolecules Applied to NMR Relaxation : 2. Interpretation of Complex Dynamics in Linear Oligosaccharides

A computational stochastic approach is applied to the description of flexible molecules. By combining (i) molecular dynamics simulations, (ii) hydrodynamics approaches, and (iii) a multidimensional diffusive description for internal and global dynamics, it is possible to build an efficient integrated approach to the interpretation of relaxation processes in flexible systems. In particular, the model is applied to the interpretation of nuclear magnetic relaxation measurements of linear oligosaccharides, namely a mannose-containing trisaccharide and the pentasaccharide LNF-1. Experimental data are reproduced with sufficient accuracy without free model parameters.