Enhanced antibody-antigen structure prediction from molecular docking using AlphaFold2.

Predicting the structure of antibody-antigen complexes has tremendous value in biomedical research but unfortunately suffers from a poor performance in real-life applications. AlphaFold2 (AF2) has provided renewed hope for improvements in the field of protein-protein docking but has shown limited success against antibody-antigen complexes due to the lack of co-evolutionary constraints. In this study, we used physics-based protein docking methods for building decoy sets consisting of low-energy docking solutions that were either geometrically close to the native structure (positives) or not (negatives).

Exploring rigid-backbone protein docking in biologics discovery: a test using the DARPin scaffold.

Accurate protein-protein docking remains challenging, especially for artificial biologics not coevolved naturally against their protein targets, like antibodies and other engineered scaffolds. We previously developed ProPOSE, an exhaustive docker with full atomistic details, which delivers cutting-edge performance by allowing side-chain rearrangements upon docking. However, extensive protein backbone flexibility limits its practical applicability as indicated by unbound docking tests.

Solvated interaction energy: from small-molecule to antibody drug design.

Scoring functions are ubiquitous in structure-based drug design as an aid to predicting binding modes and estimating binding affinities. Ideally, a scoring function should be broadly applicable, obviating the need to recalibrate and refit its parameters for every new target and class of ligands. Traditionally, drugs have been small molecules, but in recent years biologics, particularly antibodies, have become an increasingly important if not dominant class of therapeutics.

Optimizing Antibody-Antigen Binding Affinities with the ADAPT Platform.

The ADAPT (Assisted Design of Antibody and Protein Therapeutics) platform guides the selection of mutants that improve/modulate the affinity of antibodies and other biologics. Predicted affinities are based on a consensus z-score from three scoring functions. Computational predictions are interleaved with experimental validation, significantly enhancing the robustness of the design and selection of mutants.

Coevolved Canonical Loops Conformations of Single-Domain Antibodies: A Tale of Three Pockets Playing Musical Chairs.

Single-domain antibodies (sdAbs) are a promising class of biotherapeutics with unique structural traits within their paratope region. The distribution of canonical conformations explored by their complementarity determining region (CDR) loops differs to some extent from conventional two-chain Fv fragments of monoclonal antibodies (mAbs). In this study, we explored in detail the canonical structures of sdAb CDR-H1 and CDR-H2 loops and compared those with mAbs from the IGHV3 and IGHV1 gene families.

Antibody mutations favoring pH-dependent binding in solid tumor microenvironments: Insights from large-scale structure-based calculations.

Antibody-based therapeutics for treatment of various tumors have grown rapidly in recent years. Unfortunately, safety issues, attributed to off-tumor effects and cytotoxicity, are still a significant concern with the standard of care. Improvements to ensure targeted delivery of antitumor pharmaceuticals are desperately needed.

Redesigning an antibody H3 loop by virtual screening of a small library of human germline-derived sequences.

