Antibody side chain conformations are position‐dependent

Side chain prediction is an integral component of computational antibody design and structure prediction. Current antibody modelling tools use backbone‐dependent rotamer libraries with conformations taken from general proteins. Here we present our antibody‐specific rotamer library, where rotamers are binned according to their immunogenetics (IMGT) position, rather than their local backbone geometry. We find that for some amino acid types at certain positions, only a restricted number of side chain conformations are ever observed. Using this information, we are able to reduce the breadth of the rotamer sampling space. Based on our rotamer library, we built a side chain predictor, position‐dependent antibody rotamer swapper (PEARS). On a blind test set of 95 antibody model structures, PEARS had the highest average χ₁ and accuracy (78.7% and 64.8%) compared to three leading backbone‐dependent side chain predictors. Our use of IMGT position, rather than backbone ϕ/ψ, meant that PEARS was more robust to errors in the backbone of the model structure. PEARS also achieved the lowest number of side chain–side chain clashes. PEARS is freely available as a web application at http://opig.stats.ox.ac.uk/webapps/pears .     

The Dynameomics rotamer library: amino acid side chain conformations and dynamics from comprehensive molecular dynamics simulations in water

We have recently completed systematic molecular dynamics simulations of 807 different proteins representing 95% of the known autonomous protein folds in an effort we refer to as Dynameomics. Here we focus on the analysis of side chain conformations and dynamics to create a dynamic rotamer library. Overall this library is derived from 31,000 occurrences of each of 86,217 different residues, or 2.7 × 10(9) rotamers. This dynamic library has 74% overlap of rotamer distributions with rotamer libraries derived from static high-resolution crystal structures. Seventy-five percent of the residues had an assignable primary conformation, and 68% of the residues had at least one significant alternate conformation. The average correlation time for switching between rotamers ranged from 22 ps for Met to over 8 ns for Cys; this time decreased 20-fold on the surface of the protein and modestly for dihedral angles further from the main chain. Side chain S(2) axis order parameters were calculated and they correlated well with those derived from NMR relaxation experiments (R = 0.9). Relationships relating the S(2) axis order parameters to rotamer occupancy were derived. Overall the Dynameomics rotamer library offers a comprehensive depiction of side chain rotamer preferences and dynamics in solution, and more realistic distributions for dynamic proteins in solution at ambient temperature than libraries derived from crystal structures, in particular charged surface residues are better represented. Details of the rotamer library are presented here and the library itself can be downloaded at http://www.dynameomics.org.     

Preferred side‐chain conformation of arginine residues in a triple‐helical structure

The crystal structure of the triple‐helical peptide (Pro‐Hyp‐Gly)₃‐Pro‐Arg‐Gly‐(Pro‐Hyp‐Gly)₄ (POG3‐PRG‐POG4) was determined at 1.45 Å resolution. POG3‐PRG‐POG4 was designed to permit investigation of the side‐chain conformation of the Arg residues in a triple‐helical structure. Because of the alternative structure of one of three Arg residues, four side‐chain conformations were observed in an asymmetric unit. Among them, three adopt a ttg^?t conformation and the other adopts a tg^?g^?t conformation. A statistical analysis of 80 Arg residues in various triple‐helical peptides showed that, unlike those in globular proteins, they preferentially adopt a tt conformation for χ₁ and χ₂, as observed in POG3‐PRG‐POG4. This conformation permits van der Waals contacts between the side‐chain atoms of Arg and the main‐chain atoms of the adjacent strand in the same molecule. Unlike many other host–guest peptides, in which there is a significant difference between the helical twists in the guest and the host peptides, POG3‐PRG‐POG4 shows a marked difference between the helical twists in the N‐terminal peptide and those in the C‐terminal peptide, separated near the Arg residue. This suggested that the unique side‐chain conformation of the Arg residue affects not only the conformation of the guest peptide, but also the conformation of the peptide away from the Arg residue. © 2014 Wiley Periodicals, Inc. Biopolymers 101: 1000–1009, 2014.     

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.     

