Exploring avidity: understanding the potential gains in functional affinity and target residence time of bivalent and heterobivalent ligands

Bivalent ligands are increasingly important therapeutic agents. Although the naturally occurring antibodies are predominant, it is becoming more common to combine different antibody fragments or even low molecular weight compounds to generate heterobivalent ligands. Such ligands exhibit markedly increased affinity (i.e. avidity) and target residence time when both pharmacophores can bind simultaneously to their target sites. This is because binding of one pharmacophore forces the second tethered one to stay close to its corresponding site. This 'forced proximity' favours its binding and rebinding (once dissociated) to that site. However, rebinding will also take place when the diffusion of freshly dissociated ligands is merely slowed down. The present differential equation-based simulations explore the way both situations affect ligand binding. Both delay the attainment of binding equilibrium (resulting in steep saturation curves) and also increase the target residence time. Competitive ligands are able to interfere in a concentration-dependent manner, although much higher concentrations are required in the 'forced proximity' situation. Also, it is only in that situation that the ligand shows increased affinity. These simulations shed light on two practical consequences. Depending on the pharmacokinetic half-life of the bivalent ligand in the body, it may not have sufficient time to achieve equilibrium with the target. This will result in lower potency than expected, although it would have significant advantages in terms of residence time. In in vitro experiments, the manifestation of steep saturation curves and of accelerated dissociation in the presence of competitive ligands could mistakenly be interpreted as evidence for non-competitive, allosteric interactions.

The Multivalent Role of Fibronectin-Binding Proteins A and B (FnBPA and FnBPB) of Staphylococcus aureus in Host Infections

Staphylococcus aureus, one of the most important human pathogens, is the causative agent of several infectious diseases including sepsis, pneumonia, osteomyelitis, endocarditis and soft tissue infections. This pathogenicity is due to a multitude of virulence factors including several cell wall-anchored proteins (CWA). CWA proteins have modular structures with distinct domains binding different ligands. The majority of S. aureus strains express two CWA fibronectin (Fn)-binding adhesins FnBPA and FnBPB (Fn-binding proteins A and B), which are encoded by closely related genes. The N-terminus of FnBPA and FnBPB comprises an A domain which binds ligands such as fibrinogen, elastin and plasminogen. The A domain of FnBPB also interacts with histones and this binding results in the neutralization of the antimicrobial activity of these molecules. The C-terminal moiety of these adhesins comprises a long, intrinsically disordered domain composed of 11/10 fibronectin-binding repeats. These repetitive motifs of FnBPs promote invasion of cells that are not usually phagocytic via a mechanism by which they interact with integrin α₅β₁ through a Fn mediated-bridge. The FnBPA and FnBPB A domains engage in homophilic cell-cell interactions and promote biofilm formation and enhance platelet aggregation. In this review we update the current understanding of the structure and functional properties of FnBPs and emphasize the role they may have in the staphylococcal infections.

Mechanism of E-cadherin dimerization probed by NMR relaxation dispersion

Epithelial cadherin (E-cadherin), a member of the classical cadherin family, mediates calcium-dependent homophilic cell-cell adhesion. Crystal structures of classical cadherins reveal an adhesive dimer interface featuring reciprocal exchange of N-terminal β-strands between two protomers. Previous work has identified a putative intermediate (called the "X-dimer") in the dimerization pathway of wild-type E-cadherin EC1-EC2 domains, based on crystal structures of mutants not capable of strand swapping and on deceleration of binding kinetics by mutations at the X-dimer interface. In the present work, NMR relaxation dispersion spectroscopy is used to directly observe and characterize intermediate states without the need to disrupt the strand-swapped binding interface by mutagenesis. The results indicate that E-cadherin forms strand-swapped dimers predominantly by a mechanism in which formation of a weak and short-lived X-dimer-like state precedes the conformational changes required for formation of the mature strand-swapped dimeric structure. Disruption of this intermediate state through mutation reduces both association and dissociation rates by factors of ~10(4), while minimally perturbing affinity. The X-dimer interface lowers the energy barrier associated with strand swapping and enables E-cadherins to form strand-swapped dimers at a rate consistent with residence times in adherens junctions.

