Anharmonic Vibrational Analysis of Biomolecules and Solvated Molecules Using Hybrid QM/MM Computations

GENESIS 2.1: High-Performance Molecular Dynamics Software for Enhanced Sampling and Free-Energy Calculations for Atomistic, Coarse-Grained, and Quantum Mechanics/Molecular Mechanics Models

GENeralized-Ensemble SImulation System (GENESIS) is a molecular dynamics (MD) software developed to simulate the conformational dynamics of a single biomolecule, as well as molecular interactions in large biomolecular assemblies and between multiple biomolecules in cellular environments. To achieve the latter purpose, the earlier versions of GENESIS emphasized high performance in atomistic MD simulations on massively parallel supercomputers, with or without graphics processing units (GPUs). Here, we implemented multiscale MD simulations that include atomistic, coarse-grained, and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. They demonstrate high performance and are integrated with enhanced conformational sampling algorithms and free-energy calculations without using external programs except for the QM programs. In this article, we review new functions, molecular models, and other essential features in GENESIS version 2.1 and discuss ongoing developments for future releases.

Unraveling the Phosphorus Adsorption Mechanisms in Three-dimensional Reduced Graphene Oxide Materials

To prevent eutrophication, controlling the phosphate concentration levels is one of the most important issues in surface water management. One of the most utilized methods is phosphate adsorption. However, its application faces a bottleneck due to the unclear understanding of adsorption and interaction mechanisms. The present work unlocks the phosphorus adsorption mechanisms in three-dimensional reduced graphene oxide with different reduction levels and pore sizes to remove phosphate from water using experiments and multiscale simulations. Experiments were performed to evaluate the influence of pH, ionic strength, and temperature on the adsorption. Molecular Dynamics and Ab Initio simulations evaluated the influence of the pore size and oxidation degrees of the materials. We show that the adsorption capacity of the materials increases with increasing pH and ionic strength and decreasing temperature. It is observed that the more oxidized the material and the less compact the structure, the better the adsorption. These results are theoretically explained in terms of the interaction of functional groups and the clustering of phosphate ions, which results in better adsorption in materials with larger pores. The underlying mechanisms for the 3D-reduced graphene oxide performance were confirmed by spectroscopy analysis. All the results show that 3D-reduced graphene oxide can sorb phosphate in different complex water remediation systems with characteristics that can be modulated by changing the material synthesis method.

Anisotropic dynamics of an interfacial enzyme active site observed using tethered substrate analogs and ultrafast 2D IR spectroscopy

Enzymes accelerate the rates of biomolecular reactions by many orders of magnitude compared to bulk solution, and it is widely understood that this catalytic effect arises from a combination of polar pre-organization and electrostatic transition state stabilization. A number of recent reports have also implicated ultrafast (femtosecond-picosecond) timescale motions in enzymatic activity. However, complications arising from spatially-distributed disorder, the occurrence of multiple substrate binding modes, and the influence of hydration dynamics on solvent-exposed active sites still confound many experimental studies. Here we use ultrafast two-dimensional infrared (2D IR) spectroscopy and covalently-tethered substrate analogs to examine dynamical properties of the promiscuous Pyrococcus horikoshii ene-reductase (PhENR) active site in two binding configurations mimicking proposed "inactive" and "reactive" Michaelis complexes. Spectral diffusion measurements of aryl-nitrile substrate analogs reveal an end-to-end tradeoff between fast (sub-ps) and slow (>5 ps) motions. Fermi resonant aryl-azide analogs that sense interactions of coupled oscillators are described. Lineshape and quantum beat analyses of these probes reveal characteristics that correlate with aryl-nitrile frequency fluctuation correlation functions parameters, demonstrating that this anisotropy is an intrinsic property of the water-exposed active site, where countervailing gradients of fast dynamics and disorder in the reactant ground state are maintained near the hydration interface. Our results suggest several plausible factors leading to state-selective rate enhancement and promiscuity in PhENR. This study also highlights a strategy to detect perturbations to vibrational modes outside the transparent window of the mid-IR spectrum, which may be extended to other macromolecular systems.

