Conformational preferences of amino acid side chains in collagen

Conformational energy computations were carried out on collagenlike triple-stranded conformations of several poly(tripeptide)s with the general structure CH₃CO? (Gly? X? Y)₃? NHCH₃. The sequences considered had various amino acid residues in position X or Y of the central tripeptide, with either Pro or Ala as a neighbor, i.e., Gly-X-Pro, Gly-X-Ala, Gly-Pro-Y, and Gly-Ala-Y. Minimum-energy conformations were computed for the side chains, and their distributions were compared for the four sequences. The residues used were Abu (= α-aminobutyric acid), Leu, Phe, Ser, Asp, Asn, Val, Ile, and Thr. The conformational energy of a ? Ch₂? CH₃ side chain in Abu was mapped as a function of the dihedral angle χ¹. Intrastrand interactions with neighboring residues do not affect the conformations of a side chain in position Y, and they have a minor effect on it in the X-Ala sequence, but they strongly restrict the conformational freedom of the side chain in the X-Pro sequence. Conversely, interstrand interactions do not affect side chains in position X, but they strongly restrict the conformational freedom of a side chain in position Y if there is a nearby Pro residue in a neighboring strand. Hydrogen bonds with the backbone can be formed in some conformations of long polar side chains, such as Asp, Asn, or Gln. All amino acid residues can be accommodated in collagen. Because of the interactions mentioned above, steric and energetic constraints can be correlated with observed preferences of certain amino acids for positions X or Y in collagen. Hence, these preferences may be explained, in part, in terms of differences in the conformational freedom of the side chains in the triple-stranded structure.     

Computation of Conformational Coupling in Allosteric Proteins

In allosteric regulation, an effector molecule binding a protein at one site induces conformational changes, which alter structure and function at a distant active site. Two key challenges in the computational modeling of allostery are the prediction of the structure of one allosteric state starting from the structure of the other, and elucidating the mechanisms underlying the conformational coupling of the effector and active sites. Here we approach these two challenges using the Rosetta high-resolution structure prediction methodology. We find that the method can recapitulate the relaxation of effector-bound forms of single domain allosteric proteins into the corresponding ligand-free states, particularly when sampling is focused on regions known to change conformation most significantly. Analysis of the coupling between contacting pairs of residues in large ensembles of conformations spread throughout the landscape between and around the two allosteric states suggests that the transitions are built up from blocks of tightly coupled interacting sets of residues that are more loosely coupled to one another.     

Structural properties of amyloid β(1‐40) dimer explored by replica exchange molecular dynamics simulations

Replica exchange molecular dynamics simulations (300 ns) were used to study the dimerization of amyloid β(1‐40) (Aβ(1‐40)) polypeptide. Configurational entropy calculations revealed that at physiological temperature (310 K, 37°C) dynamic dimers are formed by randomly docked monomers. Free energy of binding of the two chains to each other was ?93.56 ± 6.341 kJ mol^?1. Prevalence of random coil conformations was found for both chains with the exceptions of increased β‐sheet content from residues 16‐21 and 29‐32 of chain A and residues 15‐21 and 30‐33 of chain B with β‐turn/β‐bend conformations in both chains from residues 1‐16, 21‐29 of chain A, 1‐16, and 21‐29 of chain B. There is a mixed β‐turn/β‐sheet region from residues 33‐38 of both chains. Analysis of intra‐ and interchain residue distances shows that, although the individual chains are highly flexible, the dimer system stays in a loosely packed antiparallel β‐sheet configuration with contacts between residues 17‐21 of chain A with residues 17‐21 and 31‐36 of chain B as well as residues 31‐36 of chain A with residues 17‐21 and 31‐36 of chain B. Based on dihedral principal component analysis, the antiparallel β‐sheet‐loop‐β‐sheet conformational motif is favored for many low energy sampled conformations. Our results show that Aβ(1‐40) can form dynamic dimers in aqueous solution that have significant conformational flexibility and are stabilized by collapse of the central and C‐terminal hydrophobic cores with the expected β‐sheet‐loop‐β‐sheet conformational motif. Proteins 2017; 85:1024–1045. © 2017 Wiley Periodicals, Inc.     

Charge neutralization in the active site of the catalytic trimer of aspartate transcarbamoylase promotes diverse structural changes