Analysis of a data set of paired uncomplexed protein structures: new metrics for side-chain flexibility and model evaluation

We compiled and analyzed a data set of paired protein structures containing proteins for which multiple high-quality uncomplexed atomic structures were available in the Protein Data Bank. Side-chain flexibility was quantified, yielding a set of residue- and environment-specific confidence levels describing the range of motion around chi1 and chi2 angles. As expected, buried residues were inflexible, adopting similar conformations in different crystal structure analyses. Ile, Thr, Asn, Asp, and the large aromatics also showed limited flexibility when exposed on the protein surface, whereas exposed Ser, Lys, Arg, Met, Gln, and Glu residues were very flexible. This information is different from and complementary to the information available from rotamer surveys. The confidence levels are useful for assessing the significance of observed side-chain motion and estimating the extent of side-chain motion in protein structure prediction. We compare the performance of a simple 40 degrees threshold with these quantitative confidence levels in a critical evaluation of side-chain prediction with the program SCWRL.     

The impact of protein flexibility on protein-protein docking

Accounting for protein flexibility in protein-protein docking algorithms is challenging, and most algorithms therefore treat proteins as rigid bodies or permit side-chain motion only. While the consequences are obvious when there are large conformational changes upon binding, the situation is less clear for the modest conformational changes that occur upon formation of most protein-protein complexes. We have therefore studied the impact of local protein flexibility on protein-protein association by means of rigid body and torsion angle dynamics simulation. The binding of barnase and barstar was chosen as a model system for this study, because the complexation of these 2 proteins is well-characterized experimentally, and the conformational changes accompanying binding are modest. On the side-chain level, we show that the orientation of particular residues at the interface (so-called hotspot residues) have a crucial influence on the way contacts are established during docking from short protein separations of approximately 5 A. However, side-chain torsion angle dynamics simulations did not result in satisfactory docking of the proteins when using the unbound protein structures. This can be explained by our observations that, on the backbone level, even small (2 A) local loop deformations affect the dynamics of contact formation upon docking. Complementary shape-based docking calculations confirm this result, which indicates that both side-chain and backbone levels of flexibility influence short-range protein-protein association and should be treated simultaneously for atomic-detail computational docking of proteins.     

The role of side chain conformational flexibility in surface recognition by Tenebrio molitor antifreeze protein

Two-dimensional nuclear magnetic resonance spectroscopy was used to investigate the flexibility of the threonine side chains in the beta-helical Tenebrio molitor antifreeze protein (TmAFP) at low temperatures. From measurement of the (3)J(alphabeta) (1)H-(1)H scalar coupling constants, the chi(1) angles and preferred rotamer populations can be calculated. It was determined that the threonines on the ice-binding face of the protein adopt a preferred rotameric conformation at near freezing temperatures, whereas the threonines not on the ice-binding face sample many rotameric states. This suggests that TmAFP maintains a preformed ice-binding conformation in solution, wherein the rigid array of threonines that form the AFP-ice interface matches the ice crystal lattice. A key factor in binding to the ice surface and inhibition of ice crystal growth appears to be the close surface-to-surface complementarity between the AFP and crystalline ice, and the lack of an entropic penalty associated with freezing out motions in a flexible ligand.     

Incorporation of side chain flexibility into protein binding pockets using MTflex

The plasticity of active sites plays a significant role in drug recognition and binding, but the accurate incorporation of ‘receptor flexibility’ remains a significant computational challenge. Many approaches have been put forward to address receptor flexibility in docking studies by generating relevant ensembles on the energy surface. Herein, we describe the Movable Type with flexibility (MT_flex) method that generates ensembles on the more relevant free energy surface in a computationally tractable manner. This novel approach enumerates conformational states based on side chain flexibility and then estimates their relative free energies using the MT methodology. The resultant conformational states can then be used in subsequent docking and scoring exercises. In particular, we demonstrate that using the MT_flex ensembles improves MT’s ability to predict binding free energies over docking only to the crystal structure.     

