Singular Interface Dynamics of the SARS-CoV-2 Delta Variant Explained with Contact Perturbation Analysis.

Emerging SARS-CoV-2 variants raise concerns about our ability to withstand the Covid-19 pandemic, and therefore, understanding mechanistic differences of those variants is crucial. In this study, we investigate disparities between the SARS-CoV-2 wild type and five variants that emerged in late 2020, focusing on the structure and dynamics of the spike protein interface with the human angiotensin-converting enzyme 2 (ACE2) receptor, by using crystallographic structures and extended analysis of microsecond molecular dynamics simulations. Dihedral angle principal component analysis (PCA) showed the strong similarities in the spike receptor binding domain (RBD) dynamics of the Alpha, Beta, Gamma, and Delta variants, in contrast with those of WT and Epsilon. Dynamical perturbation networks and contact PCA identified the peculiar interface dynamics of the Delta variant, which cannot be directly imputable to its specific L452R and T478K mutations since those residues are not in direct contact with the human ACE2 receptor. Our outcome shows that in the Delta variant the L452R and T478K mutations act synergistically on neighboring residues to provoke drastic changes in the spike/ACE2 interface; thus a singular mechanism of action eventually explains why it dominated over preceding variants.

Excited States Computation of Models of Phenylalanine Protein Chains: TD-DFT and Composite CC2/TD-DFT Protocols.

The present benchmark calculations testify to the validity of time-dependent density functional theory (TD-DFT) when exploring the low-lying excited states potential energy surfaces of models of phenylalanine protein chains. Among three functionals suitable for systems exhibiting charge-transfer excited states, LC-ωPBE, CAM-B3LYP, and ωB97X-D, which were tested on a reference peptide system, we selected the ωB97X-D functional, which gave the best results compared to the approximate coupled-cluster singles and doubles (CC2) method. A quantitative agreement for both the geometrical parameters and the vibrational frequencies was obtained for the lowest singlet excited state (a ππ* state) of the series of capped peptides. In contrast, only a qualitative agreement was met for the corresponding adiabatic zero-point vibrational energy (ZPVE)-corrected excitation energies. Two composite protocols combining CC2 and DFT/TD-DFT methods were then developed to improve these calculations. Both protocols substantially reduced the error compared to CC2 and experiment, and the best of both even led to results of CC2 quality at a lower cost, thus providing a reliable alternative to this method for very large systems.

Ion Pair Supramolecular Structure Identified by ATR-FTIR Spectroscopy and Simulations in Explicit Solvent*.

The present work uses ATR-FTIR spectroscopy assisted by simulations in explicit solvent and frequency calculations to investigate the supramolecular structure of carboxylate alkali-metal ion pairs in aqueous solutions. ATR-FTIR spectra in the 0.25-4.0 M concentration range displayed cation-specific behaviors, which enabled the measurement of the appearance concentration thresholds of contact ion pairs between 1.9 and 2.6 M depending on the cation. Conformational explorations performed using a non-local optimization method associated to a polarizable force-field (AMOEBA), followed by high quantum chemistry level (RI-B97-D3/dhf-TZVPP) optimizations, mode-dependent scaled harmonic frequency calculations and electron density analyses, were used to identify the main supramolecular structures contributing to the experimental spectra. A thorough analysis enables us to reveal the mechanisms responsible for the spectroscopic sensitivity of the carboxylate group and the respective role played by the cation and the water molecules, highlighting the necessity of combining advanced experimental and theoretical techniques to provide a fair and accurate description of ion pairing.

Model protein excited states: MRCI calculations with large active spaces vs CC2 method.