A classical model for allosteric regulation of enzyme activity posits an equilibrium between inactive and active conformations. An alternative view is that allosteric activation is achieved by increasing the potential for conformational changes that are essential for catalysis. In the present study, substitution of a basic residue in the active site of the catalytic (C) trimer of aspartate transcarbamoylase with a non‐polar residue results in large interdomain hinge changes in the three chains of the trimer. One conformation is more open than the chains in both the wild‐type C trimer and the catalytic chains in the holoenzyme, the second is closed similar to the bisubstrate‐analog bound conformation and the third hinge angle is intermediate to the other two. The active‐site 240s loop conformation is very different between the most open and closed chains, and is disordered in the third chain, as in the holoenzyme. We hypothesize that binding of anionic substrates may promote similar structural changes. Further, the ability of the three catalytic chains in the trimer to access the open and closed active‐site conformations simultaneously suggests a cyclic catalytic mechanism, in which at least one of the chains is in an open conformation suitable for substrate binding whereas another chain is closed for catalytic turnover. Based on the many conformations observed for the chains in the isolated catalytic trimer to date, we propose that allosteric activation of the holoenzyme occurs by release of quaternary constraint into an ensemble of active‐site conformations.     

Exploring the conformational space of protein side chains using dead-end elimination and the A* algorithm

We describe an algorithm which enables us to search the conformational space of the side chains of a protein to identify the global minimum energy combination of side chain conformations as well as all other conformations within a specified energy cutoff of the global energy minimum. The program is used to explore the side chain conformational energy surface of a number of proteins, to investigate how this surface varies with the energy model used to describe the interactions within the system and the rotamer library. Enumeration of the rotamer combinations enables us to directly evaluate the partition function, and thus calculate the side chain contribution to the conformational entropy of the folded protein. An investigation of these conformations and the relationships between them shows that most of the conformations near to the global energy minimum arise from changes in side chain conformations that are essentially independent; very few result from a concerted change in conformation of two or more residues. Some of the limitations of the approach are discussed. Proteins 33:227–239, 1998. © 1998 Wiley-Liss, Inc.     

Structure of a designed tetrahedral protein assembly variant engineered to have improved soluble expression

We recently reported the development of a computational method for the design of coassembling multicomponent protein nanomaterials. While four such materials were validated at high‐resolution by X‐ray crystallography, low yield of soluble protein prevented X‐ray structure determination of a fifth designed material, T33‐09. Here we report the design and crystal structure of T33‐31, a variant of T33‐09 with improved soluble yield resulting from redesign efforts focused on mutating solvent‐exposed side chains to charged amino acids. The structure is found to match the computational design model with atomic‐level accuracy, providing further validation of the design approach and demonstrating a simple and potentially general means of improving the yield of designed protein nanomaterials.     

A comparative analysis of the equilibrium dynamics of a designed protein inferred from NMR,X‐ray,and computations

A detailed analysis of high‐resolution structural data and computationally predicted dynamics was carried out for a designed sugar‐binding protein. The mean‐square deviations in the positions of residues derived from nuclear magnetic resonance (NMR) models and those inferred from X‐ray crystallographic B‐factors for two different crystal forms were compared with the predictions based on the Gaussian Network Model (GNM) and the results from molecular dynamics (MD) simulations. GNM systematically yielded a higher correlation than MD, with experimental data, suggesting that the lack of atomistic details in the coarse‐grained GNM is more than compensated for by the mathematically exact evaluation of fluctuations using the native contacts topology. Evidence is provided that particular loop motions are curtailed by intermolecular contacts in the crystal environment causing a discrepancy between theory and experiments. Interestingly, the information conveyed by X‐ray crystallography becomes more consistent with NMR models and computational predictions when ensembles of X‐ray models are considered. Less precise (broadly distributed) ensembles indeed appear to describe the accessible conformational space under native state conditions better than B‐factors. Our results highlight the importance of using multiple conformations obtained by alternative experimental methods, and analyzing results from both coarse‐grained models and atomic simulations, for accurate assessment of motions accessible to proteins under native state conditions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Structural properties and photophysical behavior of conformationally constrained hexapeptides functionalized with a new fluorescent analog of tryptophan and a nitroxide radical quencher

The influence of the conformational properties on the photophysics of two de novo designed hexapeptides was studied by spectroscopic measurements (ir, NMR, steady-state, and time resolved fluorescence) and molecular mechanics calculations. The peptide sequences comprise two nonproteinogenic residues: a beta-(1-azulenyl)-L-alanine (Aal) residue, obtained by formally functionalizing the Ala side chain with the azulene chromophore, and a Calpha-tetrasubstituted alpha-amino acid (TOAC), incorporating a nitroxide group in a cycloalkyl moiety. Aal represents a new fluorescent, quasi-isosteric Trp analog and TOAC a stable radical species, frequently used as a paramagnetic probe in biochemical studies. The peptide chains differ in the sequence position of the two probes and are heavily based on Aib (alpha-aminoisobutyric acid) residues to generate conformationally restricted helical structures, as confirmed by both spectroscopic and computational results. The conformationally controlled, excited state interactions, determining the photophysical relaxation of the Aal*/TOAC pair, are also discussed.     

