Stabilizing interactions between aromatic and basic side chains in alpha-helical peptides and proteins. Tyrosine effects on helix circular dichroism

Evaluating Chemical Footprinting-Induced Perturbation of Protein Higher Order Structure

Specific amino acid footprinting mass spectrometry (MS) is an increasingly utilized method for elucidating protein higher order structure (HOS). It does this by adding to certain amino acid residues a mass tag, whose reaction extent depends on solvent accessibility and microenvironment of the protein. Unlike reactive free radicals and carbenes, these specific footprinters react slower than protein unfolding. Thus, their footprinting, under certain conditions, provokes structural changes to the protein, leading to labeling on non-native structures. It is critical to establish conditions (i.e., reagent concentrations, time of reaction) to ensure that the structure of the protein following footprinting remains native. Here, we compare the efficacy of five methods in assessing protein HOS following footprinting at the intact protein level and then further localize the perturbation at the peptide level. Three are MS-based methods that provide dose-response plot analysis, evaluation of Poisson distributions of precursor and products, and determination of the average number of modifications. These MS-based methods reliably and effectively indicate HOS perturbation at the intact protein level, whereas spectroscopic methods (circular dichroism (CD) and dynamic light scattering (DLS)) are less sensitive in monitoring subtle HOS perturbation caused by footprinting. Evaluation of HOS at the peptide level indicates regions that are sensitive to localized perturbations. Peptide-level analysis also provides higher resolution of the HOS perturbation, and we recommend using it for future footprinting studies. Overall, this work shows conclusive evidence for HOS perturbation caused by footprinting. Implementation of quality control workflows can identify conditions to avoid the perturbation, for footprinting, allowing accurate and reliable identification of protein structural changes that accompany, for example, ligand interactions, mutations, and changes in solution environment.

Enzyme-Driven, Switchable Catalysis Based on Dynamic Self-Assembly of Peptides

Covalent regulatory systems of enzymes are widely used to modulate biological enzyme activities. Inspired by the regulation of reactive-site phosphorylation in organisms, we developed peptide-based catecholase mimetics with switchable catalytic activity and high selectivity through the co-assembly of nanofibers comprising peptides and copper ions (Cu2+ ). Through careful design and modification of the peptide backbone structure based on the change in the free energy of the system, we identified the peptide with the most effective reversible catalytic activity. Kinase/phosphatase switches were used to control the reversible transition of nanofiber formation and depolymerization, as well as to modulate the active-site microenvironment. Notably, the self-assembly and disassembly processes of nanofibers were simulated using coarse-grained molecular dynamics. Furthermore, theoretical calculations confirmed the coordination of the peptide and Cu2+ , forming a zipper-like four-ligand structure at the catalytically active center of the nanofibers. Additionally, we conducted a comprehensive analysis of the catalytic mechanism. This study opens novel avenues for designing biomimetic enzymes with ordered structures and dynamic catalytic activities.

Divalent metal ion binding to Staphylococcus aureus FeoB transporter regions

Transition metal ions such as iron, copper, zinc, manganese or, nickel are essential in many biological processes. Bacteria have developed a number of mechanisms for their acquisition and transport, in which numerous of proteins and smaller molecules are involved. One of the representatives of these proteins is FeoB, which belongs to the Feo (ferrous ion transporter) family. Although ferrous iron transport system is widespread among microorganisms, it is still poorly described in Gram-positive pathogens, such as Staphylococcus aureus. In this work, combined potentiometric and spectroscopic studies (UV-Vis, CD and EPR) were carried out to determine Cu(II), Fe(II) and Zn(II) binding modes to FeoB fragments (Ac-IDYHKLMK-NH₂, Ac-ETSHDKY-NH₂, and Ac-SFLHMVGS-NH₂). For the first time iron(II) complexes with peptides were characterized by potentiometry. All studied ligands are able to form a variety of thermodynamically stable complexes with transition metal ions. It was concluded that among the studied systems, the most effective metal ion binding is observed for the Ac-ETSHDKY-NH₂ peptide. Moreover, comparing preferences of all ligands towards different metal ions, copper(II) complexes are the most stable ones at physiological pH.

Ion-Triggered Hydrogels Self-Assembled from Statistical Copolypeptides

