Finding pharmaceutically relevant target conformations from an arbitrary set of protein conformations remains a challenge in structure-based virtual screening (SBVS). The growth in the number of available conformations, either experimentally determined or computationally derived, obscures the situation further. While the inflated conformation space potentially contains viable druggable targets, the increase of conformational complexity, as a consequence, poses a selection problem. To address this challenge, we took advantage of machine learning methods, namely an over-sampling and a binary classification procedure, and present a novel method to select druggable receptor conformations. Specifically, we trained a binary classifier on a set of nuclear receptor conformations, wherein each conformation was labeled with an enrichment measure for a corresponding SBVS. The classifier enabled us to formulate suggestions and identify enriching SBVS targets for six of seven nuclear receptors. Further, the classifier can be extended to other proteins of interest simply by feeding new training data sets to the classifier. Our work, thus, provides a methodology to identify pharmaceutically interesting receptor conformations for nuclear receptors and other drug targets.
Polymorphisms of the ADIPOR2 gene are frequently linked to a higher risk of developing diseases including obesity, type 2 diabetes and cardiovascular diseases. Though mutations of the ADIPOR2 gene are detrimental, there is a lack of comprehensive in silico analyses of the functional and structural impacts at the protein level. Considering the involvement of ADIPOR2 in glucose uptake and fatty acid oxidation, an in silico functional analysis was conducted to explore the possible association between genetic mutations and phenotypic variations. A genomic analysis of 82 nonsynonymous SNPs in ADIPOR2 was initiated using SIFT followed by the SNAP2, nsSNPAnalyzer, PolyPhen-2, SNPs&GO, FATHMM and PROVEAN servers. A total of 10 mutations (R126W, L160Q, L195P, F201S, L235R, L235P, L256R, Y328H, E334K and Q349H) were predicted to have deleterious effects on the ADIPOR2 protein and were therefore selected for further analysis. Theoretical models of the variants were generated by comparative modeling via MODELLER 9.16. A protein structural analysis of these amino acid variants was performed using SNPeffect, I-Mutant, ConSurf, Swiss-PDB Viewer and NetSurfP to explore their solvent accessibility, molecular dynamics and energy minimization calculations. In addition, FTSite was used to predict the ligand binding sites, while NetGlycate, NetPhos2.0, UbPerd and SUMOplot were used to predict post-translational modification sites. All of the variants showed increased free energy, though F201S exhibited the highest energy increase. The root mean square deviation values of the modeled mutants strongly indicated likely pathogenicity. Remarkably, three binding sites were detected on ADIPOR2, and two mutations at positions 328 and 201 were found in the first and second binding pockets, respectively. Interestingly, no mutations were found at the post-translational modification sites. These genetic variants can provide a better understanding of the wide range of disease susceptibility associated with ADIPOR2 and aid the development of new molecular diagnostic markers for these diseases. The findings may also facilitate the development of novel therapeutic elements for associated diseases.
In vitro Lead (Pb(2+)) binding capacity of two probiotic bacteria strains, namely Bifidobacterium longumBB79 and Lactobacillus pentosusITA23, was assessed following incubation with the intestinal contents (IC) of laying hens. Results of this study demonstrated that IC treatment significantly enhanced (P<0.01) Pb(2+) binding capacity of both bacterial strains. Fourier transform infrared analysis indicated that several functional groups (O-H or N-H, C-H, C˭O, C-O, and C-O-C) on the bacteria cell wall involved in metal ion binding were altered after IC incubation, and new groups appeared between the 3700cm(-1) and 4000cm(-1)bands. Transmission electron microscopy demonstrated that after incubation with IC, unidentified IC components created new binding sites on the bacterial cell surface. These particles also changed the mechanism of Pb(2+) binding of the two strains from intracellular accumulation to extracellular adsorption.
To characterize the binding of a widely used herbicide, 2,4-dichlorophenoxyacetic acid (2,4-D) to the major transporter in human circulation, human serum albumin (HSA), multi-spectroscopic approaches such as fluorescence, absorption and circular dichroism along with computational methods were employed. Analysis of the fluorescence and absorption spectroscopic data confirmed the 2,4-D-HSA complex formation. A static quenching mechanism was evident from the inverse temperature dependence of the KSV values. The complex was stabilized by a weak binding affinity (Ka = 5.08 × 103 M-1 at 298 K). Quantitative analysis of thermodynamic data revealed participation of hydrophobic and van der Waals interactions as well as hydrogen bonds in the binding process. Circular dichroism and three-dimensional fluorescence spectral results showed structural (secondary and tertiary) changes in HSA as well as microenvironmental perturbation around protein fluorophores (Trp and Tyr residues) upon 2,4-D binding. Addition of 2,4-D to HSA was found to improve protein's thermal stability. Competitive displacement results as well as computational analyses suggested preferred location of the 2,4-D binding site as Sudlow's site I (subdomain IIA) in HSA.
Protease is an enzyme that catalysed the hydrolysis of protein into peptide. Application of protease in industry has been linked with cost effective substrates and complex of enzyme-substrate stability. Molecular docking approach has identified casein as a preference substrates. However, lack of data on casein mode of binding to protease and enzyme stability represents a limitation for its production and structural optimization. In this study, we have used a molecular dynamic (MD) to examine the stability of complex enzyme-substrate of protease from Bacillus lehensis G1. The 3D structure of protease (BleG1_1979) was docked with substrate casein using AutoDock Vina. Structural analysis of the substrate-binding cleft revealed a binding site of casein was predominantly at the hydrophobic region of BleG1_1979. The MD of complex BleG1_1979-casein was tested with two temperatures; 298 K and 310 K using GROMACS v5.1.4. MD simulation showed a stable behaviour of BleG1_1979 over the 20 ns simulation period. The molecular docking and MD simulation suggested that the production of protease from B. lehensis G1 by utilization of casein and the stability of complex protease-casein could be a potential application to generate a cost effective enzyme to be develop for industrial use.
The presence of at least two ionizable active center groups has been detected by a study of the effect of pH upon catalysis of hydrolysis of L-alanyl-beta-naphthylamide by human liver alanine aminopeptidase and upon the inhibition of hydrolysis by inhibitors and substrate analogs. Octanoic acid, octylamine, and peptide inhibitors have been found to be competitive inhibitors and are therefore thought to bind the active center. L-Phe was previously shown to bind the active center since it was found to be a competitive inhibitor of the hydrolysis of tripeptide substrates (Garner, C. W., and Behal, F. J. (1975), Biochemistry 14, 3208). A plot of pKm vs. pH for the substrate L-Ala-beta-naphthylamide showed that binding decreased below pH 5.9 and above 7.5, the points at which the theoretical curve undergoes an integral change in slope. These points are interpreted as the pKa either of substrate ionizable groups or binding-dependent enzyme active center groups. Similar plots of pKm vs. pH for L-alanyl-p-nitroanilide (as substrate) and pKi vs. pH for L-Leu-L-Leu-L-Leu and D-Leu-L-Tyr (as inhibitors) gave pairs fo pKa values of 5.8 and 7.4, 6.0 and 7.5, and 5.7 and 7.5, respectively. All the above substrates (and D-Leu-L-Tyr) have pKa values near 7.5; therefore, the binding-dependent group with a pKa value near 7.5 is possibly this substrate group. Similar plots of pKi vs. pH for the inhibitors L-Phe, L-Met, L-Leu, octylamine, and octanoic acid had only one bending point at 7.7, 7.6, 7.4, 6.3, and 5.9, respectively. Amino acid inhibitors, octylamine, and octanoic acid have no groups with pKa values between 5 and 9. These data indicate that there are two active center ionizable groups with pKa values of approximately 6.0 and 7.5 which are involved in substrate binding or inhibitory amino acid binding but not in catalysis since Vmax was constant at all pH values tested.
Aripiprazole (ARP), a quinolinone derivative, is an atypical antipsychotic drug that is used in the treatment of schizophrenia. ARP has an extensive distribution and more than 99% of the ARP and dehydro-ARP, the main active metabolite, is bound to plasma proteins. However, information regarding the protein binding of ARP is limited. In this study, we report on a systematic study of the protein binding of ARP. The interaction of ARP and structurally related compounds with human serum albumin (HSA) was examined using equilibrium dialysis, circular dichroism (CD) spectroscopy, fluorescent probe displacement, and an X-ray crystallographic analysis. The binding affinities (nK) for ARP and its main metabolite, dehydro-ARP with HSA were found to be significantly higher than other structurally related compounds. The results of equilibrium dialysis experiments and CD spectral data indicated that the chloro-group linked to the phenylpiperazine ring in the ARP molecule plays a major role in the binding of these ligands to HSA. Furthermore, fluorescent probe displacement results indicated that ARP appears to bind at the site II pocket in subdomain III. A detailed CD spectral analysis suggests that the chloro-group linked to the phenylpiperazine ring may control the geometry of the ARP molecule when binding in the site II binding pocket. X-ray crystallographic analysis of the ARP-HSA complex revealed that the distance between the chlorine atom at the 3-positon of dichlorophenyl-piperazine on ARP and the sulfur atom of Cys392 in HSA was 3.4-3.6 Å. A similar halogen bond interaction has also been observed in the HSA structure complexed with diazepam, which also contains a chloro-group. Thus, the mechanism responsible for the binding of ARP to a protein elucidated here should be relevant for assessing the pharmacokinetics and pharmacodynamics of ARP in various clinical situations and for designing new drugs.
Liquefactions of kenaf core wood were carried out at different phenol-kenaf (P/k) ratios. Characterizations of kenaf core wood liquefied residue were carried out to measure the degree of liquefaction. This provides a new approach to understand some fundamental aspects of the liquefaction reaction. Functional groups on the raw kenaf core wood and liquefied residue were examined using Fourier transform infrared spectroscopy (FTIR). The crystallinity index of the kenaf wood liquefied residue, which represents crystallinity changes of the cellulose component after the liquefaction process, was studied using X-ray diffraction (XRD). The surface morphology of the wood residue was observed using scanning electron microscopy (SEM). The thermal behavior of the residues was analyzed using thermogravimetric analysis (TGA). Abroad peak around 3450-3400 cm-1 representing OH stretching in lignin start to disappear as P/K ratio increases. The results showed that the higher the P/K ratio the greater the liquefaction of the lignin component in the kenaf core wood. The crystallinity index (CrI) on the kenaf liquefied residues increased with the increase in P/K ratio. SEM images showed that the small fragments attached on the liquefied kenaf residue surface were gradually removed as the P/K ratio was increased from 1.5/1.0 to 2.5/1.0, which is mainly attributed to the greater chemical penetration toward reactive site of the kenaf fibres. Residue content decreased as the P/K ratio increased from 1.5/1.0 to 2.5/1.0. TGA results showed the increase of heat resistance in the residue as the P/K ratio was increased.
Despite the huge loss of lives and massive disruption of the world economy by the COVID -19 pandemic caused by SARS -CoV-2, scientists are yet to come out with an effective therapeutic against this viral disease . Several vaccines have obtained 'emergency approval ', but difficulties are being faced in the even distribution of vaccines amongst high- and low- income countries . On top of it, comorbidities associated with COVID -19 like diabetes, hypertension and malaria can seriously impede the treatment of the main disease, thus increasing the fatality rate . This is more so in the context of sub -Saharan African and south Asian countries . Our objective was to demonstrate that a single plant containing different phytoconstituents may be used for treatment of COVID -19 and comorbidities . Towards initial selection of a plant, existing scientific literature was scanned for reported relevant traditional uses, phytochemicals and pharmacological activities of a number of plants and their phytoconstituents pertaining to treatment of COVID-19 symptoms and comorbidities. Molecular docking studies were then performed with phytochemicals of the selected plant and SARS-CoV-2 components - Mpro, and spike protein receptor binding domain and hACE2 interface using AutoDock V ina. We showed that crude extracts of an indigenous African plant, Costus afer having traditional antidiabetic and antimalarial uses, has phytochemicals with high binding affinities for Mpro, and /or spike protein receptor binding domain and hACE2 interface; the various phytochemicals with predicted high binding energies include aferoside C, dibutyl phthalate, nerolidol, suginal, and ± -terpinene, making them potential therapeutics for COVID -19. The results suggest that crude extracts and phytochemicals of C. afer can function as a treatment modality for COVID -19 and comorbidities like especially diabetes and malaria .
Broadly neutralizing antibodies (bnAbs) have been considered to be potent therapeutic tools and potential vaccine candidates to enable protection against various clades of human immunodeficiency virus (HIV). The generation of bnAbs has been associated with enhanced exposure to antigen, high viral load and low CD4+ T cell counts, among other factors. However, only limited data are available on the generation of bnAbs in viraemic non-progressors that demonstrate moderate to high viraemia. Further, since HIV-1 subtype C viruses account for more than 50 % of global HIV infections, the identification of bnAbs with novel specificities is crucial to enable the development of potent tools to aid in HIV therapy and prevention. In the present study, we analysed and compared the neutralization potential of responses in 70 plasma samples isolated from ART-naïve HIV-1 subtype C-infected individuals with various disease progression profiles against a panel of 30 pseudoviruses. Among the seven samples that exhibited a neutralization breadth of ≥70 %, four were identified as 'elite neutralizers', and three of these were from viraemic non-progressors while the fourth was from a typical progressor. Analysis of the neutralization specificities revealed that none of the four elite neutralizers were reactive to epitopes in the membrane proximal external region (MPER), CD4-binding site and V1V2 or V3 glycan. However, two of the four elite neutralizers exhibited enhanced sensitivity towards viruses lacking N332 glycan, indicating high neutralization potency. Overall, our findings indicate that the identification of potent neutralization responses with distinct epitope specificities is possible from the as yet unexplored Indian population, which has a high prevalence of HIV-1 subtype C infection.
G protein-coupled receptors (GPCRs) are proteins of pharmaceutical importance, with over 30% of all drugs in clinical use targeting them. Increasing numbers of X-ray crystal (XRC) structures of GPCRs offer a wealth of data relating to ligand binding. For the β-adrenoceptors (β-ARs), XRC structures are available for human β2- and turkey β1-subtypes, in complexes with a range of ligands. While these structures provide insight into the origins of ligand structure-activity relationships (SARs), questions remain. The ligands in all published complexed XRC structures lack extensive substitution, with no obvious way the ligand-binding site can accommodate β1-AR-selective antagonists with extended side-chains para- to the common aryloxypropanolamine pharmacophore. Using standard computational docking tools with such ligands generally returns poses that fail to explain known SARs. Application of our Active Site Pressurization (ASP) modelling method to β-AR XRC structures and homology models however, reveals a dynamic area in the ligand-binding pocket that, through minor changes in amino acid side chain orientations, opens a fissure between transmembrane (TM) helices H4 and H5, exposing intra-membrane space. This fissure, which we term the 'keyhole', is ideally located to accommodate extended moieties present in many high-affinity β1-AR-selective ligands; allowing the rest of the ligand structure to adopt a canonical pose in the orthosteric binding site. We propose the keyhole may be a feature of both β1- and β2-ARs, but that subtle structural differences exist between the two, contributing to subtype-selectivity. This has consequences for the rational design of future generations of subtype-selective ligands for these therapeutically important targets.
Pazopanib (PZP) is a multi-targeting tyrosine kinase inhibitor and is currently approved by FDA for the treatment of soft tissue sarcoma and renal cancer. Molecular interaction mechanism of PZP with human serum albumin (HSA) was explored under simulated physiological conditions (pH = 7.4), using fluorescence and UV absorption spectroscopy along with computational methods. Based on the inverse correlation between the Stern-Volmer constant (Ksv) and temperature, it was concluded that PZP quenched the protein fluorescence through static quenching mechanism. This was also confirmed from the UV-vis absorption spectral results. Moderate binding affinity between PZP and HSA was evident from the Ka values (5.51 - 1.05 × 105 M-1) while PZP-HSA complex formation was driven by hydrophobic and van der Waals interactions as well as hydrogen bonds, as revealed by positive entropy change (ΔS = +98.37 J mol-1 K-1) and negative enthalpy change (ΔH = -60.31 kJ mol-1). Three-dimensional fluorescence spectral results disclosed microenvironmental perturbations around Trp and Tyr residues of the protein upon PZP binding. Interestingly, the addition of PZP to HSA significantly protected the protein against thermal stress. Competitive drug displacement results obtained with warfarin, phenylbutazone and diazepam elucidated Sudlow's Site I, positioned in subdomain IIA of HSA, as the preferred binding site of PZP which was well supported by molecular docking analysis, while molecular dynamics simulation results suggested the stability of the PZP-HSA complex.Communicated by Vsevolod Makeev.
We earlier showed that Torpedo californica acetylcholinesterase (AChE) contains a cluster of four conserved aspartates that can strongly bind divalent cations, which we named the 4D motif. Binding of the divalent metal cations greatly increases its thermal stability. Here we systematically examined all available crystallographic structures of T. californica AChE. Two additional metal-binding sites were identified, both composed of acidic and histidine residues. Relative binding to the 4D and additional sites was studied using metadynamics simulations. It was observed that in crystal structures devoid of metal ions in the 4D site, the conformation of T. californica AChE is almost identical to that in structures in which it is occupied by a divalent metal ion. Closer examination of the 4D motif reveals that three of the four acidic residues form ion pairs with conserved basic residues surrounding them. We named this new motif the 4A/3B motif. Molecular dynamics with quantum potential simulations was used to quantify the 4D motif's binding strength compared with that of the metal-binding site in the protein fXIIIa, which consists of four aspartates, but is devoid of adjacent cationic residues. Whereas fXIIIa's 4D site, in the absence of a metal cation, expanded significantly in the simulation, that of Torpedo AChE displayed only minor periodic changes in size. Furthermore, the energy of metal ion unbinding from the two sites differs by ca. 10 kcal/mol. We identified several other proteins in the PDB that contain the 4A/3B motif, whose conformations are identical in the presence or absence of a metal ion. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at https://proteopedia.org/w/Journal:Protein_Science:4.
beta-D-N-Acetylhexosaminidase, a family 20 glycosyl hydrolase, catalyzes the removal of β-1,4-linked N-acetylhexosamine residues from oligosaccharides and their conjugates. We constructed phylogenetic tree of β-hexosaminidases to analyze the evolutionary history and predicted functions of plant hexosaminidases. Phylogenetic analysis reveals the complex history of evolution of plant β-hexosaminidase that can be described by gene duplication events. The 3D structure of tomato β-hexosaminidase (β-Hex-Sl) was predicted by homology modeling using 1now as a template. Structural conformity studies of the best fit model showed that more than 98% of the residues lie inside the favoured and allowed regions where only 0.9% lie in the unfavourable region. Predicted 3D structure contains 531 amino acids residues with glycosyl hydrolase20b domain-I and glycosyl hydrolase20 superfamily domain-II including the (β/α)8 barrel in the central part. The α and β contents of the modeled structure were found to be 33.3% and 12.2%, respectively. Eleven amino acids were found to be involved in ligand-binding site; Asp(330) and Glu(331) could play important roles in enzyme-catalyzed reactions. The predicted model provides a structural framework that can act as a guide to develop a hypothesis for β-Hex-Sl mutagenesis experiments for exploring the functions of this class of enzymes in plant kingdom.
Filamentous phage random peptide libraries were used to identify the epitopes of Burkholderia pseudomallei protease by panning against IgG polyclonal sera that exhibited protease neutralizing properties. The isolated fusion peptides presented a consensus peptide sequence, TKSMALSG, which closely resembles part of the active site sequence, 435GTSMATPHVAG445, of B. pseudomallei serine metalloprotease. By comparing the consensus sequence, TKSMALSG, with the predicted three-dimensional molecular model of B. pseudomallei serine metalloprotease, it appears that the potential antibody binding epitope was buried within the molecule. This active site was conformational whereby one continuous sub-region (SMA) was located between two discontinuous sub-regions, supplied by the flanking residues in the same polypeptide. All phages selected from the biopanning with IgG polyclonal sera showed good binding towards the polyclonal antibodies when compared to the negative control. In addition, these peptide-bearing phages showed competitive inhibition of B. pseudomallei serine metalloprotease binding to the polyclonal IgG.
The renin-angiotensin system (RAS) plays an important role in the pathophysiology of cardiovascular disorders. Pharmacologic interventions targeting the RAS cascade have led to the discovery of renin inhibitors, angiotensin-converting enzyme inhibitors, and AT(1) receptor blockers (ARBs) to treat hypertension and some cardiovascular and renal disorders. Mutagenesis and modeling studies have revealed that differential functional outcomes are the results of multiple active states conformed by the AT(1) receptor upon interaction with angiotensin II (Ang II). The binding of agonist is dependent on both extracellular and intramembrane regions of the receptor molecule, and as a consequence occupies more extensive area of the receptor than a non-peptide antagonist. Both agonist and antagonist bind to the same intramembrane regions to interfere with each other's binding to exhibit competitive, surmountable interaction. The nature of interactions with the amino acids in the receptor is different for each of the ARBs given the small differences in the molecular structure between drugs. AT(1) receptors attain different conformation states after binding various Ang II analogues, resulting in variable responses through activation of multiple signaling pathways. These include both classical and non-classical pathways mediated through growth factor receptor transactivations, and provide cross-communication between downstream signaling molecules. The structural requirements for AT(1) receptors to activate extracellular signal-regulated kinases 1 and 2 through G proteins, or G protein-independently through β-arrestin, are different. We review the structural and functional characteristics of Ang II and its analogs and antagonists, and their interaction with amino acid residues in the AT(1) receptor.
Background: Andrographolide and its benzylidene derivatives, SRJ09 and SRJ23, potentially bind oncogenic K-Ras to exert anticancer activity. Their molecular interactions with K-Ras oncoproteins that lead to effective biological activity are of major interest. Methods & results: In silico docking and molecular dynamics simulation were performed using Glide and Desmond, respectively; while saturation transfer difference NMR was performed using GDP-bound K-RasG12V. SRJ23 was found to bind strongly and selectively to K-RasG12V, by anchoring to a binding pocket (namely p2) principally via hydrogen bond and hydrophobic interactions. The saturation transfer difference NMR analysis revealed the proximity of protons of functional moieties in SRJ23 to K-RasG12V, suggesting positive binding. Conclusion: SRJ23 binds strongly and interacts stably with K-RasG12V to exhibit its inhibitory activity.
Protein engineering is a very useful tool for probing structure-function relationships in proteins. Specifically, site-directed mutagenized proteins can provide useful insights into structural, binding and catalytic mechanisms of a protein, particularly when coupled with crystallization. In this chapter, we describe two protocols for performing site-directed mutagenesis of any protein-coding sequence, namely, megaprimer PCR and overlapping extension PCR (OE-PCR). We use as an example how these two SDM methods enhanced the function of a cyclodextrin glucosyltransferase (CGTase) from Bacillus lehensis strain G1.
Sixteen novel coumarin-based compounds are reported as potent acetylcholinesterase (AChE) inhibitors. The most active compound in this series, 5a (IC50 0.04 ± 0.01 µM), noncompetitively inhibited AChE with a higher potency than tacrine and galantamine. Compounds 5d, 5j, and 5 m showed a moderate antilipid peroxidation activity. The compounds showed cytotoxicity in the same range as the standard drugs in HEK-293 cells. Molecular docking demonstrated that 5a acted as a dual binding site inhibitor. The coumarin moiety occupied the peripheral anionic site and showed π-π interaction with Trp278. The tertiary amino group displayed significant cation-π interaction with Phe329. The aromatic group showed π-π interaction with Trp83 at the catalytic anionic site. The long chain of methylene lay along the gorge interacting with Phe330 via hydrophobic interaction. Molecular docking was applied to postulate the selectivity toward AChE of 5a in comparison with donepezil and tacrine. Structural insights into the selectivity of the coumarin derivatives toward huAChE were explored by molecular docking and 3D QSAR and molecular dynamics simulation for 20 ns. ADMET analysis suggested that the 2-(2-oxo-2H-chromen-4-yl)acetamides showed a good pharmacokinetic profile and no hepatotoxicity. These coumarin derivatives showed high potential for further development as anti-Alzheimer agents.