RESUMEN

The number of applications for nanobodies is steadily expanding, positioning these molecules as fast-growing biologic products in the biotechnology market. Several of their applications require protein engineering, which in turn would greatly benefit from having a reliable structural model of the nanobody of interest. However, as with antibodies, the structural modeling of nanobodies is still a challenge. With the rise of artificial intelligence (AI), several methods have been developed in recent years that attempt to solve the problem of protein modeling. In this study, we have compared the performance in nanobody modeling of several state-of-the-art AI-based programs, either designed for general protein modeling, such as AlphaFold2, OmegaFold, ESMFold, and Yang-Server, or specifically designed for antibody modeling, such as IgFold, and Nanonet. While all these programs performed rather well in constructing the nanobody framework and CDRs 1 and 2, modeling CDR3 still represents a big challenge. Interestingly, tailoring an AI method for antibody modeling does not necessarily translate into better results for nanobodies.

Asunto(s)

RESUMEN

We present NbThermo-a first-in-class database that collects melting temperatures (Tm), amino acid sequences and several other categories of useful data for hundreds of nanobodies (Nbs), compiled from an extensive literature search. This so-far unique database currently contains up-to-date, manually curated data for 564 Nbs. It represents a contribution to efforts aimed at developing new algorithms for reliable Tm prediction to assist Nb engineering for a wide range of applications of these unique biomolecules. Nbs from the two most common source organisms-llama and camel-show similar distributions of melting temperatures. A first exploratory research that takes advantage of this large data collection evidences that understanding the structural bases of Nb thermostability is a complex task, since there are no apparent differences in sequence patterns between the frameworks of Nbs with lower and higher melting temperatures, indicating that the highly variable loops play a relevant role in defining Nb thermostability. Database URL https://valdes-tresanco-ms.github.io/NbThermo.

Asunto(s)

RESUMEN

Computational alanine scanning with the molecular mechanics generalized Born surface area (MM/GBSA) method constitutes a widely used approach for identifying critical residues at protein-protein interfaces. Despite its popularity, the MM/GBSA method still has certain drawbacks due to its dependence on many factors. Here, we performed a systematical study on the impact of four different parameters, namely, the internal dielectric constant, the generalized Born model, the entropic term, and the inclusion of structural waters on the accuracy of computational alanine scanning calculations with the MM/GBSA method. Our results show that the internal dielectric constant is the most critical parameter for getting accurate predictions. The introduction of entropy and interfacial water molecules decreased the quality of the predictions, while the generalized Born model had little to no effect. Considering the significance of the internal dielectric value, we proposed a methodology based on the energetic predominance of a particular set of amino acids at the protein-protein interface for selecting an appropriate value for this variable. We hope that these results serve as a guideline for future studies of protein-protein complexes using the MM/GBSA method.

Asunto(s)

RESUMEN

Invasive fungal infections represent a public health problem that worsens over the years with the increasing resistance to current antimycotic agents. Therefore, there is a compelling medical need of widening the antifungal drug repertoire, following different methods such as drug repositioning, identification and validation of new molecular targets and developing new inhibitors against these targets. In this work we developed a structure-based strategy for drug repositioning and new drug design, which can be applied to infectious fungi and other pathogens. Instead of applying the commonly accepted off-target criterion to discard fungal proteins with close homologues in humans, the core of our approach consists in identifying fungal proteins with active sites that are structurally similar, but preferably not identical to binding sites of proteins from the so-called "human pharmacolome". Using structural information from thousands of human protein target-inhibitor complexes, we identified dozens of proteins in fungal species of the genera Histoplasma, Candida, Cryptococcus, Aspergillus and Fusarium, which might be exploited for drug repositioning and, more importantly, also for the design of new fungus-specific inhibitors. As a case study, we present the in vitro experiments performed with a set of selected inhibitors of the human mitogen-activated protein kinases 1/2 (MEK1/2), several of which showed a marked cytotoxic activity in different fungal species.

Asunto(s)

RESUMEN

Single domain antibodies from camelids, or nanobodies, are a unique class of antibody fragments with several advantageous characteristics: small monomeric size, high stability and solubility and easy tailoring for multiple applications. Nanobodies are gaining increasing acceptance as diagnostic tools and promising therapeutic agents in cancer and other diseases. While most nanobodies are obtained from immunized animals of the camelid family, a few synthetic nanobody libraries constructed in recent years have shown the capability of generating high quality nanobodies in terms of affinity and stability. Since this synthetic approach has important advantages over the use of animals, the recent advances are indeed encouraging. Here we review over a dozen synthetic nanobody libraries reported so far and discuss the different approaches followed in their construction and validation, with an emphasis on framework and hypervariable loop design as critical issues defining their potential as high-class nanobody sources.

Asunto(s)

RESUMEN

OBJECTIVE: To design and construct a new synthetic nanobody library using a structure-based approach that seeks to maintain high protein stability and increase the number of functional variants within the combinatorial space of mutations. RESULTS: Synthetic nanobody (Nb) libraries are emerging as an attractive alternative to animal immunization for the selection of stable, high affinity Nbs. Two key features define a synthetic Nb library: framework selection and CDR design. We selected the universal VHH framework from the cAbBCII10 Nb. CDR1 and CDR2 were designed with the same fixed length as in cAbBCII10, while for CDR3 we chose a 14-long loop, which creates a convex binding site topology. Based on the analysis of the cAbBCII10 crystal structure, we carefully selected the positions to be randomized and tailored the codon usage at each position, keeping at particular places amino acids that guarantee stability, favoring properties like polarity at solvent-exposed positions and avoiding destabilizing amino acids. Gene synthesis and library construction were carried out by GenScript, using our own phagemid vector. The constructed library has an estimated size of 1.75 × 108. NGS showed that the amino acid diversity and frequency at each randomized position are the expected from the codon usage.