Minimal ensembles of side chain conformers for modeling protein–protein interactions

The goal of this article is to reduce the complexity of the side chain search within docking problems. We apply six methods of generating side chain conformers to unbound protein structures and determine their ability of obtaining the bound conformation in small ensembles of conformers. Methods are evaluated in terms of the positions of side chain end groups. Results for 68 protein complexes yield two important observations. First, the end‐group positions change less than 1 Å on association for over 60% of interface side chains. Thus, the unbound protein structure carries substantial information about the side chains in the bound state, and the inclusion of the unbound conformation into the ensemble of conformers is very beneficial. Second, considering each surface side chain separately in its protein environment, small ensembles of low‐energy states include the bound conformation for a large fraction of side chains. In particular, the ensemble consisting of the unbound conformation and the two highest probability predicted conformers includes the bound conformer with an accuracy of 1 Å for 78% of interface side chains. As more than 60% of the interface side chains have only one conformer and many others only a few, these ensembles of low‐energy states substantially reduce the complexity of side chain search in docking problems. This approach was already used for finding pockets in protein–protein interfaces that can bind small molecules to potentially disrupt protein–protein interactions. Side‐chain search with the reduced search space will also be incorporated into protein docking algorithms. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Backbone conformations and side chain flexibility of two somatostatin mimics investigated by molecular dynamics simulations

Molecular dynamics simulations with two designed somatostatin mimics, SOM230 and SMS 201‐995, were performed in explicit water for a total aggregated time of 208 ns. Analysis of the runs with SOM230 revealed the presence of two clusters of conformations. Strikingly, the two sampled conformers correspond to the two main X‐ray structures in the asymmetric unit of SMS 201‐995. Structural comparison between the residues of SOM230 and SMS 201‐995 provides an explanation for the high binding affinity of SOM230 to four of five somatostatin receptors. Similarly, cluster analysis of the simulations with SMS 201‐995 shows that the backbone of the peptide interconverts between its two main crystallographic conformers. The conformations of SMS 201‐995 sampled in the two clusters violated two different sets of NOE distance constraints in agreement with a previous NMR study. Differences in side chain fluctuations between SOM230 and SMS 201‐995 observed in the simulations may contribute to the relatively higher binding affinity of SOM230 to most somatostatin receptors. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

On the chain flexibility of arabinoxylans and other beta-(1-->4) polysaccharides

The recent paper by Dervilly-Pinel and co-workers (Carbohydr. Res. 2001, 330, 365-372) presents a complete macromolecular characterisation of a series of de-esterified arabinoxylans extracted and fractionated from wheat flour. From their measurements, they were able to extract parameters related to chain flexibility, such as the Mark-Houwink exponent a and the chain persistence length. However, the estimate they obtain for the latter parameter is rather larger than would be expected, since the arabinoxylan backbone is beta-(1-->4) like cellulose and the galactomannans. By treating their data in an alternative, but well accepted manner, we are able to obtain a lower value of persistence length, which agrees well with estimates for this parameter for cellulose, in the literature, and our own recent measurements for a series of differently substituted galactomannans. These results suggest that the parameters obtained may be applicable to other beta-(1-->4) polysaccharides.     

SIDEpro: a novel machine learning approach for the fast and accurate prediction of side-chain conformations