Mechanism of β-actin mRNA Recognition by ZBP1

Zipcode binding protein 1 (ZBP1) is an oncofetal RNA-binding protein that mediates the transport and local translation of β-actin mRNA by the KH3-KH4 di-domain, which is essential for neuronal development. The high-resolution structures of KH3-KH4 with their respective target sequences show that KH4 recognizes a non-canonical GGA sequence via an enlarged and dynamic hydrophobic groove, whereas KH3 binding to a core CA sequence occurs with low specificity. A data-informed kinetic simulation of the two-step binding reaction reveals that the overall reaction is driven by the second binding event and that the moderate affinities of the individual interactions favor RNA looping. Furthermore, the concentration of ZBP1, but not of the target RNA, modulates the interaction, which explains the functional significance of enhanced ZBP1 expression during embryonic development.

•

The dynamic groove of ZBP1’s KH4 domain allows recognition of a G-rich RNA sequence

•

ZBP1’s KH3 and KH4 domains bind their target RNA sequences with similar affinities

•

RNA looping drives the ZBP1-β-actin interaction

•

The protein, rather than the RNA, concentration regulates ZBP1-β-actin mRNA binding

Rate Constants and Mechanisms of Protein-Ligand Binding

Whereas protein-ligand binding affinities have long-established prominence, binding rate constants and binding mechanisms have gained increasing attention in recent years. Both new computational methods and new experimental techniques have been developed to characterize the latter properties. It is now realized that binding mechanisms, like binding rate constants, can and should be quantitatively determined. In this review, we summarize studies and synthesize ideas on several topics in the hope of providing a coherent picture of and physical insight into binding kinetics. The topics include microscopic formulation of the kinetic problem and its reduction to simple rate equations; computation of binding rate constants; quantitative determination of binding mechanisms; and elucidation of physical factors that control binding rate constants and mechanisms.

Two Differential Binding Mechanisms of FG-Nucleoporins and Nuclear Transport Receptors

Phenylalanine-glycine-rich nucleoporins (FG-Nups) are intrinsically disordered proteins, constituting the selective barrier of the nuclear pore complex (NPC). Previous studies showed that nuclear transport receptors (NTRs) were found to interact with FG-Nups by forming an “archetypal-fuzzy” complex through the rapid formation and breakage of interactions with many individual FG motifs. Here, we use single-molecule studies combined with atomistic simulations to show that, in sharp contrast, FG-Nup214 undergoes a coupled reconfiguration-binding mechanism when interacting with the export receptor CRM1. Association and dissociation rate constants are more than an order of magnitude lower than in the archetypal-fuzzy complex between FG-Nup153 and NTRs. Unexpectedly, this behavior appears not to be encoded selectively into CRM1 but rather into the FG-Nup214 sequence. The same distinct binding mechanisms are unperturbed in O-linked β-N-acetylglucosamine-modified FG-Nups. Our results have implications for differential roles of distinctly spatially distributed FG-Nup⋅NTR interactions in the cell.

•

Identification of two differential binding mechanisms in the nuclear transport pathway

•

FG-Nup214 does not bind CRM1 via an archetypal-fuzzy complex

•

Glycosylated FG-Nups maintain their NTR-binding mechanisms

•

Linker regions of FG-Nups may have functional relevance to the binding mechanism

Cell membranes… and how long drugs may exert beneficial pharmacological activity in vivo