Statistical copolypeptides comprising lysine and tyrosine with unprecedented ion-induced gelation behavior are reported. Copolypeptides are obtained by one-step N-carboxyanhydride (NCA) ring-opening polymerization. The gelation mechanism is studied by in situ SAXS analyses, in addition to optical spectroscopy and transmission electron microscopy (TEM). It is found that the gelation of these statistically polymerized polypeptides is due to the formation of stable intermolecular β-sheet secondary structures induced by the presence of salt ions as well as the aggregation of an α-helix between the copolypeptides. This behavior is unique to the statistical lysine/tyrosine copolypeptides and was not observed in any other amino acid combination or arrangement. Furthermore, the diffusion and mechanical properties of these hydrogels can be tuned through tailoring the polypeptide chain length and ion strength.

Physics-driven coarse-grained model for biomolecular phase separation with near-quantitative accuracy

Various physics- and data-driven sequence-dependent protein coarse-grained models have been developed to study biomolecular phase separation and elucidate the dominant physicochemical driving forces. Here, we present Mpipi, a multiscale coarse-grained model that describes almost quantitatively the change in protein critical temperatures as a function of amino-acid sequence. The model is parameterised from both atomistic simulations and bioinformatics data and accounts for the dominant role of π-π and hybrid cation-π/π-π interactions and the much stronger attractive contacts established by arginines than lysines. We provide a comprehensive set of benchmarks for Mpipi and seven other residue-level coarse-grained models against experimental radii of gyration and quantitative in-vitro phase diagrams; Mpipi predictions agree well with experiment on both fronts. Moreover, it can account for protein-RNA interactions, correctly predicts the multiphase behaviour of a charge-matched poly-arginine/poly-lysine/RNA system, and recapitulates experimental LLPS trends for sequence mutations on FUS, DDX4 and LAF-1 proteins.

Advances in Application of Azobenzene as a Trigger in Biomedicine: Molecular Design and Spontaneous Assembly

Azobenzene is a well-known derivative of stimulus-responsive molecular switches and has shown superior performance as a functional material in biomedical applications. The results of multiple studies have led to the development of light/hypoxia-responsive azobenzene for biomedical use. In recent years, long-wavelength-responsive azobenzene has been developed. Matching the longer wavelength absorption and hypoxia-response characteristics of the azobenzene switch unit to the bio-optical window results in a large and effective stimulus response. In addition, azobenzene has been used as a hypoxia-sensitive connector via biological cleavage under appropriate stimulus conditions. This has resulted in on/off state switching of properties such as pharmacology and fluorescence activity. Herein, recent advances in the design and fabrication of azobenzene as a trigger in biomedicine are summarized.

Approaches to Introduce Helical Structure in Cysteine-Containing Peptides with a Bimane Group

An i-i+4 or i-i+3 bimane-containing linker was introduced into a peptide known to target Estrogen Receptor alpha (ERα), in order to stabilise an α-helical geometry. These macrocycles were studied by CD and NMR to reveal the i-i+4 constrained peptide adopts a 3₁₀ -helical structure in solution, and an α-helical conformation on interaction with the ERα coactivator recruitment surface in silico. An acyclic bimane-modified peptide is also helical, when it includes a tryptophan or tyrosine residue; but is significantly less helical with a phenylalanine or alanine residue, which indicates such a bimane modification influences peptide structure in a sequence dependent manner. The fluorescence intensity of the bimane appears influenced by peptide conformation, where helical peptides displayed a fluorescence increase when TFE was added to phosphate buffer, compared to a decrease for less helical peptides. This study presents the bimane as a useful modification to influence peptide structure as an acyclic peptide modification, or as a side-chain constraint to give a macrocycle.

Reentrant liquid condensate phase of proteins is stabilized by hydrophobic and non-ionic interactions

Liquid-liquid phase separation of proteins underpins the formation of membraneless compartments in living cells. Elucidating the molecular driving forces underlying protein phase transitions is therefore a key objective for understanding biological function and malfunction. Here we show that cellular proteins, which form condensates at low salt concentrations, including FUS, TDP-43, Brd4, Sox2, and Annexin A11, can reenter a phase-separated regime at high salt concentrations. By bringing together experiments and simulations, we demonstrate that this reentrant phase transition in the high-salt regime is driven by hydrophobic and non-ionic interactions, and is mechanistically distinct from the low-salt regime, where condensates are additionally stabilized by electrostatic forces. Our work thus sheds light on the cooperation of hydrophobic and non-ionic interactions as general driving forces in the condensation process, with important implications for aberrant function, druggability, and material properties of biomolecular condensates.