Computational analysis of the vibrational spectra and structure of aqueous cytosine

The recently developed efficient protocol for the explicit quantum mechanical modeling of the structure and IR spectra of liquids and solutions [Katsyuba et al., J. Phys. Chem. B, 2020, 124, 6664-6670] is used to describe aqueous solutions of cytosine. The same cluster model of a solute surrounded by the first solvation shell of solvent molecules was shown to be sufficient to reproduce experimental vibrational frequencies and relative IR and Raman intensities. An equally good quality of Raman spectra was provided by B3LYP-D3/def2-TZVP and B3LYP-D3/aug-cc-pVDZ simulations. Computations using the PBE functional were sufficient for modeling of the IR spectra but failed in the simulations of Raman scattering. It is shown that strong changes of frequencies and relative intensities of Raman and IR bands of cytosine, caused by its hydration, cannot be completely assigned to the influence of hydrogen bonds (HBs) with 7 or 8 closest water molecules. They are rather ascribed to the combined effect of solute-solute and solute-solvent HBs with the participation of at least 30 water molecules separating cytosine from the bulk solvent. This suggests that the vibrational modes and derivatives of the polarizability and dipole moment of the solute are mainly locally influenced by its first hydration shell, while the influence of bulk water is rather modest.

Gaussian process regression adaptive density-guided approach: Toward calculations of potential energy surfaces for larger molecules

We present a new program implementation of the Gaussian process regression adaptive density-guided approach [Schmitz et al., J. Chem. Phys. 153, 064105 (2020)] for automatic and cost-efficient potential energy surface construction in the MidasCpp program. A number of technical and methodological improvements made allowed us to extend this approach toward calculations of larger molecular systems than those previously accessible and maintain the very high accuracy of constructed potential energy surfaces. On the methodological side, improvements were made by using a Δ-learning approach, predicting the difference against a fully harmonic potential, and employing a computationally more efficient hyperparameter optimization procedure. We demonstrate the performance of this method on a test set of molecules of growing size and show that up to 80% of single point calculations could be avoided, introducing a root mean square deviation in fundamental excitations of about 3 cm-1. A much higher accuracy with errors below 1 cm-1 could be achieved with tighter convergence thresholds still reducing the number of single point computations by up to 68%. We further support our findings with a detailed analysis of wall times measured while employing different electronic structure methods. Our results demonstrate that GPR-ADGA is an effective tool, which could be applied for cost-efficient calculations of potential energy surfaces suitable for highly accurate vibrational spectra simulations.

Macrophage Cell Membrane Coating on Piperine-Loaded MIL-100(Fe) Nanoparticles for Breast Cancer Treatment

Piperine (PIP), a compound found in Piper longum, has shown promise as a potential chemotherapeutic agent for breast cancer. However, its inherent toxicity has limited its application. To overcome this challenge, researchers have developed PIP@MIL-100(Fe), an organic metal-organic framework (MOF) that encapsulates PIP for breast cancer treatment. Nanotechnology offers further treatment options, including the modification of nanostructures with macrophage membranes (MM) to enhance the evasion of the immune system. In this study, the researchers aimed to evaluate the potential of MM-coated MOFs encapsulated with PIP for breast cancer treatment. They successfully synthesized MM@PIP@MIL-100(Fe) through impregnation synthesis. The presence of MM coating on the MOF surface was confirmed through SDS-PAGE analysis, which revealed distinct protein bands. Transmission electron microscopy (TEM) images demonstrated the existence of a PIP@MIL-100(Fe) core with a diameter of around 50 nm, surrounded by an outer lipid bilayer layer measuring approximately 10 nm in thickness. Furthermore, the researchers evaluated the cytotoxicity indices of the nanoparticles against various breast cancer cell lines, including MCF-7, BT-549, SKBR-3, and MDA. The results demonstrated that the MOFs exhibited between 4 and 17 times higher cytotoxicity (IC₅₀) in all four cell lines compared to free PIP (IC₅₀ = 193.67 ± 0.30 µM). These findings suggest that MM@PIP@MIL-100(Fe) holds potential as an effective treatment for breast cancer. The study's outcomes highlight the potential of utilizing MM-coated MOFs encapsulated with PIP as an innovative approach for breast cancer therapy, offering improved cytotoxicity compared to free PIP alone. Further research and development are warranted to explore the clinical translation and optimize the efficacy and safety of this treatment strategy.

