Energy minimization method using automata network for sequence and side-chain conformation prediction from given backbone geometry

Globular proteins have high packing densities as a result of residue side chains in the core achieving a tight, complementary packing. The internal packing is considered the main determinant of native protein structure. From that point of view, we present here a method of energy minimization using an automata network to predict a set of amino acid sequences and their side-chain conformations from a desired backbone geometry for de novo design of proteins. Using discrete side-chain conformations, that is, rotamers, the sequence generation problem from a given backbone geometry becomes one of combinatorial problems. We focused on the residues composing the interior core region and predicted a set of amino acid Sequences and their side-chain conformations only from a given backbone geometry. The kinds of residues were restricted to six hydrophobic amino acids (Ala, Ile, Met, Leu, Phe, and Val) because the core regions are almost always composed of hydrophobic residues. The obtained sequences were well packed as was the native sequence. The method can be used for automated sequence generation in the de novo design of proteins. © 1994 Wiley-Liss, Inc.     

Statistical and conformational analysis of the electron density of protein side chains

Protein side chains make most of the specific contacts between proteins and other molecules, and their conformational properties have been studied for many years. These properties have been analyzed primarily in the form of rotamer libraries, which cluster the observed conformations into groups and provide frequencies and average dihedral angles for these groups. In recent years, these libraries have improved with higher resolution structures and using various criteria such as high thermal factors to eliminate side chains that may be misplaced within the crystallographic model coordinates. Many of these side chains have highly non-rotameric dihedral angles. The origin of side chains with high B-factors and/or with non-rotameric dihedral angles is of interest in the determination of protein structures and in assessing the prediction of side chain conformations. In this paper, using a statistical analysis of the electron density of a large set of proteins, it is shown that: (1) most non-rotameric side chains have low electron density compared to rotameric side chains; (2) up to 15% of chi1 non-rotameric side chains in PDB models can clearly be fit to density at a single rotameric conformation and in some cases multiple rotameric conformations; (3) a further 47% of non-rotameric side chains have highly dispersed electron density, indicating potentially interconverting rotameric conformations; (4) the entropy of these side chains is close to that of side chains annotated as having more than one chi(1) rotamer in the crystallographic model; (5) many rotameric side chains with high entropy clearly show multiple conformations that are not annotated in the crystallographic model. These results indicate that modeling of side chains alternating between rotamers in the electron density is important and needs further improvement, both in structure determination and in structure prediction.     

3d interaction homology: The structurally known rotamers of tyrosine derive from a surprisingly limited set of information‐rich hydropathic interaction environments described by maps

Sidechain rotamer libraries are obtained through exhaustive statistical analysis of existing crystallographic structures of proteins and have been applied in multiple aspects of structural biology, for example, crystallography of relatively low‐resolution structures, in homology model building and in biomolecular NMR. Little is known, however, about the driving forces that lead to the preference or suitability of one rotamer over another. Construction of 3D hydropathic interaction maps for nearly 30,000 tyrosines reveals the environment around each, in terms of hydrophobic (π–π stacking, etc.) and polar (hydrogen bonding, etc.) interactions. After partitioning the tyrosines into backbone‐dependent (?, ψ) bins, a map similarity metric based on the correlation coefficient was applied to each map‐map pair to build matrices suitable for clustering with k‐means. The first bin (?200° ≤ ? < –155°; ?205° ≤ ψ < –160°), representing 631 tyrosines, reduced to 14 unique hydropathic environments, with most diversity arising from favorable hydrophobic interactions with many different residue partner types. Polar interactions for tyrosine include surprisingly ubiquitous hydrogen bonding with the phenolic OH and a handful of unique environments surrounding the tyrosine backbone. The memberships of all but one of the 14 environments are dominated (>50%) by a single χ₁/χ₂ rotamer. The last environment has weak or no interactions with the tyrosine ring and its χ₁/χ₂ rotamer is indeterminate, which is consistent with it being composed of mostly surface residues. Each tyrosine residue attempts to fulfill its hydropathic valence and thus, structural water molecules are seen in a variety of roles throughout protein structure. Proteins 2015; 83:1118–1136. © 2015 Wiley Periodicals, Inc.     

Organization of the multiaminoacyl-tRNA synthetase complex and the cotranslational protein folding