Reentrant liquid condensate phase of proteins is stabilized by hydrophobic and non-ionic interactions

Liquid-liquid phase separation of proteins underpins the formation of membraneless compartments in living cells. Elucidating the molecular driving forces underlying protein phase transitions is therefore a key objective for understanding biological function and malfunction. Here we show that cellular proteins, which form condensates at low salt concentrations, including FUS, TDP-43, Brd4, Sox2, and Annexin A11, can reenter a phase-separated regime at high salt concentrations. By bringing together experiments and simulations, we demonstrate that this reentrant phase transition in the high-salt regime is driven by hydrophobic and non-ionic interactions, and is mechanistically distinct from the low-salt regime, where condensates are additionally stabilized by electrostatic forces. Our work thus sheds light on the cooperation of hydrophobic and non-ionic interactions as general driving forces in the condensation process, with important implications for aberrant function, druggability, and material properties of biomolecular condensates.

Structural analysis of cross α-helical nanotubes provides insight into the designability of filamentous peptide nanomaterials

The exquisite structure-function correlations observed in filamentous protein assemblies provide a paradigm for the design of synthetic peptide-based nanomaterials. However, the plasticity of quaternary structure in sequence-space and the lability of helical symmetry present significant challenges to the de novo design and structural analysis of such filaments. Here, we describe a rational approach to design self-assembling peptide nanotubes based on controlling lateral interactions between protofilaments having an unusual cross-α supramolecular architecture. Near-atomic resolution cryo-EM structural analysis of seven designed nanotubes provides insight into the designability of interfaces within these synthetic peptide assemblies and identifies a non-native structural interaction based on a pair of arginine residues. This arginine clasp motif can robustly mediate cohesive interactions between protofilaments within the cross-α nanotubes. The structure of the resultant assemblies can be controlled through the sequence and length of the peptide subunits, which generates synthetic peptide filaments of similar dimensions to flagella and pili.

Near-Ultraviolet Circular Dichroism and Two-Dimensional Spectroscopy of Polypeptides

A fully quantitative theory of the relationship between protein conformation and optical spectroscopy would facilitate deeper insights into biophysical and simulation studies of protein dynamics and folding. In contrast to intense bands in the far-ultraviolet, near-UV bands are much weaker and have been challenging to compute theoretically. We report some advances in the accuracy of calculations in the near-UV, which were realised through the consideration of the vibrational structure of the electronic transitions of aromatic side chains.

Disorder and partial folding in the regulatory subunit hinge region of Trypanosoma brucei protein kinase A: The C-linker portion inhibits the parasite's protein kinase A

Microbial pathogens, such as Trypanosoma brucei, have an enormous impact on global health and economic systems. Protein kinase A of T. brucei is an attractive drug target as it is an essential enzyme which differs significantly from its human homolog. The hinge region of this protein's regulatory domain is vital for enzymatic function, but its conformation is unknown. Here, the secondary structure of this region has been characterized using NMR and CD spectroscopies. More specifically, three overlapping peptides corresponding to residues T187-I211, G198-Y223 and V220-S245 called peptide 1, peptide 2 and peptide 3, respectively, were studied. The peptide 1 and peptide 2 are chiefly unfolded; only low populations (<10%) of α-helix were detected under the conditions studied. In contrast, the peptide 3 contains a long α-helix whose population is significantly higher; namely, 36% under the conditions studied. Utilizing the dihedral φ and ψ angles calculated on the basis of the NMR data, the conformation of the peptide 3 was calculated and revealed an α-helix spanning residues E230-N241. This α-helix showed amphiphilicity and reversible unfolding and refolding upon heating and cooling. Most fascinating, however, is its capacity to inhibit the activity of the catalytic domain of Trypanosoma equiperdum protein kinase A even though it is quite distinct from the canonical inhibitor motif. Based on this property, we advance that peptoids based on the peptide 3 α-helix could be novel leads for developing anti-trypanosomal therapeutics.

Aromatic Foldamer Helices as α-Helix Extended Surface Mimetics

