Systematic QM/MM Study for Predicting 31P NMR Chemical Shifts of Adenosine Nucleotides in Solution and Stages of ATP Hydrolysis in a Protein Environment.

NMR (nuclear magnetic resonance) spectroscopy allows for important atomistic insights into the structure and dynamics of biological macromolecules; however, reliable assignments of experimental spectra are often difficult. Herein, quantum mechanical/molecular mechanical (QM/MM) calculations can provide crucial support. A major problem for the simulations is that experimental NMR signals are time-averaged over much longer time scales, and since computed chemical shifts are highly sensitive to local changes in the electronic and structural environment, sufficiently large averages over representative structural ensembles are essential. This entails high computational demands for reliable simulations. For NMR measurements in biological systems, a nucleus of major interest is 31P since it is both highly present (e.g., in nucleic acids) and easily observable. The focus of our present study is to develop a robust and computationally cost-efficient framework for simulating 31P NMR chemical shifts of nucleotides. We apply this scheme to study the different stages of the ATP hydrolysis reaction catalyzed by p97. Our methodology is based on MM molecular dynamics (MM-MD) sampling, followed by QM/MM structure optimizations and NMR calculations. Overall, our study is one of the most comprehensive QM-based 31P studies in a protein environment and the first to provide computed NMR chemical shifts for multiple nucleotide states in a protein environment. This study sheds light on a process that is challenging to probe experimentally and aims to bridge the gap between measured and calculated NMR spectroscopic properties.

Sustained Solar H2 Evolution from a Thiazolo[5,4-d]thiazole-Bridged Covalent Organic Framework and Nickel-Thiolate Cluster in Water.

Solar hydrogen (H2) evolution from water utilizing covalent organic frameworks (COFs) as heterogeneous photosensitizers has gathered significant momentum by virtue of the COFs' predictive structural design, long-range ordering, tunable porosity, and excellent light-harvesting ability. However, most photocatalytic systems involve rare and expensive platinum as the co-catalyst for water reduction, which appears to be the bottleneck in the development of economical and environmentally benign solar H2 production systems. Herein, we report a simple, efficient, and low-cost all-in-one photocatalytic H2 evolution system composed of a thiazolo[5,4-d]thiazole-linked COF (TpDTz) as the photoabsorber and an earth-abundant, noble-metal-free nickel-thiolate hexameric cluster co-catalyst assembled in situ in water, together with triethanolamine (TEoA) as the sacrificial electron donor. The high crystallinity, porosity, photochemical stability, and light absorption ability of the TpDTz COF enables excellent long-term H2 production over 70 h with a maximum rate of 941 µmol h-1 g-1, turnover number TONNi > 103, and total projected TONNi > 443 until complete catalyst depletion. The high H2 evolution rate and TON, coupled with long-term photocatalytic operation of this hybrid system in water, surpass those of many previously known organic dyes, carbon nitride, and COF-sensitized photocatalytic H2O reduction systems. Furthermore, we gather unique insights into the reaction mechanism, enabled by a specifically designed continuous-flow system for non-invasive, direct H2 production rate monitoring, providing higher accuracy in quantification compared to the existing batch measurement methods. Overall, the results presented here open the door toward the rational design of robust and efficient earth-abundant COF-molecular co-catalyst hybrid systems for sustainable solar H2 production in water.

Calculating free energies from the vibrational density of states function: Validation and critical assessment.

We explore and show the usefulness of the density of states function for computing vibrational free energies and free energy differences between small systems. Therefore, we compare this density of states integration method (DSI) to more established schemes such as Bennett's Acceptance Ratio method (BAR), the Normal Mode Analysis (NMA), and the Quasiharmonic Analysis (QHA). The strengths and shortcomings of all methods are highlighted with three numerical examples. Furthermore, the free energy of the ionization of ammonia and the mutation from serine to cysteine are computed using extensive ab initio molecular dynamics simulations. We conclude that DSI improves upon the other frequency-based methods (NMA and QHA) regarding the treatment of anharmonicity and yielding results comparable to BAR in all cases without the need for alchemical transformations. Low-frequency modes lead to larger errors indicating that long simulation times might be required for larger systems. In addition, we introduce the use of DSI for the localization of the vibrational free energy to specific atoms or residues, leading to insights into the underlying process, a unique feature that is only offered by this method.