Characterization of (Boc‐Cys/Sec‐NHMe)2 and (Boc‐Cys/Sec‐OMe)2: Evidence of local conformational difference between disulfide and diselenide

Conformations of disulfide and diselenide were compared in (Boc‐Cys/Sec‐NHMe)₂ and (Boc‐Cys/Sec‐OMe)₂ using X‐ray crystallography, nuclear magnetic resonance (NMR) spectroscopy, density functional theory (DFT), and circular dichroism (CD) spectroscopy. Conformations of disulfide/diselenide in polypeptides are defined based on the sign of side chain torsion angle χ₃ (–CH₂–S/Se–S/Se–CH₂–); negative indicates left‐handed and positive indicates right‐handed orientation. In the crystals of (Boc‐Cys‐OMe)₂ and (Boc‐Sec‐OMe)₂, the disulfide exhibits a left‐handed and the diselenide a right‐handed orientation. Characterization of cystine and selenocystine derivatives in solution using ¹H‐NMR, natural abundant ⁷⁷Se NMR, 2D‐ROESY, and chemical shift analysis coupled to DMSO titration has indicated the symmetrical nature and antiparallel orientation of Cys/Sec residues about the disulfide/diselenide bridges. Structural calculations of cystine and selenocystine derivatives using DFT further support the antiparallel orientation of Cys/Sec residues about disulfide/diselenide. The far‐ultraviolet (UV) region CD spectra of cystine and selenocystine derivatives have exhibited the negative Cotton effect (CE) for disulfide and positive for diselenide confirming the difference in the conformational preference of disulfide and diselenide. In the previously reported polymorphic structure of (Boc‐Sec‐OMe)₂, the diselenide has right‐handed orientation. In the X‐ray structures of disulfide and diselenide analogues of Escherichia coli protein encoded by curli specific gene C (CgsC) retrieved from Protein Databank (PDB), disulfide has left‐handed and the diselenide right‐handed orientation. The current report provides the evidence for the local conformational difference between a disulfide and a diselenide group under unconstrained conditions, which may be useful for the rational replacement of disulfide by diselenide in polypeptide chains.     

New type of helix and 2(7) ribbon structure formation in poly DeltaLeu peptides: construction of a single-handed template

alpha,beta-Dehydroamino acid residues due to the presence of Calpha = Cbeta double bond influences the main chain and the side chain conformations. These residues have interesting chemical features including the increased resistance to enzymatic degradation. The chain length dependent conformational behavior of poly alpha,beta-dehydroleucine (DeltaLeu) peptides in both the pure forms Z and E and their various combinations like alternate ZE/EZ etc. have been investigated by using quantum mechanical method PCILO (perturbative configuration interaction of localized orbitals). The conformational states in alternate Z and E forms, with Phi, Psi values of -10 degrees , 105 degrees /1 degrees , -88 degrees for Z form and 35 degrees , 22 degrees /-34 degrees , -27 degrees for the E form are found to be the most stable and degenerate than the states in pure Z and E forms and the EZ form etc. The repeated Phi, Psi values give rise to altogether new types of left and right handed helices, and their stability increases with increasing chain length. These structures are stabilized by intramolecular hydrogen bonding, carbonyl-carbonyl interactions and hydrophobic interactions between the side chains of DeltaZLeu and DeltaELeu residues. The 2(7) ribbon structure (seven-membered hydrogen-bonded ring involving two consecutive amino acid residues) is found to be most stable and degenerate in the pentapeptide Ac-DeltaELeu5-NHMe, due to the formation of maximum hydrogen bonds. A right-handed template from achiral DeltaLeu peptides has been achieved by incorporating L-Leu at the C-terminal or D-Leu at the N-terminal.     

Computational studies on the backbone-dependent side-chain orientation induced by the (S,S)-CXC motif