Helically folded aromatic oligoamide foldamers have a size and geometrical parameters very distinct from those of α-helices and are not obvious candidates for α-helix mimicry. Nevertheless, they offer multiple sites for attaching side chains. It was found that some arrays of side chains at the surface of an aromatic helix make it possible to mimic extended α-helical surfaces. Synthetic methods were developed to produce quinoline monomers suitably functionalized for solid phase synthesis. A dodecamer was prepared. Its crystal structure validated the initial design and showed helix bundling involving the α-helix-like interface. These results open up new uses of aromatic helices to recognize protein surfaces and to program helix bundling in water.

Unfolding Dynamics of a Photoswitchable Helical Peptide

We present an atomistic force field for the azo-moiety of the photoswitchable FK-11-X peptide. We use the parameters to study the unfolding of the peptide through molecular dynamics simulations. The unfolded ensemble contains many different structures, ranging from a partially unfolded peptide to a fully unfolded structure. The averaged computed far-ultraviolet circular dichroism (CD) spectrum of the set of structures, which was simulated using the newly developed force field, agrees well with experiment. The rate of the simulated unfolding process was estimated to have a time constant of 5.80 ± 0.03 ns from the time evolution of the CD spectrum of the peptide, computed from the backbone conformations sampled over 40 simulated trajectories. Our estimated time constant is faster than, but not inconsistent with, previous experimental estimates from time-resolved infrared and optical rotatory dispersion spectroscopy.

Arginine/Tryptophan-Rich Cyclic α/β-Antimicrobial Peptides: The Roles of Hydrogen Bonding and Hydrophobic/Hydrophilic Solvent-Accessible Surface Areas upon Activity and Membrane Selectivity

The bacterial selectivity of an amphiphilic library of small cyclic α/β-tetra-, α/β-penta-, and α/β-hexapeptides rich in arginine/tryptophan (Arg/Trp) residues, which contains asymmetric backbone configurations and differ in hydrophobicity and alternating d,l-amino acids, was investigated against Bacillus subtilis and Escherichia coli. The structural analyses showed that the peptides tend to form assemblies of different shapes. All-l-peptides, especially the most hydrophobic pentamers, were more strongly anti-B. subtilis. With the exception to cyclo(Phe-d-Trp-β³ hArg-Arg-d-Trp) (Phe=phenylalanine), the peptides had no effects on inner membrane of E. coli, but lyzed the lipopolysaccharide layer according to their activity pattern. The activities adversely changed with a decrease in the number of amide intramolecular hydrogen bonds in assemblies of diastereomeric peptides and the ratio of hydrophobic/hydrophilic solvent-accessible surface areas. The remarkable enhanced entropic contribution for the partitioning of the least conformationally constrained cyclo(Trp-d-Phe-β³ hTrp-Arg-d-Arg) sequence into the membranes supported the strong self-assembly behavior, therefore making the peptide less penetrable through the E. coli outer layer.

Cation-π Interactions between Quaternary Ammonium Ions and Amino Acid Aromatic Groups in Aqueous Solution

Cation-π interactions play important roles in the stabilization of protein structures and protein-ligand complexes. They contribute to the binding of quaternary ammonium ligands (mainly RNH₃⁺ and RN(CH₃)₃⁺) to various protein receptors and are likely involved in the blockage of potassium channels by tetramethylammonium (TMA⁺) and tetraethylammonium (TEA⁺). Polarizable molecular models are calibrated for NH₄⁺, TMA⁺, and TEA⁺ interacting with benzene, toluene, 4-methylphenol, and 3-methylindole (representing aromatic amino acid side chains) based on the ab initio MP2(full)/6-311++G(d,p) properties of the complexes. Whereas the gas-phase affinity of the ions with a given aromatic follows the trend NH₄⁺ > TMA⁺ > TEA⁺, molecular dynamics simulations using the polarizable models show a reverse trend in water, likely due to a contribution from the hydrophobic effect. This reversed trend follows the solubility of aromatic hydrocarbons in quaternary ammonium salt solutions, which suggests a role for cation-π interactions in the salting-in of aromatic compounds in solution. Simulations in water show that the complexes possess binding free energies ranging from -1.3 to -3.3 kcal/mol (compared to gas-phase binding energies between -8.5 and -25.0 kcal/mol). Interestingly, whereas the most stable complexes involve TEA⁺ (the largest ion), the most stable solvent-separated complexes involve TMA⁺ (the intermediate-size ion).

