MOTIVATION: Predicting the structure of antibody-antigen complexes is a challenging task with significant implications for the design of better antibody therapeutics. However, the levels of success have remained dauntingly low, particularly when high standards for model quality are required, a necessity for efficient antibody design. Artificial intelligence (AI) has significantly impacted the landscape of structure prediction for antibodies, both alone and in complex with their antigens. METHODS: We utilized AI-guided antibody modeling tools to generate ensembles displaying diversity in the complementarity-determining region (CDR) and integrated those into our previously published AlphaFold2-rescored docking pipeline, a strategy called AI-augmented physics-based docking. In this study, we also compare docking performance with AlphaFold and Boltz-1, the new state-of-the-art. We distinguish between two types of success tailored to specific downstream applications: (i) criteria sufficient for epitope mapping, where gross quality is adequate and can complement experimental techniques, and (ii) criteria for producing higher-quality models suitable for engineering purposes. RESULTS: We highlight that the quality of the ensemble is crucial for docking performance, that including too many models can be detrimental, and that prioritization of models is essential for achieving good performance. In a scenario analogous to docking using a crystallized antigen, our results robustly demonstrate the advantages of AI-augmented docking over AlphaFold2, further accentuated when higher standards in quality are imposed. Docking also shows improvements over Boltz-1, but those are less pronounced. Docking performance is still noticeably lower than AlphaFold3 in both epitope mapping and antibody design use cases. We observe a strong dependence on CDR-H3 loop length for physics-based tools on their ability to successfully predict. This helps define an applicability range where physics-based docking can be competitive to the newer generation of AI tools. AVAILABILITY AND IMPLEMENTATION: The AF2 rescoring scripts are available at github.com/gaudreaultfnrc/AF2-Rescoring.

UNLABELLED: Ligand protein docking simulations play a fundamental role in understanding molecular recognition. Herein we introduce the NRGsuite, a PyMOL plugin that permits the detection of surface cavities in proteins, their refinements, calculation of volume and use, individually or jointly, as target binding-sites for docking simulations with FlexAID. The NRGsuite offers the users control over a large number of important parameters in docking simulations including the assignment of flexible side-chains and definition of geometric constraints. Furthermore, the NRGsuite permits the visualization of the docking simulation in real time. The NRGsuite give access to powerful docking simulations that can be used in structure-guided drug design as well as an educational tool. The NRGsuite is implemented in Python and C/C++ with an easy to use package installer. The NRGsuite is available for Windows, Linux and MacOS. AVAILABILITY AND IMPLEMENTATION: http://bcb.med.usherbrooke.ca/flexaid. CONTACT: [email protected]. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MOTIVATION: An increasing amount of evidence from experimental and computational analysis suggests that rare codon clusters are functionally important for protein activity. Most of the studies on rare codon clusters were performed on a limited number of proteins or protein families. In the present study, we present the Sherlocc program and how it can be used for large scale protein family analysis of evolutionarily conserved rare codon clusters and their relation to protein function and structure. This large-scale analysis was performed using the whole Pfam database covering over 70% of the known protein sequence universe. Our program Sherlocc, detects statistically relevant conserved rare codon clusters and produces a user-friendly HTML output. RESULTS: Statistically significant rare codon clusters were detected in a multitude of Pfam protein families. The most statistically significant rare codon clusters were predominantly identified in N-terminal Pfam families. Many of the longest rare codon clusters are found in membrane-related proteins which are required to interact with other proteins as part of their function, for example in targeting or insertion. We identified some cases where rare codon clusters can play a regulating role in the folding of catalytically important domains. Our results support the existence of a widespread functional role for rare codon clusters across species. Finally, we developed an online filter-based search interface that provides access to Sherlocc results for all Pfam families. AVAILABILITY: The Sherlocc program and search interface are open access and are available at http://bcb.med.usherbrooke.ca

MOTIVATION: Protein movements form a continuum from large domain rearrangements (including folding and restructuring) to side-chain rotamer changes and small rearrangements. Understanding side-chain flexibility upon binding is important to understand molecular recognition events and predict ligand binding. METHODS: In the present work, we developed a well-curated non-redundant dataset of 188 proteins in pairs of structures in the Apo (unbound) and Holo (bound) forms to study the extent and the factors that guide side-chain rotamer changes upon binding. RESULTS: Our analysis shows that side-chain rotamer changes are widespread with only 10% of binding sites displaying no conformational changes. Overall, at most five rotamer changes account for the observed movements in 90% of the cases. Furthermore, rotamer changes are essential in 32% of flexible binding sites. The different amino acids have a 11-fold difference in their probability to undergo changes. Side-chain flexibility represents an intrinsic property of amino acids as it correlates well with configurational entropy differences. Furthermore, on average b-factors and solvent accessible surface areas can discriminate flexible side-chains in the Apo form. Finally, there is a rearrangement of the hydrogen-bonding network upon binding primarily with a loss of H-bonds with water molecules and a gain of H-bonds with protein residues for flexible residues. Interestingly, only 25% of side chains capable of forming H-bonds do so with the ligand upon binding. In terms of drug design, this last result shows that there is a large number of potential interactions that may be exploited to modulate the specificity and sensitivity of inhibitors. CONTACT: [email protected].