PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “molecular docking”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33Linked to original sources

Modeled ligand-protein complexes elucidate the origin of substrate specificity and provide insight into catalytic mechanisms of phenylalanine hydroxylase and tyrosine hydroxylase.

NMR spectroscopy and X-ray crystallography have provided important insight into structural features of phenylalanine hydroxylase (PAH) and tyrosine hydroxylase (TH). Nevertheless, significant problems such as the substrate specificity of PAH and the different susceptibility of TH to feedback inhibition by l-3,4-dihydroxyphenylalanine (l-DOPA) compared with dopamine (DA) remain unresolved. Based on the crystal structures 5pah for PAH and 2toh for TH (Protein Data Bank), we have used molecular docking to model the binding of 6(R)-l-erythro-5,6,7,8-tetrahydrobiopterin (BH4) and the substrates phenylalanine and tyrosine to the catalytic domains of PAH and TH. The amino acid substrates were placed in positions common to both enzymes. The productive position of tyrosine in TH.BH4 was stabilized by a hydrogen bond with BH4. Despite favorable energy scores, tyrosine in a position trans to PAH residue His290 or TH residue His336 interferes with the access of the essential cofactor dioxygen to the catalytic center, thereby blocking the enzymatic reaction. DA and l-DOPA were directly coordinated to the active site iron via the hydroxyl residues of their catechol groups. Two alternative conformations, rotated 180 degrees around an imaginary iron-catecholamine axis, were found for DA and l-DOPA in PAH and for DA in TH. Electrostatic forces play a key role in hindering the bidentate binding of the immediate reaction product l-DOPA to TH, thereby saving the enzyme from direct feedback inhibition.

Amino Acid Sequence↗

Structure and specificity of a human valacyclovir activating enzyme: a homology model of BPHL.

Biphenyl hydrolase-like (BPHL) protein is a novel serine hydrolase which has been identified as human valacyclovirase (VACVase), catalyzing the hydrolytic activation of valine ester prodrugs of the antiviral drugs acyclovir and ganciclovir as well as other amino acid ester prodrugs of therapeutic nucleoside analogues. The broad specificity for nucleoside analogues as parent drugs suggests that BPHL may be particularly useful as a molecular target for prodrug activation. In order to develop an initial structural view of the specificity of BPHL, a homology model of BPHL based on the crystal structure of 2-hydroxy-6-oxo-7-methylocta-2,4-dienoate hydrolase was developed using the Molecular Operating Environment package (Chemical Computing Group, Montreal, Quebec), evaluated for its stereochemical quality and identification of free cysteines, and used in a molecular docking study. The BPHL model has residues S122, H255, and D227 comprising the putative catalytic triad in proximity and potential charge-charge interaction sites, M52 or D123 for the alpha-amino group. The model also suggested that the structural preference of BPHL for hydrophobic amino acyl promoieties and its limited activity for the secondary alcohol substrates may be attributed to the hydrophobic acyl-binding site formed by residues I158, G161, I162, and L229, and the spatial constraint around the catalytic site by a loop on one side, the active serine and histidine on the other side, and L53 and L179 on top. In addition, the broad specificity for nucleoside analogues may be due to the relatively less constrained nucleoside-binding site opening toward the entrance of the substrate-binding pocket. The homology model of BPHL provides a basis for further investigation of the catalytic and active site residues, can account for the observed structure activity profile of BPHL, and will be useful in the design of nucleoside prodrugs.

Acyclovir↗

Withanolides, a new class of natural cholinesterase inhibitors with calcium antagonistic properties.