Tyramine functionalization of poly(glycerol sebacate) increases the elasticity of the polymer

Poly(glycerol sebacate) (PGS) is an elastomer used widely in tissue engineering studies due to good biocompatibility. Hereby we report a tyramine functionalized PGS called PGS-TA. Tyramine adds a stronger physical bonding capability to PGS-TA. Tensile tests showed that the softness and toughness of the material were similar to PGS. However, PGS-TA demonstrated 16-folds increase of elastic deformations compared to PGS processed under identical conditions. The in vitro studies demonstrated that the viability, and metabolic activity of baboon smooth muscle cells were the same as those on tissue culture polystyrene. Porous subcutaneous implants of PGS-TA substantially degraded in vivo over two weeks, showing good biodegradability and biocompatibility. We expect PGS-TA to be useful for applications in tissues and organs that are subjected to large reversible mechanical deformations.

Molecular insight of DREAM and presenilin 1 C-terminal fragment interactions

Interactions between downstream regulatory element antagonist modulator (DREAM) and presenilin 1 (PS1) are related to numerous neuronal processes. We demonstrate that association of PS1 carboxyl peptide (residues 445-467, HL9) with DREAM is calcium dependent and stabilized by a cluster of three aromatic residues: F462 and F465 from PS1 and F252 from DREAM. Additional stabilization is provided by residues in a loop connecting α helices 7 and 8 in DREAM and residues of PS1, namely cation-π interactions between R200 in DREAM and F465 in PS1 and the salt bridges formed by R207 in DREAM and D450 and D458 in PS1.

Probing the Ca(2+)-assisted π-π interaction during Ca(2+)-dependent protein folding

Protein folding is governed by a balance of non-covalent interactions, of which cation-π and π-π play important roles. Theoretical calculations revealed a strong cooperativity between cation-π involving alkali and alkaline earth metal ions and π-π interactions, but however, no experimental evidence was provided in this regard. Here, we characterized a Ca(2+)-binding self-processing module (SPM), which mediates a highly-specific Ca(2+)-dependent autocatalytic processing of iron-regulated protein FrpC secreted by the pathogenic Gram-negative bacterium Neisseria meningitidis. The SPM undergoes a Ca(2+)-induced transition from an intrinsically unstructured conformation to the compact protein fold that is ultimately stabilized by the π-π interaction between two unique tryptophan residues arranged in the T-shaped orientation. Moreover, the pair of tryptophans is located in a close vicinity of a calcium-binding site, suggesting the involvement of a Ca(2+)-assisted π-π interaction in the stabilization of the tertiary structure of the SPM. This makes the SPM an excellent model for the investigation of the Ca(2+)-assisted π-π interaction during Ca(2+)-induced protein folding.

Protein thermostability engineering

Using structure and sequence based analysis we can engineer proteins to increase their thermal stability.

Nature of aryl-tyrosine interactions contribute to β-hairpin scaffold stability: NMR evidence for alternate ring geometry

The specific contribution of the acidic-aromatic β-sheet favouring amino acid tyrosine to the stability of short octapeptide β-hairpin structures is presented here. Solution NMR analysis in near-apolar environments suggests the energetically favourable mode of interaction to be T-shaped face-to-edge (FtE) and that a Trp-Tyr interacting pair is the most stabilizing. Alternate aryl geometries also exist in solution, which readily equilibrate between a preferred π···π conformation to an aromatic-amide conformation, without any change in the backbone structure. While the phenolic ring is readily accommodated at the "edge" of FtE aryl interactions, it exhibits an overall lowered contribution to scaffold stability in the "face" orientation. Such differential tyrosine interactions are key to its dual nature in proteins.

Implications of aromatic-aromatic interactions: From protein structures to peptide models

With increasing structural information on proteins, the opportunity to understand physical forces governing protein folding is also expanding. One of the significant non-covalent forces between the protein side chains is aromatic-aromatic interactions. Aromatic interactions have been widely exploited and thoroughly investigated in the context of folding, stability, molecular recognition, and self-assembly processes. Through this review, we discuss the contribution of aromatic interactions to the activity and stability of thermophilic, mesophilic, and psychrophilic proteins. Being hydrophobic, aromatic amino acids tend to reside in the protein hydrophobic interior or transmembrane segments of proteins. In such positions, it can play a diverse role in soluble and membrane proteins, and in α-helix and β-sheet stabilization. We also highlight here some excellent investigations made using peptide models and several approaches involving aryl-aryl interactions, as an increasingly popular strategy in protein and peptide engineering. A recent survey described the existence of aromatic clusters (trimer, tetramer, pentamer, and higher order assemblies), revealing the self-associating property of aryl groups, even in folded protein structures. The application of this self-assembly of aromatics in the generation of modern bionanomaterials is also discussed.

