Rotamer Library Research Articles

CIFDock: A novel CHARMM-based flexible receptor-flexible ligand docking protocol.

Docking studies play a critical role in the current workflow of drug discovery. However, limitations may often arise through factors including inadequate ligand sampling, a lack of protein flexibility, scoring function inadequacies (e.g., due to metals, co-factors, etc.), and difficulty in retaining explicit water molecules. Herein, we present a novel CHARMM-based induced fit docking (CIFDock) workflow that can circumvent these limitations by employing all-atom force fields coupled to enhanced sampling molecular dynamics procedures. Self-guided Langevin dynamics simulations are used to effectively sample relevant ligand conformations, side chain orientations, crystal water positions, and active site residue motion. Protein flexibility is further enhanced by dynamic sampling of side chain orientations using an expandable rotamer library. Steps in the procedure consisting of fixing individual components (e.g., the ligand) while sampling the other components (e.g., the residues in the active site of the protein) allow for the complex to adapt to conformational changes. Ultimately, all components of the complex-the protein, ligand, and waters-are sampled simultaneously and unrestrained with SGLD to capture any induced fit effects. This modular flexible docking procedure is automated using CHARMM scripting, interfaced with SLURM array processing, and parallelized to use the desired number of processors. We validated the CIFDock procedure by performing cross-docking studies using a data set comprised of 21 pharmaceutically relevant proteins. Five variants of the CHARMM-based SWISSDOCK scoring functions were created to quantify the results of the final generated poses. Results obtained were comparable to, or in some cases improved upon, commercial docking program data.

Resolving distance variations by single-molecule FRET and EPR spectroscopy using rotamer libraries

Förster resonance energy transfer (FRET) and electron paramagnetic resonance (EPR) spectroscopy are complementary techniques for quantifying distances in the nanometer range. Both approaches are commonly employed for probing the conformations and conformational changes of biological macromolecules based on site-directed fluorescent or paramagnetic labeling. FRET can be applied in solution at ambient temperature and thus provides direct access to dynamics, especially if used at the single-molecule level, whereas EPR requires immobilization or work at cryogenic temperatures but provides data that can be more reliably used to extract distance distributions. However, a combined analysis of the complementary data from the two techniques has been complicated by the lack of a common modeling framework. Here, we demonstrate a systematic analysis approach based on rotamer libraries for both FRET and EPR labels to predict distance distributions between two labels from a structural model. Dynamics of the fluorophores within these distance distributions are taken into account by diffusional averaging, which improves the agreement with experiment. Benchmarking this methodology with a series of surface-exposed pairs of sites in a structured protein domain reveals that the lowest resolved distance differences can be as small as ∼0.25 nm for both techniques, with quantitative agreement between experimental and simulated transfer efficiencies within a range of ±0.045. Rotamer library analysis thus establishes a coherent way of treating experimental data from EPR and FRET and provides a basis for integrative structural modeling, including studies of conformational distributions and dynamics of biological macromolecules using both techniques.

DEER-PREdict: Software for efficient calculation of spin-labeling EPR and NMR data from conformational ensembles.

Owing to their plasticity, intrinsically disordered and multidomain proteins require descriptions based on multiple conformations, thus calling for techniques and analysis tools that are capable of dealing with conformational ensembles rather than a single protein structure. Here, we introduce DEER-PREdict, a software program to predict Double Electron-Electron Resonance distance distributions as well as Paramagnetic Relaxation Enhancement rates from ensembles of protein conformations. DEER-PREdict uses an established rotamer library approach to describe the paramagnetic probes which are bound covalently to the protein.DEER-PREdict has been designed to operate efficiently on large conformational ensembles, such as those generated by molecular dynamics simulation, to facilitate the validation or refinement of molecular models as well as the interpretation of experimental data. The performance and accuracy of the software is demonstrated with experimentally characterized protein systems: HIV-1 protease, T4 Lysozyme and Acyl-CoA-binding protein. DEER-PREdict is open source (GPLv3) and available at github.com/KULL-Centre/DEERpredict and as a Python PyPI package pypi.org/project/DEERPREdict.

BCL::Conf: Improved Open-Source Knowledge-Based Conformation Sampling Using the Crystallography Open Database.