The time course of the beneficial pharmacological effect of a drug has long been considered to depend merely on the temporal fluctuation of its free concentration. Only in the last decade has it become widely accepted that target-binding kinetics can also affect in vivo pharmacological activity. Although current reviews still essentially focus on genuine dissociation rates, evidence is accumulating that additional micro-pharmacokinetic (PK) and -pharmacodynamic (PD) mechanisms, in which the cell membrane plays a central role, may also increase the residence time of a drug on its target. The present review provides a compilation of otherwise widely dispersed information on this topic. The cell membrane can intervene in drug binding via the following three major mechanisms: (i) by acting as a sink/repository for the drug; (ii) by modulating the conformation of the drug and even by participating in the binding process; and (iii) by facilitating the approach (and rebinding) of the drug to the target. To highlight these mechanisms, we focus on drugs that are currently used in clinical therapy, such as the antihypertensive angiotensin II type 1 receptor antagonist candesartan, the atypical antipsychotic agent clozapine and the bronchodilator salmeterol. Although the role of cell membranes in PK-PD modelling is gaining increasing interest, many issues remain unresolved. It is likely that novel biophysical and computational approaches will provide improved insights in the near future.

Scaffold rearrangement of dihydroxypyrimidine inhibitors of HIV integrase: Docking model revisited

A series of dihydroxypyrimidine (DHP) derivatives were designed as inhibitors of HIV integrase (IN) based on known homology models. Through chemical synthesis and biochemical assays it was found that the activity profile of these compounds largely deviates from predictions with existing models. With the recently disclosed IN crystal structure of prototype foamy virus (PFV), a new HIV IN homology model was constructed featuring a critical IN/DNA interface previously lacking. With this new model, docking results completely corroborated observed biological activities. This new model should provide a more accurate and improved platform for the design of new inhibitors of HIV IN.

A computational comparative analysis of the binding mechanism of molnupiravir's active metabolite to RNA-dependent RNA polymerase of wild-type and Delta subvariant AY.4 of SARS-CoV-2

The antiviral drug molnupiravir targets the SARS-CoV-2 RNA-dependent RNA polymerase (RdRP) enzyme. Early treatment with molnupiravir reduced the risk of hospitalization or death in at-risk, unvaccinated adults with COVID-19, according to phase 3 clinical trials. Many mutations have occurred within this virus as a result of its widespread distribution. The current study sought to determine whether mutations in the RdRP of Delta subvariant AY.4 (D-AY.4 RdRP) influence the interaction of the enzyme with molnupiravir triphosphate (MTP), the active metabolite of molnupiravir. The interactions between the wild-type (WT) RdRP and D-AY.4 RdRP with MTP were evaluated based on molecular docking and dynamic simulation (MD) studies. The results show that the MTP interaction is stronger and more stable with D-AY.4 RdRP than with WT RdRP. This study provides insight into the potential significance of administering MTP to patients infected with D-AY.4 RdRP, which may have a more favorable chance of alleviating the illness. According to the findings of this study, MTP has a high likelihood of becoming widely used as an anti-SARS-CoV-2 agent. The fact that MTP is not only cytotoxic but also mutagenic to mammalian cells, as well as the possibility that it may cause DNA damage in the host, have all been raised as potential concerns.

Investigations of FAK inhibitors: a combination of 3D-QSAR, docking, and molecular dynamics simulations studies

Focal adhesion kinase (FAK) is one kind of tyrosine kinases that modulates integrin and growth factor signaling pathways, which is a promising therapeutic target because of involving in cancer cell migration, proliferation, and survival. To investigate the mechanism between FAK and triazinic inhibitors and design high activity inhibitors, a molecular modeling integrated with 3D-QSAR, molecular docking, molecular dynamics simulations, and binding free energy calculations was performed. The optimum CoMFA and CoMSIA models showed good reliability and satisfactory predictability (with Q² = 0.663, R² = 0.987, [Formula: see text] = 0.921 and Q² = 0.670, R² = 0.981, [Formula: see text] = 0.953). Its contour maps could provide structural features to improve inhibitory activity. Furthermore, a good consistency between contour maps, docking, and molecular dynamics simulations strongly demonstrates that the molecular modeling is reliable. Based on it, we designed several new compounds and their inhibitory activities were validated by the molecular models. We expect our studies could bring new ideas to promote the development of novel inhibitors with higher inhibitory activity for FAK.

