Combinations Of Dihedral Angles Research Articles

Equilibrium transitions between side-chain conformations in leucine and isoleucine.

Despite recent improvements in computational methods for protein design, we still lack a quantitative, predictive understanding of the intrinsic probabilities for amino acids to adopt particular side-chain conformations. Surprisingly, this question has remained unsettled for many years, in part because of inconsistent results from different experimental approaches. To explicitly determine the relative populations of different side-chain dihedral angles, we performed all-atom hard-sphere Langevin Dynamics simulations of leucine (Leu) and isoleucine (Ile) dipeptide mimetics with stereo-chemical constraints and repulsive-only steric interactions between non-bonded atoms. We determine the relative populations of the different χ(1) and χ(2) dihedral angle combinations as a function of the backbone dihedral angles ϕ and ψ. We also propose, and test, a mechanism for inter-conversion between the different side-chain conformations. Specifically, we discover that some of the transitions between side-chain dihedral angle combinations are very frequent, whereas others are orders of magnitude less frequent, because they require rare coordinated motions to avoid steric clashes. For example, to transition between different values of χ(2), the Leu side-chain bond angles κ(1) and κ(2) must increase, whereas to transition in χ(1), the Ile bond angles λ(1) and λ(2) must increase. These results emphasize the importance of computational approaches in stimulating further experimental studies of the conformations of side-chains in proteins. Moreover, our studies emphasize the power of simple steric models to inform our understanding of protein structure, dynamics, and design.

New Insights into the Interdependence between Amino Acid Stereochemistry and Protein Structure

To successfully design new proteins and understand the effects of mutations in natural proteins, we must understand the geometric and physicochemical principles underlying protein structure. The side chains of amino acids in peptides and proteins adopt specific dihedral angle combinations; however, we still do not have a fundamental quantitative understanding of why some side-chain dihedral angle combinations are highly populated and others are not. Here we employ a hard-sphere plus stereochemical constraint model of dipeptide mimetics to enumerate the side-chain dihedral angles of leucine (Leu) and isoleucine (Ile), and identify those conformations that are sterically allowed versus those that are not as a function of the backbone dihedral angles ϕ and ψ. We compare our results with the observed distributions of side-chain dihedral angles in proteins of known structure. With the hard-sphere plus stereochemical constraint model, we obtain agreement between the model predictions and the observed side-chain dihedral angle distributions for Leu and Ile. These results quantify the extent to which local, geometrical constraints determine protein side-chain conformations.

163 Enhanced sampling of peptides and proteins with a new biasing replica exchange method

Classical MD simulations (cMD) are limited by the sampling of relevant states of the peptides. Replica exchange (REMD) methods aim to search the conformational space of proteins more efficiently (reviewed in Ostermeir & Zacharias, 2013). We have developed a Hamiltonian REMD method that takes advantage of an intrinsic property of proteins, the specific Φ ϕ dihedral angle combinations along the polymer backbone. By employing a coupled two-dimensional biasing potential the energy barriers along the polymer backbone are reduced more effectively than by a previous approach based on a one-D biasing potential (Kannan & Zacharias, 2007). Thus, adjacent amino acids along the polymers backbone can easily switch between favourable regions in the Ramachandran plot. Additionally, energy barriers of rotameric states of amino acid side chains of proteins are also biased in the replica runs. The method improves the sampling of conformational substates of proteins at a modest number of replicas (nine replicas in the standard set-up with one replica running without biasing potential) compared to much larger numbers necessary in the case of standard temperature (T)-REMD simulations. A further improvement is achieved by a dynamical adjustment of the penalty potential levels in the replicas such that high exchange rates and improved mixing of conformations between different replicas are guaranteed. The biasing potential (BP)-REMD method turns out to be suitable to speed up both the folding of spaghetti-like test peptides and the refinement of loop decoy structures. Starting from extended structures, an α-helical oligo-alanine and β-hairpin chignolin and the Trp-cage protein fold more rapidly in near-native structures than in cMD simulations. The BP-REMD simulations not only accelerate the folding process of test proteins but also enlarge the variety of sampled configurations in conformational space. Since flexible parts of the protein can be penalized selectively, this method provides a precise tool to investigate regions of interest of the protein.

Motivated Proteins: A web application for studying small three-dimensional protein motifs

BackgroundSmall loop-shaped motifs are common constituents of the three-dimensional structure of proteins. Typically they comprise between three and seven amino acid residues, and are defined by a combination of dihedral angles and hydrogen bonding partners. The most abundant of these are αβ-motifs, asx-motifs, asx-turns, β-bulges, β-bulge loops, β-turns, nests, niches, Schellmann loops, ST-motifs, ST-staples and ST-turns.We have constructed a database of such motifs from a range of high-quality protein structures and built a web application as a visual interface to this.DescriptionThe web application, Motivated Proteins, provides access to these 12 motifs (with 48 sub-categories) in a database of over 400 representative proteins. Queries can be made for specific categories or sub-categories of motif, motifs in the vicinity of ligands, motifs which include part of an enzyme active site, overlapping motifs, or motifs which include a particular amino acid sequence. Individual proteins can be specified, or, where appropriate, motifs for all proteins listed. The results of queries are presented in textual form as an (X)HTML table, and may be saved as parsable plain text or XML. Motifs can be viewed and manipulated either individually or in the context of the protein in the Jmol applet structural viewer. Cartoons of the motifs imposed on a linear representation of protein secondary structure are also provided. Summary information for the motifs is available, as are histograms of amino acid distribution, and graphs of dihedral angles at individual positions in the motifs.ConclusionMotivated Proteins is a publicly and freely accessible web application that enables protein scientists to study small three-dimensional motifs without requiring knowledge of either Structured Query Language or the underlying database schema.

