Probing binding mechanism of interleukin-6 and olokizumab: in silico design of potential lead antibodies for autoimmune and inflammatory diseases

Computer-aided antibody engineering has been successful in the design of new biologics for disease diagnosis and therapeutic interventions. Interleukin-6 (IL-6), a well-recognized drug target for various autoimmune and inflammatory diseases such as rheumatoid arthritis, multiple sclerosis, and psoriasis, was investigated in silico to design potential lead antibodies. Here, crystal structure of IL-6 along with monoclonal antibody olokizumab was explored to predict antigen–antibody (Ag???Ab)-interacting residues using DiscoTope, Paratome, and PyMOL. Tyr56, Tyr103 in heavy chain and Gly30, Ile31 in light chain of olokizumab were mutated with residues Ser, Thr, Tyr, Trp, and Phe. A set of 899 mutant macromolecules were designed, and binding affinity of these macromolecules to IL-6 was evaluated through Ag???Ab docking (ZDOCK, ClusPro, and Rosetta server), binding free-energy calculations using Molecular Mechanics/Poisson Boltzman Surface Area (MM/PBSA) method, and interaction energy estimation. In comparison to olokizumab, eight newly designed theoretical antibodies demonstrated better result in all assessments. Therefore, these newly designed macromolecules were proposed as potential lead antibodies to serve as a therapeutics option for IL-6-mediated diseases.     

Dynamics and modulation studies of human voltage gated Kv1.5 channel

The voltage gated Kv1.5 channels conduct the ultrarapid delayed rectifier current (IK_ur) and play critical role in repolarization of action potential duration. It is the most rapidly activated channel and has very little or no inactivated states. In human cardiac cells, these channels are expressed more extensively in atrial myocytes than ventricle. From the evidences of its localization and functions, Kv1.5 has been declared a selective drug target for the treatment of atrial fibrillation (AF). In this present study, we have tried to identify the rapidly activating property of Kv1.5 and studied its mode of inhibition using molecular modeling, docking, and simulation techniques. Channel in open conformation is found to be stabilized quickly within the dipalmitoylphosphatidylcholine membrane, whereas most of the secondary structure elements were lost in closed state conformation. The obvious reason behind its ultra-rapid property is possibly due to the amino acid alteration in S4–S5 linker; the replacement of Lysine by Glutamine and vice versa. The popular published drugs as well as newly identified lead molecules were able to inhibit the Kv1.5 in a very similar pattern, mainly through the nonpolar interactions, and formed sable complexes. V512 is found as the main contributor for the interaction along with the other important residues such as V505, I508, A509, V512, P513, and V516. Furthermore, two screened novel compounds show surprisingly better inhibitory potency and can be considered for the future perspective of antiarrhythmic survey.     

Comparison of different methods of determining cell viability after exposure to cytotoxic compounds

Cell-survival (of DON and L1210 cells) after treatment with cytotoxic compounds was assessed by measuring cloning efficiency, exclusion of trypan blue and erythrosin B, [⁵¹Cr] release, and attachment of DON cells to glass. Cell survival as measured by cloning efficiency did not correlate with survival measured by any of the other methods. We found that the stainability of cells after drug exposure depended on the cell line used. For example, after 3 h exposure to tubercidin although 100% of both DON and L1210 cells were killed (on basis of cloning efficiency), only 11% of DON cells and 68% of L1210 cells were dead as indicated by staining with erythrosin B. The stainability of cells also depended on the particular drug used. For example, after 24 h exposure of L1210 cells to adriamycin and tubercidin (both killed >99% of cells on basis of cloning efficiency) 21% of cells exposed to adriamycin and 99% of cells exposed to tubercidin were stained. The results obtained with several other cytotoxic compounds are discussed.     

Trimethyloxonium modification of single batrachotoxin-activated sodium channels in planar bilayers. Changes in unit conductance and in block by saxitoxin and calcium

Single batrachotoxin-activated sodium channels from rat brain were modified by trimethyloxonium (TMO) after incorporation in planar lipid bilayers. TMO modification eliminated saxitoxin (STX) sensitivity, reduced the single channel conductance by 37%, and reduced calcium block of inward sodium currents. These effects always occurred concomitantly, in an all-or-none fashion. Calcium and STX protected sodium channels from TMO modification with potencies similar to their affinities for block. Calcium inhibited STX binding to rat brain membrane vesicles and relieved toxin block of channels in bilayers, apparently by competing with STX for the toxin binding site. These results suggest that toxins, permeant cations, and blocking cations can interact with a common site on the sodium channel near the extracellular surface. It is likely that permeant cations transiently bind to this superficial site, as the first of several steps in passing inward through the channel.     