Accurate protein side-chain conformation prediction is crucial for protein modeling and existing methods for the task are widely used; however, faster and more accurate methods are still required. Here we present a new machine learning approach to the problem where an energy function for each rotamer in a structure is computed additively over pairs of contacting atoms. A family of 156 neural networks indexed by amino acid and contacting atom types is used to compute these rotamer energies as a function of atomic contact distances. Although direct energy targets are not available for training, the neural networks can still be optimized by converting the energies to probabilities and optimizing these probabilities using Markov Chain Monte Carlo methods. The resulting predictor SIDEpro makes predictions by initially setting the rotamer probabilities for each residue from a backbone-dependent rotamer library, then iteratively updating these probabilities using the trained neural networks. After convergences of the probabilities, the side-chains are set to the highest probability rotamer. Finally, a post processing clash reduction step is applied to the models. SIDEpro represents a significant improvement in speed and a modest, but statistically significant, improvement in accuracy when compared with the state-of-the-art for rapid side-chain prediction method SCWRL4 on the following datasets: (1) 379 protein test set of SCWRL4; (2) 94 proteins from CASP9; (3) a set of seven large protein-only complexes; and (4) a ribosome with and without the RNA. Using the SCWRL4 test set, SIDEpro's accuracy (χ(1) 86.14%, χ(1+2) 74.15%) is slightly better than SCWRL4-FRM (χ(1) 85.43%, χ(1+2) 73.47%) and it is 7.0 times faster. On the same test set SIDEpro is clearly more accurate than SCWRL4-rigid rotamer model (RRM) (χ(1) 84.15%, χ(1+2) 71.24%) and 2.4 times faster. Evaluation on the additional test sets yield similar accuracy results with SIDEpro being slightly more accurate than SCWRL4-flexible rotamer model (FRM) and clearly more accurate than SCWRL4-RRM; however, the gap in CPU time is much more significant when the methods are applied to large protein complexes. SIDEpro is part of the SCRATCH suite of predictors and available from: http://scratch.proteomics.ics.uci.edu/.     

Fast and accurate side-chain topology and energy refinement (FASTER) as a new method for protein structure optimization

We have developed an original method for global optimization of protein side-chain conformations, called the Fast and Accurate Side-Chain Topology and Energy Refinement (FASTER) method. The method operates by systematically overcoming local minima of increasing order. Comparison of the FASTER results with those of the dead-end elimination (DEE) algorithm showed that both methods produce nearly identical results, but the FASTER algorithm is 100-1000 times faster than the DEE method and scales in a stable and favorable way as a function of protein size. We also show that low-order local minima may be almost as accurate as the global minimum when evaluated against experimentally determined structures. In addition, the new algorithm provides significant information about the conformational flexibility of individual side-chains. We observed that strictly rigid side-chains are concentrated mainly in the core of the protein, whereas highly flexible side-chains are found almost exclusively among solvent-oriented residues.     

The molecular structure and solution conformation of an acidic heteropolysaccharide from Auricularia auricula-judae

A water soluble acidic heteropolysaccharide named WAF was isolated from Auricularia auricula‐judae by extracting with 0.9% NaCl solution. By using gas chromatography, gas chromatography‐mass spectrometry, and NMR, its chemical structure was determined to be composed of a backbone of α‐(1→3)‐linked D ‐mannopyranose residues with pendant side groups of β‐D ‐xylose, β‐D ‐glucose, or β‐D ‐glucuronic acid at position O6 or O2. Six fractions prepared from WAF with a weight‐average molecular mass (M_w) between 5.9 × 10⁴ and 64.7 × 10⁴ g/mol were characterized with laser light scattering and viscometry in 0.1M NaCl at 25°C. The dependence of intrinsic viscosity ([η]) and radius of gyration (R_g) on M_w for this polysaccharide were found to be [η] = 1.79 × 10^?3M_w^0.96 cm³ g^?1 and R_g = 6.99 × 10^?2 M_w^0.54 nm. The molar mass per unit contour length (M_L) and the persistence length (L_p) were estimated to be 1124 nm^?1 and 11 nm, respectively. The WAF exhibited a semirigid character typical of linear polysaccharides. Molecular modeling was then used to predict the ordered and disordered states of WAF; the simulated M_L and L_p were however much smaller than the experimental values. Taken altogether, the results suggested that WAF formed a duplex in solution. © 2010 Wiley Periodicals, Inc. Biopolymers 95: 217–227, 2011.     

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 χ₁ and χ₂ 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 χ₂, the Leu side‐chain bond angles κ₁ and κ₂ must increase, whereas to transition in χ₁, the Ile bond angles λ₁ and λ₂ 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. Proteins 2015; 83:1488–1499. © 2015 Wiley Periodicals, Inc.     