Disulfide cyclization is a well-known procedure to impose conformational restriction to peptides undergoing backbone flexibility. Rigid conformations are induced only for small rings with a specific combination of amino acids. In this work, we present a computational search of the backbone and backbone-dependent side-chain orientation of two series of linear and cyclic peptide analogs. The -C[XY]C- scaffold (where X,Y is arginine, aspartic acid or alanine residue) in its open and (S,S) cyclic form was used for the design of the studied analogs. Thirty-six compounds, resulting from the extension with one residue at either the N- or the C-terminus were studied with classical MD. The local backbone conformation and the relative orientation of the X and Y side chains induced by either cyclization and/or the presence of the charged residues are discussed. From the present study it is concluded that cyclization has a great impact on the synplanar orientation of the X and Y side chains in the (S,S)Ac-XCYC-NH2 series of compounds while charge-charge interaction has only a weak synergic effect. On the contrary, the antiplanar orientation is favored in the case of (S,S)Ac-CXCY-NH2.     

Preferred conformations of RGDX tetrapeptides to inhibit the binding of fibrinogen to platelets

The conformational study on Arg-Gly-Asp (RGD)-containing tetrapeptides in the unhydrated and hydrated states has been carried out using the force field ECEPP/3 and the hydration shell model. The tetrapeptides studied here are H-RGDX-OH (X = Trp, Tyr, Phe, Leu, Val, Cys, Gln, and Ser), which show the inhibitory activity for binding of fibrinogen to platelets in the order of RGDW approximately equal to RGDY approximately equal to RGDF approximately equal to RGDL > RGDV > or = RGDC > or = RGDQ > or = RGDS. The backbone conformations with two C(7) backbone-to-backbone hydrogen bonds between Asp and Arg residues and between Xaa and Gly residues are in common most probable for the RGD sequence of RGDX tetrapeptides in the hydrated state. The dominant beta-turns for RGDX are found to be the types V' and IV at Gly-Asp and Asp-Xaa sequences, respectively, which are quite similar to the types II' and I (or II), respectively. However, it cannot be ruled out that the extended conformations are also remarkably feasible for RGDX tetrapeptides in water by peering the distributions of backbone conformations. These calculated results are consistent with the experimental results on RGD-containing proteins and conformationally constrained RGD-containing peptides. The reason why the RGDX becomes more potent as the side chain of the X residue is more hydrophobic may be ascribed to that the more hydrophobic is the residue X, the more populated are beta-turn structures for the Gly-Asp sequence. The hydrophobic side chain of X residue exposed to water is likely to interact with the hydrophobic region of receptor easily.     

NMR characterization of the conformational fluctuations of the human lymphocyte function‐associated antigen‐1 I‐domain

Lymphocyte function‐associated antigen‐1 (LFA‐1) is an integrin protein that transmits information across the plasma membrane through the so‐called inside‐out and outside‐in signaling mechanisms. To investigate these mechanisms, we carried out an NMR analysis of the dynamics of the LFA‐1 I‐domain, which has enabled us to characterize the motions of this domain on a broad range of timescales. We studied first the internal motions on the nanosecond timescale by spin relaxation measurements and model‐free analysis. We then extended this analysis to the millisecond timescale motions by measuring ¹⁵N‐¹H residual dipolar couplings of the backbone amide groups. We analyzed these results in the context of the three major conformational states of the I‐domain using their corresponding X‐ray crystallographic structures. Our results highlight the importance of the low‐frequency motions of the LFA‐1 I‐domain in the inactive apo‐state. We found in particular that α‐helix 7 is in a position in the apo‐closed state that cannot be fully described by any of the existing X‐ray structures, as it appears to be in dynamic exchange between different conformations. This type of motion seems to represent an inherent property of the LFA‐1 I‐domain and might be relevant for controlling the access to the allosteric binding pocket, as well as for the downward displacement of α‐helix 7 that is required for the activation of LFA‐1.     

Effect of subdomain interactions on methyl group dynamics in the hydrophobic core of villin headpiece protein

Thermostable villin headpiece protein (HP67) consists of the N‐terminal subdomain (residues 10–41) and the autonomously folding C‐terminal subdomain (residues 42–76) which pack against each other to form a structure with a unified hydrophobic core. The X‐ray structures of the isolated C‐terminal subdomain (HP36) and its counterpart in HP67 are very similar for the hydrophobic core residues. However, fine rearrangements of the free energy landscape are expected to occur because of the interactions between the two subdomains. We detect and characterize these changes by comparing the µs‐ms time scale dynamics of the methyl‐bearing side chains in isolated HP36 and in HP67. Specifically, we probe three hydrophobic side chains at the interface of the two subdomains (L42, V50, and L75) as well as at two residues far from the interface (L61 and L69). Solid‐state deuteron NMR techniques are combined with computational modeling for the detailed characterization of motional modes in terms of their kinetic and thermodynamic parameters. The effect of interdomain interactions on side chain dynamics is seen for all residues but L75. Thus, changes in dynamics because of subdomain interactions are not confined to the site of perturbation. One of the main results is a two‐ to threefold increase in the value of the activation energies for the rotameric mode of motions in HP67 compared with HP36. Detailed analysis of configurational entropies and heat capacities complement the kinetic view of the degree of the disorder in the folded state.     