C-H···O Hydrogen Bonding. The Prototypical Methane-Formaldehyde System: A Critical Assessment.

Distinguishing the functionality of C-H···O hydrogen bonds (HBs) remains challenging, because their properties are difficult to quantify reliably. Herein, we present a study of the model methane-formaldehyde complex (MFC). Six stationary points on the MFC potential energy surface (PES) were obtained at the CCSD(T)/ANO2 level. The CCSDT(Q)/CBS interaction energies of the conformers range from only -1.12 kcal mol-1 to -0.33 kcal mol-1, denoting a very flat PES. Notably, only the lowest energy stationary point (MFC1) corresponds to a genuine minimum, whereas all other stationary points-including the previously studied ideal case of ae(C-H···O) = 180°-exhibit some degree of freedom that leads to MFC1. Despite the flat PES, we clearly see that the HB properties of MFC1 align with those of the prototypical water dimer O-H···O HB. Each HB property generally becomes less prominent in the higher-energy conformers. Only the MFC1 conformer prominently exhibits (1) elongated C-H donor bonds, (2) attractive C-H···O═C interactions, (3) n(O) → σ*(C-H) hyperconjugation, (4) critical points in the electron density from Bader's method and from the noncovalent interactions method, (5) positively charged donor hydrogen, and (6) downfield NMR chemical shifts and nonzero 2J(CM-HM···OF) coupling constants. Based on this research, some issues merit further study. The flat PES hinders reliable determinations of the HB-induced shifts of the C-H stretches; a similarly difficult challenge is observed for the experiment. The role of charge transfer in HBs remains an intriguing open question, although our BLW and NBO computations suggest that it is relevant to the C-H···O HB geometries. These issues notwithstanding, the prominence of the HB properties in MFC1 serves as clear evidence that the MFC is predominantly bound by a C-H···O HB.

Exploiting Noncovalent Interactions in an Imine-Based Covalent Organic Framework for Quercetin Delivery.

Covalent organic frameworks (COFs) are a new class of nanoporous polymeric vector showing promise as drug-delivery vehicles with high loading capacity and biocompatibility. The interaction between the carrier and the cargo is specifically tailored on a molecular level by H-bonding. Cell-proliferation studies indicate higher efficacy of the drug in cancer cells by nanocarrier delivery mediated by the COF.

Structural, Biochemical, and Computational Studies Reveal the Mechanism of Selective Aldehyde Dehydrogenase 1A1 Inhibition by Cytotoxic Duocarmycin Analogues.

Analogues of the natural product duocarmycin bearing an indole moiety were shown to bind aldehyde dehydrogenase 1A1 (ALDH1A1) in addition to DNA, while derivatives without the indole solely addressed the ALDH1A1 protein. The molecular mechanism of selective ALDH1A1 inhibition by duocarmycin analogues was unraveled through cocrystallization, mutational studies, and molecular dynamics simulations. The structure of the complex shows the compound embedded in a hydrophobic pocket, where it is stabilized by several crucial π-stacking and van der Waals interactions. This binding mode positions the cyclopropyl electrophile for nucleophilic attack by the noncatalytic residue Cys302, thereby resulting in covalent attachment, steric occlusion of the active site, and inhibition of catalysis. The selectivity of duocarmycin analogues for ALDH1A1 is unique, since only minor alterations in the sequence of closely related protein isoforms restrict compound accessibility.

A base-independent repair mechanism for DNA glycosylase--no discrimination within the active site.

The ubiquitous occurrence of DNA damages renders its repair machinery a crucial requirement for the genomic stability and the survival of living organisms. Deficiencies in DNA repair can lead to carcinogenesis, Alzheimer, or Diabetes II, where increased amounts of oxidized DNA bases have been found in patients. Despite the highest mutation frequency among oxidized DNA bases, the base-excision repair process of oxidized and ring-opened guanine, FapydG (2,6-diamino-4-hydroxy-5-formamidopyrimidine), remained unclear since it is difficult to study experimentally. We use newly-developed linear-scaling quantum-chemical methods (QM) allowing us to include up to 700 QM-atoms and achieving size convergence. Instead of the widely assumed base-protonated pathway we find a ribose-protonated repair mechanism which explains experimental observations and shows strong evidence for a base-independent repair process. Our results also imply that discrimination must occur during recognition, prior to the binding within the active site.

