Role of CSF1R 550th-tryptophan in kusunokinin and CSF1R inhibitor binding and ligand-induced structural effect

Binding affinity is an important factor in drug design to improve drug-target selectivity and specificity. In this study, in silico techniques based on molecular docking followed by molecular dynamics (MD) simulations were utilized to identify the key residue(s) for CSF1R binding affinity among 14 pan-tyrosine kinase inhibitors and 15 CSF1R-specific inhibitors. We found tryptophan at position 550 (W550) on the CSF1R binding site interacted with the inhibitors' aromatic ring in a π-π way that made the ligands better at binding. Upon W550-Alanine substitution (W550A), the binding affinity of trans-(-)-kusunokinin and imatinib to CSF1R was significantly decreased. However, in terms of structural features, W550 did not significantly affect overall CSF1R structure, but provided destabilizing effect upon mutation. The W550A also did not either cause ligand to change its binding site or conformational changes due to ligand binding. As a result of our findings, the π-π interaction with W550's aromatic ring could be still the choice for increasing binding affinity to CSF1R. Nevertheless, our study showed that the increasing binding to W550 of the design ligand may not ensure CSF1R specificity and inhibition since W550-ligand bound state did not induce significantly conformational change into inactive state.

Conformational features and interaction mechanisms of VH H antibodies with β-hairpin CDR3: A case of Nb8-HigB2 interaction

The β-hairpin conformation is regarded as an important basic motif to form and regulate protein-protein interactions. Single-domain VH H antibodies are potential therapeutic and diagnostic tools, and the third complementarity-determining regions of the heavy chains (CDR3s) of these antibodies are critical for antigen recognition. Although the sequences and conformations of the CDR3s are diverse, CDR3s sometimes adopt β-hairpin conformations. However, characteristic features and interaction mechanisms of β-hairpin CDR3s remain to be fully elucidated. In this study, we investigated the molecular recognition of the anti-HigB2 VH H antibody Nb8, which has a CDR3 that forms a β-hairpin conformation. The interaction was analyzed by evaluation of alanine-scanning mutants, molecular dynamics simulations, and hydrogen/deuterium exchange mass spectrometry. These experiments demonstrated that positions 93 and 94 (Chothia numbering) in framework region 3, which is right outside CDR3 by definition, play pivotal roles in maintaining structural stability and binding properties of Nb8. These findings will facilitate the design and optimization of single-domain antibodies.

Degradation of ochratoxins A and B by lipases: A kinetic study unraveled by molecular modeling

Mycotoxins are toxic substances produced by fungi and, frequently, different mycotoxins cooccur in food commodities. Ochratoxin A (OTA) and Ochratoxin B (OTB) may co-occur in a variety of foods, like red wines and wheat, presenting a significant risk of population exposure. In this study, we investigated the potential of five lipases (Candida rugosa Lipase, Candida antarctica B Lipase, Thermomyces lanuginosus Lipase, Amano Lipase A from Aspergillus niger (ANL) and Porcine Pancreas Lipase (PPL)) to hydrolyze OTA and OTB into non-hazardous products. Only ANL and PPL degraded both substrates, however, with varying degrees of efficiency. PPL completely degraded OTB (9 h), but only 43% of OTA (25 h). Molecular simulations indicated a high binding energy of OTA to PPL, that can be explained by the impact of the chlorine group, impairing hydrolysis. ANL was able to completely degrade both mycotoxins, OTA in 3 h and OTB in 10 h. The ANL enzyme showed also high specificity to OTA, however, the activity of this enzyme is not affected by chlorine and hydrolyzes OTA faster than OTB. These two enzymes were found to be able to detoxify co-occurring ochratoxins A and B, making isolated enzymes an alternative to the direct use of microorganisms for mycotoxin mitigation in food.

Investigations on the role of cation-π interactions in active centers of superoxide dismutase

In this study, we have analysed the influence of cation-? interactions on stability and properties of superoxide dismutase (SOD) active centres. The number of interactions formed by arginine is higher than lysine in the cationic group, while histidine is comparatively higher in the ? group. The energy contribution resulting from most frequent cation-? interactions was in the lower range of strong hydrogen bonds. The cation-? interactions involving transition metal ions as cation have energy more negative than -418.4 kJ mol-1. Stabilization centres for these proteins showed that all residues involved in cation-? interactions were important in locating one or more of such centres. The majority of the residues involved in cation-p interactions were evolutionarily conserved and might have a significant contribution towards the stability of SOD proteins. The results presented in this work can be very useful for understanding the contribution of cation-? interactions to the stability of SOD active canters.

Structural complexes of the agonist, inverse agonist and antagonist bound C5a receptor: insights into pharmacology and signaling

The C5a receptor (C5aR) is a pharmacologically important G-protein coupled receptor (GPCR) that interacts with (h)C5a, by recruiting both the "orthosteric" sites (site1 at the N-terminus and site2 at the ECS, extra cellular surface) on C5aR in a two site-binding model. However, the complex pharmacological landscape and the distinguishing chemistry operating either at the "orthosteric" site1 or at the functionally important "orthosteric" site2 of C5aR are still not clear, which greatly limits the understanding of C5aR pharmacology. One of the major bottlenecks is the lack of an experimental structure or a refined model structure of C5aR with appropriately defined active sites. The study attempts to understand the pharmacology at the "orthosteric" site2 of C5aR rationally by generating a highly refined full-blown model structure of C5aR through advanced molecular modeling techniques, and further subjecting it to automated docking and molecular dynamics (MD) studies in the POPC bilayer. The first series of structural complexes of C5aR respectively bound to a linear native peptide agonist ((h)C5a-CT), a small molecule inverse agonist (NDT) and a cyclic peptide antagonist (PMX53) are reported, apparently establishing the unique pharmacological landscape of the "orthosteric" site2, which also illustrates an energetically distinct but coherent competitive chemistry ("cation-π" vs. "π-π" interactions) involved in distinguishing the established ligands known for targeting the "orthosteric" site2 of C5aR. Over a total of 1 μs molecular dynamics (MD) simulation in the POPC bilayer, it is evidenced that while the agonist prefers a "cation-π" interaction, the inverse agonist prefers a "cogwheel/L-shaped" interaction in contrast to the "edge-to-face/T-shaped" type π-π interactions demonstrated by the antagonist by engaging the F275(7.28) of the C5aR. In the absence of a NMR or crystallographically guided model structure of C5aR, the computational model complexes not only provide valuable insights for understanding the C5aR pharmacology, but also emerge as a promising platform for the design and discovery of future potential drug candidates targeting the (h)C5a-C5aR signaling axes.

