Electronic and vibrational polarizabilities of the twenty naturally occurring amino acids

Ultrafast terahertz Stark spectroscopy reveals the excited-state dipole moments of retinal in bacteriorhodopsin

The photoinduced all-trans to 13-cis isomerization of the retinal Schiff base represents the ultrafast first step in the reaction cycle of bacteriorhodopsin (BR). Extensive experimental and theoretical work has addressed excited-state dynamics and isomerization via a conical intersection with the ground state. In conflicting molecular pictures, the excited state potential energy surface has been modeled as a pure S[Formula: see text] state that intersects with the ground state, or in a 3-state picture involving the S[Formula: see text] and S[Formula: see text] states. Here, the photoexcited system passes two crossing regions to return to the ground state. The electric dipole moment of the Schiff base in the S[Formula: see text] and S[Formula: see text] state differs strongly and, thus, its measurement allows for assessing the character of the excited-state potential. We apply the method of ultrafast terahertz (THz) Stark spectroscopy to measure electric dipole changes of wild-type BR and a BR D85T mutant upon electronic excitation. A fully reversible transient broadening and spectral shift of electronic absorption is induced by a picosecond THz field of several megavolts/cm and mapped by a 120-fs optical probe pulse. For both BR variants, we derive a moderate electric dipole change of 5 [Formula: see text] 1 Debye, which is markedly smaller than predicted for a neat S[Formula: see text]-character of the excited state. In contrast, S[Formula: see text]-admixture and temporal averaging of excited-state dynamics over the probe pulse duration gives a dipole change in line with experiment. Our results support a picture of electronic and nuclear dynamics governed by the interaction of S[Formula: see text] and S[Formula: see text] states in a 3-state model.

Hybrid quantum-classical polarizability model for single molecule biosensing

Optical whispering gallery mode biosensors are able to detect single molecules through effects of their polarizability. We address the factors that affect the polarizability of amino acids, which are the building blocks of life, via electronic structure theory. Amino acids are detected in aqueous environments, where their polarizability is different compared to the gasphase due to solvent effects. Solvent effects include structural changes, protonation and the local field enhancement through the solvent (water). We analyse the impact of these effects and find that all contribute to an increased effective polarizability in the solvent. We also address the excess polarizability relative to the displaced water cavity and develop a hybrid quantum-classical model that is in good agreement with self-consistent calculations. We apply our model to calculate the excess polarizability of 20 proteinogenic amino acids and determine the minimum resolution required to distinguish the different molecules and their ionised conformers based on their polarizability.

Molecular properties of linear amino acids in water

Four linear amino acids of increased separation of the carboxyl and amino groups, namely glycine (aminoacetic acid), β-alanine (3-aminopropanoic acid), GABA (4-aminobutanoic acid) and DAVA (5-aminopentanoic acid), have been studied by quantum chemical ab initio and DFT methods including the solvent effect in order to get electronic structure and molecular descriptors, such as ionisation energy, electron affinity, molecular electronegativity, chemical hardness, electrophilicity index, dipole moment, quadrupole moment and dipole polarizability. Thermodynamic functions (zero-point energy, inner energy, enthalpy, entropy, and the Gibbs energy) were evaluated after the complete vibrational analysis at the true energy minimum provided by the full geometry optimization. Reaction Gibbs energy allows evaluating the absolute redox potentials on reduction and/or oxidation. The non-local non-additive molecular descriptors were compared along the series showing which of them behave as extensive, varying in match with the molar mass and/or separation of the carboxyl and amino groups. Amino acidic forms and zwitterionic forms of the substances were studied in parallel in order to compare their relative stability and redox properties. In total, 24 species were investigated by B3LYP/def2-TZVPD method (M1) including neutral molecules, molecular cations and molecular anions. For comparison, MP2/def2-TZVPD method (M2) with full geometry optimization and vibrational analysis in water has been applied for 12 species; analogously, for 24 substances, DLPNO-CCSD(T)/aug-cc-pVTZ method (M3) has been applied in the geometry obtained by MP2 and/or B3LYP. It was found that the absolute oxidation potential correlates with the adiabatic ionisation energy; the absolute reduction potential correlates with the adiabatic electron affinity and the electrophilicity index. In order to validate the used methodology with experimental vertical ionisation energies and vibrational spectrum obtained in gas phase, calculations were done also in vacuo.

