Tapta Kanchan Roy

Vibrational self-consistent field calculations for spectroscopy of biological molecules: new algorithmic developments and applications

This review describes the vibrational self-consistent field (VSCF) method and its other variants for computing anharmonic vibrational spectroscopy of biological molecules. The superiority and limitations of this algorithm are discussed with examples. The spectroscopic accuracy of the VSCF method is compared with experimental results and other available state-of-the-art algorithms for various biologically important systems. For large biological molecules with many vibrational modes, the scaling of computational effort is investigated. The accuracy of the vibrational spectra of biological molecules using the VSCF approach for different electronic structure methods is also assessed. Finally, a few open problems and challenges in this field are discussed.

Approximate first-principles anharmonic calculations of polyatomic spectra using MP2 and B3LYP potentials: comparisons with experiment.

Anharmonic vibrational spectroscopy calculations using MP2 and B3LYP computed potential surfaces are carried out for a series of molecules, and frequencies and intensities are compared with those from experiment. The vibrational self-consistent field with second-order perturbation correction (VSCF-PT2) is used in computing the spectra. The test calculations have been performed for the molecules HNO3, C2H4, C2H4O, H2SO4, CH3COOH, glycine, and alanine. Both MP2 and B3LYP give results in good accord with experimental frequencies, though, on the whole, MP2 gives very slightly better agreement. A statistical analysis of deviations in frequencies from experiment is carried out that gives interesting insights. The most probable percentage deviation from experimental frequencies is about -2% (to the red of the experiment) for B3LYP and +2% (to the blue of the experiment) for MP2. There is a higher probability for relatively large percentage deviations when B3LYP is used. The calculated intensities are also found to be in good accord with experiment, but the percentage deviations are much larger than those for frequencies. The results show that both MP2 and B3LYP potentials, used in VSCF-PT2 calculations, account well for anharmonic effects in the spectroscopy of molecules of the types considered.

Conformational Structures of a Decapeptide Validated by First Principles Calculations and Cold Ion Spectroscopy

Calculated structures of the two most stable conformers of a protonated decapeptide gramicidin S in the gas phase have been validated by comparing the vibrational spectra, calculated from first- principles and measured in a wide spectral range using infrared (IR)-UV double resonance cold ion spectroscopy. All the 522 vibrational modes of each conformer were calculated quantum mechanically and compared with the experiment without any recourse to an empirical scaling. The study demonstrates that first-principles calculations, when accounting for vibrational anharmonicity, can reproduce high-resolution experimental spectra well enough for validating structures of molecules as large as of 200 atoms. The validated accurate structures of the peptide may serve as templates for in silico drug design and absolute calibration of ion mobility measurements.

A Tandem In Situ Peptide Cyclization through Trifluoroacetic Acid Cleavage

We present a new approach for peptide cyclization during solid phase synthesis under highly acidic conditions. Our approach involves simultaneous in situ deprotection, cyclization and trifluoroacetic acid (TFA) cleavage of the peptide, which is achieved by forming an amide bond between a lysine side chain and a succinic acid linker at the peptide N-terminus. The reaction proceeds via a highly active succinimide intermediate, which was isolated and characterized. The structure of a model cyclic peptide was solved by NMR spectroscopy. Theoretical calculations support the proposed mechanism of cyclization. Our new methodology is applicable for the formation of macrocycles in solid-phase synthesis of peptides and organic molecules.

A Decapeptide Hydrated by Two Waters: Conformers Determined by Theory and Validated by Cold Ion Spectroscopy

The intrinsic structures of biomolecules in the gas phase may not reflect their native solution geometries. Microsolvation of the molecules bridges the two environments, enabling a tracking of molecular structural changes upon hydration at the atomistic level. We employ density functional calculations to compute a large pool of structures and vibrational spectra for a gas-phase complex, in which a doubly protonated decapeptide, gramicidin S, is solvated by two water molecules. Though most vibrations of this large complex are treated in a harmonic approximation, the water molecules and the vibrations of the host ion coupled to them are locally described by a quantum mechanical vibrational self-consistent field theory with second-order perturbation correction (VSCF-PT2). Guided and validated by the available cold ion spectroscopy data, the computational analysis identifies structures of the three experimentally observed conformers of the complex. They, mainly, differ by the hydration sites, of which the one at the Orn side chain is the most important for reshaping the peptide toward its native structure. The study demonstrates the ability of a quantum chemistry approach that intelligently combines the semiempirical and ab initio computations to disentangle a complex interplay of intra- and intermolecular hydrogen bonds in large molecular systems.

Mechanistic studies of malonic acid‐mediated in situ acylation

We have previously introduced an easy to perform, cost‐effective and highly efficient acetylation technique for solid phase synthesis (SPPS). Malonic acid is used as a precursor and the reaction proceeds via a reactive ketene that acetylates the target amine. Here we present a detailed mechanistic study of the malonic acid‐mediated acylation. The influence of reaction conditions, peptide sequence and reagents was systematically studied. Our results show that the methodology can be successfully applied to different types of peptides and nonpeptidic molecules irrespective of their structure, sequence, or conformation. Using alkyl, phenyl, and benzyl malonic acid, we synthesized various acyl peptides with almost quantitative yields. The ketenes obtained from the different malonic acid derived precursors were characterized by in situ 1H‐NMR. The reaction proceeded in short reaction times and resulted in excellent yields when using uronium‐based coupling agents, DIPEA as a base, DMF/DMSO/NMP as solvents, Rink amide/Wang/Merrifield resins, temperature of 20°C, pH 8–12 and 5 min preactivation at inert atmosphere. The reaction was unaffected by Lewis acids, transition metal ions, surfactants, or salt. DFT studies support the kinetically favorable concerted mechanism for CO2 and ketene formation that leads to the thermodynamically stable acylated products. We conclude that the malonic acid‐mediated acylation is a general method applicable to various target molecules.