Unusual binding interactions in PDZ domain crystal structures help explain binding mechanisms

PDZ domains most commonly bind the C-terminus of their protein targets. Typically the C-terminal four residues of the protein target are considered as the binding motif, particularly the C-terminal residue (P0) and third-last residue (P-2) that form the major contacts with the PDZ domain's "binding groove". We solved crystal structures of seven human PDZ domains, including five of the seven PDLIM family members. The structures of GRASP, PDLIM2, PDLIM5, and PDLIM7 show a binding mode with only the C-terminal P0 residue bound in the binding groove. Importantly, in some cases, the P-2 residue formed interactions outside of the binding groove, providing insight into the influence of residues remote from the binding groove on selectivity. In the GRASP structure, we observed both canonical and noncanonical binding in the two molecules present in the asymmetric unit making a direct comparison of these binding modes possible. In addition, structures of the PDZ domains from PDLIM1 and PDLIM4 also presented here allow comparison with canonical binding for the PDLIM PDZ domain family. Although influenced by crystal packing arrangements, the structures nevertheless show that changes in the positions of PDZ domain side-chains and the alpha B helix allow noncanonical binding interactions. These interactions may be indicative of intermediate states between unbound and fully bound PDZ domain and target protein. The noncanonical "perpendicular" binding observed potentially represents the general form of a kinetic intermediate. Comparison with canonical binding suggests that the rearrangement during binding involves both the PDZ domain and its ligand.

A conserved mechanism for binding of p53 DNA-binding domain and anti-apoptotic Bcl-2 family proteins

The molecular interaction between tumor suppressor p53 and the anti-apoptotic Bcl-2 family proteins plays an essential role in the transcription-independent apoptotic pathway of p53. In this study, we investigated the binding of p53 DNA-binding domain (p53DBD) with the anti-apoptotic Bcl-2 family proteins, Bcl-w, Mcl-1, and Bcl-2, using GST pull-down assay and NMR spectroscopy. The GST pull-down assays and NMR experiments demonstrated the direct binding of the p53DBD with Bcl-w, Mcl-1, and Bcl-2. Further, NMR chemical shift perturbation data showed that Bcl-w and Mcl-1 bind to the positively charged DNA-binding surface of p53DBD. Noticeably, the refined structural models of the complexes between p53DBD and Bcl-w, Mcl-1, and Bcl-2 showed that the binding mode of p53DBD is highly conserved among the anti-apoptotic Bcl-2 family proteins. Furthermore, the chemical shift perturbations on Bcl-w, Mcl-1, and Bcl-2 induced by p53DBD binding occurred not only at the p53DBD-binding acidic region but also at the BH3 peptide-binding pocket, which suggests an allosteric conformational change similar to that observed in Bcl-XL. Taken altogether, our results revealed a structural basis for a conserved binding mechanism between p53DBD and the anti-apoptotic Bcl-2 family proteins, which shed light on to the molecular understanding of the transcription-independent apoptosis pathway of p53.

Probing the correlation between ligand efficacy and conformational diversity at the α1A-adrenoreceptor reveals allosteric coupling of its microswitches

G protein-coupled receptors (GPCRs) use a series of conserved microswitches to transmit signals across the cell membrane via an allosteric network encompassing the ligand-binding site and the G protein-binding site. Crystal structures of GPCRs provide snapshots of their inactive and active states, but poorly describe the conformational dynamics of the allosteric network that underlies GPCR activation. Here, we analyzed the correlation between ligand binding and receptor conformation of the α1A-adrenoreceptor, a GPCR that stimulates smooth muscle contraction in response to binding noradrenaline. NMR of [13CϵH3]methionine-labeled α1A-adrenoreceptor variants, each exhibiting differing signaling capacities, revealed how different classes of ligands modulate the conformational equilibria of this receptor. [13CϵH3]Methionine residues near the microswitches exhibited distinct states that correlated with ligand efficacies, supporting a conformational selection mechanism. We propose that allosteric coupling among the microswitches controls the conformation of the α1A-adrenoreceptor and underlies the mechanism of ligand modulation of GPCR signaling in cells.