We previously described BCL::Conf, a knowledge-based conformation sampling algorithm utilizing a small molecule fragment rotamer library derived from the Cambridge Structural Database (CSD, license required), as a component of the BioChemical Library (BCL) cheminformatics toolkit. This paper describes substantial improvements made to the BCL::Conf algorithm and a transition to a rotamer library derived from molecules in the Crystallography Open Database (COD, no license required). We demonstrate the performance of the new BCL::Conf on native conformer recovery in the Platinum dataset of high-quality protein-ligand complexes. This set of 2859 structures has previously been used to assess the performance of over a dozen conformer generation algorithms, including the Conformator, Balloon, RDKit DG, ETKDG, Confab, Frog2, MultiConf-DOCK, CSD conformer generator, ConfGenX-OPSL3 force field, Omega, excalc, iCon, and MOE. These benchmarks suggest that the CSD conformer generator is at the state of the art of reported conformer generators. Our results indicate that the improved BCL::Conf significantly outperforms the CSD conformer generation algorithm at binding conformer recovery across a range of ensemble sizes and with similarly fast rates of conformer generation. BCL::Conf is now distributed with the COD-derived rotamer library and is free for academic use. The BCL can be downloaded at http://meilerlab.org/bclcommons for Windows, Linux, or Apple operating systems. BCL::Conf can now also be accessed via webserver at http://meilerlab.org/bclconf.

OPUS-Rota3: Improving Protein Side-Chain Modeling by Deep Neural Networks and Ensemble Methods.

Side-chain modeling is critical for protein structure prediction since the uniqueness of the protein structure is largely determined by its side-chain packing conformation. In this paper, differing from most approaches that rely on rotamer library sampling, we first propose a novel side-chain rotamer prediction method based on deep neural networks, named OPUS-RotaNN. Then, on the basis of our previous work OPUS-Rota2, we propose an open-source side-chain modeling framework, OPUS-Rota3, which integrates the results of different methods into its rotamer library as the sampling candidates. By including OPUS-RotaNN into OPUS-Rota3, we conduct our experiments on three native backbone test sets and one non-native backbone test set. On the native backbone test set, CAMEO-Hard61 for example, OPUS-Rota3 successfully predicts 51.14% of all side-chain dihedral angles with a tolerance criterion of 20° and outperforms OSCAR-star (50.87%), SCWRL4 (50.40%), and FASPR (49.85%). On the non-native backbone test set DB379-ITASSER, the accuracy of OPUS-Rota3 is 52.49%, better than OSCAR-star (48.95%), FASPR (48.69%), and SCWRL4 (48.29%). All the source codes including the training codes and the data we used are available at https://github.com/thuxugang/opus_rota3.

GWOVina: A grey wolf optimization approach to rigid and flexible receptor docking.

Protein–ligand docking programs are indispensable tools for predicting the binding pose of a ligand to the receptor protein. In this paper, we introduce an efficient flexible docking method, gwovina, which is a variant of the Vina implementation using the grey wolf optimizer (GWO) and random walk for the global search, and the Dunbrack rotamer library for side‐chain sampling. The new method was validated for rigid and flexible‐receptor docking using four independent datasets. In rigid docking, gwovina showed comparable docking performance to Vina in terms of ligand pose RMSD, success rate, and affinity prediction. In flexible‐receptor docking, gwovina has improved success rate compared to Vina and AutoDockFR. It ran 2 to 7 times faster than Vina and 40 to 100 times faster than AutoDockFR. Therefore, gwovina can play a role in solving the complex flexible‐receptor docking cases and is suitable for virtual screening of compound libraries. gwovina is freely available at https://cbbio.cis.um.edu.mo/software/gwovina for testing.

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 (CP3) consists of atom-centered charges and isotropic dipole-dipole polarizabilities of the standard amino acids. In terms of reproduction of electric potentials, the CP3 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 CP3 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 CP3 is used only to describe amino acids beyond a threshold distance from the central quantum part.

Toward the Accuracy and Speed of Protein Side-Chain Packing: ASystematic Study on Rotamer Libraries.