Molecular tweezers with varying anions: a comparative study.

Selective binding of the phosphate-substituted molecular tweezer 1a to protein lysine residues was suggested to explain the inhibition of certain enzymes and the aberrant aggregation of amyloid petide Aß42 or α-synuclein, which are assumed to be responsible for Alzheimer's and Parkinson's disease, respectively. In this work we systematically investigated the binding of four water-soluble tweezers 1a-d (substituted by phosphate, methanephosphonate, sulfate, or O-methylenecarboxylate groups) to amino acids and peptides containing lysine or arginine residues by using fluorescence spectroscopy, NMR spectroscopy, and isothermal titration calorimetry (ITC). The comparison of the experimental results with theoretical data obtained by a combination of QM/MM and ab initio(1)H NMR shift calculations provides clear evidence that the tweezers 1a-c bind the amino acid or peptide guest molecules by threading the lysine or arginine side chain through the tweezers' cavity, whereas in the case of 1d the guest molecule is preferentially positioned outside the tweezer's cavity. Attractive ionic, CH-π, and hydrophobic interactions are here the major binding forces. The combination of experiment and theory provides deep insight into the host-guest binding modes, a prerequisite to understanding the exciting influence of these tweezers on the aggregation of proteins and the activity of enzymes.

Effects of terminal functional groups on the stability of the polyproline II structure: a combined experimental and theoretical study.

The conformational stability of the polyproline II (PPII) helix with respect to the functional groups at the C- and N-termini was examined both experimentally and theoretically. Oligoprolines AcN-[Pro](12)-CONH(2) (1), HN-[Pro](12)-CONH(2) (2), AcN-[Pro](12)-CO(2)H (3), and HN-[Pro](12)-CO(2)H (4) with charged and capped termini served as model compounds, and the relative ease with which they switch from the PPII to the polyproline I (PPI) helix was used as a measure to analyze their conformational stabilities. CD spectroscopic studies demonstrate that a positively charged N-terminus and a negatively charged C-terminus destabilize the PPII helix and favor the PPI helix, whereas capped termini favor the PPII over the PPI helix. These experimental findings are supported by the energy differences between the PPII and PPI helices of oligoprolines 1-4 computed by ab initio methods including electron-correlation effects (second-order Møller-Plesset perturbation theory, MP2). Furthermore, these quantum-chemical calculations show that differences in charge-dipole interactions are responsible for the experimentally and computationally observed relative stabilities. Although these electrostatic interactions between the terminal charges and the amide dipoles stabilize both helices, they are significantly stronger in the PPI helix where the amide bonds are oriented almost linear to the helix axis as compared to the PPII helix in which the amides are nearly perpendicular to the axis. Moreover, we demonstrate that a negative charge at the C-terminus has a more pronounced effect on the relative stability as compared to a positive charge at the N-terminus due to destabilization of the PPII helix by repulsive interaction between the C-terminal carboxylate with the neighboring amide bond. Studies at different pH values verified the electrostatic nature of the observed effects and demonstrate how changes in the protonation state can be used to deliberately stabilize the PPII helix over the PPI helix or vice versa.

A convergence study of QM/MM isomerization energies with the selected size of the QM region for peptidic systems.

A systematic study of the convergence of QM/MM results with respect to the chosen size of the QM region is presented for two examples of peptidic systems. For this purpose, we increased the QM region to up to 1637 atoms at the HF/SVP and 383 atoms at the SOS-AO-MP2/6-31G** level. While the convergence behavior is almost independent of the chosen method and basis set, the study clearly shows that for the considered proton-transfer energy the QM/MM treatment leads to a significantly faster convergence than the pure QM treatment. This behavior can be rationalized by the fair description of the surrounding of the active center using MM methods, even though the MM description of the active center is not adequate in our present case. At the same time, the observed convergence is quite insensitive to a variation of charge surroundings in the chosen model peptides. Although the QM/MM results do converge much quicker with the system size than the pure QM ones, the data show that even for the chosen simple model systems about 150-300 QM atoms are needed to achieve accuracies in the order of 10 kJ/mol and about 300-1000 atoms for an accuracy of 2 kJ/mol with respect to a convergence with the QM-region size.

