Amide isomerization pathways: Electronic and structural background of protonation- and deprotonation-mediated cis-trans interconversions.

The cis-trans isomerization of amide bonds leads to wide range of structural and functional changes in proteins and can easily be the rate-limiting step in folding. The trans isomer is thermodynamically more stable than the cis, nevertheless the cis form can play a role in biopolymers' function. The molecular system of N-methylacetamide · 2H2O is complex enough to reveal energetics of the cis-trans isomerization at coupled cluster single-double and coupled cluster single-double and perturbative triple [CCSD(T)] levels of theory. The cis-trans isomerization cannot be oversimplified by a rotation along ω, since this rotation is coupled with the N-atom pyramidal inversion, requesting the introduction of a second dihedral angle "α." Full f(ω,α) potential energy surfaces of the different amide protonation states, critical points and isomerization reaction paths were determined, and the barriers of the neutral, O-protonated and N-deprotonated amides were found too high to allow cis-trans interconversion at room temperature: ∼85, ∼140, and ∼110 kJ mol-1, respectively. For the N-protonated amide bond, the cis form (ω = 0°) is a maximum rather than a minimum, and each ω state is accessible for less than ∼10 kJ mol-1. Here we outline a cis-trans isomerization pathway with a previously undescribed low energy transition state, which suggests that the proton is transferred from the more favorable O- to the N-protonation site with the aid of nearby water molecules, allowing the trans → cis transition to occur at an energy cost of ≤11.6 kJ mol-1. Our results help to explain why isomerase enzymes operate via protonated amide bonds and how N-protonation of the peptide bond occurs via O-protonation.

Off-pathway 3D-structure provides protection against spontaneous Asn/Asp isomerization: shielding proteins Achilles heel.

Spontaneous deamidation prompted backbone isomerization of Asn/Asp residues resulting in - most cases - the insertion of an extra methylene group into the backbone poses a threat to the structural integrity of proteins. Here we present a systematical analysis of how temperature, pH, presence of charged residues, but most importantly backbone conformation and dynamics affect isomerization rates as determined by nuclear magnetic resonance in the case of designed peptide-models. We demonstrate that restricted mobility (such as being part of a secondary structural element) may safeguard against isomerization, but this protective factor is most effective in the case of off-pathway folds which can slow the reaction by several magnitudes compared to their on-pathway counterparts. We show that the geometric descriptors of the initial nucleophilic attack of the isomerization can be used to classify local conformation and contribute to the design of stable protein drugs, antibodies or the assessment of the severity of mutations.At any ­Asn/AspGly­ sites in proteins a spontaneous backbone isomerization occurs within days under physiological conditions leading to various forms of proteopathy. This unwanted transformation especially harmful to long-lived proteins (e.g. hemoglobin and crystallins), can be slowed down, though never stopped, by a rigid three-dimensional protein fold, if it can delay in the conformational maze, on-pathway intermediates from occurring.

Four faces of the interaction between ions and aromatic rings.

Non-covalent interactions between ions and aromatic rings play an important role in the stabilization of macromolecular complexes; of particular interest are peptides and proteins containing aromatic side chains (Phe, Trp, and Tyr) interacting with negatively (Asp and Glu) and positively (Arg and Lys) charged amino acid residues. The structures of the ion-aromatic-ring complexes are the result of an interaction between the large quadrupole moment of the ring and the charge of the ion. Four attractive interaction types are proposed to be distinguished based on the position of the ion with respect to the plane of the ring: perpendicular cation-π (CP⊥ ), co-planar cation-π (CP∥ ), perpendicular anion-π (AP⊥ ), and co-planar anion-π (AP∥ ). To understand more than the basic features of these four interaction types, a systematic, high-level quantum chemical study is performed, using the X- + C6 H6 , M+ + C6 H6 , X- + C6 F6 , and M+ + C6 F6 model systems with X- = H- , F- , Cl- , HCOO- , CH3 COO- and M+ = H+ , Li+ , Na+ , NH4+, CH3 NH3+, whereby C6 H6 and C6 F6 represent an electron-rich and an electron-deficient π system, respectively. Benchmark-quality interaction energies with small uncertainties, obtained via the so-called focal-point analysis (FPA) technique, are reported for the four interaction types. The computations reveal that the interactions lead to significant stabilization, and that the interaction energy order, given in kcal mol-1 in parentheses, is CP⊥ (23-37) > AP⊥ (14-21) > CP∥ (9-22) > AP∥ (6-16). A natural bond orbital analysis performed leads to a deeper qualitative understanding of the four interaction types. To facilitate the future quantum chemical characterization of ion-aromatic-ring interactions in large biomolecules, the performance of three density functional theory methods, B3LYP, BHandHLYP, and M06-2X, is tested against the FPA benchmarks, with the result that the M06-2X functional performs best. © 2017 Wiley Periodicals, Inc.