Benchmarking calculations on excited states of models of phenylalanine protein chains are presented to assess the ability of alternative methods to the standard and most commonly used multiconfigurational wave function-based method, the complete active space self-consistent field (CASSCF), in recovering the non-dynamical correlation for systems that become not affordable by the CASSCF. The exploration of larger active spaces beyond the CASSCF limit is benchmarked through three strategies based on the reduction in the number of determinants: the restricted active space self-consistent field, the generalized active space self-consistent field (GASSCF), and the occupation-restricted multiple active space (ORMAS) schemes. The remaining dynamic correlation effects are then added by the complete active space second-order perturbation theory and by the multireference difference dedicated configuration interaction methods. In parallel, the approximate second-order coupled cluster (CC2), already proven to be successful for small building blocks of model proteins in one of our previous works [Ben Amor et al., J. Chem. Phys. 148, 184105 (2018)], is investigated to assess its performances for larger systems. Among the different alternative strategies to CASSCF, our results highlight the greatest efficiency of the GASSCF and ORMAS schemes in the systematic reduction of the configuration interaction expansion without loss of accuracy in both nature and excitation energies of both singlet ππ* and nπ* CO excited states with respect to the equivalent CASSCF calculations. Guidelines for an optimum applicability of this scheme to systems requiring active spaces beyond the complete active space limit are then proposed. Finally, the extension of the CC2 method to such large systems without loss of accuracy is demonstrated, highlighting the great potential of this method to treat accurately excited states, mainly single reference, of very large systems.

Tinker-HP: Accelerating Molecular Dynamics Simulations of Large Complex Systems with Advanced Point Dipole Polarizable Force Fields Using GPUs and Multi-GPU Systems.