Modelling family 2 cystatins and their interaction with papain

Cystatins are extensively studied cysteine protease inhibitors, found in wide range of organisms with highly conserved structural folds. S-type of cystatins is well known for their abundance in saliva, high selectivity and poorer activity towards host cysteine proteases in comparison to their immediate ancestor cystatin C. Despite more than 90% sequence similarity, the members of this group show highly dissimilar binding affinity towards papain. Cystatin M/E is a potent inhibitor of legumain and papain like cysteine proteases and recognized for its involvement in skin barrier formation and potential role as a tumor suppressor gene. However, the structures of these proteins and their complexes with papain or legumain are still unknown. In the present study, we have employed computational methods to get insight into the interactions between papain and cystatins. Three-dimensional structures of the cystatins are generated by homology modelling, refined with molecular dynamics simulation, validated through numerous web servers and finally complexed with papain using ZDOCK algorithm in Discovery Studio. A high degree of shape complementarity is observed within the complexes, stabilized by numerous hydrogen bonds (HB) and hydrophobic interactions. Using interaction energy, HB and solvent accessible surface area analyses, we have identified a series of key residues that may be involved in papain–cystatin interaction. Differential approaches of cystatins towards papain are also noticed which are possibly responsible for diverse inhibitory activity within the group. These findings will improve our understanding of fundamental inhibitory mechanisms of cystatin and provide clues for further research.     

Lens-specific regulation of the thioredoxin-1 gene, but not thioredoxin-2, upon in vivo photochemical oxidative stress in the Emory mouse

The angiotensin II (Ang II) type 1 receptor mediates various actions of Ang II, whereas the function of the type 2 (AT2) receptor is not well understood. In the mice lacking the gene encoding the AT2 receptor, the pressor response to Ang II was increased although the underlying mechanism is unknown. We tested the hypothesis that vasoconstrictor response is exaggerated in the AT2 receptor null mice. We measured hemodynamic parameters and evaluated systemic vascular resistance (SVR) in the anesthetized open-chest wild-type and AT2 receptor null mice. Ang II infusion caused dose-dependent increases in SVR in both strains, while the response was significantly higher at 0.5 microgram/kg Ang II in the AT2 receptor null mice (305 +/- 53% of baseline) than in the wild-type mice (179 +/- 27% of baseline). To investigate further the vascular contractility, we examined contraction of aortic rings in vitro. The contraction induced by 1 microM Ang II was increased in the AT2 receptor null mice compared with that in the wild-type mice (0.82 +/- 0.11 vs. 0.54 +/- 0.12 g). Ang II-induced contraction was still greater in the AT2 receptor null mice when calibrated by the maximum tension induced by 90 mM KCl. These data suggest that the AT2 receptor modulates vascular contractility, which may influence blood pressure.     

Backbone dynamics of free barnase and its complex with barstar determined by 15N NMR relaxation study

Backbone dynamics of uniformly ¹⁵N-labeled free barnase and its complex with unlabelled barstar have been studied at 40°C, pH 6.6, using ¹⁵N relaxation data obtained from proton-detected 2D {¹H}-¹⁵N NMR spectroscopy. ¹⁵N spin-lattice relaxation rate constants (R₁), spin-spin relaxation rate constants (R₂), and steady-state heteronuclear {¹H}-¹⁵N NOEs have been measured at a magnetic field strength of 14.1 Tesla for 91 residues of free barnase and for 90 residues out of a total of 106 in the complex (excluding three prolines and the N-terminal residue) backbone amide ¹⁵N sites of barnase. The primary relaxation data for both the cases have been analyzed in the framework of the model-free formalism using both isotropic and axially symmetric models of the rotational diffusion tensor. As per the latter, the overall rotational correlation times (_m) are 5.0 and 9.5 ns for the free and complexed barnase, respectively. The average order parameter is found to be 0.80 for free barnase and 0.86 for the complex. However, the changes are not uniform along the backbone and for about 5 residues near the binding interface there is actually a significant decrease in the order parameters on complex formation. These residues are not involved in the actual binding. For the residues where the order parameter increases, the magnitudes vary significantly. It is observed that the complex has much less internal mobility, compared to free barnase. From the changes in the order parameters, the entropic contribution of NH bond vector motion to the free energy of complex formation has been calculated. It is apparent that these motions cause significant unfavorable contributions and therefore must be compensated by many other favorable contributions to effect tight complex formation. The observed variations in the motion and their different locations with regard to the binding interface may have important implications for remote effects and regulation of the enzyme action.