The design of superior biologic therapeutics, including antibodies and engineered proteins, involves optimizing their specific ability to bind to disease-related molecular targets. Previously, we developed and applied the Assisted Design of Antibody and Protein Therapeutics (ADAPT) platform for virtual affinity maturation of antibodies (Vivcharuk et al. in PLoS One 12(7):e0181490, https://doi.

A Novel Actinobacterial Cutinase Containing a Noncatalytic Polymer-Binding Domain.

The single putative cutinase-encoding gene from the genome of Kineococcus radiotolerans SRS30216 was cloned and expressed in Escherichia coli as a secreted fusion protein, designated YebF-KrCUT, where YebF is the extracellular carrier protein. The 294-amino-acid sequence of KrCUT is unique among currently characterized cutinases by having a C-terminal extension that consists of a short (Pro-Thr)-rich linker and a 55-amino-acid region resembling the substrate binding domain of poly(hydroxybutyrate) (PHB) depolymerases. Phylogenetically, KrCUT takes a unique position among known cutinases and cutinase-like proteins of bacterial and fungal origins.

Structure-based engineering of pH-dependent antibody binding for selective targeting of solid-tumor microenvironment.

Recent development of monoclonal antibodies as mainstream anticancer agents demands further optimization of their safety for use in humans. Potent targeting and/or effector activities on normal tissues is an obvious toxicity concern. Optimization of specific tumor targeting could be achieved by taking advantage of the extracellular acidity of solid tumors relative to normal tissues.

An accurate TMT-based approach to quantify and model lysine susceptibility to conjugation via N-hydroxysuccinimide esters in a monoclonal antibody.

Conjugation of small molecules to proteins through N-hydroxysuccinimide (NHS) esters results in a random distribution of small molecules on lysine residues and the protein N-terminus. While mass spectrometry methods have improved characterization of these protein conjugates, it remains a challenge to quantify the occupancy at individual sites of conjugation. Here, we present a method using Tandem Mass Tags (TMT) that enabled the accurate and sensitive quantification of occupancy at individual conjugation sites in the NIST monoclonal antibody.

ProPOSE: Direct Exhaustive Protein-Protein Docking with Side Chain Flexibility.

Despite decades of development, protein-protein docking remains a largely unsolved problem. The main difficulties are the immense space spanned by the translational and rotational degrees of freedom and the prediction of the conformational changes of proteins upon binding. FFT is generally the preferred method to exhaustively explore the translation-rotation space at a fine grid resolution, albeit with the trade-off of approximating force fields with correlation functions.

Binding pose and affinity prediction in the 2016 D3R Grand Challenge 2 using the Wilma-SIE method.

The Farnesoid X receptor (FXR) exhibits significant backbone movement in response to the binding of various ligands and can be a challenge for pose prediction algorithms. As part of the D3R Grand Challenge 2, we tested Wilma-SIE, a rigid-protein docking method, on a set of 36 FXR ligands for which the crystal structures had originally been blinded. These ligands covered several classes of compounds.

Assisted Design of Antibody and Protein Therapeutics (ADAPT).

Effective biologic therapeutics require binding affinities that are fine-tuned to their disease-related molecular target. The ADAPT (Assisted Design of Antibody and Protein Therapeutics) platform aids in the selection of mutants that improve/modulate the affinity of antibodies and other biologics. It uses a consensus z-score from three scoring functions and interleaves computational predictions with experimental validation, significantly enhancing the robustness of the design and selection of mutants.

Medicinal Chemistry Projects Requiring Imaginative Structure-Based Drug Design Methods.

Computational methods for docking small molecules to proteins are prominent in drug discovery. There are hundreds, if not thousands, of documented examples-and several pertinent cases within our research program. Fifteen years ago, our first docking-guided drug design project yielded nanomolar metalloproteinase inhibitors and illustrated the potential of structure-based drug design.

Assessment of Solvated Interaction Energy Function for Ranking Antibody-Antigen Binding Affinities.

Affinity modulation of antibodies and antibody fragments of therapeutic value is often required in order to improve their clinical efficacies. Virtual affinity maturation has the potential to quickly focus on the critical hotspot residues without the combinatorial explosion problem of conventional display and library approaches. However, this requires a binding affinity scoring function that is capable of ranking single-point mutations of a starting antibody.

Docking ligands into flexible and solvated macromolecules. 7. Impact of protein flexibility and water molecules on docking-based virtual screening accuracy.

The use of predictive computational methods in the drug discovery process is in a state of continual growth. Over the last two decades, an increasingly large number of docking tools have been developed to identify hits or optimize lead molecules through in-silico screening of chemical libraries to proteins. In recent years, the focus has been on implementing protein flexibility and water molecules.

Raman and ab initio investigation of aqueous Cu(I) chloride complexes from 25 to 80 °C.

Temperature-dependent Raman studies of aqueous copper(I) chloride complexes have been carried out up to 80 °C, along with supporting ab initio calculations for the species [CuCl(n)(H2O)m](1-n), n = 0-4 and hydration numbers m = 0-6. Normalized reduced isotropic Raman spectra were obtained from perpendicular and parallel polarization measurements, with perchlorate anion, ClO4(-), as an internal standard. Although the Raman spectra were not intense, spectra could be corrected by solvent baseline subtraction, to yield quantitative reduced molar scattering coefficients for the symmetric vibrational bands at 297 ± 3 and 247 ± 3 cm(-1).

Variability in docking success rates due to dataset preparation.

The results of cognate docking with the prepared Astex dataset provided by the organizers of the "Docking and Scoring: A Review of Docking Programs" session at the 241st ACS national meeting are presented. The MOE software with the newly developed GBVI/WSA dG scoring function is used throughout the study. For 80 % of the Astex targets, the MOE docker produces a top-scoring pose within 2 Å of the X-ray structure.

Integrating medicinal chemistry, organic/combinatorial chemistry, and computational chemistry for the discovery of selective estrogen receptor modulators with Forecaster, a novel platform for drug discovery.

As part of a large medicinal chemistry program, we wish to develop novel selective estrogen receptor modulators (SERMs) as potential breast cancer treatments using a combination of experimental and computational approaches. However, one of the remaining difficulties nowadays is to fully integrate computational (i.e.

Toward a computational tool predicting the stereochemical outcome of asymmetric reactions: development of the molecular mechanics-based program ACE and application to asymmetric epoxidation reactions.

The development and application of ACE, a program that predicts the stereochemical outcome of asymmetric reactions is presented. As major implementations, ACE includes a genetic algorithm to carry out an efficient global conformational search combined with a conjugate gradient minimization routine for local optimization and a corner flap algorithm to search ring conformations. Further improvements have been made that enable ACE to generate Boltzmann populations of conformations, to investigate highly asynchronous reactions, to compute fluctuating partial atomic charges and solvation energy and to automatically construct reactants and products from libraries of catalysts and substrates.

Rapid Prediction of Solvation Free Energy. 2. The First-Shell Hydration (FiSH) Continuum Model.

Local ordering of water in the first hydration shell around a solute is different from isotropic bulk water. This leads to effects that are captured by explicit solvation models and missed by continuum solvation models which replace the explicit waters with a continuous medium. In this paper, we introduce the First-Shell Hydration (FiSH) model as a first attempt to introduce first-shell effects within a continuum solvation framework.

Rapid Prediction of Solvation Free Energy. 1. An Extensive Test of Linear Interaction Energy (LIE).

The present study provides a comprehensive systematic analysis on the applicability of the linear interaction energy (LIE) approximation to the prediction of gas-to-water transfer (hydration) free energy. The study is based on molecular dynamics simulations in explicit solvent for an extensive and diverse hydration data set comprising 564 neutral compounds with measured hydration free energies, including a "traditional" data set and the more challenging drug-like SAMPL1 data set. A highly correlative LIE model was achieved without empirical scaling of the solute-solvent interaction energy terms along with a cavity term calibrated to the experiment.

Rapid prediction of solvation free energy. 3. Application to the SAMPL2 challenge.

The SAMPL2 hydration free energy blind prediction challenge consisted of a data set of 41 molecules divided into three subsets: explanatory, obscure and investigatory, where experimental hydration free energies were given for the explanatory, withheld for the obscure, and not known for the investigatory molecules. We employed two solvation models for this challenge, a linear interaction energy (LIE) model based on explicit-water molecular dynamics simulations, and the first-shell hydration (FiSH) continuum model previously calibrated to mimic LIE data. On the 23 compounds from the obscure (blind) dataset, the prospectively submitted LIE and FiSH models provided predictions highly correlated with experimental hydration free energy data, with mean-unsigned-errors of 1.

Docking ligands into flexible and solvated macromolecules. 3. Impact of input ligand conformation, protein flexibility, and water molecules on the accuracy of docking programs.

Several modifications and additions to Fitted1.5 led to the development of Fitted2.6.

Functional characterization and in silico docking of full and partial GluK2 kainate receptor agonists.

Two structural models have been developed to explain how agonist binding leads to ionotropic glutamate receptor (iGluR) activation. At alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) iGluRs, full and partial agonists close the agonist-binding domain (ABD) to different degrees whereas agonist-induced domain closure is apparently fixed at N-methyl-D-aspartate receptors. Although kainate (KA) iGluRs are thought to behave like AMPA receptors, the issue has not been formally tested because of the paucity of available receptor agonists.