Conformational effects on the circular dichroism of Human Carbonic Anhydrase II: a multilevel computational study

Circular Dichroism (CD) spectroscopy is a powerful method for investigating conformational changes in proteins and therefore has numerous applications in structural and molecular biology. Here a computational investigation of the CD spectrum of the Human Carbonic Anhydrase II (HCAII), with main focus on the near-UV CD spectra of the wild-type enzyme and it seven tryptophan mutant forms, is presented and compared to experimental studies. Multilevel computational methods (Molecular Dynamics, Semiempirical Quantum Mechanics, Time-Dependent Density Functional Theory) were applied in order to gain insight into the mechanisms of interaction between the aromatic chromophores within the protein environment and understand how the conformational flexibility of the protein influences these mechanisms. The analysis suggests that combining CD semi empirical calculations, crystal structures and molecular dynamics (MD) could help in achieving a better agreement between the computed and experimental protein spectra and provide some unique insight into the dynamic nature of the mechanisms of chromophore interactions.

Molecular dynamics of β-hairpin models of epigenetic recognition motifs

The conformations and stabilities of the β-hairpin model peptides of Waters (Riemen, A. J.; Waters, M. L. Biochemistry 2009, 48, 1525; Hughes, R. M.; Benshoff, M. L.; Waters, M. L. Chemistry 2007, 13, 5753) have been experimentally characterized as a function of lysine ε-methylation. These models were developed to explore molecular recognition of known epigenetic recognition motifs. This system offered an opportunity to computationally examine the role of cation-π interactions, desolvation of the ε-methylated ammonium groups, and aromatic/aromatic interactions on the observed differences in NMR spectra. AMOEBA, a second-generation force field (Ponder, J. W.; Wu, C.; Ren, P.; Pande, V. S.; Chodera, J. D.; Schnieders, M. J.; Haque, I.; Mobley, D. L.; Lambrecht, D. S.; DiStasio, R. A., Jr.; Head-Gordon, M.; Clark, G. N.; Johnson, M. E.; Head-Gordon, T. J. Phys. Chem. B 2010, 114, 2549), was chosen as it includes both multipole electrostatics and polarizability thought to be essential to accurately characterize such interactions. Independent parametrization of ε-methylated amines was required from which aqueous solvation free energies were estimated and shown to agree with literature values. Molecular dynamics simulations (100 ns) using the derived parameters with model peptides, such as Ac-R-W-V-W-V-N-G-Orn-K(Me)(n)-I-L-Q-NH(2), where n = 0, 1, 2, or 3, were conducted in explicit solvent. Distances between the centers of the indole rings of the two-tryptophan residues, 2 and 4, and the ε-methylated ammonium group on Lys-9 as well as the distance between the N- and C-termini were monitored to estimate the strength and orientation of the cation-π and aromatic/aromatic interactions. In agreement with the experimental data, the stability of the β-hairpin increased significantly with lysine ε-methylation. The ability of MD simulations to reproduce the observed NOEs for the four peptides was further estimated for the monopole-based force fields, AMBER, CHARMM, and OPLSAA. AMOEBA correctly predicted over 80% of the observed NOEs for all 4 peptides, while the three-monopole force fields were 40-50% predictive in only 2 cases and approximately 10% in the other 10 examples. Preliminary analysis suggests that the decreased cost of desolvation of the substituted ammonium group significantly compensated for the reduced cation-π interaction resulting from the increased separation due to steric bulk of the ε-methylated amines.

Amino acid pair- and triplet-wise groupings in the interior of α-helical segments in proteins

