SDRL: a sequence-dependent protein side-chain rotamer library

Exploiting Sequence-Dependent Rotamer Information in Global Optimization of Proteins

Rotamers, namely amino acid side chain conformations common to many different peptides, can be compiled into libraries. These rotamer libraries are used in protein modeling, where the limited conformational space occupied by amino acid side chains is exploited. Here, we construct a sequence-dependent rotamer library from simulations of all possible tripeptides, which provides rotameric states dependent on adjacent amino acids. We observe significant sensitivity of rotamer populations to sequence and find that the library is successful in locating side chain conformations present in crystal structures. The library is designed for applications with basin-hopping global optimization, where we use it to propose moves in conformational space. The addition of rotamer moves significantly increases the efficiency of protein structure prediction within this framework, and we determine parameters to optimize efficiency.

Cost-Effective Potential for Accurate Polarizable Embedding Calculations in Protein Environments

The fragment-based polarizable embedding (PE) model combined with an appropriate electronic structure method constitutes a highly efficient and accurate multiscale approach for computing spectroscopic properties of a central moiety including effects from its molecular environment through an embedding potential. There is, however, a comparatively high computational overhead associated with the computation of the embedding potential, which is derived from first-principles calculations on individual fragments of the environment. To reduce the computational cost associated with the calculation of embedding potential parameters, we developed a set of amino acid-specific transferable parameters tailored for large-scale PE-based calculations that include proteins. The amino acid-based parameters are obtained by simultaneously fitting to a set of reference electric potentials based on structures derived from a backbone-dependent rotamer library. The developed cost-effective polarizable protein potential (CP³) consists of atom-centered charges and isotropic dipole-dipole polarizabilities of the standard amino acids. In terms of reproduction of electric potentials, the CP³ is shown to perform consistently and with acceptable accuracy across both small tripeptide test systems and larger proteins. We show, through applications on realistic protein systems, that acceptable accuracy can be obtained by using a pure CP³ representation of the protein environment, thus altogether omitting the cost associated with the calculation of embedding potential parameters. High accuracy comparable to that of the full fragment-based approach can be achieved through a mixed description where the CP³ is used only to describe amino acids beyond a threshold distance from the central quantum part.

Side-Chain Conformational Preferences Govern Protein-Protein Interactions

Protein secondary structures serve as geometrically constrained scaffolds for the display of key interacting residues at protein interfaces. Given the critical role of secondary structures in protein folding and the dependence of folding propensities on backbone dihedrals, secondary structure is expected to influence the identity of residues that are important for complex formation. Counter to this expectation, we find that a narrow set of residues dominates the binding energy in protein-protein complexes independent of backbone conformation. This finding suggests that the binding epitope may instead be substantially influenced by the side-chain conformations adopted. We analyzed side-chain conformational preferences in residues that contribute significantly to binding. This analysis suggests that preferred rotamers contribute directly to specificity in protein complex formation and provides guidelines for peptidomimetic inhibitor design.

Effect of introducing aib in a designed helical inhibitor of hdm2-p53 interaction: A molecular dynamics study

Although p53 is an intrinsically disordered protein, upon binding to Hdm2, a short stretch (residues 19-25) comprising the binding epitope assumes a helical backbone. Because the allowed conformational space of α-aminoisobutyric acid (Aib) is restricted to only the helical basin, Aib-containing helical mimics of p53 (binding epitope) are expected to inhibit interaction between p53 and Hdm2 with a much stronger affinity than the wild type p53 peptide (binding epitope), due to the entropic advantage associated with Aib. However, the IC50 values for the disruption of p53-Hdm2 interaction by Aib-p53 peptides and wild type p53 peptide were found to be comparable (J. Peptide Res. 2002, 60:88-94). To understand why incorporation of Aib didn't substantially increase Hdm2 affinity of Aib-p53 peptides, a series of molecular dynamics simulations were performed. It was found that despite stabilizing a helical backbone in the unbound state, the Aib residues in Aib-p53 peptide arrested two functionally important side-chains (F19 and W23) in non-productive conformations, resulting in relative side-chain orientations of the binding triad F19-W23-L26 incompatible with the bound conformation. Therefore, although a Aib-induced pre-formed helical peptide backbone in the unbound state is expected to favor binding, the locked side-chain orientations of the binding triad in non-productive modes would disfavor binding. This study shows that when using Aib to design functionally important helical peptides, care must be taken to consider potential interactions between side-chains of neighboring residues and Aib in the unbound state.