Asunto(s)

RESUMEN

Vps34 is the only isoform of the PI3K family in fungi, making this protein an attractive target to develop new treatments against pathogenic fungi. The high structural similarity between the active sites of the human and fungal Vps34 makes repurposing of human Vps34 inhibitors an appealing strategy. Nonetheless, while some of the cross-reactive inhibitors might have the potential to treat fungal infections, a safer approach to prevent undesired side effects would be to identify molecules that specifically inhibit the fungal Vps34. This study presents the parameterization of four LIE models for estimating the binding free energy of Vps34-inhibitor complexes. Two models are parameterized using a multiparametric linear regression leaving one or more free parameters, while the other two are based on the LIE-D model. All of the models show good predictive capacity (R 2 > 0.7, r > 0.85) and a low mean absolute error (MAE < 0.71 kcal/mol). The current study highlights the advantages of LIE-D-derived models when predicting the weight of the different contributions to the binding free energy. It is expected that this study will provide researchers with a valuable tool to identify new Vps34 inhibitors for relevant applications such as cancer treatment and the development of new antimicrobial agents.

RESUMEN

Molecular mechanics/Poisson-Boltzmann (Generalized-Born) surface area is one of the most popular methods to estimate binding free energies. This method has been proven to balance accuracy and computational efficiency, especially when dealing with large systems. As a result of its popularity, several programs have been developed for performing MM/PB(GB)SA calculations within the GROMACS community. These programs, however, present several limitations. Here we present gmx_MMPBSA, a new tool to perform end-state free energy calculations from GROMACS molecular dynamics trajectories. gmx_MMPBSA provides the user with several options, including binding free energy calculations with different solvation models (PB, GB, or 3D-RISM), stability calculations, computational alanine scanning, entropy corrections, and binding free energy decomposition. Noteworthy, several promising methodologies to calculate relative binding free energies such as alanine scanning with variable dielectric constant and interaction entropy have also been implemented in gmx_MMPBSA. Two additional tools-gmx_MMPBSA_test and gmx_MMPBSA_ana-have been integrated within gmx_MMPBSA to improve its usability. Multiple illustrating examples can be accessed through gmx_MMPBSA_test, while gmx_MMPBSA_ana provides fast, easy, and efficient access to different graphics plotted from gmx_MMPBSA output files. The latest version (v1.4.3, 26/05/2021) is available free of charge (documentation, test files, and tutorials included) at https://github.com/Valdes-Tresanco-MS/gmx_MMPBSA.

RESUMEN

AMDock (Assisted Molecular Docking) is a user-friendly graphical tool to assist in the docking of protein-ligand complexes using Autodock Vina and AutoDock4, including the option of using the Autodock4Zn force field for metalloproteins. AMDock integrates several external programs (Open Babel, PDB2PQR, AutoLigand, ADT scripts) to accurately prepare the input structure files and to optimally define the search space, offering several alternatives and different degrees of user supervision. For visualization of molecular structures, AMDock uses PyMOL, starting it automatically with several predefined visualization schemes to aid in setting up the box defining the search space and to visualize and analyze the docking results. One particularly useful feature implemented in AMDock is the off-target docking procedure that allows to conduct ligand selectivity studies easily. In summary, AMDock's functional versatility makes it a very useful tool to conduct different docking studies, especially for beginners. The program is available, either for Windows or Linux, at https://github.com/Valdes-Tresanco-MS . REVIEWERS: This article was reviewed by Alexander Krah and Thomas Gaillard.

Asunto(s)

RESUMEN

The calculation of absolute binding affinities for protein-inhibitor complexes remains as one of the main challenges in computational structure-based ligand design. The present work explored the calculations of surface fractal dimension (as a measure of surface roughness) and the relationship with experimental binding free energies of Plasmepsin II complexes. Plasmepsin II is an attractive target for novel therapeutic compounds to treat malaria. However, the structural flexibility of this enzyme is a drawback when searching for specific inhibitors. Concerning that, we performed separate explicitly solvated molecular dynamics simulations using the available high-resolution crystal structures of different Plasmepsin II complexes. Molecular dynamics simulations allowed a better approximation to systems dynamics and, therefore, a more reliable estimation of surface roughness. This constitutes a novel approximation in order to obtain more realistic values of fractal dimension, because previous works considered only x-ray structures. Binding site fractal dimension was calculated considering the ensemble of structures generated at different simulation times. A linear relationship between binding site fractal dimension and experimental binding free energies of the complexes was observed within 20 ns. Previous studies of the subject did not uncover this relationship. Regression model, coined FD model, was built to estimate binding free energies from binding site fractal dimension values. Leave-one-out cross-validation showed that our model reproduced accurately the absolute binding free energies for our training set (R2  = 0.76; <|error|> =0.55 kcal/mol; SDerror  = 0.19 kcal/mol). The fact that such a simple model may be applied raises some questions that are addressed in the article.

Asunto(s)