Structural and Vibrational Properties of Amino Acids from Composite Schemes and Double-Hybrid DFT: Hydrogen Bonding in Serine as a Test Case

The structures, relative stabilities, and vibrational wavenumbers of the two most stable conformers of serine, stabilized by the O-H···N, O-H···O═C and N-H···O-H intramolecular hydrogen bonds, have been evaluated by means of state-of-the-art composite schemes based on coupled-cluster (CC) theory. The so-called "cheap" composite approach (CCSD(T)/(CBS+CV)^MP2) allowed determination of accurate equilibrium structures and harmonic vibrational wavenumbers, also pointing out significant corrections beyond the CCSD(T)/cc-pVTZ level. These accurate results stand as a reference for benchmarking selected hybrid and double-hybrid, dispersion-corrected DFT functionals. B2PLYP-D3 and DSDPBEP86 in conjunction with a triple-ζ basis set have been confirmed as effective methodologies for structural and spectroscopic studies of medium-sized flexible biomolecules, also showing intramolecular hydrogen bonding. These best performing double-hybrid functionals have been employed to simulate IR spectra by means of vibrational perturbation theory, also considering hybrid CC/DFT schemes. The best overall agreement with experiment, with mean absolute error of 8 cm^-1, has been obtained by combining CCSD(T)/(CBS+CV)^MP2 harmonic wavenumbers with B2PLYP-D3/maug-cc-pVTZ anharmonic corrections. Finally, a composite scheme entirely based on CCSD(T) calculations (CCSD(T)/CBS+CV) has been employed for energetics, further confirming that serine II is the most stable conformer, also when zero-point vibrational energy corrections are included.

TOAC spin-labeled peptides tailored for DNP-NMR studies in lipid membrane environments

The benefit of combining in-cell solid-state dynamic nuclear polarization (DNP) NMR and cryogenic temperatures is providing sufficient signal/noise and preservation of bacterial integrity via cryoprotection to enable in situ biophysical studies of antimicrobial peptides. The radical source required for DNP was delivered into cells by adding a nitroxide-tagged peptide based on the antimicrobial peptide maculatin 1.1 (Mac1). In this study, the structure, localization, and signal enhancement properties of a single (T-MacW) and double (T-T-MacW) TOAC (2,2,6,6-tetramethylpiperidine-N-oxyl-4-amino-4-carboxylic acid) spin-labeled Mac1 analogs were determined within micelles or lipid vesicles. The solution NMR and circular dichroism results showed that the spin-labeled peptides adopted helical structures in contact with micelles. The peptides behaved as an isolated radical source in the presence of multilamellar vesicles, and the electron paramagnetic resonance (EPR) electron-electron distance for the doubly spin-labeled peptide was ∼1 nm. The strongest paramagnetic relaxation enhancement (PRE) was observed for the lipid NMR signals near the glycerol-carbonyl backbone and was stronger for the doubly spin-labeled peptide. Molecular dynamics simulation of the T-T-MacW radical source in phospholipid bilayers supported the EPR and PRE observations while providing further structural insights. Overall, the T-T-MacW peptide achieved better ¹³C and ¹⁵N signal NMR enhancements and ¹H spin-lattice T₁ relaxation than T-MacW.

Side-Chain Polarity Modulates the Intrinsic Conformational Landscape of Model Dipeptides

The intrinsic conformational preferences of small peptides may provide additional insight into the thermodynamics and kinetics of protein folding. In this study, we explore the underlying energy landscapes of two model peptides, namely, Ac-Ala-NH₂ and Ac-Ser-NH₂, using geometry-optimization-based tools developed within the context of energy landscape theory. We analyze not only how side-chain polarity influences the structural preferences of the dipeptides, but also other emergent properties of the landscape, including heat capacity profiles, and kinetics of conformational rearrangements. The contrasting topographies of the free energy landscape agree with recent results from Fourier transform microwave spectroscopy experiments, where Ac-Ala-NH₂ was found to exist as a mixture of two conformers, while Ac-Ser-NH₂ remained structurally locked, despite exhibiting an apparently rich conformational landscape.

Backbone-Constrained Peptides: Temperature and Secondary Structure Affect Solid-State Electron Transport