First-Principles Simulation of Dielectric Function in Biomolecules

The dielectric spectra of complex biomolecules reflect the molecular heterogeneity of the proteins and are particularly important for the calculations of electrostatic (Coulomb) and electrodynamic (van der Waals) interactions in protein physics. The dielectric response of the proteins can be decomposed into different components depending on the size, structure, composition, locality, and environment of the protein in general. We present a new robust simulation method anchored in rigorous ab initio quantum mechanical calculations of explicit atomistic models, without any indeterminate parameters to compute and gain insight into the dielectric spectra of small proteins under different conditions. We implement this methodology to a polypeptide RGD-4C (1FUV) in different environments, and the SD1 domain in the spike protein of SARS-COV-2. Two peaks at 5.2-5.7 eV and 14.4-15.2 eV in the dielectric absorption spectra are observed for 1FUV and SD1 in vacuum as well as in their solvated and salted models.

Light-induced effects in glycine aqueous solution studied by Fourier transform infrared-emission spectroscopy and ultraviolet-visible spectroscopy

Although amino acids are insensitive to visible light, as is generally accepted, we show that particular light-matter interaction can break this obviousness. Using sensitive (FT-IR)-technique in a combination of a broadband visible light source, we registered emission spectra of glycine in the range 2500-500 cm^-1. Sensitivity of the infrared emission spectrum to the exciting power -induced changes in the glycine structure was demonstrated experimentally.Vibrational spectra of glycine displayed the prominent spectral features of CH₂, COO-, COOH, NH⁺₃ groups in the "fingerprint region". Simultaneous appearance of ionised COO- and unionised COOH forms of glycine in solution at neutral pH clearly indicated that visible light induces the partial protonation of COO- groups; if so, visible light irradiation should lead to occurrence of dimers or dimeric hydrogen - bonded structures. Spectroscopic and microscopic evidence of visible light-mediated formation of aggregates and nucleus in aqueous solution was presented.Electronic absorption/emission spectra of glycine in aqueous solution were primarily characterized in the near ultraviolet-visible region (240-600 nm). Negligible absorption near 270 nm was observed for a 1.0 M solution and dramatically enhanced with its "aging". Moreover, an extension of the absorption edge into the region above 400 nm could be seen. Due to the visible light irradiation, we observed modification of electronic structure or occurrence of additional species causing changes in absorption of glycine amino acid. For "aged" solution, it was shown that excitation spectra corresponding to the different emission wavelengths were entirely different, at that each excitation-spectral band had a characteristic emission band.Communicated by Ramaswamy H. Sarma.

Enhanced Dispersion and Polarization Interactions Achieved through Dithiophosphate Group Incorporation Yield a Dramatic Binding Affinity Increase for an RNA Aptamer-Thrombin Complex