Allosteric regulation of β-reaction stage I in tryptophan synthase upon the α-ligand binding

Tryptophan synthase (TRPS) is a bifunctional enzyme consisting of α- and β-subunits that catalyzes the last two steps of L-tryptophan (L-Trp) biosynthesis. The first stage of the reaction at the β-subunit is called β-reaction stage I, which converts the β-ligand from an internal aldimine [E(Ain)] to an α-aminoacrylate [E(A-A)] intermediate. The activity is known to increase 3-10-fold upon the binding of 3-indole-D-glycerol-3'-phosphate (IGP) at the α-subunit. The effect of α-ligand binding on β-reaction stage I at the distal β-active site is not well understood despite the abundant structural information available for TRPS. Here, we investigate the β-reaction stage I by carrying out minimum-energy pathway searches based on a hybrid quantum mechanics/molecular mechanics (QM/MM) model. The free-energy differences along the pathway are also examined using QM/MM umbrella sampling simulations with QM calculations at the B3LYP-D3/aug-cc-pVDZ level of theory. Our simulations suggest that the sidechain orientation of βD305 near the β-ligand likely plays an essential role in the allosteric regulation: a hydrogen bond is formed between βD305 and the β-ligand in the absence of the α-ligand, prohibiting a smooth rotation of the hydroxyl group in the quinonoid intermediate, whereas the dihedral angle rotates smoothly after the hydrogen bond is switched from βD305-β-ligand to βD305-βR141. This switch could occur upon the IGP-binding at the α-subunit, as evidenced by the existing TRPS crystal structures.

How does multiple substrate binding lead to substrate inhibition of CYP2D6 metabolizing dextromethorphan? A theoretical study

CYP2D6 is one of the most important metalloenzymes involved in the biodegradation of many drug molecules in the human body. It has been found that multiple substrate binding can lead to substrate inhibition of CYP2D6 metabolizing dextromethorphan (DM), but the corresponding theoretical mechanism is rarely reported. Therefore, we chose DM as the probe and performed molecular dynamics simulations and quantum mechanical calculations on CYP2D6-DM systems to investigate the mechanism of how the multiple substrate binding leads to the substrate inhibition of CYP2D6 metabolizing substrates. According to our results, three gate residues (Arg221, Val374, and Phe483) for the catalytic pocket are determined. We also found that the multiple substrate binding can lead to substrate inhibition by reducing the stability of CYP2D6 binding DM and increasing the reactive activation energy of the rate-determining step. Our findings would help to understand the substrate inhibition of CYP2D6 metabolizing the DM and enrich the knowledge of the drug-drug interactions for the cytochrome P450 superfamily.

General Design Strategy to Precisely Control the Emission of Fluorophores via a Twisted Intramolecular Charge Transfer (TICT) Process

Fluorogenic probes for bioimaging have become essential tools for life science and medicine, and the key to their development is a precise understanding of the mechanisms available for fluorescence off/on control, such as photoinduced electron transfer (PeT) and Förster resonance energy transfer (FRET). Here we establish a new molecular design strategy to rationally develop activatable fluorescent probes, which exhibit a fluorescence off/on change in response to target biomolecules, by controlling the twisted intramolecular charge transfer (TICT) process. This approach was developed on the basis of a thorough investigation of the fluorescence quenching mechanism of N-phenyl rhodamine dyes (commercially available as the QSY series) by means of time-dependent density functional theory (TD-DFT) calculations and photophysical evaluation of their derivatives. To illustrate and validate this TICT-based design strategy, we employed it to develop practical fluorogenic probes for HaloTag and SNAP-tag. We further show that the TICT-controlled fluorescence off/on mechanism is generalizable by synthesizing a Si-rhodamine-based fluorogenic probe for HaloTag, thus providing a palette of chemical dyes that spans the visible and near-infrared range.