Predictable Conformational Diversity in Foldamers of Sugar Amino Acids.

A systematic conformational search was carried out for monomers and homohexamers of furanoid ß-amino acids: cis-(S,R) and trans-(S,S) stereoisomers of aminocyclopentane carboxylic acid (ACPC), two different aminofuranuronic acids (AFUα and AFUß), their isopropylidene derivatives (AFU(ip)), and the key intermediate ß-aminotetrahydrofurancarboxylic acid (ATFC). The stereochemistry of the building blocks was chosen to match that of the natural sugar amino acid (xylose and ribose) precursors (XylAFU and RibAFU). The results show that hexamers of cis-furanoid ß-amino acids show great variability: while hydrophobic cyclopentane (cis-ACPC)6 and hydrophilic (XylAFUα/ß)6 foldamers favor two different zigzagged conformation as hexamers, the backbone fold turns into a helix in the case of (cis-ATFC)6 (10-helix) and (XylAFU(ip))6 (14-helix). Trans stereochemistry resulted in hexamers exclusively with the right-handed helix conformation, (H12P)6, regardless of their polarity. We found that the preferred oligomeric structure of XylAFUα/ß is conformationally compatible with ß-pleated sheets, while that of the trans/(S,S) units matches with α-helices of proteins.

Aromatic Cluster Sensor of Protein Folding: Near-UV Electronic Circular Dichroism Bands Assigned to Fold Compactness.

Both far- and near-UV electronic circular dichroism (ECD) spectra have bands sensitive to thermal unfolding of Trp and Tyr residues containing proteins. Beside spectral changes at 222 nm reporting secondary structural variations (far-UV range), Lb bands (near-UV range) are applicable as 3D-fold sensors of protein's core structure. In this study we show that both Lb (Tyr) and Lb (Trp) ECD bands could be used as sensors of fold compactness. ECD is a relative method and thus requires NMR referencing and cross-validation, also provided here. The ensemble of 204 ECD spectra of Trp-cage miniproteins is analysed as a training set for "calibrating" Trp↔Tyr folded systems of known NMR structure. While in the far-UV ECD spectra changes are linear as a function of the temperature, near-UV ECD data indicate a non-linear and thus, cooperative unfolding mechanism of these proteins. Ensemble of ECD spectra deconvoluted gives both conformational weights and insight to a protein folding↔unfolding mechanism. We found that the Lb293 band is reporting on the 3D-structure compactness. In addition, the pure near-UV ECD spectrum of the unfolded state is described here for the first time. Thus, ECD folding information now validated can be applied with confidence in a large thermal window (5≤T≤85 °C) compared to NMR for studying the unfolding of Trp↔Tyr residue pairs. In conclusion, folding propensities of important proteins (RNA polymerase II, ubiquitin protein ligase, tryptase-inhibitor etc.) can now be analysed with higher confidence.

The role of entropy in initializing the aggregation of peptides: a first principle study on oligopeptide oligomerization.

The initiation and progression of Alzheimer's disease is coupled to the oligo- and polymerization of amyloid peptides in the brain. Amyloid like aggregates of protein domains were found practically independent of their primary sequences. Thus, the driving force of the transformation from the original to a disordered amyloid fold is expected to lie in the protein backbone common to all proteins. In order to investigate the thermodynamics of oligomerization, full geometry optimizations and frequency calculations were performed both on parallel and antiparallel ß-pleated sheet model structures of [HCO-(Ala)(1-6)-NH(2)](2) and (For-Ala(1-2)-NH(2))(1-6) peptides, both at the B3LYP and M05-2X/6-311++G(d,p)//M05-2X/6-31G(d) levels of theory, both in vacuum and in water. Our results show that relative entropy and enthalpy both show a hyperbolic decrease with increasing residue number and with increasing number of strands as well. Thus, di- and oligomerization are always thermodynamically favored. Antiparallel arrangements were found to have greater stability than parallel arrangements of the polypeptide backbones. During our study the relative changes in thermodynamic functions are found to be constant for long enough peptides, indicating that stability and entropy terms are predictable. All thermodynamic functions of antiparallel di- and oligomers show a staggered nature along the increasing residue number. By identifying and analyzing the 6 newly emerging dimer vibrational modes of the 10- and 14-membered building units, the staggered nature of the entropy function can be rationalized. Thus, the vanishing rotational and translational modes with respect to single strands are converted into entropy terms "holding tight" the dimers and oligomers formed, rationalizing the intrinsic adherence of natural polypeptide backbones to aggregate.