An arginine switch in the species B adenovirus knob determines high-affinity engagement of cellular receptor CD46.

Adenoviruses (Ads) are icosahedral, nonenveloped viruses with a double-stranded DNA genome. The 51 known Ad serotypes exhibit profound variations in cell tropism and disease types. The number of observed Ad infections is steadily increasing, sometimes leading to fatal outcomes even in healthy individuals. Species B Ads can cause kidney infections, hemorrhagic cystitis, and severe respiratory infections, and most of them use the membrane cofactor protein CD46 as a cellular receptor. The crystal structure of the human Ad type 11 (Ad11) knob complexed with CD46 is known; however, the determinants of CD46 binding in related species B Ads remain unclear. We report here a structural and functional analysis of the Ad11 knob, as well as the Ad7 and Ad14 knobs, which are closely related in sequence to the Ad11 knob but have altered CD46-binding properties. The comparison of the structures of the three knobs, which we determined at very high resolution, provides a platform for understanding these differences and allows us to propose a mechanism for productive high-affinity engagement of CD46. At the center of this mechanism is an Ad knob arginine that needs to switch its orientation in order to engage CD46 with high affinity. Quantum chemical calculations showed that the CD46-binding affinity of Ad11 is significantly higher than that of Ad7. Thus, while Ad7 and Ad14 also bind CD46, the affinity and kinetics of these interactions suggest that these Ads are unlikely to use CD46 productively. The proposed mechanism is likely to determine the receptor usage of all CD46-binding Ads.

Hierarchical self-assembly of aminopyrazole peptides into nanorosettes in water.

Self-association of aminopyrazole peptide hybrid 1 leads to stacked nanorosettes. This remarkable, well-ordered structure obeys the laws of nucleic acid self-assembly. In a strictly hierarchical process, formation of aminopyrazole "base" triplets via a hydrogen bond network is accompanied by pi-stacking with a second rosette and final dimerization of two double rosettes to a four-layer superstructure, stabilized by a six-fold half-crown alkylammonium lock. The final complex is soluble in organic as well as in aqueous solution. It was characterized in the solid state by X-ray crystallography, in water by NMR spectroscopy, and in silico by quantum chemical shift calculation. All these methods provide strong evidence for the same hexameric complex geometry. Its structural features bear striking similarity to nucleic acid architectures and their peptidic counterparts, especially alanyl-PNA. The whole self-assembly process is highly solvent- and temperature-dependent and occurs with a high degree of cooperativity--no intermediates are observed. Formation and dissociation of the nanorosette, however, are kinetically slow. The limitation to a hexameric aggregate can be explained by six sterically demanding valine residues, whose replacement by alanines may result in formation of infinite fibers.

Selective complexation of N-alkylpyridinium salts: binding of NAD+ in water.

A new class of receptor molecules is presented that is highly selective for N-alkylpyridinium ions and electron-poor aromatics. Its key feature is the combination of a well-preorganized molecular clip with an electron-rich inner cavity and strategically placed, flanking bis-phosphonate monoester anions. This shape and arrangement of binding sites attracts predominantly flat electron-poor aromatics in water, binds them mainly by pi-cation, pi-pi, CH-pi, and hydrophobic interactions, and leads to their highly efficient desolvation. NAD(+) and NADP, the important cofactors of many redox enzymes, are recognized by the new receptor molecule, which embraces the catalytically active nicotinamide site and the adenine unit. Even nucleosides such as adenosine are likewise drawn into the clip's cavity. Complex formation and structures were examined by one- and two-dimensional NMR spectroscopy, Job plot analyses, and isothermal titration microcalorimetric (ITC) measurements, as well as quantum chemical calculations of (1)H NMR shifts. The new receptor molecule is a promising tool for controlling enzymatic oxidation processes and for DNA chemistry.