Multistructural microiteration combined with QM/MM-ONIOM electrostatic embedding

Multistructural microiteration (MSM) is a method to take account of contributions of multiple surrounding structures in a geometrical optimization or reaction path calculation using the quantum mechanics/molecular mechanics (QM/MM) ONIOM method. In this study, we combined MSM with the electrostatic embedding (EE) scheme of the QM/MM-ONIOM method by extending its original formulation for mechanical embedding (ME). MSM-EE takes account of the polarization in the QM region induced by point charges assigned to atoms in the multiple surrounding structures, where the point charges are scaled by the weight factor of each surrounding structure determined through MSM. The performance of MSM-EE was compared with that of the other methods, i.e., ONIOM-ME, ONIOM-EE, and MSM-ME, by applying them to three chemical processes: (1) chorismate-to-prephenate transformation in aqueous solution, (2) the same transformation as (1) in an enzyme, and (3) hydroxylation in p-hydroxybenzoate hydroxylase. These numerical tests of MSM-EE yielded barriers and reaction energies close to experimental values with computational costs comparable to those of the other three methods.

Retinal Vibrations in Bacteriorhodopsin are Mechanically Harmonic but Electrically Anharmonic: Evidence From Overtone and Combination Bands

Fundamental vibrations of the chromophore in the membrane protein bacteriorhodopsin (BR), a protonated Schiff base retinal, have been studied for decades, both by resonance Raman and by infrared (IR) difference spectroscopy. Such studies started comparing vibrational changes between the initial BR state (all-trans retinal) and the K intermediate (13-cis retinal), being later extended to the rest of intermediates. They contributed to our understanding of the proton-pumping mechanism of BR by exploiting the sensitivity of fundamental vibrational transitions of the retinal to its conformation. Here, we report on new bands in the 2,500 to 1,800 cm⁻¹ region of the K-BR difference FT-IR spectrum. We show that the bands between 2,500 and 2,300 cm⁻¹ originate from overtone and combination transitions from C-C stretches of the retinal. We assigned bands below 2,300 cm⁻¹ to the combination of retinal C-C stretches with methyl rocks and with hydrogen-out-of-plane vibrations. Remarkably, experimental C-C overtone bands appeared at roughly twice the wavenumber of their fundamentals, with anharmonic mechanical constants ≤3.5 cm⁻¹, and in some cases of ∼1 cm⁻¹. Comparison of combination and fundamental bands indicates that most of the mechanical coupling constants are also very small. Despite the mechanical quasi-harmonicity of the C-C stretches, the area of their overtone bands was only ∼50 to ∼100 times smaller than of their fundamental bands. We concluded that electrical anharmonicity, the second mechanism giving intensity to overtone bands, must be particularly high for the retinal C-C stretches. We corroborated the assignments of negative bands in the K-BR difference FT-IR spectrum by ab initio anharmonic vibrational calculations of all-trans retinal in BR using a quantum-mechanics/molecular mechanics approach, reproducing reasonably well the small experimental anharmonic and coupling mechanical constants. Yet, and in spite accounting for both mechanical and electrical anharmonicities, the intensity of overtone C-C transitions was underestimated by a factor of 4–20, indicating room for improvement in state-of-the-art anharmonic vibrational calculations. The relatively intense overtone and combination bands of the retinal might open the possibility to detect retinal conformational changes too subtle to significantly affect fundamental transitions but leaving a footprint in overtone and combination transitions.

Weight average approaches for predicting dynamical properties of biomolecules