A statistical approach has been applied to analyse primary structure patterns at inner positions of α-helices in proteins. A systematic survey was carried out in a recent sample of non-redundant proteins selected from the Protein Data Bank, which were used to analyse α-helix structures for amino acid pairing patterns. Only residues more than three positions apart from both termini of the α-helix were considered as inner. Amino acid pairings i, i+k (k=1, 2, 3, 4, 5), were analysed and the corresponding 20×20 matrices of relative global propensities were constructed. An analysis of (i, i+4, i+8) and (i, i+3, i+4) triplet patterns was also performed. These analysis yielded information on a series of amino acid patterns (pairings and triplets) showing either high or low preference for α-helical motifs and suggested a novel approach to protein alphabet reduction. In addition, it has been shown that the individual amino acid propensities are not enough to define the statistical distribution of these patterns. Global pair propensities also depend on the type of pattern, its composition and orientation in the protein sequence. The data presented should prove useful to obtain and refine useful predictive rules which can further the development and fine-tuning of protein structure prediction algorithms and tools.

Amino acid pairing at the N- and C-termini of helical segments in proteins

A systematic survey was carried out in an unbiased sample of 815 protein chains with a maximum of 20% homology selected from the Protein Data Bank, whose structures were solved at a resolution higher than 1.6 A and with a R-factor lower than 25%. A set of 5556 subsequences with alpha-helix or 3(10)-helix motifs was extracted from the protein chains considered. Global and local propensities were then calculated for all possible amino acid pairs of the type (i, i + 1), (i, i + 2), (i, i + 3), and (i, i + 4), starting at the relevant helical positions N1, N2, N3, C3, C2, C1, and N-int (interior positions), and also at the first nonhelical positions in both termini of the helices, namely, N-cap and C-cap. The statistical analysis of the propensity values has shown that pairing is significantly dependent on the type of the amino acids and on the position of the pair. A few sequences of three and four amino acids were selected and their high prevalence in helices is outlined in this work. The Glu-Lys-Tyr-Pro sequence shows a peculiar distribution in proteins, which may suggest a relevant structural role in alpha-helices when Pro is located at the C-cap position. A bioinformatics tool was developed, which updates automatically and periodically the results and makes them available in a web site.

A mixed-α,β miniprotein stereochemically reprogrammed to high-binding affinity for acetylcholine

The protein-structure space is limited to L configuration in the asymmetric alpha-amino acid structures; the function space, on other hand, seems limitless because of the chemical diversity in the amino acid side chain structures. The chemical diversity in side chain structure may be multiplied beneficially with the stereochemical diversity in main chain structure; thus, de novo protein design may be explored for customizing molecular structures stereochemically and molecular functions chemically. Illustrating de novo design in the structure space of L and D alphabet, canonical all-beta folds of poly-L structure were reprogrammed as bracelet, boat, and canoe-shaped molecules-the "boat" as a receptor-like pocket and the "canoe" as a metal-ion receptor-simply by mutating specific L-amino acid residues to the corresponding D stereochemical structure. Demonstrating customization of molecular shape stereochemically and function chemically, a 15-residue mixed-alpha, beta miniprotein of canonical poly-L structure is now reprogrammed stereochemically as a cup-shaped receptor for acetylcholine via cation-pi interaction with a triad of aromatic side chains placed in mimicry of the acetylcholine-receptor sites both natural and artificial. Evidence from CD, fluorescence, NMR, DSC, ITC, MD, and molecular-docking studies is presented to show that a rationally designed 15-residue mixed-L, D peptide is a cooperatively ordered molecular fold in the stereochemically specified molecular morphology, submicromolar in affinity of acetylcholine and thus an acetylcholine receptor exceptionally small and simple. .

Relative strength of cation-pi vs salt-bridge interactions: the Gtalpha(340-350) peptide/rhodopsin system

Interactions between cationic and aromatic side chains of amino acid residues, the so-called cation-pi interaction, are thought to contribute to the overall stability of the folded structure of peptides and proteins. The transferred NOE NMR structure of the G(t)alpha(340-350) peptide bound to photoactivated rhodopsin (R*) geometrically suggests a cation-pi interaction stabilizing the structure between the epsilon-amine of Lys341 and the aromatic ring of the C-terminal residue, Phe350. This interaction has been explored by varying substituents on the phenyl ring to alter the electron density of the aromatic ring of Phe350 and observing the impact on binding of the peptide to R*. The results suggest that while a cation-pi interaction geometrically exists in the G(t)alpha(340-350) peptide when bound to R*, its energetic contribution to the stability of the receptor-bound structure is relatively insignificant, as it was not observed experimentally. The presence of an adjacent and competing salt-bridge interaction between the epsilon-amine of Lys341 and the C-terminal carboxylate of Phe350 effectively shields the charge of the ammonium group. Experimental data supporting a significant cation-pi interaction can be regained through a series of Phe350 analogues where the C-terminal carboxyl has been converted to the neutral carboxamide, thus eliminating the shielding salt-bridge. TrNOE NMR experiments confirmed the existence of the cation-pi interaction in the carboxamide analogues. Various literature estimates of the strength of cation-pi interactions, including some that estimate strengths in excess of salt-bridges, are compromised by omission of the relevant anion in the calculations.