The primary sequence and secondary structure of a peptide are crucial to charge migration, not only in solution (electron transfer, ET), but also in the solid-state (electron transport, ETp). Hence, understanding the charge migration mechanisms is fundamental to the development of biomolecular devices and sensors. We report studies on four Aib-containing helical peptide analogues: two acyclic linear peptides with one and two electron-rich alkene-based side chains, respectively, and two peptides that are further rigidified into a macrocycle by a side bridge constraint, containing one or no alkene. ETp was investigated across Au/peptide/Au junctions, between 80 and 340 K in combination with the molecular dynamic (MD) simulations. The results reveal that the helical structure of the peptide and electron-rich side chain both facilitate the ETp. As temperature increases, the loss of helical structure, change of monolayer tilt angle, and increase of thermally activated fluctuations affect the conductance of peptides. Specifically, room temperature conductance across the peptide monolayers correlates well with previously observed ET rate constants, where an interplay between backbone rigidity and electron-rich side chains was revealed. Our findings provide new means to manipulate electronic transport across solid-state peptide junctions.

A spin-1 representation for dual-funnel energy landscapes

The interconversion between the left- and right-handed helical folds of a polypeptide defines a dual-funneled free energy landscape. In this context, the funnel minima are connected through a continuum of unfolded conformations, evocative of the classical helix-coil transition. Physical intuition and recent conjectures suggest that this landscape can be mapped by assigning a left- or right-handed helical state to each residue. We explore this possibility using all-atom replica exchange molecular dynamics and an Ising-like model, demonstrating that the energy landscape architecture is at odds with a two-state picture. A three-state model-left, right, and unstructured-can account for most key intermediates during chiral interconversion. Competing folds and excited conformational states still impose limitations on the scope of this approach. However, the improvement is stark: Moving from a two-state to a three-state model decreases the fit error from 1.6 kBT to 0.3 kBT along the left-to-right interconversion pathway.

Chain length, temperature and solvent effects on the structural properties of α-aminoisobutyric acid homooligopeptides

Non-coded α-amino acids, originally exploited by nature, have been successfully reproduced by recent synthetic strategies to confer special structural and functional properties to small peptides. The most known and well-studied atypical residue is α-aminoisobutyric acid (Aib), which is contained in a fairly large number of peptides with known antibiotic effects. Here, we report on a molecular dynamics (MD) study of a series of homooligopeptides based on α-aminoisobutyric acid (Aib) with increasing length (Ac-(Aib)n-NMe, n = 5, 6, 7 and 10) and at various temperatures, employing a recent extension of the AMBER force field tailored for the Aib residue. Solvent effects have been analyzed by comparative MD simulations of a heptapeptide in water and dimethylsulfoxide at different temperatures. Our results show that the preference for the 310- and/or α-helix structures, which typically characterize Aib based peptides, is finely tuned by several factors including the chain length, temperature and solvent nature. While the transitions between intra-molecular i → i + 3 and i → i + 4 hydrogen bonds characterizing 310 and α-helices, respectively, are rather fast in small peptides (in the picosecond timescale), our analysis shows that the above physical and chemical factors modulate the relative equilibrium populations of the two helical structures. The obtained results nicely agree with available experimental data and support the use of the new force field for modeling Aib containing peptides.

The influence of solvent on conformational properties of peptides with Aib residue-a DFT study

The conformational propensities of the Aib residue on the example of two model peptides Ac-Aib-NHMe (1) and Ac-Aib-NMe2 (2), were studied by B3LYP and M06-2X functionals, in the gas phase and in the polar solvents. To verify the reliability of selected functionals, we also performed MP2 calculations for the tested molecules in vacuum. Polarizable continuum models (PCM and SMD) were used to estimate the solvent effect. Ramachandran maps were calculated to find all energy minima. Noncovalent intramolecular interactions due to hydrogen-bonds and dipole attractions between carbonyl groups are responsible for the relative stabilities of the conformers. In order to verify the theoretical results, the available conformations of similar X-ray structures from the Cambridge Crystallographic Data Center (CCDC) were analyzed. The results of the calculations show that both derivatives with the Aib residue in the gas phase prefer structures stabilized by intramolecular N-H⋯O hydrogen bonds, i.e., C5 and C7 conformations, while polar solvent promotes helical conformation with φ, ψ values equal to +/-60°, +/-40°. In addition, in the case of molecule 2, the helical conformation is the only one available in the polar environment. This result is fully consistent with the X-ray data. Graphical abstract Effect of solvent on the Ramachandran maps of the model peptides with Aib residue.