Aminoacyl-tRNA synthetases (ARSs) play an essential role in the protein synthesis by catalyzing an attachment of their cognate amino acids to tRNAs. Unlike their prokaryotic counterparts, ARSs in higher eukaryotes form a multiaminoacyl-tRNA synthetase complex (MARS), consisting of the subset of ARS polypeptides and three auxiliary proteins. The intriguing feature of MARS complex is the presence of only nine out of twenty ARSs, specific for Arg, Asp, Gln, Glu, Ile, Leu, Lys, Met, and Pro, regardless of the organism, cell, or tissue types. Although existence of MARSs complex in higher eukaryotes has been already known for more than four decades, its functional significance remains elusive. We found that seven of the nine corresponding amino acids (Arg, Gln, Glu, Ile, Leu, Lys, and Met) together with Ala form a predictor of the protein α-helicity. Remarkably, all amino acids (besides Ala) in the predictor have the highest possible number of side-chain rotamers. Therefore, compositional bias of a typical α-helix can contribute to the helix''s stability by increasing the entropy of the folded state. It also appears that position-specific α-helical propensity, specifically periodic alternation of charged and hydrophobic residues in the helices, may well be provided by the structural organization of the complex. Considering characteristics of MARS complex from the perspective of the α-helicity, we hypothesize that specific composition and structure of the complex represents a functional mechanism for coordination of translation with the fast and correct folding of amphiphilic α-helices.     

Dead-End Based Modeling Tools to Explore the Sequence Space That Is Compatible with a Given Scaffold

The dead-end elimination algorithm has proven to be a powerful tool in protein homology modeling since it allows one to determine rapidly the global minimum-energy conformation (GMEC) of an arbitrarily large collection of side chains, given fixed backbone coordinates. After introducing briefly the necessary background, we focus on logic arguments that increase the efficacy of the dead-end elimination process. Second, we present new theoretical considerations on the use of the dead-end elimination method as a tool to identify sequences that are compatible with a given scaffold structure. Third, we initiate a search for properties derived from the computed GMEC structure to predict whether a given sequence can be well packed in the core of a protein. Three properties will be considered: the nonbonded energy, the accessible surface area, and the extent by which the GMEC side-chain conformations deviate from a locally optimal conformation.     

Free‐energy profiles for ions in the influenza M2‐TMD channel

M₂ transmembrane domain channel (M₂‐TMD) permeation properties are studied using molecular dynamics simulations of M₂‐TMD (1NYJ) embedded in a lipid bilayer (DMPC) with 1 mol/kg NaCl or KCl saline solution. This study allows examination of spontaneous cation and anion entry into the selectivity filter. Three titration states of the M₂‐TMD tetramer are modeled for which the four His³⁷ residues, forming the selectivity filter, are net uncharged, +2 charged, or +3 charged. M₂‐TMD structural properties from our simulations are compared with the properties of other models extracted from NMR and X‐ray studies. During 10 ns simulations, chloride ions occasionally occupy the positively‐charged selectivity filter region, and from umbrella sampling simulations, Cl^? has a lower free‐energy barrier in the selectivity‐filter region than either Na⁺ or NH, and NH has a lower free‐energy barrier than Na⁺. For Na⁺ and Cl^?, the free‐energy barriers are less than 5 kcal/mol, suggesting that the 1NYJ conformation would probably not be exquisitely proton selective. We also point out a rotameric configuration of Trp⁴¹ that could fully occlude the channel. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Bayesian statistical analysis of protein side-chain rotamer preferences.

We present a Bayesian statistical analysis of the conformations of side chains in proteins from the Protein Data Bank. This is an extension of the backbone-dependent rotamer library, and includes rotamer populations and average chi angles for a full range of phi, psi values. The Bayesian analysis used here provides a rigorous statistical method for taking account of varying amounts of data. Bayesian statistics requires the assumption of a prior distribution for parameters over their range of possible values. This prior distribution can be derived from previous data or from pooling some of the present data. The prior distribution is combined with the data to form the posterior distribution, which is a compromise between the prior distribution and the data. For the chi 2, chi 3, and chi 4 rotamer prior distributions, we assume that the probability of each rotamer type is dependent only on the previous chi rotamer in the chain. For the backbone-dependence of the chi 1 rotamers, we derive prior distributions from the product of the phi-dependent and psi-dependent probabilities. Molecular mechanics calculations with the CHARMM22 potential show a strong similarity with the experimental distributions, indicating that proteins attain their lowest energy rotamers with respect to local backbone-side-chain interactions. The new library is suitable for use in homology modeling, protein folding simulations, and the refinement of X-ray and NMR structures.     