Influence of EGCG on α-synuclein (αS) aggregation and identification of their possible binding mode: A computational study using molecular dynamics simulation

The accumulation of intrinsically disordered α-synuclein (αS) protein that can form β-sheet-rich fibrils is linked to Parkinson's disease. (-)-Epigallocatechin-3-gallate (EGCG) is the most abundant active component in green tea and can inhibit the fibrillation of αS. The elucidation of this molecular mechanism will be helpful to understand the inhibition mechanism of EGCG to the fibrillation of αS and also to find more potential small molecules that can inhibit the aggregation of αS. In this work, to study the influence of EGCG on the structure of β-sheet-rich fibrils of αS and identification of their possible binding mode, molecular dynamics simulations of pentamer and decamer aggregates of αS in complex with EGCG were performed. The obtained results indicate that EGCG can remodel the αS fibrils and break the initial ordered pattern by reducing the β-sheet content. EGCG can also break the Greek conformation of αS by the disappeared H-bond in the secondary structure of turn. The results from our study can not only reveal the specific interaction between EGCG and β-sheet-rich fibrils of αS, but also provide the useful guidance for the discovery of other potential inhibitors.

Impact of the K24N mutation on the transactivation domain of p53 and its binding to murine double-minute clone 2

The level of the p53 transcription factor is negatively regulated by the E3 ubiquitin ligase murine double-minute clone 2 (MDM2). The interaction between p53 and MDM2 is essential for the maintenance of genomic integrity for most eukaryotes. Previous structural studies revealed that MDM2 binds to p53 transactivation domain (p53TAD) from residues 17 to 29. The K24N mutation of p53TAD changes a lysine at position 24 to an asparagine. This mutation occurs naturally in the bovine family and is also found in a rare form of human gestational cancer called choriocarcinoma. In this study, we have investigated how the K24N mutation affects the affinity, structure, and dynamics of p53TAD binding to MDM2. Nuclear magnetic resonance studies of p53TAD show that the K24N mutant is more flexible and has less transient helical secondary structure than the wild type. Isothermal titration calorimetry measurements demonstrate that these changes in structure and dynamics do not significantly change the binding affinity for p53TAD-MDM2. Finally, free-energy perturbation and standard molecular dynamic simulations suggest the negligible affinity change is due to a compensating interaction energy between the K24N mutant and the MDM2 when it is bound. Overall, the data suggest that the K24N-MDM2 complex is able to, at least partly, compensate for an increase in the conformational entropy in unbound K24N with an increase in the bound-state electrostatic interaction energy.

Extreme Fuzziness: Direct Interactions between Two IDPs

Protein interactions involving intrinsically disordered proteins (IDPs) greatly extend the range of binding mechanisms available to proteins. In interactions employing coupled folding and binding, IDPs undergo disorder-to-order transitions to form a complex with a well-defined structure. In many other cases, IDPs retain structural plasticity in the final complexes, which have been defined as the fuzzy complexes. While a large number of fuzzy complexes have been characterized with variety of fuzzy patterns, many of the interactions are between an IDP and a structured protein. Thus, whether two IDPs can interact directly to form a fuzzy complex without disorder-to-order transition remains an open question. Recently, two studies of interactions between IDPs (4.1G-CTD/NuMA and H1/ProTα) have found a definite answer to this question. Detailed characterizations combined with nuclear magnetic resonance (NMR), single-molecule Förster resonance energy transfer (smFRET) and molecular dynamics (MD) simulation demonstrate that direct interactions between these two pairs of IDPs do form fuzzy complexes while retaining the conformational dynamics of the isolated proteins, which we name as the extremely fuzzy complexes. Extreme fuzziness completes the full spectrum of protein-protein interaction modes, suggesting that a more generalized model beyond existing binding mechanisms is required. Previous models of protein interaction could be applicable to some aspects of the extremely fuzzy interactions, but in more general sense, the distinction between native and nonnative contacts, which was used to understand protein folding and binding, becomes obscure. Exploring the phenomenon of extreme fuzziness may shed new light on molecular recognition and drug design.