Cryo-EM structure of acylpeptide hydrolase reveals substrate selection by multimerization and a multi-state serine-protease triad.

The first structure of tetrameric mammalian acylaminoacyl peptidase, an enzyme that functions as an upstream regulator of the proteasome through the removal of terminal N-acetylated residues from its protein substrates, was determined by cryo-EM and further elucidated by MD simulations. Self-association results in a toroid-shaped quaternary structure, guided by an amyloidogenic ß-edge and unique inserts. With a Pro introduced into its central ß-sheet, sufficient conformational freedom is awarded to the segment containing the catalytic Ser587 that the serine protease catalytic triad alternates between active and latent states. Active site flexibility suggests that the dual function of catalysis and substrate selection are fulfilled by a novel mechanism: substrate entrance is regulated by flexible loops creating a double-gated channel system, while binding of the substrate to the active site is required for stabilization of the catalytic apparatus - as a second filter before hydrolysis. The structure not only underlines that within the family of S9 proteases homo-multimerization acts as a crucial tool for substrate selection, but it will also allow drug design targeting of the ubiquitin-proteasome system.

A carbapenem antibiotic inhibiting a mammalian serine protease: structure of the acylaminoacyl peptidase-meropenem complex.

The structure of porcine AAP (pAAP) in a covalently bound complex with meropenem was determined by cryo-EM to 2.1 Å resolution, showing the mammalian serine-protease inhibited by a carbapenem antibiotic. AAP is a modulator of the ubiquitin-proteasome degradation system and the site of a drug-drug interaction between the widely used antipsychotic, valproate and carbapenems. The active form of pAAP - a toroidal tetramer - binds four meropenem molecules covalently linked to the catalytic Ser587 of the serine-protease triad, in an acyl-enzyme state. AAP is hindered from fully processing the antibiotic by the displacement and protonation of His707 of the catalytic triad. We show that AAP is made susceptible to the association by its unusually sheltered active pockets and flexible catalytic triads, while the carbapenems possess sufficiently small substituents on their ß-lactam rings to fit into the shallow substrate-specificity pocket of the enzyme.

The inherent flexibility of peptides and protein fragments quantitized by CD in conjunction with CCA+.

