Hydrogen bonding in the urea dimers and adenine–thymine DNA base pair: anharmonic effects in the intermolecular H-bond and intramolecular H-stretching vibrations

Influence of donor point modifications on the assembly of chalcogen-bonded organic frameworks

Incremental, single-atom substitutions of Se-based chalcogen bond (Ch-bond) donors with stronger donating Te centers were implemented in two new triptycene tris(1,2,5-chalcogenadiazole) tectons. The appreciably more favorable Ch-bonding ability of the Te-based donors promotes assembly of low-density networks and more stable Ch-bonded organic frameworks (ChOFs).

Sequential flipping: the donor-acceptor exchange mechanism in water trimers

The donor-acceptor exchange (DAE) is a significant hydrogen bond network rearrangement (HBNR) mechanism because it can lead to the change of the hydrogen bond direction. In this work, we report a new DAE mechanism found in water trimers that is realized by sequential flipping (SF) of all molecules rather than the well-known proton transfer (PT) process. Meanwhile, the SF process has a much smaller potential barrier (0.262 eV) than the previously predicted collective rotation process (about 1.7 eV), implying that the SF process is the main flipping process that can lead to DAE. Importantly, high-precision ab initio calculations show that SF-DAE can make the water ring to show a clear chiral difference from PT-DAE, which brings the prospect of distinguishing the two confusing processes based on circular dichroism spectra. The reaction rate analysis including quantum tunneling indicates an obvious temperature-dependent competitive relationship between the SF and PT processes; specifically, the SF process dominates above 65 K, while the PT process dominates below 65 K. Therefore, in most cases, the contribution for DAE mainly comes from the flipping process, rather than the PT process as previously thought. Our work enriches the understanding of the DAE mechanism in water trimers and provides a piece of the jigsaw that has been sought for the HBNR mechanism.

SARS-CoV-2 Spike Protein Unlikely to Bind to Integrins via the Arg-Gly-Asp (RGD) Motif of the Receptor Binding Domain: Evidence From Structural Analysis and Microscale Accelerated Molecular Dynamics

The Receptor Binding Domain (RBD) of SARS-CoV-2 virus harbors a sequence of Arg-Gly-Asp tripeptide named RGD motif, which has also been identified in extracellular matrix proteins that bind integrins as well as other disintegrins and viruses. Accordingly, integrins have been proposed as host receptors for SARS-CoV-2. However, given that the microenvironment of the RGD motif imposes a structural hindrance to the protein-protein association, the validity of this hypothesis is still uncertain. Here, we used normal mode analysis, accelerated molecular dynamics microscale simulation, and protein-protein docking to investigate the putative role of RGD motif of SARS-CoV-2 RBD for interacting with integrins. We found, that neither RGD motif nor its microenvironment showed any significant conformational shift in the RBD structure. Highly populated clusters of RBD showed no capability to interact with the RGD binding site in integrins. The free energy landscape revealed that the RGD conformation within RBD could not acquire an optimal geometry to allow the interaction with integrins. In light of these results, and in the event where integrins are confirmed to be host receptors for SARS-CoV-2, we suggest a possible involvement of other residues to stabilize the interaction.

Supramolecular co-assembly of water-soluble nucleobase-containing copolymers: bioinspired synthetic platforms towards new biomimetic materials

This study presents the development of new co-assembled copolymers architectures at physiological pH (pH=7.4) formed via H-bonds between complementary nucleobase-containing copolymers. Well-defined hydrophilic copolymers were synthesised by RAFT polymerisation: statistical...

Quantum and classical effects in DNA point mutations: Watson-Crick tautomerism in AT and GC base pairs

Proton transfer along the hydrogen bonds of DNA can lead to the creation of short-lived, but biologically relevant point mutations that can further lead to gene mutation and, potentially, cancer. In this work, the energy landscape of the canonical A-T and G-C base pairs (standard, amino-keto) to tautomeric A*-T* and G*-C* (non-standard, imino-enol) Watson-Crick DNA base pairs is modelled with density functional theory and machine-learning nudge-elastic band methods. We calculate the energy barriers and tunnelling rates of hydrogen transfer between and within each base monomer (A, T, G and C). We show that the role of tunnelling in A-T tautomerisation is statistically unlikely due to the presence of a small reverse reaction barrier. On the contrary, the thermal populations of the G*-C* point mutation could be non-trivial and propagate through the replisome. For the direct intramolecular transfer, the reaction is hindered by a substantial energy barrier. However, our calculations indicate that tautomeric bases in their monomeric form have remarkably long lifetimes.

