QSLE-v1.0: A New Software Package for the Calculation of Coupled Quantum-Classical Dynamics in Condensed Phases Based on a Stochastic Approach.

Recently, a stochastic version of the quantum-classical Liouville equation has been proposed [Campeggio, J.; Cortivo, R.; Zerbetto, M.

Analysis of Velocity Autocorrelation Functions from Molecular Dynamics Simulations of a Small Peptide by the Generalized Langevin Equation with a Power-Law Kernel.

Internal motions play an essential role in the biological functions of proteins and have been the subject of numerous theoretical and spectroscopic studies. Such complex environments are associated with anomalous diffusion where, in contrast to the classical Brownian motion, the relevant correlation functions have power law decays with time. In this work, we investigate the presence of long memory stochastic processes through the analysis of atomic velocity autocorrelation functions.

Intrinsic Flexibility beyond the Highly Ordered DNA Tetrahedron: An Integrative Spectroscopic and Molecular Dynamics Approach.

Since the introduction of DNA-based architectures, in the past decade, DNA tetrahedrons have aroused great interest. Applications of such nanostructures require structural control, especially in the perspective of their possible functionalities. In this work, an integrated approach for structural characterization of a tetrahedron structure is proposed with a focus on the fundamental biophysical aspects driving the assembly process.

The roto-conformational diffusion tensor as a tool to interpret molecular flexibility.

Stochastic modeling approaches can be used to rationalize complex molecular dynamical behaviours in solution, helping to interpret the coupling mechanisms among internal and external degrees of freedom, providing insight into reaction mechanisms and extracting structural and dynamical data from spectroscopic observables. However, the definition of comprehensive models is usually limited by (i) the difficulty in defining - without resorting to phenomenological assumptions - a representative reduced ensemble of molecular coordinates able to capture essential dynamical properties and (ii) the complexity of numerical or approximate treatments of the resulting equations. In this paper, we address the first of these two issues.

Parameter free evaluation of S2 reaction rates for halide substitution in halomethane.

We estimate the kinetic constants of a series of archetypal S2 reactions, , the nucleophilic substitutions of halides in halomethane. A parameter free, multiscale approach recently developed [Campeggio , , 2020, , 3455] is employed. The protocol relies on quantum mechanical calculations for the description of the energy profile along the intrinsic reaction coordinate, which is then mapped onto a reaction coordinate conveniently built for the reactive process.

Stochastic Modelling of C NMR Spin Relaxation Experiments in Oligosaccharides.

A framework for the stochastic description of relaxation processes in flexible macromolecules including dissipative effects has been recently introduced, starting from an atomistic view, describing the joint relaxation of internal coordinates and global degrees of freedom, and depending on parameters recoverable from classic force fields (energetics) and medium modelling at the continuum level (friction tensors). The new approach provides a rational context for the interpretation of magnetic resonance relaxation experiments. In its simplest formulation, the semi-flexible Brownian (SFB) model has been until now shown to reproduce correctly correlation functions and spectral densities related to orientational properties obtained by direct molecular dynamics simulations of peptides.

Enabling Circular Economy: The Overlooked Role of Inorganic Materials Chemistry.

Circular economy is considered a new chance to build a more sustainable world from both the social and the economic point of view. In this Essay, the possible contribution of inorganic chemistry towards a smooth transition to circularity in inorganic materials design and production is discussed by adopting an interdisciplinary approach.

Multiscale modeling of reaction rates: application to archetypal S2 nucleophilic substitutions.

We propose an approach to the evaluation of kinetic rates of elementary chemical reactions within Kramers' theory based on the definition of the reaction coordinate as a linear combination of natural, pseudo Z-matrix, internal coordinates of the system. The element of novelty is the possibility to evaluate the friction along the reaction coordinate, within a hydrodynamic framework developed recently [J. Campeggio et al.

Glycosidic linkage flexibility: The ψ torsion angle has a bimodal distribution in α-L-Rhap-(1→2)-α-L-Rhap-OMe as deduced from C NMR spin relaxation.

The molecular dynamics (MD) computer simulation technique is powerful for the investigation of conformational equilibrium properties of biomolecules. In particular, free energy surfaces of the torsion angles (those degrees of freedom from which the geometry mostly depends) allow one to access conformational states, as well as kinetic information, i.e.

Similarity and Specificity of Chlorophyll Triplet State in Comparison to Chlorophyll as Revealed by EPR/ENDOR and DFT Calculations.

An investigation of the photoexcited triplet state of chlorophyll (Chl) has been carried out by means of electron nuclear double resonance, both in a frozen organic solvent and in a protein environment provided by the water-soluble chlorophyll protein of . Density functional theory calculations have allowed the complete assignment of the observed hyperfine couplings corresponding to the methine protons and the methyl groups, leading to a complete picture of the spin density distribution of the triplet state in the tetrapyrrole macrocycle. The triplet-state properties of Chl are found to be similar, in many respects, to those previously reported for Chl , although some specificities have been highlighted.

Stochastic modeling of macromolecules in solution. I. Relaxation processes.

A framework for the stochastic description of relaxation processes in flexible macromolecules, including dissipative effects, is introduced from an atomistic point of view. Projection-operator techniques are employed to obtain multidimensional Fokker-Planck operators governing the relaxation of internal coordinates and global degrees of freedom and depending upon parameters fully recoverable from classic force fields (energetics) and continuum models (friction tensors). A hierarchy of approaches of different complexity is proposed in this unified context, aimed primarily at the interpretation of magnetic resonance relaxation experiments.

Stochastic modeling of macromolecules in solution. II. Spectral densities.

In Paper I [Polimeno et al., J. Chem.

Changes in the fraction of strongly attached cross bridges in mouse atrophic and hypertrophic muscles as revealed by continuous wave electron paramagnetic resonance.

Electron paramagnetic resonance (EPR), coupled with site-directed spin labeling, has been proven to be a particularly suitable technique to extract information on the fraction of myosin heads strongly bound to actin upon muscle contraction. The approach can be used to investigate possible structural changes occurring in myosin of fiber s altered by diseases and aging. In this work, we labeled myosin at position Cys707, located in the SH1-SH2 helix in the myosin head cleft, with iodoacetamide spin label, a spin label that is sensitive to the reorientational motion of this protein during the ATPase cycle and characterized the biochemical states of the labeled myosin head by means of continuous wave EPR.

DiTe2: Calculating the diffusion tensor for flexible molecules.

We report on an extended hydrodynamic modeling of the friction tensorial properties of flexible molecules including all types of natural, Z-Matrix like, internal coordinates. We implement the new methodology by extending and updating the software DiTe [Barone et al. J.

Evaluating rotation diffusion properties of molecules from short trajectories.

We show that under proper assumptions it is possible to estimate with good precision the principal values of the rotational diffusion tensor of proteins from the analysis of short (up to 2-3 ns) molecular dynamics trajectories. We apply this analysis to a few model cases: three polyalanine peptides (2, 5, and 10 aminoacids), the fragment B3 of protein G (GB3), the bovine pancreatic trypsin inhibitor (BPTI), the hen egg-white lysozyme (LYS), the B1 domain of plexin (PB1), and thrombin. The protocol is based on the analysis of the global angular momentum autocorrelation functions, complementing the standard approach based on rotational autocorrelation functions, which requires much longer trajectories.

Differential Dynamics at Glycosidic Linkages of an Oligosaccharide as Revealed by C NMR Spin Relaxation and Stochastic Modeling.

Among biomolecules, carbohydrates are unique in that not only can linkages be formed through different positions, but the structures may also be branched. The trisaccharide β-d-Glcp-(1→3)[β-d-Glcp-(1→2)]-α-d-Manp-OMe represents a model of a branched vicinally disubstituted structure. A C site-specific isotopologue, with labeling in each of the two terminal glucosyl residues, enabled the acquisition of high-quality C NMR relaxation parameters, T and T, and heteronuclear NOE, with standard deviations of ≤0.

Integrated Computational Approach to the Electron Paramagnetic Resonance Characterization of Rigid 3-Helical Peptides with TOAC Nitroxide Spin Labels.

We address the interpretation, via an integrated computational approach, of the experimental continuous-wave electron paramagnetic resonance (cw-EPR) spectra of a complete set of conformationally highly restricted, stable 3-helical peptides from hexa- to nonamers, each bis-labeled with nitroxide radical-containing TOAC (4-amino-1-oxyl-2,2,6,6-tetramethylpiperidine-4-carboxylic acid) residues. The usefulness of TOAC for this type of analysis has been shown already to be due to its cyclic piperidine side chain, which is rigidly connected to the peptide backbone α-carbon. The TOAC α-amino acids are separated by two, three, four, and five intervening residues.

Decomposition of Proteins into Dynamic Units from Atomic Cross-Correlation Functions.

In this article, we present a clustering method of atoms in proteins based on the analysis of the correlation times of interatomic distance correlation functions computed from MD simulations. The goal is to provide a coarse-grained description of the protein in terms of fewer elements that can be treated as dynamically independent subunits. Importantly, this domain decomposition method does not take into account structural properties of the protein.

Light-Induced Porphyrin-Based Spectroscopic Ruler for Nanometer Distance Measurements.

We present a novel pulsed electron paramagnetic resonance (EPR) spectroscopic ruler to test the performance of a recently developed spin-labeling method based on the photoexcited triplet state (S=1). Four-pulse electron double resonance (PELDOR) experiments are carried out on a series of helical peptides, labeled at the N-terminal end with the porphyrin moiety, which can be excited to the triplet state, and with the nitroxide at various sequence positions, spanning distances in the range 1.8-8 nm.

Loop Electrostatics Asymmetry Modulates the Preexisting Conformational Equilibrium in Thrombin.

Thrombin exists as an ensemble of active (E) and inactive (E*) conformations that differ in their accessibility to the active site. Here we show that redistribution of the E*-E equilibrium can be achieved by perturbing the electrostatic properties of the enzyme. Removal of the negative charge of the catalytic Asp102 or Asp189 in the primary specificity site destabilizes the E form and causes a shift in the 215-217 segment that compromises substrate entrance.

General AMBER Force Field Parameters for Diphenyl Diselenides and Diphenyl Ditellurides.

The General AMBER Force Field (GAFF) has been extended to describe a series of selenium and tellurium diphenyl dichalcogenides. These compounds, besides being eco-friendly catalysts for numerous oxidations in organic chemistry, display peroxidase activity, i.e.

Flexibility at a glycosidic linkage revealed by molecular dynamics, stochastic modeling, and (13)C NMR spin relaxation: conformational preferences of α-L-Rhap-α-(1 → 2)-α-L-Rhap-OMe in water and dimethyl sulfoxide solutions.

The monosaccharide L-rhamnose is common in bacterial polysaccharides and the disaccharide α-L-Rhap-α-(1 → 2)-α-L-Rhap-OMe represents a structural model for a part of Shigella flexneri O-antigen polysaccharides. Utilization of [1'-(13)C]-site-specific labeling in the anomeric position at the glycosidic linkage between the two sugar residues facilitated the determination of transglycosidic NMR (3)JCH and (3)JCC coupling constants. Based on these spin-spin couplings the major state and the conformational distribution could be determined with respect to the ψ torsion angle, which changed between water and dimethyl sulfoxide (DMSO) as solvents, a finding mirrored by molecular dynamics (MD) simulations with explicit solvent molecules.

Selenocysteine oxidation in glutathione peroxidase catalysis: an MS-supported quantum mechanics study.

Glutathione peroxidases (GPxs) are enzymes working with either selenium or sulfur catalysis. They adopted diverse functions ranging from detoxification of H(2)O(2) to redox signaling and differentiation. The relative stability of the selenoenzymes, however, remained enigmatic in view of the postulated involvement of a highly unstable selenenic acid form during catalysis.

On the reliability of NMR relaxation data analyses: a Markov Chain Monte Carlo approach.

The analysis of NMR relaxation data is revisited along the lines of a Bayesian approach. Using a Markov Chain Monte Carlo strategy of data fitting, we investigate conditions under which relaxation data can be effectively interpreted in terms of internal dynamics. The limitations to the extraction of kinetic parameters that characterize internal dynamics are analyzed, and we show that extracting characteristic time scales shorter than a few tens of ps is very unlikely.

Energetics of oxo- and thio-dipeptide formation via amino acid condensation: a systematic computational analysis.

Oxo-dipeptides and thio-dipeptides are built via condensation between couples of amino acids and amino thioacids, the latter with the carbonyl oxygen replaced by an sp(2) sulfur. We explored via in silico methods (PBE0/6-31G(d,p) and PBE0/6-311G(d,p)) all the possible combinations and built 800 dipeptides, whose structures were fully optimized. Maps of condensation energies are presented to highlight optimal partners leading to stable dipeptides and critical situations for which lower stability or instability is predicted in terms of Gibbs reaction free energies.