Regiospecific replacement of a single phosphate (PO2) by a dithiophosphate (PS2) group in an RNA can dramatically increase its binding affinity for a target protein. Thus, complexes between antithrombin and anti-VEGF RNA aptamers with single dithiophosphate moieties and thrombin and VEGF, respectively, display equilibrium dissociation constants K_D of ca. 1 pM, 1000-fold tighter than the native RNA complexes (ca. 1 nM). Inspection of crystal structures of the native and PS2-RNA aptamer:thrombin complexes reveals an RNA-induced fit in the latter. This leads to a close approach between PS2 and the phenyl ring edge of Phe-232 that is surrounded by pairs of lysines and arginines. To better understand the origins of the tighter binding and individual contributions to the interaction energy, we carried out QM calculations with phosphate- and dithiophosphate-benzene and dimethyl phosphate- and dimethyl dithiophosphate-benzene model systems. These calculations demonstrate that the dithiophosphate-benzene interaction is much stronger than the corresponding interaction with phosphate. QM/MM calculations with the full complexes confirmed this finding and support the hypothesis that the electric field generated by basic residues surrounding Phe-232 is key to the polarization of the PS2 moiety. Thus, disparate polarization and dispersion energies between the PO2 and PS2 complexes contribute critically to the difference in binding affinity. By comparison, easier desolvation of the dithiophosphate group compared to phosphate does not contribute decisively to the observed difference in binding affinity. Favorable polarization and dispersion energies may be a general feature of the dramatic affinity gains seen for complexes between RNAs carrying dithiophosphate groups and their binding proteins.

UV Resonance Raman Spectroscopy as a Tool to Probe Membrane Protein Structure and Dynamics

Ultraviolet resonance Raman (UVRR) spectroscopy is a vibrational technique that reveals structures and dynamics of biological macromolecules without the use of extrinsic labels. By tuning the Raman excitation wavelength to the deep UV region (e.g., 228 nm), Raman signal from tryptophan and tyrosine residues are selectively enhanced, allowing for the study of these functionally relevant amino acids in lipid and aqueous environments. In this chapter, we present methods on the UVRR data acquisition and analysis of the tryptophan vibrational modes of a model β-barrel membrane protein, OmpA, in folded and unfolded conformations.

Insights into Protein Structure and Dynamics by Ultraviolet and Visible Resonance Raman Spectroscopy

Raman spectroscopy is a form of vibrational spectroscopy based on inelastic scattering of light. In resonance Raman spectroscopy, the wavelength of the incident light falls within an absorption band of a chromophore, and this overlap of excitation and absorption energy greatly enhances the Raman scattering efficiency of the absorbing species. The ability to probe vibrational spectra of select chromophores within a complex mixture of molecules makes resonance Raman spectroscopy an excellent tool for studies of biomolecules. In this Current Topic, we discuss the type of molecular insights obtained from steady-state and time-resolved resonance Raman studies of a prototypical photoactive protein, rhodopsin. We also review recent efforts in ultraviolet resonance Raman investigations of soluble and membrane-associated biomolecules, including integral membrane proteins and antimicrobial peptides. These examples illustrate that resonance Raman is a sensitive, selective, and practical method for studying the structures of biological molecules, and the molecular bonding, geometry, and environments of protein cofactors, the backbone, and side chains.

Anharmonic IR and Raman spectra and electronic and vibrational (hyper)polarizabilities of barbituric, 2-thiobarbituric and 2-selenobarbituric acids

Infrared, Raman and electronic absorption spectra, electronic and vibrational (hyper)polarizabilities, of barbituric, 2-thiobarbituric and 2-selenobarbituric acids were studied in gas using ab initio and density functional theory levels. The vibrational spectra were computed using harmonic and anharmonic methods. Anharmonic contributions improve the agreement between calculated and available experimental wavenumbers, especially in the highest-energy spectral region (wavenumbers >1700 cm(-1)). Vibrational and electronic transitions potentially useful to identify the investigated compounds were explored. The electronic and vibrational hyperpolarizabilities for the IDRI nonlinear optical (NLO) process at the λ value of 790 nm were computed. Supported by spectroscopic results, electronic and vibrational polarizabilities and second-order hyperpolarizabilities increase progressively in the order barbituric acid<2-thiobarbituric acid<2-selenobarbituric acid. The seleno-derivative is predicted to be ca. three/four times more hyperpolarizable than the barbituric acid. The Se→O or Se→S substitutions can be practical strategies to enhances the NLO properties of barbituric and thiobarbituric acid-based materials.

Electron correlation effects and density analysis of the first-order hyperpolarizability of neutral guanine tautomers