Recent advances in atomistic molecular dynamics (MD) simulations of biomolecules allow us to explore their conformational spaces widely, observing large-scale conformational fluctuations or transitions between distinct structures. To reproduce or refine experimental data using MD simulations, structure ensembles, which are characterized by multiple structures and their statistical weights on the rugged free-energy landscapes, are often used. Here, we summarize weight average approaches for various experimental measurements. Weight average approaches are now applied to hybrid quantum mechanics/molecular mechanics MD simulations to predict fast vibrational motions in a protein with a high accuracy for better understanding of molecular functions from atomic structures.

Segregation on the nanoscale coupled to liquid water polyamorphism in supercooled aqueous ionic-liquid solution

The most intriguing hypothesis explaining many water anomalies is a metastable liquid-liquid phase transition (LLPT) at high pressure and low temperatures, experimentally hidden by homogeneous nucleation. Recent infrared spectroscopic experiments showed that upon addition of hydrazinium trifluoroacetate to water, the supercooled ionic solution undergoes a sharp, reversible LLPT at ambient pressure, possible offspring of that in pure water. Here, we calculate the temperature-dependent signature of the OH-stretching band, reporting on the low/high density phase of water, in neat water and in the same experimentally investigated ionic solution. The comparison between the infrared signature of the pure liquid and that of the ionic solution can be achieved only computationally, providing insight into the nature of the experimentally observed phase transition and allowing us to investigate the effects of ionic compounds on the high to low density supercooled liquid water transition. We show that the experimentally observed crossover behavior in the ionic solution can be reproduced only if the phase transition between the low- and high-density liquid states of water is coupled to a mixing-unmixing transition between the water component and the ions: at low temperatures, water and ions are separated and the water component is a low density liquid. At high temperatures, water and ions get mixed and the water component is a high-density liquid. The separation at low temperatures into ion-rich and ion-poor regions allows unveiling the polyamorphic nature of liquid water, leading to a crossover behavior resembling that observed in supercooled neat water under high pressure.

Advances, challenges and perspectives of quantum chemical approaches in molecular spectroscopy of the condensed phase