Thermodynamic entropy of organic oxidation in the water environment: experimental evaluation compared to semi-empirical calculation

Residual organic matters in the secondary effluent are usually less biodegradable in terms of the total organic carbon content, and when discharged into a receiving water body, their further decomposition most likely mainly occurs due to chemical oxidation. Using this scenario, a semi-empirical method was previously developed to calculate the thermodynamic entropy of organic oxidation to quantitatively evaluate the impact of organic discharge on the water environment. In this study, the relationship between the entropy increase (ΔS_C) and excess organic mass (ΔTOC) was experimentally verified via combustion heat measurement using typical organic chemicals and mixtures. For individual organic chemicals, a linear relationship was detected between ΔS_C and ΔTOC with the same proportionality coefficient, 54.0 kJ/g, determined in the previous semi-empirical relationship. For the organic mixtures, a linear relationship was also identified; however, the proportionality coefficient was 69.2 kJ/g, indicating an approximately 28 % increase in the oxidation heat required to decompose the same organic mass. This increase in energy can likely be attributed to the synergistic effects of hydrogen bonding, hydrophobic interactions, π-π interactions, and van der Waals interactions between functional groups of different organic compounds. Intermolecular interactions may result in 17-32 % more dissociation energy for organic mixtures compared to the organic components' chemical structures. Because organics discharged into a water body are always a mixture of organic compounds, the proportionality coefficient obtained using organic mixtures should be adopted to modify the previously proposed semi-empirical equation.

Contribution of cation-π interactions to the stability of Sm/LSm oligomeric assemblies

In this work, we have analyzed the influence of cation-π interactions to the stability of Sm/LSm assemblies and their environmental preferences. The number of interactions formed by arginine is higher than lysine in the cationic group, while histidine is comparatively higher than phenylalanine and tyrosine in the π group. Arg-Tyr interactions are predominant among the various pairs analyzed. The furcation level of multiple cation-π interactions is much higher than that of single cation-π interactions in Sm/LSm interfaces. We have found hot spot residues forming cation-π interactions, and hot spot composition is similar for all aromatic residues. The Arg-Phe pair has the strongest interaction energy of -8.81 kcal mol(-1) among all the possible pairs of amino acids. The extent of burial of the residue side-chain correlates with the ΔΔG of binding for residues in the core and also for hot spot residues cation-π bonded across the interface. Secondary structure of the cation-π residues shows that Arg and Lys preferred to be in strand. Among the π residues, His prefers to be in helix, Phe prefers to be in turn, and Tyr prefers to be in strand. Stabilization centers for these proteins showed that all the five residues found in cation-π interactions are important in locating one or more of such centers. More than 50 % of the cation-π interacting residues are highly conserved. It is likely that the cation-π interactions contribute significantly to the overall stability of Sm/LSm proteins.

Molecular dynamics simulations of wild type and mutants of botulinum neurotoxin A complexed with synaptic vesicle protein 2C

Using molecular dynamics simulations, we investigate the relationship between the conformational changes of BoNT/A-RBD:SV2C-LD and the interfacial interactions.

Thermodynamic, energetic, and topological properties of crystal packing of pyrazolo[1,5-a]pyrimidines governed by weak electrostatic intermolecular interactions

The relationship between energetic and topological properties of crystals with weak electrostatic intermolecular interactions and thermodynamic data are presented.

Anion–π interactions in protein–porphyrin complexes

In this work, we have analyzed the influence of anion–π interactions on the stability of high resolution protein–porphyrin complex crystal structures.

Contribution of anion-π interactions to the stability of Sm/LSm proteins

We have analyzed the influence of anion-π interactions to the stability of Sm/LSm assemblies. The side chain of Glu is more likely to be in anion-π interactions than Asp. Phe has the highest occurrence in these interactions than the other two π residues. Among the anion-π residue pairs, Glu-Phe residue pair showed the maximum number of anion-π. We have found hot-spot residues forming anion-π interactions, and Glu-Phe is the most common hot-spot interacting pair. The significant numbers of anion-π interacting residues identified in the dataset were involved in the formation of multiple anion-π interactions. More than half of the residues involved in these interactions are evolutionarily conserved. The anion-π interaction energies are distance and orientation dependent. It was found that anion-π interactions showed energy less than -15 kcal mol(-1), and most of them have energy in the range -2 to -9 kcal mol(-1). Solvent accessibility pattern of Sm/LSm proteins reveals that all of the interacting residues are preferred to be in buried regions. Most of the interacting residues preferred to be in strand. A significant percentage of anion-π interacting residues are located as stabilization centers and thus might provide additional stability to these proteins. The simultaneous interaction of anions and cations on different faces of the same π-system has been observed. On the whole, the results presented in this work will be very useful for understanding the contribution of anion-π interaction to the stability of Sm/LSm proteins.