Dipole moments (μ), charge distributions, and static electronic first-order hyperpolarizabilities (β(μ)) of the two lowest-energy keto tautomers of guanine (7H and 9H) were determined in the gas phase using Hartree-Fock, Møller-Plesset perturbation theory (MP2 and MP4), and DFT (PBE1PBE, B97-1, B3LYP, CAM-B3LYP) methods with Dunning's correlation-consistent aug-cc-pVDZ and d-aug-cc-pVDZ basis sets. The most stable isomer 7H exhibits a μ value smaller than that of the 9H form by a factor of ca. 3.5. The β μ value of the 9H tautomer is strongly dependent on the computational method employed, as it dramatically influences the β(μ) (9H)/β(μ) (7H) ratio, which at the highest correlated MP4/aug-cc-pVDZ level is predicted to be ca. 5. The Coulomb-attenuating hybrid exchange-correlation CAM-B3LYP method is superior to the conventional PBE1PBE, B3LYP, and B97-1 functionals in predicting the β(μ) values. Differences between the largest diagonal hyperpolarizability components were clarified through hyperpolarizability density analyses. Dipole moment and first-order hyperpolarizability are molecular properties that are potentially useful for distinguishing the 7H from the 9H tautomer.

A THEORETICAL STUDY ON SECOND HARMONIC GENERATION HYPERPOLARIZABILITIES OF PHENYLALANINE POLYPEPTIDES

The second harmonic generation (SHG) hyperpolarizabilities of phenylalanine and homopolypeptides are investigated by configuration interaction among singly excited configurations (CIS) technique combined with the sum-over-states (SOS) method. The geometries of peptides containing phenylalanine ( Phe )n(n = 1–8) are optimized by B3LYP/6-31g(d) method, and they form the special structures like β-sheet (a common protein secondary structure). It is found that the energy gaps of various peptides are reduced and the hyperpolarizabilities are increased with the peptide chains lengthened. We discuss the origin of the second-order nonlinear optical response in phenylalanine homopolypeptides and confirm that the π → π* transitions in the aromatic residue of phenylalanine make the most important contributions to the second-order polarizability. Our results strongly suggest that the hyperpolarizabilities are dominated from the propagation direction of peptide chains.

The molecular refractive function of lens γ-Crystallins

γ-Crystallins constitute the major protein component in the nucleus of the vertebrate eye lens. Present at very high concentrations, they exhibit extreme solubility and thermodynamic stability to prevent scattering of light and formation of cataracts. However, functions beyond this structural role have remained mostly unclear. Here, we calculate molecular refractive index increments of crystallins. We show that all lens γ-crystallins have evolved a significantly elevated molecular refractive index increment, which is far above those of most proteins, including nonlens members of the βγ-crystallin family from different species. The same trait has evolved in parallel in crystallins of different phyla, including S-crystallins of cephalopods. A high refractive index increment can lower the crystallin concentration required to achieve a suitable refractive power of the lens and thereby reduce their propensity to aggregate and form cataracts. To produce a significant increase in the refractive index increment, a substantial global shift in amino acid composition is required, which can naturally explain the highly unusual amino acid composition of γ-crystallins and their functional homologues. This function provides a new perspective for interpreting their molecular structure.

On the distribution of protein refractive index increments

The protein refractive index increment, dn/dc, is an important parameter underlying the concentration determination and the biophysical characterization of proteins and protein complexes in many techniques. In this study, we examine the widely used assumption that most proteins have dn/dc values in a very narrow range, and reappraise the prediction of dn/dc of unmodified proteins based on their amino acid composition. Applying this approach in large scale to the entire set of known and predicted human proteins, we obtain, for the first time, to our knowledge, an estimate of the full distribution of protein dn/dc values. The distribution is close to Gaussian with a mean of 0.190 ml/g (for unmodified proteins at 589 nm) and a standard deviation of 0.003 ml/g. However, small proteins <10 kDa exhibit a larger spread, and almost 3000 proteins have values deviating by more than two standard deviations from the mean. Due to the widespread availability of protein sequences and the potential for outliers, the compositional prediction should be convenient and provide greater accuracy than an average consensus value for all proteins. We discuss how this approach should be particularly valuable for certain protein classes where a high dn/dc is coincidental to structural features, or may be functionally relevant such as in proteins of the eye.