Total synthesis,structural, and biological evaluation of stylissatin A and related analogs

The natural product cyclic peptide stylissatin A ( 1a ) was reported to inhibit nitric oxide production in LPS‐stimulated murine macrophage RAW 264.7 cells. In the current study, solid‐phase total synthesis of stylissatin A was performed by using a safety‐catch linker and yielded the peptide with a trans‐Phe⁷‐Pro⁶ linkage, whereas the natural product is the cis rotamer at this position as evidenced by a marked difference in NMR chemical shifts. In order to preclude the possibility of 1b being an epimer of the natural product, we repeated the synthesis using d ‐allo‐Ile in place of l ‐Ile and a different site for macrocyclization. The resulting product (d ‐allo‐Ile²)‐stylissatin A ( 1c ) was also found to have the trans‐Phe⁷‐Pro⁶ peptide conformations like rotamer 1b . Applying the second route to the synthesis of stylissatin A itself, we obtained stylissatin A natural rotamer 1a accompanied by rotamer 1b as the major product. Rotamers 1a , 1b , and the epimer 1c were separable by HPLC, and 1a was found to match the natural product in structure and biological activity. Six related analogs 2–7 of stylissatin A were synthesized on Wang resin and characterized by spectral analysis. The natural product ( 1a ), the rotamer ( 1b ), and (d ‐allo‐Ile²)‐stylissatin A ( 1c ) exhibited significant inhibition of NO^.. Further investigations were focused on 1b , which also inhibited proliferation of T‐cells and inflammatory cytokine IL‐2 production. The analogs 2–7 weakly inhibited NO^. production, but strongly inhibited IL‐2 cytokine production compared with synthetic peptide 1b . All analogs inhibited the proliferation of T‐cells, with analog 7 having the strongest effect. In the analogs, the Pro⁶ residue was replaced by Glu/Ala, and the SAR indicates that the nature of this residue plays a role in the biological function of these peptides. Copyright © 2016 European Peptide Society and John Wiley & Sons, Ltd.     

Metadynamics modelling of the solvent effect on primary hydroxyl rotamer equilibria in hexopyranosides

Accurate modelling of rotamer equilibria for the primary hydroxyl groups of monosaccharides continues to be a great challenge of computational glycochemistry. The metadynamics technique was applied to study the conformational free energy surfaces of methyl α-d-glucopyranoside and methyl α-d-galactopyranoside, employing the glycam06 force field. For both molecules, seven to eight conformational free-energy minima, differing in the ω (O-5–C-5–C-6–O-6) and χ (C-3–C-4–O-4–HO-4) dihedral angles, were identified in vacuum or in a water environment. The calculated rotamer equilibrium of the primary hydroxyl group is significantly different in vacuum than in water. The major effect of a water environment is the destabilisation of a hydrogen bond between O-4–HO-4 and O-6–HO-6 groups. It was possible to calculate the free-energy differences of individual rotamers with an accuracy of better than 2 kJ/mol. The calculated gg, gt and tg rotamer populations in water are in close agreement with experimental measurements, and therefore support the theoretical background of metadynamics.     

The role of carbon-donor hydrogen bonds in stabilizing tryptophan conformations

Although most side-chain torsion angles correspond to low-energy rotameric positions, deviations occur with significant frequency. One striking example arises in Trp residues, which have an important role in stabilizing protein structures because of their size and mixed hydrophobic/hydrophilic character. Ten percent of Trp side-chains have unexplained conformations with chi(2) near 0 degrees instead of the expected 90 degrees. The current work is a structural and energetic analysis of these conformations. It is shown that many Trp residues with these orientations are stabilized by three-center carbon-donor hydrogen bonds of the form C-H...X...H-C, where X is a polar hydrogen-bond acceptor in the environment of the side-chain. The bridging hydrogen bonds occur both within the Trp side-chain and between the side-chain and the local protein backbone. Free energy maps of an isolated Trp residue in an explicit water environment show a minimum corresponding to the off-rotamer peak observed in the crystallographic data. Bridging carbon-donor hydrogen bonds are also shown to stabilize on-rotamer Trp conformations, and similar bridging hydrogen bonds also stabilize some off-rotamer Asp conformations. The present results suggest a previously unrecognized role for three-center carbon-donor hydrogen bonds in protein structures and support the view that the off-rotamer Trp side-chain orientations are real rather than artifacts of crystallographic refinements. Certain of the off-rotamer Trp conformations appear to have a functional role.