We present the extension of the Tinker-HP package (Lagardère, Chem. Sci. 2018, 9, 956-972) to the use of Graphics Processing Unit (GPU) cards to accelerate molecular dynamics simulations using polarizable many-body force fields. The new high-performance module allows for an efficient use of single- and multiple-GPU architectures ranging from research laboratories to modern supercomputer centers. After detailing an analysis of our general scalable strategy that relies on OpenACC and CUDA, we discuss the various capabilities of the package. Among them, the multiprecision possibilities of the code are discussed. If an efficient double precision implementation is provided to preserve the possibility of fast reference computations, we show that a lower precision arithmetic is preferred providing a similar accuracy for molecular dynamics while exhibiting superior performances. As Tinker-HP is mainly dedicated to accelerate simulations using new generation point dipole polarizable force field, we focus our study on the implementation of the AMOEBA model. Testing various NVIDIA platforms including 2080Ti, 3090, V100, and A100 cards, we provide illustrative benchmarks of the code for single- and multicards simulations on large biosystems encompassing up to millions of atoms. The new code strongly reduces time to solution and offers the best performances to date obtained using the AMOEBA polarizable force field. Perspectives toward the strong-scaling performance of our multinode massive parallelization strategy, unsupervised adaptive sampling and large scale applicability of the Tinker-HP code in biophysics are discussed. The present software has been released in phase advance on GitHub in link with the High Performance Computing community COVID-19 research efforts and is free for Academics (see https://github.com/TinkerTools/tinker-hp).

Tinker-HP : Accelerating Molecular Dynamics Simulations of Large Complex Systems with Advanced Point Dipole Polarizable Force Fields using GPUs and Multi-GPUs systems.

We present the extension of the Tinker-HP package (Lagard\`ere et al., Chem. Sci., 2018,9, 956-972) to the use of Graphics Processing Unit (GPU) cards to accelerate molecular dynamics simulations using polarizable many-body force fields. The new high-performance module allows for an efficient use of single- and multi-GPU architectures ranging from research laboratories to modern supercomputer centers. After detailing an analysis of our general scalable strategy that relies on OpenACC and CUDA, we discuss the various capabilities of the package. Among them, the multi-precision possibilities of the code are discussed. If an efficient double precision implementation is provided to preserve the possibility of fast reference computations, we show that a lower precision arithmetic is preferred providing a similar accuracy for molecular dynamics while exhibiting superior performances. As Tinker-HP is mainly dedicated to accelerate simulations using new generation point dipole polarizable force field, we focus our study on the implementation of the AMOEBA model. Testing various NVIDIA platforms including 2080Ti, 3090, V100 and A100 cards, we provide illustrative benchmarks of the code for single- and multi-cards simulations on large biosystems encompassing up to millions of atoms. The new code strongly reduces time to solution and offers the best performances to date obtained using the AMOEBA polarizable force field. Perspectives toward the strong-scaling performance of our multi-node massive parallelization strategy, unsupervised adaptive sampling and large scale applicability of the Tinker-HP code in biophysics are discussed. The present software has been released in phase advance on GitHub in link with the High Performance Computing community COVID-19 research efforts and is free for Academics (see https://github.com/TinkerTools/tinker-hp).

Identification of ion pairs in solution by IR spectroscopy: crucial contributions of gas phase data and simulations.

In a context where structure elucidation of ion pairs in solution remains a contemporary challenge, this work explores an original approach where accurate gas phase spectroscopic data are used to refine high level quantum chemistry calculations of ion pairs in solution, resulting in an unprecedented level of accuracy in vibrational frequency prediction. First, gas phase studies focus on a series of isolated contact ion pairs (M+, Ph-CH2-COO-, with M = Li, Na, K, Rb, Cs) for which conformer-selective IR spectra in the CO2- stretch region are recorded. These experiments reveal the interactions at play in isolated contact ion pairs, and provide vibrational frequencies enabling us to assess the accuracy of the theoretical approach used, i.e., mode-dependent scaled harmonic frequency calculations at the RI-B97-D3/dhf-TZVPP level. This level of calculation is then employed on large water clusters embedding either a free acetate ion or its contact or solvent-shared pairs with a sodium cation in order to simulate the individual vibrational spectra of these species in solution. This study shows that the stretching modes of carboxylate are sensitive to both solvent-shared and contact ion pair formation. FTIR spectra of solutions of increasing concentrations indeed reveal several spectral changes consistent with the presence of specific types of solvent-shared and contact ion pairs. By providing relevant guidelines for the interpretation of solution phase IR spectra, this work illustrates the potential of the approach for the elucidation of supramolecular structures in electrolyte solutions.

Photophysical properties of ruthenium(II) polypyridyl DNA intercalators: effects of the molecular surroundings investigated by theory.

The environmental effects on the structural and photophysical properties of [Ru(L)2 (dppz)](2+) complexes (L=bpy=2,2'-bipyridine, phen=1,10-phenanthroline, tap=1,4,5,8-tetraazaphenanthrene; dppz=dipyrido[3,3-a:2',3'-c]phenazine), used as DNA intercalators, have been studied by means of DFT, time-dependent DFT, and quantum mechanics/molecular mechanics calculations. The electronic characteristics of the low-lying triplet excited states in water, acetonitrile, and DNA have been investigated to decipher the influence of the environment on the luminescent behavior of this class of molecules. The lowest triplet intra-ligand (IL) excited state calculated at λ≈800 nm for the three complexes and localized on the dppz ligand is not very sensitive to the environment and is available for electron transfer from a guanine nucleobase. Whereas the lowest triplet metal-to-ligand charge-transfer ((3)MLCT) states remain localized on the ancillary ligand (tap) in [Ru(tap)2 (dppz)](2+), regardless of the environment, their character is drastically modified in the other complexes [Ru(phen)2 (dppz)](2+) and [Ru(bpy)2 (dppz)](2+) upon going from acetonitrile (MLCTdppz/phen or MLCTdppz/bpy) to water (MLCTdppz) and DNA (MLCTphen and MLCTbpy). The change in the character of the low-lying (3) MLCT states accompanying nuclear relaxation in the excited state controls the emissive properties of the complexes in water, acetonitrile, and DNA. The light-switching effect has been rationalized on the basis of environment-induced control of the electronic density distributed in the lowest triplet excited states.

Toward the understanding of the environmental effects on core ionizations.

Experimental X-ray absorption spectra are extensively used to determine electronic structure of small molecules but remain difficult to exploit for proteins due to the large number of peaks within their spectra. For such complex systems, theoretical tools like quantum mechanics/molecular mechanics methodology can greatly ease the assignment of the spectra. This study presents a systematic methodology to evaluate core-ionization energies (E(ion)) in proteins with the help of the asymptotic projection approach (Glushkov and Tsaune, Z. Vichislit. Matem. Mat. Fiz. 1985, 25, 298; Glushkov, Chem. Phys. Lett. 1997, 273, 122; Glushkov, Chem. Phys. Lett. 1998, 287, 189; Glushkov, J. Math. Chem. 2002, 31, 91; Glushkov, Opt. Spectrosc. 2002, 93, 15). An in-depth inspection of E(ion) of systems of increasing complexity is considered, going from amino acids to polyglycine and to glycine in human serum albumin (HSA). Computational analysis can help to better understand experimental data and to discriminate environmental effects by tracing them back to individual and collective electrostatic contributions. In the present work, it was found that E(ion) of alpha carbon of glycine residues in HSA ranges from 285 to 295 eV depending on their surroundings.

UV-vis absorption spectrum of a novel Ru(II) complex intercalated in DNA: [Ru(2,2'-bipy)(dppz)(2,2'-ArPy)]⁺.

The synthesis of a new Ru(II) complex is reported. Its absorption spectrum when interacting with DNA in water was calculated at the hybrid quantum mechanics molecular mechanics level of theory and compared with experimental data. The vertical transitions were computed using time-dependent density functional theory in the linear response approximation. The complex and its environment were treated at the quantum mechanical and molecular mechanical levels, respectively. The effects of the environment were investigated in detail and conveniently classified into electrostatic and polarization effects. The latter were modeled using the computationally inexpensive "electronic response of the surroundings" method. It was found that the main features of the experimental spectrum are nicely reproduced by the theoretical calculations. Moreover, analysis of the most intense transitions utilizing the natural transition orbital formalism revealed important insights into their nature and their potential role in the irreversible oxidation of DNA, a phenomenon that could be relevant in the field of cancer therapy.

A QM/MM study of the absorption spectrum of harmane in water solution and interacting with DNA: the crucial role of dynamic effects.

We present a time-dependent density functional theory computation of the absorption spectra of one ß-carboline system: the harmane molecule in its neutral and cationic forms. The spectra are computed in aqueous solution. The interaction of cationic harmane with DNA is also studied. In particular, the use of hybrid quantum mechanics/molecular mechanics methods is discussed, together with its coupling to a molecular dynamics strategy to take into account dynamic effects of the environment and the vibrational degrees of freedom of the chromophore. Different levels of treatment of the environment are addressed starting from purely mechanical embedding to electrostatic and polarizable embedding. We show that a static description of the spectrum based on equilibrium geometry only is unable to give a correct agreement with experimental results, and dynamic effects need to be taken into account. The presence of two stable noncovalent interaction modes between harmane and DNA is also presented, as well as the associated absorption spectrum of harmane cation.

Insight on the interaction of polychlorobiphenyl with nucleic acid-base.

The interaction between one polychlorobiphenyl (3,3',4,4',-tetrachlorobiphenyl, coded PCB77) and the four DNA nucleic acid-base is studied by means of quantum mechanics calculations in stacked conformations. It is shown that even if the intermolecular dispersion energy is the largest component of the total interaction energy, some other contributions play a non negligible role. In particular the electrostatic dipole-dipole interaction and the charge transfer from the nucleobase to the PCB are responsible for the relative orientation of the monomers in the complexes. In addition, the charge transfer tends to flatten the PCB, which could therefore intercalate more easily between DNA base pairs. From these seminal results, we predict that PCB could intercalate completely between two base pairs, preferably between Guanine:Cytosine pairs.

Spectral properties of polypyridyl ruthenium complexes intercalated in DNA: theoretical insights into the surrounding effects of [Ru(dppz)(bpy)2]2+.

The UV/Visible absorption properties of a polypyridyl ruthenium complex upon intercalation on DNA are studied at the mixed quantum mechanics molecular mechanics level of theory. Vertical excitation transitions are computed by time dependent density functional theory. Particular emphasis is put on the different levels at which the macromolecular environment is treated, and in particular on the analysis of the effect of mechanical, electrostatic and polarizable embedding. We show that with the highest level of theory the experimental absorption wavelengths are reproduced with a difference of only 2 or 3 nm for the low energy bands. The systematic analysis of the individual vertical transitions allows us to get much more insights into the role played by the environment, in particular, in metal to ligand and intra ligand charge transfer transitions that can lead to the production of DNA oxidative lesions exploitable in phototherapy.

Improved Treatment of Surrounding Effects: UV/vis Absorption Properties of a Solvated Ru(II) Complex.

The UV/vis and circular-dichroism spectra of a bis-bipyridinyl ruthenium complex are computed at the density functional theory level and the time dependent density functional level of theory. The effects of the solvent, here water, have been taken into account, by polarizable continuum methods and by a hybrid quantum-mechanics/molecular-mechanics approach combined with molecular dynamics. The effects of the solvent have been decomposed in geometric, electrostatic, and polarization of the environment. The principal transitions have been analyzed by means of natural transition orbitals.