The withanolides 1-3 and 4-5 isolated from Ajuga bracteosa and Withania somnifera, respectively, inhibited acetylcholinesterase (AChE, EC 3.1.1.7) and butyrylcholinesterase (BChE, EC 3.1.1.8) enzymes in a concentration-dependent fashion with IC50 values ranging between 20.5 and 49,2 microm and 29.0 and 85.2 microm for AChE and BChE, respectively. Lineweaver-Burk as well as Dixon plots and their secondary replots indicated that compounds 1, 3, and 5 are the linear mixed-type inhibitors of AChE, while 2 and 4 are non-competitive inhibitors of AChE with K(i) values ranging between 20.0 and 45.0 microm. All compounds were found to be non-competitive inhibitors of BChE with K(i) values ranging between 27.7 and 90.6 microm. Molecular docking study revealed that all the ligands are completely buried inside the aromatic gorge of AChE, while compounds 1, 3, and 5 extend up to the catalytic triad. A comparison of the docking results showed that all ligands generally adopt the same binding mode and lie parallel to the surface of the gorge. The superposition of the docked structures demonstrated that the non-flexible skeleton of the ligands always penetrates the aromatic gorge through the six-membered ring A, allowing their simultaneous interaction with more than one subsite of the active center. The affinity of ligands with AChE was found to be the cumulative effects of number of hydrophobic contacts and hydrogen bonding. Furthermore, all compounds also displayed dose-dependent (0.005-1.0 mg/mL) spasmolytic and Ca2+ antagonistic potentials in isolated rabbit jejunum preparations, compound 4 being the most active with an ED50 value of 0.09 +/- 0.001 mg/mL and 0.22 +/- 0.01 microg/mL on spontaneous and K+ -induced contractions, respectively. The cholinesterase inhibitory potential along with calcium antagonistic ability and safe profile in human neutrophil viability assay could make compounds 1-5 possible drug candidates for further study to treat Alzheimer's disease and associated problems.

Acetylcholinesterase↗

Microscopic evidence for a minus-end-directed power stroke in the kinesin motor ncd.

We used cryo-electron microscopy and image reconstruction to investigate the structure and microtubule-binding configurations of dimeric non-claret disjunctional (ncd) motor domains under various nucleotide conditions, and applied molecular docking using ncd's dimeric X-ray structure to generate a mechanistic model for force transduction. To visualize the alpha-helical coiled-coil neck better, we engineered an SH3 domain to the N-terminal end of our ncd construct (296-700). Ncd exhibits strikingly different nucleotide-dependent three-dimensional conformations and microtubule-binding patterns from those of conventional kinesin. In the absence of nucleotide, the neck adapts a configuration close to that found in the X-ray structure with stable interactions between the neck and motor core domain. Minus-end-directed movement is based mainly on two key events: (i) the stable neck-core interactions in ncd generate a binding geometry between motor and microtubule which places the motor ahead of its cargo in the minus-end direction; and (ii) after the uptake of ATP, the two heads rearrange their position relative to each other in a way that promotes a swing of the neck in the minus-end direction.

Adenosine Diphosphate↗

Electron transfer complexes of cytochrome c peroxidase from Paracoccus denitrificans containing more than one cytochrome.

According to the model proposed in previous papers [Pettigrew, G. W., Prazeres, S., Costa, C., Palma, N., Krippahl, L., and Moura, J. J. (1999) The structure of an electron-transfer complex containing a cytochrome c and a peroxidase, J. Biol. Chem. 274, 11383-11389; Pettigrew, G. W., Goodhew, C. F., Cooper, A., Nutley, M., Jumel, K., and Harding, S. E. (2003) Electron transfer complexes of cytochrome c peroxidase from Paracoccus denitrificans, Biochemistry 42, 2046-2055], cytochrome c peroxidase of Paracoccus denitrificans can accommodate horse cytochrome c and Paracoccus cytochrome c(550) at different sites on its molecular surface. Here we use (1)H NMR spectroscopy, analytical ultracentrifugation, molecular docking simulation, and microcalorimetry to investigate whether these small cytochromes can be accommodated simultaneously in the formation of a ternary complex. The pattern of perturbation of heme methyl and methionine methyl resonances in binary and ternary solutions shows that a ternary complex can be formed, and this is confirmed by the increase in the sedimentation coefficient upon addition of horse cytochrome c to a solution in which cytochrome c(550) fully occupies its binding site on cytochrome c peroxidase. Docking experiments in which favored binary solutions of cytochrome c(550) bound to cytochrome c peroxidase act as targets for horse cytochrome c and the reciprocal experiments in which favored binary solutions of horse cytochrome c bound to cytochrome c peroxidase act as targets for cytochrome c(550) show that the enzyme can accommodate both cytochromes at the same time on adjacent sites. Microcalorimetric titrations are difficult to interpret but are consistent with a weakened binding of horse cytochrome c to a binary complex of cytochrome c peroxidase and cytochrome c(550) and binding of cytochrome c(550) to the cytochrome c peroxidase that is affected little by the presence of horse cytochrome c in the other site. The presence of a substantial capture surface for small cytochromes on the cytochrome c peroxidase has implications for rate enhancement mechanisms which ensure that the two electrons required for re-reduction of the enzyme after reaction with hydrogen peroxide are delivered efficiently.