Material-aware differential mesh deformation using sketching interface

In this paper, we present a material-aware mesh deformation method using a sketching interface. Guided by user-specified material properties, our method can deform the surface mesh in a non-uniform way, while previous deformation techniques are mainly designed for uniform materials. The non-uniform deformation is achieved by material-dependent gradient field manipulation and Poisson-based reconstruction. Compared with previous material-oblivious deformation techniques, our method supplies better control of the deformation process and can generate more realistic results. We propose a novel detail representation that transforms geometric details between successive surface levels as a combination of dihedral angles and barycentric coordinates. This detail representation is similarity-invariant and fully compatible with material properties. Based on these two methods, we implement a multi-resolution deformation tool, allowing the user to edit a mesh inside a hierarchy in a material-aware manner. We demonstrate the effectiveness and robustness of our methods by several examples with real-world data.

Efficient Search on Energy Minima for Structure Prediction of Nucleic Acid Motifs

Structure prediction of non-canonical motifs such as mismatches, extra unmatched nucleotides or internal and hairpin loop structures in nucleic acids is of great importance for understanding the function and design of nucleic acid structures. Systematic conformational analysis of such motifs typically involves the generation of many possible combinations of backbone dihedral torsion angles for a given motif and subsequent energy minimization (EM) and evaluation. Such approach is limited due to the number of dihedral angle combinations that grows very rapidly with the size of the motif. Two conformational search approaches have been developed that allow both an effective crossing of barriers during con-formational searches and the computational demand grows much less with system size then search methods that explore all combinations of backbone dihedral torsion angles. In the first search protocol single torsion angles are flipped into favorable states using constraint EM and subsequent relaxation without constraints. The approach is repeated in an iterative manner along the backbone of the structural motif until no further energy improvement is obtained. In case of two test systems, a DNA-trinucleotide loop (sequence: GCA) and a RNA tetraloop (sequence: UUCG), the approach successfully identified low energy states close to experiment for two out of five start structures. In the second method randomly selected combinations of up to six backbone torsion angles are simultaneously flipped into preset ranges by a short constraint EM followed by unconstraint EM and acceptance according to a Metropolis acceptance criterion. This combined stochastic/EM search was even more effective than the single torsion flip approach and selected low energy states for the two test cases in between two and four cases out of five start structures.

Conformational Analysis of Alanine Dipeptide from Dipolar Couplings in a Water-Based Liquid Crystal

The proton NMR spectra of unlabeled alanine dipeptide (Ac-L-Ala-NHMe) at 300 MHz and of alanine dipeptide with a single 13C label at 500 MHz are obtained in the lyotropic liquid-crystalline solvent cesium pentadecafluorooctanoate in water (CsPFO/H2O). Simulations of the spectra yield 9 and 13 dipolar couplings Dij, respectively, many with absolute sign determined. We fit the set of dipolar couplings by systematically varying the flexible dihedral angles φ and ψ while freezing local geometric details from the electronic structure calculations of Suhai and co-workers (Han, W. G.; Jalkanen, K. J.; Elstner, M.; Suhai, S. J. Phys. Chem. B 1998, 102, 2587). The orientation tensor is optimized at each combination of dihedral angles. Remarkably, a single conformer PII (φ ≈ −85°, ψ ≈ +160°) fits both sets of couplings within experimental error. The orientation tensor can be understood in terms of a simple rocking motion that dips the central methyl group into the fluorocarbon core of the CsPFO bicelle while alternately exposing both hydrogen-bonding pockets of PII to interfacial or bulk water. The search for a minority conformer such as αR (right-handed alpha helix, often favored by theory) using the larger data set was inconclusive. The data support localization of the peptide within the PII well rather than the broad sampling of φ in the range from −60° to −180° (a “β/PII minimum”) found by certain models. We suggest that the PII geometry is stable primarily because it maximizes the opportunity for peptide−water cooperative hydrogen bonding, whereas the αR geometry is stable primarily because of its large dipole moment. Our result corroborates recent work on short polypeptides suggesting that they preferentially sample configurations that fluctuate about PII-like structures, in contrast to the usual random-coil models.

New strategy for the conformational analysis of carbohydrates based on NOE and 13C NMR coupling constants. Application to the flexible polysaccharide of Streptococcus mitis J22.

For complex oligosaccharides, which are relatively rigid with modest excursions from a single minimum energy conformation, it is straightforward to build conformational models from NOE data. Other oligosaccharides are more flexible with transitions between distinct minima separated by substantial energy barriers. We show that modeling based on scalar coupling data is superior to NOE-based modeling for the latter case. Long range 13C-13C and 13C-1H coupling constants measured for the heptasaccharide repeating subunit of the cell wall polysaccharide from Streptococcus mitis J22 are correlated with individual glycosidic dihedral angles, effectively uncoupling the degrees of freedom of the oligosaccharide and allowing a search for combinations of dihedral angles which are energetically reasonable, i.e., with no bad van der Waals contacts, and which can be combined to satisfy all the measured J values. Allowed values of the individual angles can then be combined to search for overall oligosaccharide conformations which contribute to the ensemble. We show that while the polysaccharide from S. mitis J22 is flexible, requiring multiple conformations, most of the flexibility is localized to a few bonds and only a rather small number of conformations is required to reproduce the experimental NOE and scalar coupling data.