Pinacolone-Alcohol Gas-Phase Solvation Balances as Experimental Dispersion Benchmarks

The influence of distant London dispersion forces on the docking preference of alcohols of different size between the two lone electron pairs of the carbonyl group in pinacolone was explored by infrared spectroscopy of the OH stretching fundamental in supersonic jet expansions of 1:1 solvate complexes. Experimentally, no pronounced tendency of the alcohol to switch from the methyl to the bulkier tert-butyl side with increasing size was found. In all cases, methyl docking dominates by at least a factor of two, whereas DFT-optimized structures suggest a very close balance for the larger alcohols, once corrected by CCSD(T) relative electronic energies. Together with inconsistencies when switching from a C4 to a C5 alcohol, this points at deficiencies of the investigated B3LYP and in particular TPSS functionals even after dispersion correction, which cannot be blamed on zero point energy effects. The search for density functionals which describe the harmonic frequency shift, the structural change and the energy difference between the docking isomers of larger alcohols to unsymmetric ketones in a satisfactory way is open.

Proton transfer in guanine-cytosine base pair analogues studied by NMR spectroscopy and PIMD simulations

It has been hypothesised that proton tunnelling between paired nucleobases significantly enhances the formation of rare tautomeric forms and hence leads to errors in DNA replication. Here, we study nuclear quantum effects (NQEs) using deuterium isotope-induced changes of nitrogen NMR chemical shifts in a model base pair consisting of two tautomers of isocytosine, which form hydrogen-bonded dimers in the same way as the guanine-cytosine base pair. Isotope effects in NMR are consequences of NQEs, because ro-vibrational averaging of different isotopologues gives rise to different magnetic shielding of the nuclei. The experimental deuterium-induced chemical shift changes are compared with those calculated by a combination of path integral molecular dynamics (PIMD) simulations with DFT calculations of nuclear shielding. These calculations can directly link the observable isotope-induced shifts with NQEs. A comparison of the deuterium-induced changes of 15N chemical shifts with those predicted by PIMD simulations shows that inter-base proton transfer reactions do not take place in this system. We demonstrate, however, that NMR isotope shifts provide a unique possibility to study NQEs and to evaluate the accuracy of the computational methods used for modelling quantum effects in molecules. Calculations based on the PBE functional from the general-gradient-approximation family provided significantly worse predictions of deuterium isotope shifts than those with the hybrid B3LYP functional.

Modelling proton tunnelling in the adenine–thymine base pair

Environment-assisted quantum tunnelling is unlikely to play a significant role in adenine–thymine tautomerisation.

Evaluation of coupling terms between intra- and intermolecular vibrations in coarse-grained normal-mode analysis: does a stronger acid make a stiffer hydrogen bond?

Using theory of harmonic normal-mode vibration analysis, we developed a procedure for evaluating the anisotropic stiffness of intermolecular forces. Our scheme for coarse-graining of molecular motions is modified so as to account for intramolecular vibrations in addition to relative translational/rotational displacement. We applied this new analytical scheme to four carboxylic acid dimers, for which coupling between intra- and intermolecular vibrations is crucial for determining the apparent stiffness of the intermolecular double hydrogen bond. The apparent stiffness constant was analyzed on the basis of a conjunct spring model, which defines contributions from true intermolecular stiffness and molecular internal stiffness. Consequently, the true intermolecular stiffness was in the range of 43-48 N m(-1) for all carboxylic acids studied, regardless of the molecules' acidity. We concluded that the difference in the apparent stiffness can be attributed to differences in the internal stiffness of the respective molecules.

Localization and anharmonicity of the vibrational modes for GC Watson-Crick and Hoogsteen base pairs

We present an ab initio study of the vibrational properties of cytosine and guanine in the Watson-Crick and Hoogsteen base pair configurations. The results are obtained by using two different implementations of the DFT method. We assign the vibrational frequencies to cytosine or to guanine using the vibrational density of states. Next, we investigate the importance of anharmonic corrections for the vibrational modes. In particular, the unusual anharmonic effect of the H(+) vibration in the case of the Hoogsteen base pair configuration is discussed.