Tryptophan side chain conformers monitored by NMR and time-resolved fluorescence spectroscopies

We have inserted a tryptophan (F77W) in the core of the regulatory domain of cardiac troponin C (cNTnC), and previously determined the structure of this mutant with and without the cosolvent trifluoroethanol (TFE). Interestingly, the orientations of the indole side chain of the Trp are in opposite directions in the two structures (Julien et al., Protein Sci 2009; 18:1165-1174). Fluorescence decay experiments for single Trp-containing proteins often show several lifetimes, which have been interpreted as reflecting conformational heterogeneity of the Trp side chain resulting from different rotamers. To test this interpretation, we monitored the effect of TFE on wild type, F77W and F77W-V82A calcium-saturated cNTnC using 2D (13)C-HSQC NMR and time-correlated single photon counting fluorescence spectroscopies. The time dependence of the Trp fluorescence decay was fit with three lifetimes. Addition of TFE caused a gradual, but limited decrease of the lifetimes due to dynamic quenching. For F77W cNTnC, the amplitude fractions of the lifetimes also changed upon addition of TFE-the long lifetime increased from 13 to 29%, while the middle lifetime decreased from 63 to 50% and the short lifetime remained relatively unchanged. For F77W-V82A cNTnC, comparable NMR changes are observed, confirming the switch in rotamer conformation, but only much smaller changes in fluorescence decay parameters were detected. These data indicate that the balance between the rotamer states can be changed without changing the lifetime amplitude fractions appreciably, and suggest that the environment(s) of the indole ring, responsible for the different lifetimes, can result from factors other than the intrinsic rotamer state of the tryptophan.     

Assessment of protein side‐chain conformation prediction methods in different residue environments

Computational prediction of side‐chain conformation is an important component of protein structure prediction. Accurate side‐chain prediction is crucial for practical applications of protein structure models that need atomic‐detailed resolution such as protein and ligand design. We evaluated the accuracy of eight side‐chain prediction methods in reproducing the side‐chain conformations of experimentally solved structures deposited to the Protein Data Bank. Prediction accuracy was evaluated for a total of four different structural environments (buried, surface, interface, and membrane‐spanning) in three different protein types (monomeric, multimeric, and membrane). Overall, the highest accuracy was observed for buried residues in monomeric and multimeric proteins. Notably, side‐chains at protein interfaces and membrane‐spanning regions were better predicted than surface residues even though the methods did not all use multimeric and membrane proteins for training. Thus, we conclude that the current methods are as practically useful for modeling protein docking interfaces and membrane‐spanning regions as for modeling monomers. Proteins 2014; 82:1971–1984. © 2014 Wiley Periodicals, Inc.     

On the flexibility of poly(γ-benzyl L-glutamate) in helicogenic solvents

The molecular-weight dependence of the rms radius of gyration of poly(γ-benzyl L -glutamate) (PBLG) in helicogenic solvents shows negative and positive deviations from expectations for an intact and rigid α-helix in the higher and lower molecular-weight ranges, respectively. In order to study the reason for both deviations, we compare the extant experimental data of with those computed for wormlike chain, freely jointed rod, and a rigid rod having random-coil portions at both ends. The computation for the freely jointed rod and the rigid rod having frayed ends is carried out by a simulation method of Muroga. From the Zimm and Bragg theory and the above comparisons, it is concluded that both deviations can be self-consistently explained if PBLG in helicogenic solvents has an essentially intact α-helical structure with some flexibility arising from random fluctuations in hydrogen bond length. This flexibility explains the negative deviations in the high molecular weight region. The positive deviations in the low molecular weight region result from the tendency of helices to unwind at the ends. © 1998 John Wiley & Sons, Inc. Biopoly 45: 281–288, 1998