Protein rotamers refer to the conformational isomers taken by the side-chains of amino acids to accommodate specific structural folding environments. Since accurate modeling of atomic interactions is difficult, rotamer information collected from experimentally solved protein structures is often used to guide side-chain packing in protein folding and sequence design studies. Many rotamer libraries have been built in the literature but there is little quantitative guidance on which libraries should be chosen for different structural modeling studies. Here, we performed a comparative study of six widely used rotamer libraries and systematically examined their suitability for protein folding and sequence design in four aspects: (1) side-chain match accuracy, (2) side-chain conformation prediction, (3) de novo protein sequence design, and (4) computational time cost. We demonstrated that, compared to the backbone-dependent rotamer libraries (BBDRLs), the backbone-independent rotamer libraries (BBIRLs) generated conformations that more closely matched the native conformations due to the larger number of rotamers in the local rotamer search spaces. However, more practically, using an optimized physical energy function incorporated into a simulated annealing Monte Carlo searching scheme, we showed that utilization of the BBDRLs could result in higher accuracies in side-chain prediction and higher sequence recapitulation rates in protein design experiments. Detailed data analyses showed that the major advantage of BBDRLs lies in the energy term derived from the rotamer probabilities that are associated with the individual backbone torsion angle subspaces. This term is important for distinguishing between amino acid identities as well as the rotamer conformations of an amino acid. Meanwhile, the backbone torsion angle subspace-specific rotamer search drastically speeds up the searching time, despite the significantly larger number of total rotamers in the BBDRLs. These results should provide important guidance for the development and selection of rotamer libraries for practical protein design and structure prediction studies.

Isoindoline-Based Nitroxides as Bioresistant Spin Labels for Protein Labeling through Cysteines and Alkyne-Bearing Noncanonical Amino Acids.

Electron paramagnetic resonance (EPR) spectroscopy in combination with site‐directed spin labeling (SDSL) is a powerful tool in protein structural research. Nitroxides are highly suitable spin labeling reagents, but suffer from limited stability, particularly in the cellular environment. Herein we present the synthesis of a maleimide‐ and an azide‐modified tetraethyl‐shielded isoindoline‐based nitroxide (M‐ and Az‐TEIO) for labeling of cysteines or the noncanonical amino acid para‐ethynyl‐l‐phenylalanine (pENF). We demonstrate the high stability of TEIO site‐specifically attached to the protein thioredoxin (TRX) against reduction in prokaryotic and eukaryotic environments, and conduct double electron–electron resonance (DEER) measurements. We further generate a rotamer library for the new residue pENF‐Az‐TEIO that affords a distance distribution that is in agreement with the measured distribution.

Gem-Diethyl Pyrroline Nitroxide Spin Labels: Synthesis, EPR Characterization, Rotamer Libraries and Biocompatibility.

Invited for this month's cover picture is the group of Professor Enrica Bordignon at the Ruhr University Bochum. The cover picture shows an artistic view of E. coli cells and two spin-labeled recombinantly produced proteins, which can be inserted into the cells for EPR studies. The primary sequence of the proteins is schematically shown with the one-letter amino acid code, and cysteine residues are functionalized with the two new gem diethyl nitroxide spin labels designed to better sustain the reducing cellular environment. Read the full text of their Full Paper at 10.1002/open.201900119.

Gem-Diethyl Pyrroline Nitroxide Spin Labels: Synthesis, EPR Characterization, Rotamer Libraries and Biocompatibility.

The availability of bioresistant spin labels is crucial for the optimization of site‐directed spin labeling protocols for EPR structural studies of biomolecules in a cellular context. As labeling can affect proteins’ fold and/or function, having the possibility to choose between different spin labels will increase the probability to produce spin‐labeled functional proteins. Here, we report the synthesis and characterization of iodoacetamide‐ and maleimide‐functionalized spin labels based on the gem‐diethyl pyrroline structure. The two nitroxide labels are compared to conventional gem‐dimethyl analogs by site‐directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy, using two water soluble proteins: T4 lysozyme and Bid. To foster their use for structural studies, we also present rotamer libraries for these labels, compatible with the MMM software. Finally, we investigate the “true” biocompatibility of the gem‐diethyl probes comparing the resistance towards chemical reduction of the NO group in ascorbate solutions and E. coli cytosol at different spin concentrations.

Variable Neighborhood Search with Cost Function Networks To Solve Large Computational Protein Design Problems.

Computational protein design (CPD) aims to predict amino acid sequences that fold to specific structures and perform desired functions. CPD depends on a rotamer library, an energy function, and an algorithm to search the sequence/conformation space. Variable neighborhood search (VNS) with cost function networks is a powerful framework that can provide tight upper bounds on the global minimum energy. We propose a new CPD heuristic based on VNS in which a subset of the solution space (a "neighborhood") is explored, whose size is gradually increased with a dedicated probabilistic heuristic. The algorithm was tested on 99 protein designs with fixed backbones involving nine proteins from the SH2, SH3, and PDZ families. The number of mutating positions was 20, 30, or all of the amino acids, while the rest of the protein explored side-chain rotamers. VNS was more successful than Monte Carlo (MC), replica-exchange MC, and a heuristic steepest-descent energy minimization, providing solutions with equal or lower best energies in most cases. For complete protein redesign, it gave solutions that were 2.5 to 11.2 kcal/mol lower in energy than those obtained with the other approaches. VNS is implemented in the toulbar2 software. It could be very helpful for large and/or complex design problems.

