Atomic Details of Carbon-Based Nanomolecules Interacting with Proteins

Exploring Scoring Function Space: Developing Computational Models for Drug Discovery

BACKGROUND

The idea of scoring function space established a systems-level approach to address the development of models to predict the affinity of drug molecules by those interested in drug discovery.

OBJECTIVE

Our goal here is to review the concept of scoring function space and how to explore it to develop machine learning models to address protein-ligand binding affinity.

METHODS

We searched the articles available in PubMed related to the scoring function space. We also utilized crystallographic structures found in the protein data bank (PDB) to represent the protein space.

RESULTS

The application of systems-level approaches to address receptor-drug interactions allows us to have a holistic view of the process of drug discovery. The scoring function space adds flexibility to the process since it makes it possible to see drug discovery as a relationship involving mathematical spaces.

CONCLUSION

The application of the concept of scoring function space has provided us with an integrated view of drug discovery methods. This concept is useful during drug discovery, where we see the process as a computational search of the scoring function space to find an adequate model to predict receptor-drug binding affinity.

Topology of electrostatic potential and electron density reveals a covalent to non-covalent carbon-carbon bond continuum

The covalent and non-covalent nature of carbon-carbon (CC) interactions in a wide range of molecular systems can be characterized using various methods, including the analysis of molecular electrostatic potential (MESP), represented as V(r), and the molecular electron density (MED), represented as ρ(r). These techniques provide valuable insights into the bonding between carbon atoms in different molecular environments. By uncovering a fundamental exponential relationship between the distance of the CC bond and the highest eigenvalue (λv1) of V(r) at the bond critical point (BCP), this study establishes the continuum model for all types of CC interactions, including transition states. The continuum model is further delineated into three distinct regions, namely covalent, borderline cases, and non-covalent, based on the gradient, , with the bond distance of the CC interaction. For covalent interactions, this parameter exhibits a more negative value than -5.0 a.u. Å-1, while for non-covalent interactions, it is less negative than -1.0 a.u. Å-1. Borderline cases, which encompass transition state structures, fall within the range of -1.0 to -5.0 a.u. Å-1. Furthermore, this study expands upon Popelier's analysis of the Laplacian of the MED, denoted as ∇2ρ, to encompass the entire spectrum of covalent, non-covalent, and borderline cases of CC interactions. Therefore, the present study presents compelling evidence supporting the concept of a continuum model for CC bonds in chemistry. Additionally, this continuum model is further explored within the context of C-N, C-O, C-S, N-N, O-O, and S-S interactions, albeit with a limited dataset.

Domain-agnostic predictions of nanoscale interactions in proteins and nanoparticles

Although challenging, the accurate and rapid prediction of nanoscale interactions has broad applications for numerous biological processes and material properties. While several models have been developed to predict the interaction of specific biological components, they use system-specific information that hinders their application to more general materials. Here we present NeCLAS, a general and efficient machine learning pipeline that predicts the location of nanoscale interactions, providing human-intelligible predictions. NeCLAS outperforms current nanoscale prediction models for generic nanoparticles up to 10-20 nm, reproducing interactions for biological and non-biological systems. Two aspects contribute to these results: a low-dimensional representation of nanoparticles and molecules (to reduce the effect of data uncertainty), and environmental features (to encode the physicochemical neighborhood at multiple scales). This framework has several applications, from basic research to rapid prototyping and design in nanobiotechnology.

Quenching of Protein Fluorescence by Fullerenol C60(OH)36 Nanoparticles

The effect of the interaction between fullerenol C₆₀(OH)₃₆ (FUL) and alcohol dehydrogenase (ADH) from Saccharomyces cerevisiae and human serum albumin (HSA) was studied by absorption spectroscopy, fluorescence spectroscopy, and time-resolved fluorescence spectroscopy. As shown in the study, the fluorescence intensities of ADH and HSA at excitation wavelengths λ_ex = 280 nm (Trp, Tyr) and λ_ex = 295 nm (Trp) are decreased with the increase in the FUL concentration. The results of time-resolved measurements indicate that both quenching mechanisms, dynamic and static, are present. The binding constant K_b and the number of binding sites were obtained for HSA and ADH. Thus, the results indicated the formation of FUL complexes and proteins. However, the binding of FUL to HSA is much stronger than that of ADH. The transfer of energy from the protein to FUL was also proved.

Enhanced Uptake and Phototoxicity of C60@albumin Hybrids by Folate Bioconjugation

Fullerenes are considered excellent photosensitizers, being highly suitable for photodynamic therapy (PDT). A lack of water solubility and low biocompatibility are, in many instances, still hampering the full exploitation of their potential in nanomedicine. Here, we used human serum albumin (HSA) to disperse fullerenes by binding up to five fullerene cages inside the hydrophobic cavities. Albumin was bioconjugated with folic acid to specifically address the folate receptors that are usually overexpressed in several solid tumors. Concurrently, tetramethylrhodamine isothiocyanate, TRITC, a tag for imaging, was conjugated to C₆₀@HSA in order to build an effective phototheranostic platform. The in vitro experiments demonstrated that: (i) HSA disperses C₆₀ molecules in a physiological environment, (ii) HSA, upon C₆₀ binding, maintains its biological identity and biocompatibility, (iii) the C₆₀@HSA complex shows a significant visible-light-induced production of reactive oxygen species, and (iv) folate bioconjugation improves both the internalization and the PDT-induced phototoxicity of the C₆₀@HSA complex in HeLa cells.