Tryptophan-lipid interactions in membrane protein folding probed by ultraviolet resonance Raman and fluorescence spectroscopy

Aromatic amino acids of membrane proteins are enriched at the lipid-water interface. The role of tryptophan on the folding and stability of an integral membrane protein is investigated with ultraviolet resonance Raman and fluorescence spectroscopy. We investigate a model system, the β-barrel outer membrane protein A (OmpA), and focus on interfacial tryptophan residues oriented toward the lipid bilayer (trp-7, trp-170, or trp-15) or the interior of the β-barrel pore (trp-102). OmpA mutants with a single tryptophan residue at a nonnative position 170 (Trp-170) or a native position 7 (Trp-7) exhibit the greatest stability, with Gibbs free energies of unfolding in the absence of denaturant of 9.4 and 6.7 kcal/mol, respectively. These mutants are more stable than the tryptophan-free OmpA mutant, which exhibits a free energy of unfolding of 2.6 kcal/mol. Ultraviolet resonance Raman spectra of Trp-170 and Trp-7 reveal evolution of a hydrogen bond in a nonpolar environment during the folding reaction, evidenced by systematic shifts in hydrophobicity and hydrogen bond markers. These observations suggest that the hydrogen bond acceptor is the lipid acyl carbonyl group, and this interaction contributes significantly to membrane protein stabilization. Other spectral changes are observed for a tryptophan residue at position 15, and these modifications are attributed to development of a tryptophan-lipid cation-π interaction that is more stabilizing than an intraprotein hydrogen bond by ∼2 kcal/mol. As expected, there is no evidence for lipid-protein interactions for the tryptophan residue oriented toward the interior of the β-barrel pore. These results highlight the significance of lipid-protein interactions, and indicate that the bilayer provides more than a hydrophobic environment for membrane protein folding. Instead, a paradigm of lipid-assisted membrane protein folding and stabilization must be adopted.

Formation of peptide radical cations (m+·) in electron capture dissociation of peptides adducted with group IIB metal ions

Peptides adducted with different divalent Group IIB metal ions (Zn(2+), Cd(2+), and Hg(2+)) were found to give very different ECD mass spectra. ECD of Zn(2+) adducted peptides gave series of c-/z-type fragment ions with and without metal ions. ECD of Cd(2+) and Hg(2+) adducted model peptides gave mostly a-type fragment ions with M(+•) and fragment ions corresponding to losses of neutral side chain from M(+•). No detectable a-ions could be observed in ECD spectra of Zn(2+) adducted peptides. We rationalized the present findings by invoking both proton-electron recombination and metal-ion reduction processes. As previously postulated, divalent metal-ions adducted peptides could adopt several forms, including (a) [M + Cat](2+), (b) [(M + Cat - H) + H](2+), and (c) [(M + Cat - 2H) + 2H](2+). The relative population of these precursor ions depends largely on the acidity of the metal-ion peptide complexes. Peptides adducted with divalent metal-ions of small ionic radii (i.e., Zn(2+)) would form predominantly species (b) and (c); whereas peptides adducted with metal ions of larger ionic radii (i.e., Hg(2+)) would adopt predominantly species (a). Species (b) and (c) are believed to be essential for proton-electron recombination process to give c-/z-type fragments via the labile ketylamino radical intermediates. Species (c) is particularly important for the formation of non-metalated c-/z-type fragments. Without any mobile protons, species (a) are believed to undergo metal ion reduction and subsequently induce spontaneous electron transfer from the peptide moiety to the charge-reduced metal ions. Depending on the exothermicity of the electron transfer reaction, the peptide radical cations might be formed with substantial internal energy and might undergo further dissociation to give structural related fragment ions.