The mechanism of high affinity pentasaccharide binding to antithrombin, insights from Gaussian accelerated molecular dynamics simulations

The activity of antithrombin (AT), a serpin protease inhibitor, is enhanced by heparin and heparin analogs against its target proteases, mainly thrombin, factors Xa and IXa. Considerable amount of information is available on the multistep mechanism of the heparin pentasaccharide binding and conformational activation. However, much of the details were inferred from 'static' structures obtained by X-ray diffraction. Moreover, limited information is available for the early steps of binding mechanism other than kinetic studies with various ligands. To gain insights into these processes, we performed enhanced sampling molecular dynamics (MD) simulations using the Gaussian Accelerated Molecular Dynamics (GAMD) method, applied previously in drug binding studies. We were able to observe the binding of the pentasaccharide idraparinux to a 'non-activated' AT conformation in two separate trajectories with low root mean square deviation (RMSD) values compared to X-ray structures of the bound state. These trajectories along with further simulations of the AT-pentasaccharide complex provided insights into the mechanisms of multiple conformational transitions, including the expulsion of the hinge region, the extension of helix D and the conformational behavior of the reactive center loop (RCL). We could also confirm the high stability of helix P in non-activated AT conformations, such states might play an important role in heparin binding. 'Generalized correlation' matrices revealed possible paths of allosteric signal propagation to the binding sites for the target proteases, factors Xa and IXa. Enhanced MD simulations of ligand binding to AT may assist the design of new anticoagulant drugs.Communicated by Ramaswamy H. Sarma.

The Observation of Ligand-Binding-Relevant Open States of Fatty Acid Binding Protein by Molecular Dynamics Simulations and a Markov State Model

As a member of the fatty acids transporter family, the heart fatty acid binding proteins (HFABPs) are responsible for many important biological activities. The binding mechanism of fatty acid with FABP is critical to the understanding of FABP functions. The uncovering of binding-relevant intermediate states and interactions would greatly increase our knowledge of the binding process. In this work, all-atom molecular dynamics (MD) simulations were performed to characterize the structural properties of nativelike intermediate states. Based on multiple 6 μs MD simulations and Markov state model (MSM) analysis, several "open" intermediate states were observed. The transition rates between these states and the native closed state are in good agreement with the experimental measurements, which indicates that these intermediate states are binding relevant. As a common property in the open states, the partially unfolded α2 helix generates a larger portal and provides the driving force to facilitate ligand binding. On the other side, there are two kinds of open states for the ligand-binding HFABP: one has the partially unfolded α2 helix, and the other has the looser β-barrel with disjointing βD-βE strands. Our results provide atomic-level descriptions of the binding-relevant intermediate states and could improve our understanding of the binding mechanism.

Importance of Incorporating Protein Flexibility in Molecule Modeling: A Theoretical Study on Type I1/2 NIK Inhibitors

NF-κB inducing kinase (NIK), which is considered as the central component of the non-canonical NF-κB pathway, has been proved to be an important target for the regulation of the immune system. In the past few years, NIK inhibitors with various scaffolds have been successively reported, among which type I^1/2 inhibitors that can not only bind in the ATP-binding pocket at the DFG-in state but also extend into an additional back pocket, make up the largest proportion of the NIK inhibitors, and are worthy of more attention. In this study, an integration protocol that combines molecule docking, MD simulations, ensemble docking, MM/GB(PB)SA binding free energy calculations, and decomposition was employed to understand the binding mechanism of 21 tricyclic type I^1/2 NIK inhibitors. It is found that the docking accuracy is largely dependent on the selection of docking protocols as well as the crystal structures. The predictions given by the ensemble docking based on multiple receptor conformations (MRCs) and the MM/GB(PB)SA calculations based on MD simulations showed higher linear correlations with the experimental data than those given by conventional rigid receptor docking (RRD) methods (Glide, GOLD, and Autodock Vina), highlighting the importance of incorporating protein flexibility in predicting protein–ligand interactions. Further analysis based on MM/GBSA demonstrates that the hydrophobic interactions play the most essential role in the ligand binding to NIK, and the polar interactions also make an important contribution to the NIK-ligand recognition. A deeper comparison of several pairs of representative derivatives reveals that the hydrophobic interactions are vitally important in the structural optimization of analogs as well. Besides, the H-bond interactions with some key residues and the large desolvation effect in the back pocket devote to the affinity distinction. It is expected that our study could provide valuable insights into the design of novel and potent type I^1/2 NIK inhibitors.