Precision DEER Distances from Spin-Label Ensemble Refinement.

Double electron-electron resonance (DEER) experiments probe nanometer-scale distances in spin-labeled proteins and nucleic acids. Rotamer libraries of the covalently attached spin-labels help reduce position uncertainties. Here we show that rotamer reweighting is essential for precision distance measurements, making it possible to resolve Ångstrom-scale domain motions. We analyze extensive DEER measurements on the three N-terminal polypeptide transport-associated (POTRA) domains of the outer membrane protein Omp85. Using the "Bayesian inference of ensembles" maximum-entropy method, we extract rotamer weights from the DEER measurements. Small weight changes suffice to eliminate otherwise significant discrepancies between experiments and model and unmask 1-3 Å domain motions relative to the crystal structure. Rotamer-weight refinement is a simple yet powerful tool for precision distance measurements that should be broadly applicable to label-based measurements including DEER, paramagnetic relaxation enhancement, and fluorescence resonance energy transfer (FRET).

Molecular dynamics-derived rotamer libraries for d-amino acids within homochiral and heterochiral polypeptides.

Computational resources have contributed to the design and engineering of novel proteins by integrating genomic, structural and dynamic aspects of proteins. Non-canonical amino acids, such as d-amino acids, expand the available sequence space for designing and engineering proteins; however, the rotamer libraries for d-amino acids are usually constructed as the mirror images of l-amino acid rotamer libraries, an assumption that has not been tested. To this end, we have performed molecular dynamics (MD) simulations of model host-guest peptide systems containing d-amino acids. Our simulations systematically address the applicability of the mirror image convention as well as the effects of neighboring residue chirality. Rotamer libraries derived from these systems provide realistic rotamer distributions suitable for use in both rational and computational design workflows. Our simulations also address the impact of chirality on the intrinsic conformational preferences of amino acids, providing fundamental insights into the relationship between chirality and biomolecular dynamics. While d-amino acids are rare in naturally occurring proteins, they are used in designed proteins to stabilize a desired conformation, increase bioavailability or confer favorable biochemical and physical attributes. Here, we present d-amino acid rotamer libraries derived from MD simulations of alanine-based host-guest pentapeptides and show how certain residues can deviate from mirror image symmetry. Our simulations directly model d-amino acids as guest residues within the chiral l-Ala and d-Ala pentapeptide series to explicitly incorporate any contributions resulting from the chiralities of neighboring residues.

Cost function network-based design of protein-protein interactions: predicting changes in binding affinity.

Accurate and economic methods to predict change in protein binding free energy upon mutation are imperative to accelerate the design of proteins for a wide range of applications. Free energy is defined by enthalpic and entropic contributions. Following the recent progresses of Artificial Intelligence-based algorithms for guaranteed NP-hard energy optimization and partition function computation, it becomes possible to quickly compute minimum energy conformations and to reliably estimate the entropic contribution of side-chains in the change of free energy of large protein interfaces. Using guaranteed Cost Function Network algorithms, Rosetta energy functions and Dunbrack's rotamer library, we developed and assessed EasyE and JayZ, two methods for binding affinity estimation that ignore or include conformational entropic contributions on a large benchmark of binding affinity experimental measures. If both approaches outperform most established tools, we observe that side-chain conformational entropy brings little or no improvement on most systems but becomes crucial in some rare cases. as open-source Python/C++ code at sourcesup.renater.fr/projects/easy-jayz. Supplementary data are available at Bioinformatics online.

Scoring of Side-Chain Packings: An Analysis of Weight Factors and Molecular Dynamics Structures.