Animals↗

Can we separate active from inactive conformations?

Molecular modeling methodologies such as molecular docking, pharmacophore modeling, and 3D-QSAR, rely on conformational searches of small molecules as a starting point. All of these methodologies seek conformations of the small molecules as they bind to target proteins, i.e., their active conformations. Thus the question as to whether active conformations can be separated from inactive conformations is extremely relevant. In this paper, 3D-descriptors that separate random conformations from active conformations of small molecules are sought. To select appropriate descriptors, 65 protein-ligand complexes were taken from the protein data bank. For each ligand the active conformation was compared to randomly generated low energy conformations. Descriptors such as solvent accessible surface area, number of internal interactions and radius of gyration appear to be useful for separating the active conformations from the random conformations. The results with all these descriptors indicate that active conformations are less compact that random conformations, i.e., they have more solvent accessible surface area, fewer internal interactions and a larger radius of gyration than random conformations. Thus these descriptors could be useful as weights to bias conformational search procedures to conformations more likely to bind to proteins or as filters to eliminate conformations unlikely to bind to any protein.

Computer Simulation↗

The use of in vitro peptide binding profiles and in silico ligand-receptor interaction profiles to describe ligand-induced conformations of the retinoid X receptor alpha ligand-binding domain.

It is hypothesized that different ligand-induced conformational changes can explain the different interactions of nuclear receptors with regulatory proteins, resulting in specific biological activities. Understanding the mechanism of how ligands regulate cofactor interaction facilitates drug design. To investigate these ligand-induced conformational changes at the surface of proteins, we performed a time-resolved fluorescence resonance energy transfer assay with 52 different cofactor peptides measuring the ligand-induced cofactor recruitment to the retinoid X receptor-alpha (RXRalpha) in the presence of 11 compounds. Simultaneously we analyzed the binding modes of these compounds by molecular docking. An automated method converted the complex three-dimensional data of ligand-protein interactions into two-dimensional fingerprints, the so-called ligand-receptor interaction profiles. For a subset of compounds the conformational changes at the surface, as measured by peptide recruitment, correlate well with the calculated binding modes, suggesting that clustering of ligand-receptor interaction profiles is a very useful tool to discriminate compounds that may induce different conformations and possibly different effects in a cellular environment. In addition, we successfully combined ligand-receptor interaction profiles and peptide recruitment data to reveal structural elements that are possibly involved in the ligand-induced conformations. Interestingly, we could predict a possible binding mode of LG100754, a homodimer antagonist that showed no effect on peptide recruitment. Finally, the extensive analysis of the peptide recruitment profiles provided novel insight in the potential cellular effect of the compound; for the first time, we showed that in addition to the induction of coactivator peptide binding, all well-known RXRalpha agonists also induce binding of corepressor peptides to RXRalpha.

Amino Acid Sequence↗

Molecular dynamics simulations on SDF-1alpha: binding with CXCR4 receptor.