The purpose of this review is to demonstrate advances, challenges and perspectives of quantum chemical approaches in molecular spectroscopy of the condensed phase. Molecular spectroscopy, particularly vibrational spectroscopy and electronic spectroscopy, has been used extensively for a wide range of areas of chemical sciences and materials science as well as nano- and biosciences because it provides valuable information about structure, functions, and reactions of molecules. In the meantime, quantum chemical approaches play crucial roles in the spectral analysis. They also yield important knowledge about molecular and electronic structures as well as electronic transitions. The combination of spectroscopic approaches and quantum chemical calculations is a powerful tool for science, in general. Thus, our article, which treats various spectroscopy and quantum chemical approaches, should have strong implications in the wider scientific community. This review covers a wide area of molecular spectroscopy from far-ultraviolet (FUV, 120-200 nm) to far-infrared (FIR, 400-10 cm-1)/terahertz and Raman spectroscopy. As quantum chemical approaches, we introduce several anharmonic approaches such as vibrational self-consistent field (VSCF) and the combination of periodic harmonic calculations with anharmonic corrections based on finite models, grid-based techniques like the Numerov approach, the Cartesian coordinate tensor transfer (CCT) method, Symmetry-Adapted Cluster Configuration-Interaction (SAC-CI), and the ZINDO (Semi-empirical calculations at Zerner's Intermediate Neglect of Differential Overlap). One can use anharmonic approaches and grid-based approaches for both infrared (IR) and near-infrared (NIR) spectroscopy, while CCT methods are employed for Raman, Raman optical activity (ROA), FIR/terahertz and low-frequency Raman spectroscopy. Therefore, this review overviews cross relations between molecular spectroscopy and quantum chemical approaches, and provides various kinds of close-reality advanced spectral simulation for condensed phases.

Anharmonic Vibrational Calculations Based on Group-Localized Coordinates: Applications to Internal Water Molecules in Bacteriorhodopsin

An efficient anharmonic vibrational method is developed exploiting the locality of molecular vibration. Vibrational coordinates localized to a group of atoms are employed to divide the potential energy surface (PES) of a system into intra- and inter-group contributions. Then, the vibrational Schrödinger equation is solved based on a PES, in which the inter-group coupling is truncated at the harmonic level while accounting for the intra-group anharmonicity. The method is applied to a pentagonal hydrogen bond network (HBN) composed of internal water molecules and charged residues in a membrane protein, bacteriorhodopsin. The PES is calculated by the quantum mechanics/molecular mechanics (QM/MM) calculation at the level of B3LYP-D3/aug-cc-pVDZ. The infrared (IR) spectrum is computed using a set of coordinates localized to each water molecule and amino acid residue by second-order vibrational quasi-degenerate perturbation theory (VQDPT2). Benchmark calculations show that the proposed method yields the N-D/O-D stretching frequencies with an error of 7 cm-1 at the cost reduced by more than five times. In contrast, the harmonic approximation results in a severe error of 150 cm-1. Furthermore, the size of QM regions is carefully assessed to find that the QM regions should include not only the pentagonal HBN itself but also its HB partners. VQDPT2 calculations starting from transient structures obtained by molecular dynamics simulations have shown that the structural sampling has a significant impact on the calculated IR spectrum. The incorporation of anharmonicity, sufficiently large QM regions, and structural samplings are of essential importance to reproduce the experimental IR spectrum. The computational spectrum paves the way for decoding the IR signal of strong HBNs and helps elucidate their functional roles in biomolecules.

A computational insight into the relationship between side chain IR line shapes and local environment in fibril-like structures

Infrared spectroscopy is a widely used technique to characterize protein structures and protein mediated processes. While the amide I band provides information on proteins' secondary structure, amino acid side chains are used as infrared probes for the investigation of protein reactions and local properties. In this paper, we use a hybrid quantum mechanical/classical molecular dynamical approach based on the perturbed matrix method to compute the infrared band due to the C=O stretching mode of amide-containing side chains. We calculate, at first, the infrared band of zwitterionic glutamine in water and obtain results in very good agreement with the experimental data. Then, we compute the signal arising from glutamine side chains in a microcrystal of the yeast prion Sup35-derived peptide, GNNQQNY, with a fibrillar structure. The infrared bands obtained by selective isotopic labeling of the two glutamine residues, Q4 and Q5, of each peptide were experimentally used to investigate the local hydration in the fibrillar microcrystal. The experimental spectra of the two glutamine residues, which experience different hydration environments, feature different spectral signals that are well reproduced by the corresponding calculated spectra. In addition, the analysis of the simulated spectra clarifies the molecular origin of the experimentally observed spectroscopic differences that arise from the different local electric field experienced by the two glutamine residues, which is, in turn, determined by a different hydrogen bonding pattern.

Current and future research directions in computer-aided near-infrared spectroscopy: A perspective

The present review aims to draw a perspective on the vibrational spectroscopy combined with the tools of computational chemistry. This includes an overview of the accomplishments made so far, the assessment of the present development trends and the prospects for continuing these advances. State-of-the-art methods, current challenges and the expected future advances are evaluated from the point-of-view of the practical application in vibrational spectroscopy. A special attention is given to near-infrared (NIR) spectroscopy, which occupies a distinct position among the techniques of vibrational spectroscopy. As the result of intrinsically complex spectra, reliance on the anharmonicity as well as keen interest given to complex materials, NIR spectroscopy may particularly benefit from computational chemistry. The present key limitations hindering development of NIR spectroscopy are identified; these constitute primarily the limit in the treatable system size and the inability to effectively include chemical matrix effects. Given the expanding role of NIR spectroscopy in science and industry, lifting these limitations would directly enhance the general potential of this technique.

Exploring the Minimum-Energy Pathways and Free-Energy Profiles of Enzymatic Reactions with QM/MM Calculations

Understanding molecular mechanisms of enzymatic reactions is of vital importance in biochemistry and biophysics. Here, we introduce new functions of hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in the GENESIS program to compute the minimum-energy pathways (MEPs) and free-energy profiles of enzymatic reactions. For this purpose, an interface in GENESIS is developed to utilize a highly parallel electronic structure program, QSimulate-QM (https://qsimulate.com), calling it as a shared library from GENESIS. Second, algorithms to search the MEP are implemented, combining the string method (E et al. J. Chem. Phys. 2007, 126, 164103) with the energy minimization of the buffer MM region. The method implemented in GENESIS is applied to an enzyme, triosephosphate isomerase, which converts dihyroxyacetone phosphate to glyceraldehyde 3-phosphate in four proton-transfer processes. QM/MM-molecular dynamics simulations show performances of greater than 1 ns/day with the density functional tight binding (DFTB), and 10-30 ps/day with the hybrid density functional theory, B3LYP-D3. These performances allow us to compute not only MEP but also the potential of mean force (PMF) of the enzymatic reactions using the QM/MM calculations. The barrier height obtained as 13 kcal mol-1 with B3LYP-D3 in the QM/MM calculation is in agreement with the experimental results. The impact of conformational sampling in PMF calculations and the level of electronic structure calculations (DFTB vs B3LYP-D3) suggests reliable computational protocols for enzymatic reactions without high computational costs.

Adaptive density-guided approach to double incremental potential energy surface construction

We present a combination of the recently developed double incremental expansion of potential energy surfaces with the well-established adaptive density-guided approach to grid construction. This unique methodology is based on the use of an incremental expansion for potential energy surfaces, known as n-mode expansion; an incremental many-body representation of the electronic energy; and an efficient vibrational density-guided approach to automated determination of grid dimensions and granularity. The reliability of the method is validated calculating potential energy surfaces and obtaining fundamental excitation energies for three moderate-size chain-like molecular systems. The use of our methodology leads to considerable computational savings for potential energy surface construction compared to standard approaches while maintaining a high level of accuracy in the resulting potential energy surfaces. Additional investigations indicate that our method can be applied to covalently bound and strongly interacting molecular systems, even though these cases are known to be very unfavorable for fragmentation schemes. We therefore conclude that the presented methodology is a robust and flexible approach to potential energy surface construction, which introduces considerable computational savings without compromising the accuracy of vibrational spectra calculations.

Modeling amino-acid side chain infrared spectra: the case of carboxylic residues

Infrared (IR) spectroscopy is commonly utilized for the investigation of protein structures and protein-mediated processes.

Spectra-Structure Correlations in Isotopomers of Ethanol (CX3CX2OX; X = H, D): Combined Near-Infrared and Anharmonic Computational Study

The effect of isotopic substitution on near-infrared (NIR) spectra has not been studied in detail. With an exception of few major bands, it is difficult to follow the spectral changes due to complexity of NIR spectra. Recent progress in anharmonic quantum mechanical calculations allows for accurate reconstruction of NIR spectra. Taking this opportunity, we carried out a systematic study of NIR spectra of six isotopomers of ethanol (CX₃CX₂OX; X = H, D). Besides, we calculated the theoretical spectra of two other isotopomers (CH₃CD₂OD and CD₃CH₂OD) for which the experimental spectra are not available. The anharmonic calculations were based on generalized vibrational second-order perturbation theory (GVPT2) at DFT and MP2 levels with several basis sets. We compared the accuracy and efficiency of various computational methods. It appears that the best results were obtained with B2PLYP-GD3BJ/def2-TZVP//CPCM approach. Our simulations included the first and second overtones, as well as binary and ternary combinations bands. This way, we reliably reproduced even minor bands in the spectra of diluted samples (0.1 M in CCl₄). On this basis, the effect of isotopic substitution on NIR spectra of ethanol was accurately reproduced and comprehensively explained.