The protein side-chain packing problem (PSCPP) is a central task in computational protein design. The problem is usually modeled as a combinatorial optimization problem, which consists of searching for a set of rotamers, from a given rotamer library, that minimizes a scoring function (SF). The SF is a weighted sum of terms, that can be decomposed in physics-based and knowledge-based terms. Although there are many methods to obtain approximate solutions for this problem, all of them have similar performances and there has not been a significant improvement in recent years. Studies on protein structure prediction and protein design revealed the limitations of current SFs to achieve further improvements for these two problems. In the same line, a recent work reported a similar result for the PSCPP. In this work, we ask whether or not this negative result regarding further improvements in performance is due to (i) an incorrect weighting of the SFs terms or (ii) the constrained conformation resulting from the protein crystallization process. To analyze these questions, we (i) model the PSCPP as a bi-objective combinatorial optimization problem, optimizing, at the same time, the two most important terms of two SFs of state-of-the-art algorithms and (ii) performed a preprocessing relaxation of the crystal structure through molecular dynamics to simulate the protein in the solvent and evaluated the performance of these two state-of-the-art SFs under these conditions. Our results indicate that (i) no matter what combination of weight factors we use the current SFs will not lead to better performances and (ii) the evaluated SFs will not be able to improve performance on relaxed structures. Furthermore, the experiments revealed that the SFs and the methods are biased toward crystallized structures.

Relative Contribution of the Refinement Steps to the Protein-Protein Docking Success Rate

Structural modeling of protein-protein complexes is a challenging task as proteins may undergo significant conformational changes upon binding. Flexible docking of proteins typically consists of low-resolution rigid-body global docking search, followed by the structural refinement that takes into account conformational flexibility. The refinement typically involves consecutive steps associated with different weakly coupled sets of degrees of freedom (DOF). Previously, we developed a docking refinement pipeline consisting of three steps: (i) rigid-body adjustment of the protein position, which minimizes a semi-empirical atom-atom contact potential (rigid-body rotational/translational DOF), (ii) conformer selection procedure applied to the backbone loops at the interface, to maximize the Ramachandran probability (backbone loops torsional DOF), and (iii) repacking of the interface side chains using a library of rotamers, to minimize the steric clashes (side-chain torsional DOF). In this study, we investigate relative contributions of these steps to the refinement success, in terms of the fraction of the native contacts, ligand RMSD, and interface RMSD, for bound and unbound proteins in the Dockground X-ray docking benchmark set 4.0 (http://dockground.compbio.ku.edu).

Identifying the Optimal Level of Regularization in Deer Data Analysis

Double Electron-Electron Resonance (DEER) spectroscopy is an Electron Paramagnetic Resonance (EPR) experiment that accurately reports the distribution of distances, P(r), between two or more unpaired electrons over the range of about 10 to 160 ångströms. DEER is becoming increasingly important as a tool for structural biology due to its ability to elucidate protein conformational landscapes, as well as its amenability to complex samples such as membrane-bound proteins and large protein complexes. The spin-spin distance distribution, P(r), is most often inferred from the time-domain experimental data with a regularized least-squares method known as Tikhonov regularization. Although this method typically provides reliable estimates of P(r), it requires the input of a regularization parameter, α, which scales the magnitude of a penalty term with respect to the fitting error. This regularization parameter significantly affects the estimated P(r), so selecting its optimal value is crucial for accurate analysis of DEER data. The most popular procedure for selecting α in DEER data analysis is a graphical method called the L-curve criterion, which has empirical justification. However, to date, a systematic exploration of the many extant regularization parameter selection methods has not been carried out in the context of DEER data analysis. Here, we evaluate 16 methods, including the L-curve criterion, against a large, physically-derived test set of simulated DEER data. The test set of P(r)s was generated using a side-chain rotamer library in conjunction with a crystal structure (2LZM). This test set was then used to generate a large set of experimentally realistic simulated time-domain data.

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 χ1 and χ1+2 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.

Assessing the capability of in silico mutation protocols for predicting the finite temperature conformation of amino acids.

Mutation protocols are a key tool in computational biophysics for modelling unknown side chain conformations. In particular, these protocols are used to generate the starting structures for molecular dynamics simulations. The accuracy of the initial side chain and backbone placement is crucial to obtain a stable and quickly converging simulation. In this work, we assessed the performance of several mutation protocols in predicting the most probable conformer observed in finite temperature molecular dynamics simulations for a set of protein-peptide crystals differing only by single-point mutations in the peptide sequence. Our results show that several programs which predict well the crystal conformations fail to predict the most probable finite temperature configuration. Methods relying on backbone-dependent rotamer libraries have, in general, a better performance, but even the best protocol fails in predicting approximately 30% of the mutations.