Insights into the interacting mode of CXCR4 with SDF-1alpha are crucial in understanding the structural and functional characteristics of CXCR4 receptor. In this paper a computational pipeline, integrating protein structure prediction, molecular dynamics simulations, automated molecular docking, and Brownian dynamics simulations were employed to investigate the dynamic and energetic aspects of CXCR4 associating with SDF-1alpha. The entire simulation revealed the surface distribution feature of electrostatic potentials and conformational "open-close" process of the receptor. The possible binding conformation of CXCR4 was identified, and the CXCR4-SDF-1alpha binding complex was generated. Arg188-Glu277 salt bridge plays an important role for both the extracellular domain conformational change and SDF-1alpha binding. Two binding sites were mapped at the extracellular domain (Site 1) and inside the transmembrane domain (Site 2), which are composed of conserved residues. Sites 1 and 2 contribute approximately 60% and 40% to the binding affinity with SDF-1alpha, respectively. The binding model is in agreement with most of the experimental data. Transmembrane VI has more significant motion in the harmonious conformational transition of CXCR4 during SDF-1alpha binding, which may be possibly associated with signal transduction. Based on the modeling and simulation, a binding mechanism hypothesis between CXCR4 and SDF-1alpha and its relationship to the signal transduction has been proposed.

Algorithms↗

Structural and functional characterization of human microsomal prostaglandin E synthase-1 by computational modeling and site-directed mutagenesis.

Microsomal prostaglandin (PG) E synthase-1 (mPGES-1) has recently been recognized as a novel, promising drug target for inflammation-related diseases. Functional and pathological studies on this enzyme further stimulate to understand its structure and the structure-function relationships. Using an approach of the combined structure prediction, molecular docking, site-directed mutagenesis, and enzymatic activity assay, we have developed the first three-dimensional (3D) model of the substrate-binding domain (SBD) of mPGES-1 and its binding with substrates prostaglandin H2 (PGH2) and glutathione (GSH). In light of the 3D model, key amino acid residues have been identified for the substrate binding and the obtained experimental activity data have confirmed the computationally determined substrate-enzyme binding mode. Both the computational and experimental results show that Y130 plays a vital role in the binding with PGH2 and, probably, in the catalytic reaction process. R110 and T114 interact intensively with the carboxyl tail of PGH2, whereas Q36 and Q134 only enhance the PGH2-binding affinity. The modeled binding structure indicates that substrate PGH2 interacts with GSH through hydrogen binding between the peroxy group of PGH2 and the -SH group of GSH. The -SH group of GSH is expected to attack the peroxy group of PGH2, initializing the catalytic reaction transforming PGH2 to prostaglandin E2 (PGE2). The overall agreement between the calculated and experimental results demonstrates that the predicted 3D model could be valuable in future rational design of potent inhibitors of mPGES-1 as the next-generation inflammation-related therapeutic.

Algorithms↗

A novel approach to the design of inhibitors of human secreted phospholipase A2 based on native peptide inhibition.

Human Type IIA secreted phospholipase A(2) (sPLA(2)-IIA) is an important modulator of cytokine-dependent inflammatory responses and a member of a growing superfamily of structurally related phospholipases. We have previously shown that sPLA(2)-IIA is inhibited by a pentapeptide sequence comprising residues 70-74 of the native sPLA(2)-IIA protein and that peptides derived from the equivalent region of different sPLA(2)-IIA species specifically inhibit the enzyme from which they are derived. We have now used an analogue screen of the human pentapeptide (70)FLSYK(74) in which side-chain residues were substituted, together with molecular docking approaches that modeled low-energy conformations of (70)FLSYK(74) bound to human sPLA(2)-IIA, to generate inhibitors with improved potency. Importantly, the modeling studies showed a close association between the NH(2) and COOH termini of the peptide, predicting significant enhancement of the potency of inhibition by cyclization. Cyclic compounds were synthesized and indeed showed 5-50-fold increased potency over the linear peptide in an Escherichia coli membrane assay. Furthermore, the potency of inhibition correlated with steady-state binding of the cyclic peptides to sPLA(2)-IIA as determined by surface plasmon resonance studies. Two potential peptide interaction sites were identified on sPLA(2)-IIA from the modeling studies, one in the NH(2)-terminal helix and the other in the beta-wing region, and in vitro association assays support the potential for interaction of the peptides with these sites. The inhibitors were effective at nanomolar concentrations in blocking sPLA(2)-IIA-mediated amplification of cytokine-induced prostaglandin synthesis in human rheumatoid synoviocytes in culture. These studies provide an example where native peptide sequences can be used for the development of potent and selective inhibitors of enzyme function.