ECD spectroscopy is traditionally used for rapid, non-atomic level structure analysis of natural products such as peptides and proteins. Unlike globular proteins, peptides less frequently adopt a single 3D-fold in a time average manner. Moreover, they exhibit an ensemble of conformers composed of a multitude of substantially different structures. In principle, both ECD- and vibrational circular dichroism (VCD)-spectroscopy are sensitive enough to pick up structural information on these dynamic ensembles. However, the interpretation of the raw spectral data of these highly dynamic molecular systems can be cumbersome. The herein presented Convex Constraint Analysis Plus method, or CCA+ for short (http://www.chem.elte.hu/departments/protnmr/cca/), provides a unique opportunity for spectral ensemble analysis of peptides, glycopeptides, peptidomimetics, and other foldamers. The precision and accuracy of the approach is presented here through different peptide model systems. An interesting temperature and pH dependent folding and unfolding of a miniprotein (e.g. Tc5b variant) is also described. Analysis of CD spectra sets strongly affected by solvent and ion type is also introduced to account for severe environmental-induced structure influencing effect(s). The deconvolution makes always possible the quantitative data analysis even when the interpretation of the deconvolution resulted in pure CD curves is complex.

Atropisomerism of the Asn α radicals revealed by Ramachandran surface topology.

C radicals are typically trigonal planar and thus achiral, regardless of whether they originate from a chiral or an achiral C-atom (e.g., C-H + (•)OH → C• + H2O). Oxidative stress could initiate radical formation in proteins when, for example, the H-atom is abstracted from the Cα-carbon of an amino acid residue. Electronic structure calculations show that such a radical remains achiral when formed from the achiral Gly, or the chiral but small Ala residues. However, when longer side-chain containing proteogenic amino acid residues are studied (e.g., Asn), they provide radicals of axis chirality, which in turn leads to atropisomerism observed for the first time for peptides. The two enantiomeric extended backbone structures, •ßL and •ßD, interconvert via a pair of enantiotopic reaction paths, monitored on a 4D Ramachandran surface, with two distinct transition states of very different Gibbs-free energies: 37.4 and 67.7 kJ/mol, respectively. This discovery requires the reassessment of our understanding on radical formation and their conformational and stereochemical behavior. Furthermore, the atropisomerism of proteogenic amino acid residues should affect our understanding on radicals in biological systems and, thus, reframes the role of the D-residues as markers of molecular aging.

Intrinsically stable secondary structure elements of proteins: a comprehensive study of folding units of proteins by computation and by analysis of data determined by X-ray crystallography.

Different protein architectures show strong similarities regardless of their amino acid composition: the backbone folds of the different secondary structural elements exhibit nearly identical geometries. To investigate the principles of folding and stability properties, oligopeptide models (that is, HCO-(NH-L-CHR-CO)(n)-NH(2)) have been studied. Previously, ab initio structure determinations have provided a small amount of information on the conformational building units of di- and tripeptides. A maximum of nine differently folded backbone types is available for any natural alpha-amino acid residue, with the exception of proline. All of these conformers have different relative energies. The present study compiles an ab inito database of optimized HCO-(L-Xxx)(n)-NH(2) structures, where 1

Peptide models. XXXIII. Extrapolation of low-level Hartree-Fock data of peptide conformation to large basis set SCF, MP2, DFT, and CCSD(T) results. The Ramachandran surface of alanine dipeptide computed at various levels of theory.

At the dawn of the new millenium, new concepts are required for a more profound understanding of protein structures. Together with NMR and X-ray-based 3D-structure determinations in silico methods are now widely accepted. Homology-based modeling studies, molecular dynamics methods, and quantum mechanical approaches are more commonly used. Despite the steady and exponential increase in computational power, high level ab initio methods will not be in common use for studying the structure and dynamics of large peptides and proteins in the near future. We are presenting here a novel approach, in which low- and medium-level ab initio energy results are scaled, thus extrapolating to a higher level of information. This scaling is of special significance, because we observed previously on molecular properties such as energy, chemical shielding data, etc., determined at a higher theoretical level, do correlate better with experimental data, than those originating from lower theoretical treatments. The Ramachandran surface of an alanine dipeptide now determined at six different levels of theory [RHF and B3LYP 3-21G, 6-31+G(d) and 6-311++G(d,p)] serves as a suitable test. Minima, first-order critical points and partially optimized structures, determined at different levels of theory (SCF, DFT), were completed with high level energy calculations such as MP2, MP4D, and CCSD(T). For the first time three different CCSD(T) sets of energies were determined for all stable B3LYP/6-311++G(d,p) minima of an alanine dipeptide. From the simplest ab initio data (e.g., RHF/3-21G) to more complex results [CCSD(T)/6-311+G(d,p)//B3LYP/6-311++G(d,p)] all data sets were compared, analyzed in a comprehensive manner, and evaluated by means of statistics.

Relative stability of major types of beta-turns as a function of amino acid composition: a study based on Ab initio energetic and natural abundance data.

Folding properties of small globular proteins are determined by their amino acid sequence (primary structure). This holds both for local (secondary structure) and for global conformational features of linear polypeptides and proteins composed from natural amino acid derivatives. It thus provides the rational basis of structure prediction algorithms. The shortest secondary structure element, the beta-turn, most typically adopts either a type I or a type II form, depending on the amino acid composition. Herein we investigate the sequence-dependent folding stability of both major types of beta-turns using simple dipeptide models (-Xxx-Yyy-). Gas-phase ab initio properties of 16 carefully selected and suitably protected dipeptide models (for example Val-Ser, Ala-Gly, Ser-Ser) were studied. For each backbone fold most probable side-chain conformers were considered. Fully optimized 321G RHF molecular structures were employed in medium level [B3LYP/6-311++G(d,p)//RHF/3-21G] energy calculations to estimate relative populations of the different backbone conformers. Our results show that the preference for beta-turn forms as calculated by quantum mechanics and observed in Xray determined proteins correlates significantly.