Mechanistic Insights into Binding of Ligands with Thiazolidinedione Warhead to Human Histone Deacetylase 4

Recently, we have reported that non-hydroxamate thiazolidinedione (TZD) analogs are capable of inhibiting human deacetylase 4 (HDAC4). This study aims at the dissection of the molecular determinants and kinetics of the molecular recognition of TZD ligands by HDAC4. For this purpose, a structure activity relationship analysis of 225 analogs was combined with a comprehensive study of the enzyme and binding kinetics of a variety of HDAC4 mutant variants. The experimental data were rationalized by docking to the two major conformations of HDAC4. TZD ligands are competitive inhibitors and bind via a two-step mechanism involving principal molecular recognition and induced fit. The residence time of 24 g is (34 ± 3) min and thus much larger than that of the canonical pan-HDAC inhibitor SAHA ((5 ± 2) min). Importantly, the binding kinetics can be tuned by varying the structure of the CAP group.

An Insight Based on Computational Analysis of the Interaction between the Receptor-Binding Domain of the Omicron Variants and Human Angiotensin-Converting Enzyme 2

Simple Summary

The Omicron variant has recently been divided into BA.1, BA.2, and BA.3 subvariants. In the present study, we focused on comparing the interaction between the receptor-binding domain (RBD) of BA.1 and BA.2 spike proteins with human angiotensin-converting enzyme 2 (hACE2) using a computational approach. The RBD BA.2 was modeled after the BA.1. The results from molecular docking and molecular dynamics studies showed that RBD BA.2 has a higher and more stable affinity for hACE2 compared to RBD BA.1.

Abstract

Concerns have been raised about the high number of mutations in the spike protein of the new emergence of the highly transmissible Omicron variant (B.1.1529 lineage) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This variant’s extraordinary ability to evade antibodies would significantly impair the current vaccination program. This present study aimed to computationally analyze the interaction between the receptor-binding domain (RBD) in the spike protein of Omicron variants and human angiotensin-converting enzyme 2 (hACE2). The docking results indicated that Omicron BA.2 has exceptionally strong interactions with hACE2 in comparison to Omicron BA.1, Delta, and wild-type, as indicated by various parameters such as salt bridge, hydrogen bond, and non-bonded interactions. The results of the molecular dynamics simulation study corroborate these findings, indicating that Omicron BA.2 has a strong and stable interaction with hACE2. This study provides insight into the development of an effective intervention against this variant.

The Effect of Carbon Nanotubes on the Mechanical Properties of Wood Plastic Composites by Selective Laser Sintering

Wood-plastic composites (WPCs) made by selective laser sintering (SLS) approach of 3D printing offer many advantages over single polymer materials, such as low cost, sustainability, and better sintering accuracy. However, WPCs made via SLS are too weak to have widespread applications. In order to increase the mechanical properties of WPCs, a novel type of WPCs containing 0, 0.05, 0.1 and 0.15 wt % carbon nanotubes (CNT), 14 wt % wood fibers, 86 wt % polyether sulfone (PES) was manufactured via SLS. The experimental results showed that the addition of small amount of CNTs can significantly increase the mechanical properties of the wood/PES composite material. The tensile strength, bending strength, and elasticity modulus were 76.3%, 227.9%, and 128.7% higher with 0.1 wt % CNTs than those without CNTs. The mechanical properties of specimens first increased and then decreased with the addition of CNTs. The SEM results of the specimens' fracture morphology indicate that the preferable bonding interfaces between wood flour grains and PES grains were achieved by adding CNTs to the composites. There are two reasons why the composites possessed superior mechanical properties: CNTs facilitate the laser sintering process of WPCs due to their thermal conductivities, and CNTs directly reinforce WPCs.