Amino Acid Sequence↗

Structure-based design of selective inhibitors of dihydrofolate reductase: synthesis and antiparasitic activity of 2, 4-diaminopteridine analogues with a bridged diarylamine side chain.

As part of a larger search for potent as well as selective inhibitors of dihydrofolate reductase (DHFR) enzymes from opportunistic pathogens found in patients with AIDS and other immune disorders, N-[(2,4-diaminopteridin-6-yl)methyl]dibenz[b,f]azepine (4a) and the corresponding dihydrodibenz[b,f]azepine, dihydroacridine, phenoxazine, phenothiazine, carbazole, and diphenylamine analogues were synthesized from 2, 4-diamino-6-(bromomethyl)pteridine in 50-75% yield by reaction with the sodium salts of the amines in dry tetrahydrofuran at room temperature. The products were tested for the ability to inhibit DHFR from Pneumocystis carinii (pcDHFR), Toxoplasma gondii (tgDHFR), Mycobacterium avium (maDHFR), and rat liver (rlDHFR). The member of the series with the best combination of potency and species selectivity was 4a, with IC(50) values against the four enzymes of 0. 21, 0.043, 0.012, and 4.4 microM, respectively. The dihydroacridine, phenothiazine, and carbazole analogues were also potent, but nonselective. Of the compounds tested, 4a was the only one to successfully combine the potency of trimetrexate with the selectivity of trimethoprim. Molecular docking simulations using published 3D structural coordinates for the crystalline ternary complexes of pcDHFR and hDHFR suggested a possible structural interpretation for the binding selectivity of 4a and the lack of selectivity of the other compounds. According to this model, 4a is selective because of a unique propensity of the seven-membered ring in the dibenz[b,f]azepine moiety to adopt a puckered orientation that allows it to fit more comfortably into the active site of the P. carinii enzyme than into the active site of the human enzyme. Compound 4a was also evaluated for the ability to be taken up into, and retard the growth of, P. carinii and T. gondii in culture. The IC(50) of 4a against P. carinii trophozoites after 7 days of continuous drug treatment was 1.9 microM as compared with previously observed IC(50) values of >340 microM for trimethoprim and 0.27 microM for trimetrexate. In an assay involving [(3)H]uracil incorporation into the nuclear DNA of T. gondii tachyzoites as the surrogate endpoint for growth, the IC(50) of 4a after 5 h of drug exposure was 0.077 microM. The favorable combination of potency and enzyme selectivity shown by 4a suggests that this novel structure may be an interesting lead for structure-activity optimization.

Animals↗

Molecular modeling of a putative antagonist binding site on helix III of the beta-adrenoceptor.