Calculations of the second virial coefficients of protein solutions with an extended fast multipole method

The osmotic second virial coefficients B(2) are directly related to the solubility of protein molecules in electrolyte solutions and can be useful to narrow down the search parameter space of protein crystallization conditions. Using a residue level model of protein-protein interaction in electrolyte solutions B(2) of bovine pancreatic trypsin inhibitor and lysozyme in various solution conditions such as salt concentration, pH and temperature are calculated using an extended fast multipole method in combination with the boundary element formulation. Overall, the calculated B(2) are well correlated with the experimental observations for various solution conditions. In combination with our previous work on the binding affinity calculations it is reasonable to expect that our residue level model can be used as a reliable model to describe protein-protein interaction in solutions.

Hydrogen bonding and solvent polarity markers in the uv resonance raman spectrum of tryptophan: application to membrane proteins

Ultraviolet resonance Raman (UVRR) spectra of tryptophan compounds in various solvents and a model peptide are presented and reveal systematic changes that reflect solvent polarity, hydrogen bond strength, and cation-pi interaction. The commonly utilized UVRR spectral marker for environment polarity that has been based on off-resonance Raman data, the tryptophan Fermi doublet ratio I1360/I1340, exhibits different values in on- and off-resonance Raman spectra as well as for different tryptophan derivatives. Specifically, the UVRR Fermi doublet ratio for indole ranges from 0.3 in polar solvents to 0.8 in nonpolar solvents, whereas the respective values reported here and previously for off-resonance Raman spectra are 0.5-1.3. UVRR Fermi doublet ratios for the more biologically relevant molecule, N-acetyl tryptophan ethyl ester (NATEE), are in a smaller range of 1.1 (polar solvent) to 1.7 (nonpolar solvent) and correlate to the solvent polarity/polarization parameters pi* and ETN. As has been reported previously, several UVRR modes are also sensitive to the hydrogen bond strength of the indole N-H moiety. Here, we report a new unambiguous marker for H-bonding: the ratio of the W10 (approximately 1237 cm-1) intensity to that of the W9 (approximately 1254 cm-1) mode (RW10). This ratio is 0.7 for NATEE in the absence of hydrogen bond acceptors and increases to 3.1 in the presence of strong hydrogen bond acceptors, with a value of 2.3 in water. The W8 and W17 modes shift more than +10 and approximately -5 cm-1 upon increase in hydrogen bond strength; this range for W17 is smaller than that reported previously and reflects a more realistic range for proteins and peptides in solution. Finally, our data provide evidence for change in the W18 and W16 relative intensity in the presence of cation-pi interactions. These UVRR markers are utilized to interpret spectra of model membrane-bound systems tryptophan octyl ester and the peptide toxin melittin. These spectra reveal the importance of intra- and intermolecular hydrogen bonding and cation-pi interactions that likely influence the partitioning of membrane-associated biomolecules to lipid bilayers or self-associated soluble oligomers. The UVRR analysis presented here modifies and augments prior reports and provides an unambiguous set of spectral makers that can be applied to elucidate the molecular microenvironment and structure of a wide range of complex systems, including anchoring tryptophan residues in membrane proteins and peptides.

Toward understanding amino acid adsorption at metallic interfaces: a density functional theory study

In examining adsorption of a few selected single amino acids on Au and Pd cluster models by density functional theory calculations, we have shown that specific side-chain binding affinity to the surface may occur because of a combination of effects, including charge transfer. Larger binding was calculated at the Pd interface. In addition, the interplay between amino acid solvation and adsorption at the interface was considered from first principles. This analysis serves as the first step toward gaining a more accurate understanding of specific interactions at the interface of biological-metal nanostructures than has been attempted in the past.