Geometry of nonbonded interactions involving planar groups in proteins

Although hydrophobic interaction is the main contributing factor to the stability of the protein fold, the specificity of the folding process depends on many directional interactions. An analysis has been carried out on the geometry of interaction between planar moieties of ten side chains (Phe, Tyr, Trp, His, Arg, Pro, Asp, Glu, Asn and Gln), the aromatic residues and the sulfide planes (of Met and cystine), and the aromatic residues and the peptide planes within the protein tertiary structures available in the Protein Data Bank. The occurrence of hydrogen bonds and other nonconventional interactions such as C-H...pi, C-H...O, electrophile-nucleophile interactions involving the planar moieties has been elucidated. The specific nature of the interactions constraints many of the residue pairs to occur with a fixed sequence difference, maintaining a sequential order, when located in secondary structural elements, such as alpha-helices and beta-turns. The importance of many of these interactions (for example, aromatic residues interacting with Pro or cystine sulfur atom) is revealed by the higher degree of conservation observed for them in protein structures and binding regions. The planar residues are well represented in the active sites, and the geometry of their interactions does not deviate from the general distribution. The geometrical relationship between interacting residues provides valuable insights into the process of protein folding and would be useful for the design of protein molecules and modulation of their binding properties.

Study on synthesis, structure, and DNA-binding of lanthanide complexes with 2-carboxylbenzaldehyde thiosemicarbazone

2-Carboxylbenzaldehyde thiosemicarbazone (HL), and its three lanthanide (III) complexes, LnL(3) x 4H(2)O [Ln(III)=La, Sm, Eu], have been synthesized in water. The complexes were characterized by elemental analyses, molar conductivity and IR spectra. The crystal structure of [Sm(2)L(6)(CH(3)OH)(4)] x 7.5CH(3)OH x 0.5H(2)O obtained from methanol solution was determined by X-ray diffraction analysis, crystallized in the triclinic system, space group P-1, Z=1, a=12.217 (2)A, b=14.706 (2)A, c=15.035 (2)A, alpha=111.84(1) degrees , beta=103.47(1) degrees , gamma=104.24(1) degrees , R(1)=0.0290. It has symmetrical (mu-OCO)(2), (mu-O)(2) and disamarium(III) units. The coordination geometry of each Sm(III) ion is a distorted tetradecahedron with nine oxygen atoms. In addition, the DNA-binding properties of the ligand and its complexes have been investigated by absorption, fluorescence, and viscosity measurements. The experimental results indicate that the ligand and the Sm-complex can bind to DNA, but the other two complexes cannot; the binding affinity of the Sm-complex is higher than that of the ligand and the intrinsic binding constant K(b) of the complex is 3.22 x 10(5)M(-1).

Circular and linear dichroism of proteins

Circular dichroism (CD) is an important technique in the structural characterisation of proteins, and especially for secondary structure determination. The CD of proteins can be calculated from first principles using the so-called matrix method, with an accuracy which is almost quantitative for helical proteins. Thus, for proteins of unknown structure, CD calculations and experimental data can be used in conjunction to aid structure analysis. Linear dichroism (LD) can be calculated using analogous methodology and has been used to establish the relative orientations of subunits in proteins and protein orientation in an environment such as a membrane. However, simple analysis of LD data is not possible, due to overlapping transitions. So coupling the calculations and experiment is an important strategy. In this paper, the use of LD for the determination of protein orientation and how these data can be interpreted with the aid of calculations, are discussed. We review methods for the calculation of CD spectra, focusing on semiempirical and ab initio parameter sets used in the matrix method. Lastly, a new web interface for online CD and LD calculation is presented.