Empirical evaluation of the influence of side chains on the conformational entropy of the polypeptide backbone

Changes in amino acid side chains have long been recognized to alterthe range and distribution of ?, ψ angles found in the main chain of polypeptides. Altering the range and distribution of ?, ψ angles also alters the conformational entropy of the flexible denatured state and may thus stabilize or destabilize it relative to the comparatively conformationally rigid native state. A database of 12,320 residues from 61 nonhomologous, high resolution crystal structures was examined to determine the ?, ψ conformational preferences of each of the 20 amino acids. These observed distributions in the native state of proteins are assumed to also reflect the distributions found in the denatured state. The distributionswere used to approximate the energy surface for each residue, allowing the calculation of relative conformational entropies for each residue relative to glycine. In the most extreme case, replacement of glycine by proline, conformational entropy changes will stabilize the native state relative to the denatured state by ?0.82 ± 0.08 kcal/mol at 20°C. Surprisingly, alanine is found to be the most ordered residue other than proline. This unexpected result is a result of the high percentage of alanines found in helical conformations. This either indicates that the observed distributions in the native state do not reflect the distributions in the denatured state, or that alanine is much more likely to adopt a helical conformation in the denatured state than residues with longer side chains. Among those residues with ?, ψ angles compatible with helix incorporation the percentage of alanines actually in helices is very similar to other residues. This and the consistent ordering of alanine relative to other residues regardless of secondary structure are evidence that ?, ψ distributions in native states reflect those in the denatured states. © 1995 Wiley-Liss, Inc.     

Structural basis of hierarchical multiple substates of a protein. III: Side chain and main chain local conformations

An analysis is carried out of differences in the minimum energy conformations obtained in the previous paper by energy minimization starting from conformations sampled by a Monte Carlo simulation of conformational fluctuations in the native state of a globular protein, bovine pancreatic trypsin inhibitor. Main conformational differences in each pair of energy minima are found usually localized in several side chains and in a few local main chain segments. Such side chains and local main chain segments are found to take a few distinct local conformations in the minimum energy conformations. Energy minimum conformations can thus be described in terms of combinations of these multiple local conformations.     

Peptide conformations--49(1): synthesis and structure-activity relationships of side chain modified peptides of cyclo(-D-Pro-Phe-Thr-Lys-Trp-Phe.)

Cyclic hexapeptide analogues representing the modified retro sequence of the amino acid residues 7-11 of natural somatostatin are known to protect liver cells from phalloidin poisoning. To determine the influence of steric, lipophilic, and charge effects on (a) the conformation of the backbone and the aromatic side chains and (b) the biological response, the side chains of Phe2, Lys4, and Phe6 of cyclo(-D-Pro1-Phe2-Thr3-Lys(Z)4-Trp5-Phe6-), 1a, one of the most active peptides found so far, were modified by various residues. The discussion of conformationally relevant parameters proves that neither backbone conformations nor populations of aromatic side chain rotamers were altered by these substitutions. The potency of these derivatives in a cytoprotection assay varies by at most one order of magnitude (more or less active than the parent peptide 1a). A qualitative evaluation of lipophilic, steric, and charge effects reveals the dominance of lipophilic effects of aromatic residues; the most potent compounds contain aromatic substructures in the side chain of Lys4.     

Conformation analysis of the N-glycosylation site Asn-X-Thr/Ser in glycoproteins

Theoretical conformational analysis of oligopeptides CH3CO-Asn-X-Thr-NHCH3 (X = Gly, Ala, Pro), modelling N-glycosylation site, and their glycosylated derivatives CH3CO-(GlcNAc beta 1-4GlcNAc beta 1) Asn-X-Thr-NHCH3 has been carried out. Active conformations of the site are found, corresponding to structural prerequisities of N-glycosylation: Asn residue's position in beta-turn and hydrogen bond formation between side chains of Asn and Thr/Ser residues. In this case the L conformation of the central residue X is most probable. Since Pro residue does not possess this conformation, sequences with X = Pro are not glycosylated. It is shown that glycosylation of the above-mentioned sites is accompanied by reorientation of the Asn residue's side chains.