Non-histone binding functions of PHD fingers

Plant homeodomain (PHD) fingers comprise a large and well-established family of epigenetic readers that recognize histone H3. A typical PHD finger binds to the unmodified or methylated amino-terminal tail of H3. This interaction is highly specific and can be regulated by post-translational modifications (PTMs) in H3 and other domains present in the protein. However, a set of PHD fingers has recently been shown to bind non-histone proteins, H3 mimetics, and DNA. In this review, we highlight the molecular mechanisms by which PHD fingers interact with ligands other than the amino terminus of H3 and discuss similarities and differences in engagement with histone and non-histone binding partners.

Deciphering the promiscuous interactions between intrinsically disordered transactivation domains and the KIX domain

The kinase-inducible domain interacting (KIX) domain of the transcriptional coactivator CBP protein carries 2 isolated binding sites (designated as the c-Myb site and the MLL site) and is capable of binding numerous intrinsically disordered transactivation domains (TADs), including c-Myb and pKID via the c-Myb site, and MLL, E2A and c-Jun via the MLL site. In this study we compared the kinetics for binding of various disordered TADs to the KIX domain via computational biophysical analyses. We found that the binding rates are heavily affected by long-range electrostatic interactions. The basal rate constants for forming the encounter complexes are similar for different KIX binding peptides, favorable electrostatic interactions between the MLL site and the peptides result in greater association rates when peptides bind to the MLL site. FOXO3a and p53 TAD each contains 2 copies of KIX binding motif and each motif interacts with both the MLL site and the c-Myb site. Our kinetics studies suggest that binding of FOXO3a or p53 TAD to the KIX domain is via a sequential mechanism, where one KIX binding motif binds to the MLL site first and then the other KIX binding motif binds to the c-Myb site. Considering the promiscuous interactions between FOXO3a and KIX, and p53 TAD and KIX, electrostatic steering simplifies the binding mechanism. This study highlights the importance of long-range electrostatic interactions in molecular recognition process involving multi-motif intrinsically disordered proteins and promiscuous interactions.

Molecular mechanism of interspecies differences in the binding affinity of TD139 to Galectin-3

Galectin-3 (Gal-3), a β-galactoside-binding lectin, has been implicated in a plethora of pathological disorders including fibrosis, inflammation, cancer and metabolic diseases. TD139-a thio-digalactoside inhibitor developed by Galecto Biotech as a potential therapeutic for idiopathic pulmonary fibrosis-is the most advanced small-molecule Gal-3 inhibitor in clinical studies. It binds to human Gal-3 with high affinity but has lower affinity towards mouse and rat homologs, which is also manifested in the differential inhibition of Gal-3 function. Using biophysical methods and high-resolution X-ray co-crystal structures of TD139 and Gal-3 proteins, we demonstrate that a single amino acid change corresponding to A146 in human Gal-3 is sufficient for the observed reduction in the binding affinity of TD139 in rodents. Site-directed mutagenesis of A146V (in human Gal-3) and V160A (in mouse Gal-3) was sufficient to interchange the affinities, mainly by affecting the off rates of the inhibitor binding. In addition, molecular dynamics simulations of both wild-type and mutant structures revealed the sustained favorable noncovalent interactions between the fluorophenyl ring and the active site A146 (human Gal-3 and mouse V160A) that corroborate the finding from biophysical studies. Current findings have ramifications in the context of optimization of drug candidates against Gal-3.