In recent biochemical studies it was demonstrated that residue Asp113 of the beta-adrenoceptor (beta-AR) is an indispensable amino acid for the binding of beta-AR antagonists. Earlier fluorescence studies showed that a tryptophan-rich region of the beta-AR is involved in the binding of propranolol, the prototype beta-AR antagonist. Bearing these two biochemical findings in mind, we explored the beta-AR part containing Asp113, for an energetically favorable antagonist binding site. This was done by performing molecular docking studies with the antagonist propranolol and a specific beta-AR peptide which included, besides Asp113, two possibly relevant tryptophan residues. In the docking calculations, the propranolol molecule was allowed to vary all its internal torsional angles. The receptor peptide was kept in an alpha-helix conformation, while side chains relevant to ligand binding were flexible to enable optimal adaptations to the ligand's binding conformation. By means of force-field calculations the total energy was minimized, consisting of the intramolecular energies of both ligand and receptor peptide, and the intermolecular energy. We found an antagonist binding site, consisting of amino acids Asp113 and Trp109, which enabled energetically favorable interactions with the receptor-binding groups of propranolol. According to these results, binding involves three main interaction points: (i) a reinforced ionic bond; (ii) a hydrogen bond; and (iii) a hydrophobic/charge transfer interaction. The deduced binding site shows a difference in affinity between the levo- and dextrorotatory isomers of propranolol caused by a difference in ability to form a hydrogen bond, which is in conformity with the experimentally observed stereoselectivity. Moreover, it also provides an explanation for the beta 1-selectivity of p-phenyl substituted phenoxypropanolamines like betaxolol. The p-phenyl substituent of betaxolol was shown to be sterically hindered upon binding to the beta 2-AR peptide, whereas this hindrance is very likely to be much less with the beta 1-AR peptide. Finally, the proposed antagonist binding site is discussed in the light of some recent biochemical findings and theories.

Adrenergic beta-Antagonists↗

Discovery of a daunorubicin analogue that exhibits potent antitumor activity and overcomes P-gp-mediated drug resistance.