Fullerenes against COVID-19: Repurposing C60 and C70 to Clog the Active Site of SARS-CoV-2 Protease

The persistency of COVID-19 in the world and the continuous rise of its variants demand new treatments to complement vaccines. Computational chemistry can assist in the identification of moieties able to lead to new drugs to fight the disease. Fullerenes and carbon nanomaterials can interact with proteins and are considered promising antiviral agents. Here, we propose the possibility to repurpose fullerenes to clog the active site of the SARS-CoV-2 protease, M^pro. Through the use of docking, molecular dynamics, and energy decomposition techniques, it is shown that C₆₀ has a substantial binding energy to the main protease of the SARS-CoV-2 virus, M^pro, higher than masitinib, a known inhibitor of the protein. Furthermore, we suggest the use of C₇₀ as an innovative scaffold for the inhibition of SARS-CoV-2 M^pro. At odds with masitinib, both C₆₀ and C₇₀ interact more strongly with SARS-CoV-2 M^pro when different protonation states of the catalytic dyad are considered. The binding of fullerenes to M^pro is due to shape complementarity, i.e., vdW interactions, and is aspecific. As such, it is not sensitive to mutations that can eliminate or invert the charges of the amino acids composing the binding pocket. Fullerenic cages should therefore be more effective against the SARS-CoV-2 virus than the available inhibitors such as masinitib, where the electrostatic term plays a crucial role in the binding.

Dissecting the Supramolecular Dispersion of Fullerenes by Proteins/Peptides: Amino Acid Ranking and Driving Forces for Binding to C60

Molecular dynamics simulations were used to quantitatively investigate the interactions between the twenty proteinogenic amino acids and C₆₀. The conserved amino acid backbone gave a constant energetic interaction ~5.4 kcal mol⁻¹, while the contribution to the binding due to the amino acid side chains was found to be up to ~5 kcal mol⁻¹ for tryptophan but lower, to a point where it was slightly destabilizing, for glutamic acid. The effects of the interplay between van der Waals, hydrophobic, and polar solvation interactions on the various aspects of the binding of the amino acids, which were grouped as aromatic, charged, polar and hydrophobic, are discussed. Although π–π interactions were dominant, surfactant-like and hydrophobic effects were also observed. In the molecular dynamics simulations, the interacting residues displayed a tendency to visit configurations (i.e., regions of the Ramachandran plot) that were absent when C₆₀ was not present. The amino acid backbone assumed a “tepee-like” geometrical structure to maximize interactions with the fullerene cage. Well-defined conformations of the most interactive amino acids (Trp, Arg, Met) side chains were identified upon C₆₀ binding.

Porous assembly of an antifungal protein mediated by zinc and sulfonato-calix[8]arene

Controlled protein assembly holds great potential in the fabrication of biohybrid materials. However, the structural diversity and complexity of proteins present an obstacle to controlled assembly. Supramolecular chemistry is a possible solution as it offers tools to mediate self-assembly with molecular precision. This paper deals with the calixarene- and zinc-mediated assembly and crystallization of the histidine-rich Penicillium chrysogenum antifungal protein B (PAFB). We report crystal structures of pure PAFB, PAFB in complex with Zn²⁺, and the ternary complex of PAFB, Zn²⁺ and sulfonato-calix[8]arene (sclx₈). A comparison of the three crystal structures revealed the structural plasticity of PAFB. While the flexible and highly anionic sclx₈ resulted in large molecular weight aggregates of PAFB in solution, diffraction-quality crystals of PAFB-sclx₈ were not obtained. We report crystals of PAFB-Zn²⁺-sclx₈ in which a trinuclear zinc cluster occurred adjacent to a calixarene binding site. Interestingly, the combination of sclx₈ complexation and zinc coordination resulted in a porous framework with channels of circa 2 nm diameter. Detailed analysis of the crystal structure highlighted novel molecular recognition features. This research enriches the set of supramolecular interactions available to promote protein assembly.

NMR Spectroscopy of supramolecular chemistry on protein surfaces

As one of the few analytical methods that offer atomic resolution, NMR spectroscopy is a valuable tool to study the interaction of proteins with their interaction partners, both biomolecules and synthetic ligands. In recent years, the focus in chemistry has kept expanding from targeting small binding pockets in proteins to recognizing patches on protein surfaces, mostly via supramolecular chemistry, with the goal to modulate protein–protein interactions. Here we present NMR methods that have been applied to characterize these molecular interactions and discuss the challenges of this endeavor.