Anthracyclines are considered to be some of the most effective anticancer drugs for cancer therapy. However, drug resistance and cardiotoxicity of anthracyclines limit their clinical application. We hypothesize that direct modifications of the sugar moiety of anthracyclines avert P-glycoprotein (P-gp) recognition and efflux, increase drug intracellular concentration in cancer cells, and thus overcome P-gp-mediated drug resistance. Daunorubicin (DNR) analogues with sugar modifications were synthesized by directly transforming the amino group of DNR to an azido group or triazole group. Molecular docking showed that the lead compound (3'-azidodaunorubicin, ADNR) averts P-gp binding, while daunorubicin (DNR) extensively interacts with multidrug-resistance (MDR) protein through H-bonds and electrostatic interactions. FACS assay demonstrated that these new compounds abolished P-gp drug efflux and accumulated high intracellular concentration in the drug-resistant leukemia K562/Dox. P-gp inhibition by CsA confirmed that these new analogues are no longer P-gp substrates. ADNR exhibited potent anticancer activity in both drug-sensitive (K562) and drug-resistant leukemia cells (K562/Dox), with a 25-fold lower drug resistance index than DNR. An in vivo xenograft model demonstrated that ADNR showed more than 2.5-fold higher maximum growth inhibition rate against drug-resistant cancers and significant improvement for animal survival rate versus DNR. No significant body weight reduction in mice was observed for ADNR at the maximum tolerable dose, as compared to more than 70% body weight reduction for DNR. These data suggest that sugar modifications of anthracyclines avert P-gp binding, abolish P-gp-mediated drug efflux, increase intracellular drug concentration, and thus overcome P-gp-mediated drug resistance in cancer therapy.

ATP Binding Cassette Transporter, Subfamily B, Mem↗

Juliflorine: a potent natural peripheral anionic-site-binding inhibitor of acetylcholinesterase with calcium-channel blocking potential, a leading candidate for Alzheimer's disease therapy.

The alkaloid juliflorine (1) from Prosopis juliflora inhibited acetylcholinesterase (AChE, EC 3.1.1.7) and butyrylcholinesterase (BChE, EC 3.1.1.8) enzymes in a concentration-dependent fashion with IC50 values 0.42 and 0.12 microM, respectively. Lineweaver-Burk as well as Dixon plots and their secondary replots indicated that the nature of inhibition was purely of non-competitive type with Ki values 0.4 and 0.1 microM, against AChE and BChE, respectively. By molecular docking studies compound 1 was found to be ideally spaced inside the aromatic gorge of AChE with rings A/B remaining at the top and rings C/D penetrating deep into the gorge, that might be due to the greater hydrophobicity of rings C/D as compared to rings A/B, allowing their simultaneous interaction with the peripheral anionic and quaternary ammonium-binding sites. The 1-AChE complex was found to be stabilized by hydrophobic contacts, hydrogen bonding, and pi-pi stacking between the compound 1 and amino acid residues of the aromatic gorge of AChE. Amino acid residues Tyr70, Asp72, Tyr121, Trp279, and Tyr334 of the peripheral anionic site (PAS) of AChE were found to be exclusively involved in the hydrophobic contacts with compound 1 that might be responsible for the competitive mode of inhibition. Compound 1 also showed dose-dependent (30-500 microg/mL) spasmolytic and Ca2+-channel blocking activities in isolated rabbit jejunum preparations. The cholinesterase inhibitory potential along with calcium-channel blocking activity of compound 1 and safe profile in human neutrophils viable assay could make it a possible drug candidate for Alzheimer's disease.

Alkaloids↗

Insights into the functional architecture of the catalytic center of a maize beta-glucosidase Zm-p60.1.

The maize (Zea mays) beta-glucosidase Zm-p60.1 has been implicated in regulation of plant development by the targeted release of free cytokinins from cytokinin-O-glucosides, their inactive storage forms. The crystal structure of the wild-type enzyme was solved at 2.05-A resolution, allowing molecular docking analysis to be conducted. This indicated that the enzyme specificity toward substrates with aryl aglycones is determined by aglycone aromatic system stacking with W373, and interactions with edges of F193, F200, and F461 located opposite W373 in a slot-like aglycone-binding site. These aglycone-active site interactions recently were hypothesized to determine substrate specificity in inactive enzyme substrate complexes of ZM-Glu1, an allozyme of Zm-p60.1. Here, we test this hypothesis by kinetic analysis of F193I/Y/W mutants. The decreased K(m) of all mutants confirmed the involvement of F193 in determining enzyme affinity toward substrates with an aromatic aglycone. It was unexpected that a 30-fold decrease in k(cat) was found in F193I mutant compared with the wild type. Kinetic analysis and computer modeling demonstrated that the F193-aglycone-W373 interaction not only contributes to aglycone recognition as hypothesized previously but also codetermines catalytic rate by fixing the glucosidic bond in an orientation favorable for attack by the catalytic pair, E186 and E401. The catalytic pair, assigned initially by their location in the structure, was confirmed by kinetic analysis of E186D/Q and E401D/Q mutants. It was unexpected that the E401D as well as C205S and C211S mutations dramatically impaired the assembly of a catalysis-competent homodimer, suggesting novel links between the active site structure and dimer formation.

Amino Acid Sequence↗

Modeling subtype-selective agonists binding with alpha4beta2 and alpha7 nicotinic acetylcholine receptors: effects of local binding and long-range electrostatic interactions.

The subtype-selective binding of 14 representative agonists with alpha4beta2 and alpha7 nicotinic acetylcholine receptors (nAChRs) has been studied by performing homology modeling, molecular docking, geometry optimizations, and microscopic and phenomenological binding free energy calculations. All of the computational results demonstrate that the subtype selectivity of the agonists binding with alpha4beta2 and alpha7 7 nAChRs is affected by both local binding and long-range electrostatic interactions between the receptors and the protonated structures of the agonists. The effects of the long-range electrostatic interactions are mainly due to the distinct difference in the net charge of the ligand-binding domain between the two nAChR subtypes. For the alpha4beta2-selective agonists examined, the microscopic binding modes with the alpha4beta2 nAChR are very similar to the corresponding modes with the alpha7 nAChR, and therefore, the subtype selectivity of these agonists binding with alpha4beta2 and alpha7 nAChRs is dominated by the long-range electrostatic interactions. For the alpha7-selective agonists, their microscopic binding modes with the alpha7 nAChR are remarkably different from those with the alpha4beta2 nAChR so that the local binding (including the hydrogen bonding and cation-pi interactions) with the alpha7 nAChR is much stronger than that with the alpha4beta2 nAChR. The calculated phenomenological binding free energies are in good agreement with available experimental data for the relative binding free energies concerning the subtype selectivity of agonists binding with the two different nAChR subtypes. The fundamental insights obtained in the present study should be valuable for future rational design of potential therapeutic agents targeted to specific nAChR subtypes.

Amino Acid Sequence↗

Deguanylation of human immunodeficiency virus (HIV-1) RNA by recombinant pokeweed antiviral protein.

Modeling studies, combined with the molecular docking of the trinucleotide GGG into the active site of the deadenylating RNA N-glycosidase pokeweed antiviral protein (PAP), indicated that a guanine base can fit into the active site pocket of PAP without disturbing its unique geometry and is sandwiched between residues Tyr(72) and Tyr(123) very much like an adenine base. The guanine base can form two specific hydrogen bonds with the active site residues Ser(121) and Val(73) and the attached negatively charged phosphate groups can entertain stabilizing electrostatic interactions with two clusters of positively charged patches on the PAP surface formed by Lys(210) and Arg(179) from one side and Arg(122) and Arg(135) from the other side of the active site. These observations prompted the hypothesis that the RNA depurinating activity of PAP may not be restricted to adenine residues and PAP should be capable of deguanylating ribosomal and viral RNA as well. This hypothesis was experimentally confirmed by direct demonstration that guanine base is released from both ribosomal and HIV-1 RNA after treatment with purified recombinant PAP using quantitative high performance liquid chromatography. Recombinant PAP released adenine and guanine residues at a 1:1 ratio from HIV-1 RNA and at an approximately 3:1 (adenine:guanine) ratio from Escherichia coli ribosomal RNA. At a concentration of 5 microM, recombinant PAP released 263 +/- 10 pmol of adenine and 100 +/- 11 pmol of guanine from 1 microgram of E. coli ribosomal RNA (16S + 23S) within 4 h of treatment. By comparison, 138 +/- 12 pmol of adenine and 143 +/- 10 pmol of guanine were released from 1 microgram of HIV-1 RNA under identical treatment conditions (5 microM recombinant PAP, 4 h treatment). The deguanylation of the ribosomal and viral RNA targets by recombinant PAP was concentration-dependent and is abolished by alanine substitutions of the catalytic active site residues Tyr(72) and Tyr(123). To our knowledge, these findings provide the first evidence that PAP can deguanylate both ribosomal and viral RNA.

Adenine Nucleotides↗

Ligand-Mediated Reprogramming Redirects Liver-Tropic Ionizable Lipid Nanoparticles for Lung-Selective mRNA Delivery.

Systemic delivery of messenger RNA (mRNA) to target tissues and cells using lipid nanoparticles (LNPs) holds transformative potential for gene therapy. However, most clinically validated LNP exhibit strong liver tropism, and redirecting their organ specificity without redesigning entirely new chemistries remains challenging. Here we present a ligand-mediated lipid reprogramming approach that repurposes chemically defined, liver-tropic, ionizable lipids (lipidoids) for mRNA delivery beyond the liver. From a library of 90 degradable lipidoids, we identified 2-t6b as a potent liver-targeting platform. By site-specific displaying of small molecule ligands onto 2-t6b headgroup, we engineered a series of reconfigured lipidoids that achieve lung-specific targeting while retaining the parent delivery scaffold. Ligand7-2-t6b-lipid-functionalized LNP achieved over 200-fold higher mRNA translation in the lungs compared to the parent liver-tropic LNP. Proteomics and molecular docking analysis revealed enhanced binding of the modified lipid to vitronectin, a serum glycoprotein that improves integrin binding and thus promotes cellular uptake and translation efficiency. Ligand-mediated 2-t6b/ligand7 LNPs achieved outperformed efficacy and therapeutic potential in lung-specific genome editing relative to SORT-constructed 2-t6b LNP system. Our modular reprogramming strategy provides a generalizable framework to upgrade existing liver-biased LNPs into lung-selective mRNA carriers, advancing next-generation tissue-specific mRNA therapies for gene editing, protein replacement therapy, and regenerative medicine.

RNA, Messenger↗