PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “molecular dynamic simulation assay”

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 127 records · Page 7Linked to original sources

Conformational analyses of somatostatin-related cyclic hexapeptides containing peptoid residues.

We report the conformational analysis by 1H NMR in DMSO and computer simulations involving distance geometry and molecular dynamics simulations of a series of peptoid analogues of the cyclic hexapeptide c[Phe11-Pro6-Phe7-d-Trp8-Lys9-Thr10] (1). The proline residue in compound 1 is replaced with the peptoid residues N-benzylglycine (Nphe) (compound 2), N-(S)-alpha-methylbenzylglycine [(S)-beta-MeNphe] (compound 3), and N-(R)-alpha-methylbenzylglycine [(R)-beta-MeNphe] (compound 4). The peptoid analogues 2 and 4 exhibit potent binding activities to the hsst2 receptor, while the binding affinities to the hsst5 and to the hsst3 receptors are reduced compared to that of the parent compound 1. Compound 3 shows reduced binding activities to the hsst2, hsst3, and hsst5 receptors compared to compound 1. The results of in vivo assays indicate that these compounds inhibit the growth hormone release but do not affect the insulin release. These peptoid-containing analogues show two sets of NMR signals corresponding to cis and trans conformations of the peptide bond between Phe11 and Nxaa6. We demonstrate that the backbone conformation and the orientation of the relevant side chains of compound 1 are maintained in the cis isomers of the peptoid analogues which adopt a type VI beta-turn centered around residues 11 and 6 and a type II' beta-turn with d-Trp in the i+1 position. The enhanced selectivity of the peptoid-containing analogues compared to compound 1 and the results of the conformational analysis suggest that the presence of a conformationally constrained hydrophobic group in position 6 in complementary topology to the Phe11 side chain enhances selective binding to the hsst2 receptor.

Computer Simulation↗

Modeling of flap endonuclease interactions with DNA substrate.

Structure-specific 5' nucleases play an important role in DNA replication and repair uniquely recognizing an overlap flap DNA substrate and processing it into a DNA nick. However, in the absence of a high-resolution structure of the enzyme/DNA complex, the mechanism underlying this recognition and substrate specificity, which is key to the enzyme's function, remains unclear. Here, we propose a three-dimensional model of the structure-specific 5' flap endonuclease from Pyrococcus furiosus in its complex with DNA. The model is based on the known X-ray structure of the enzyme and a variety of biochemical and molecular dynamics (MD) data utilized in the form of distance restraints between the enzyme and the DNA. Contacts between the 5' flap endonuclease and the sugar-phosphate backbone of the overlap flap substrate were identified using enzyme activity assays on substrates with methylphosphonate or 2'-O-methyl substitutions. The enzyme footprint extends two to four base-pairs upstream and eight to nine base-pairs downstream of the cleavage site, thus covering 10-13 base-pairs of duplex DNA. The footprint data are consistent with a model in which the substrate is bound in the DNA-binding groove such that the downstream duplex interacts with the helix-hairpin-helix motif of the enzyme. MD simulations to identify the substrate orientation in this model are consistent with the results of the enzyme activity assays on the methylphosphonate and 2'-O-methyl-modified substrates. To further refine the model, 5' flap endonuclease variants with alanine point substitutions at amino acid residues expected to contact phosphates in the substrate and one deletion mutant were tested in enzyme activity assays on the methylphosphonate-modified substrates. Changes in the enzyme footprint observed for two point mutants, R64A and R94A, and for the deletion mutant in the enzyme's beta(A)/beta(B) region, were interpreted as being the result of specific interactions in the enzyme/DNA complex and were used as distance restraints in MD simulations. The final structure suggests that the substrate's 5' flap interacts with the enzyme's helical arch and that the helix-hairpin-helix motif interacts with the template strand in the downstream duplex eight base-pairs from the cleavage site. This model suggests specific interactions between the 3' end of the upstream oligonucleotide and the enzyme. The proposed structure presents the first detailed description of substrate recognition by structure-specific 5' nucleases.

Base Sequence↗

Molecular identification of canine podocalyxin-like protein 1 as a renal tubulogenic regulator.

GP135 is an apical membrane protein expressed in polarized MDCK epithelial cells. When cultured in three-dimensional collagen gel, MDCK cells form branching tubules in response to hepatocyte growth factor stimulation in a manner that simulates the embryonic renal development. During this process, GP135 displays transient loss of membranous localization but reappears at the cell surface when nascent lumen emerges from the developing tubules. Despite being used for decades as the canonical hallmark of apical surface, the molecular identity and the significance of the dynamic expression of GP135 during the tubulogenic process remain elusive. For exploring the function of GP135, the full-length cDNA encoding GP135 was obtained. Sequence alignments and features analysis confirm GP135 as a canine homolog of podocalyxin, confirming the finding of an earlier independent study. Immunohistochemical assays on canine kidney sections identified both glomerular and tubular distribution of GP135 along the nephron. Mutant MDCK cells expressing siRNA targeted at two regions of GP135 show defects in hepatocyte growth factor-induced tubulogenesis. Re-expression of full-length and an O-linked glycosylation abbreviated construct of GP135 could recapitulate the tubulogenesis process lacking in siRNA knockdown cells; however, a deletion construct devoid of the cytoplasmic domain failed to rescue the phenotype. In summary, the data identify the MDCK apical domain marker GP135 as a tubular form of podocalyxin and provide evidence for its importance in renal tubulogenesis.

Animals↗

A bioactive somatostatin analog without a type II' beta-turn: synthesis and conformational analysis in solution.

A cyclic somatostatin analog [structure: see text] (1) has been synthesized. Biological assays show that this compound has strong binding affinities to somatostatin hsst2 and hsst5 receptor subtypes (5.2 and 1.2 nM, respectively, and modest affinity to hsst4 (41.1 nM)). Our conformational analysis carried out in DMSO-d6 indicates that this compound exists as two structures arising from the trans and cis configurations of the peptide bond between Phe7 and N-alkylated Gly8. However, neither conformer exhibits a type II' beta-turn. This is the first report of a potent bioactive somatostatin analog that does not exhibit a type II' beta-turn in solution. Molecular dynamics simulations (500 ps) carried out at 300 K indicate that the backbone of compound 1 is more flexible than other cyclic somatostatin analogs formed by disulfide bonds.

Amino Acid Sequence↗

Antagonist profile and molecular dynamic simulation of a Drosophila melanogaster muscarinic acetylcholine receptor.

A stably-transfected, Drosophila cell line (S2-DMl-1) expressing the Drosophila DMl muscarinic acetylcholine receptor (mAChR) exhibits high-affinity, saturable, specific binding of the radiolabelled muscarinic antagonist [3H]-N-methyl scopolamine ([3H]-NMS) with an equilibrium dissociation constant (Kd) of 0.67 +/- 0.02 and a Bmax of 1.53 +/- 0.3 pmol/mg protein. Displacement of [3H]-NMS by mAChR antagonists results in the pharmacological profile: 4-diphenylacetoxy-N-methylpiperidine methiodide (4-DAMP) > hexahydrosiladifenidol > p-fluorohexahydrosiladifenidol > nitrocaramiphen > pirenzepine > methoctramine > AFDX-116. This antagonist profile most closely resembles that of the vertebrate M3 mAChR subtype. In this study, however, we have demonstrated that the antagonist profile of DM1 is distinct from those of vertebrate mAChR subtypes. Molecular dynamic simulations of the Drosophila muscarinic receptor are presented in the free, carbamylcholine-bound and NMS-bound forms. Theoretical, quantitative structure-activity relationship models have been developed; a good correlation is observed between the interaction energies of the minimized ligand-receptor complexes and the pharmacological affinities of the antagonists tested.

Animals↗

Extending the understanding of mutagenicity: structural insights into primer-extension past a benzo[a]pyrene diol epoxide-DNA adduct.

DNA polymerase enzymes employ a number of innate fidelity mechanisms to ensure the faithful replication of the genome. However, when confronted with DNA damage, their fidelity mechanisms can be evaded, resulting in a mutation that may contribute to the carcinogenic process. The environmental carcinogen benzo[a]pyrene is metabolically activated to reactive intermediates, including the tumorigenic (+)-anti-benzo[a]pyrene diol epoxide, which can attack DNA at the exocyclic amino group of guanine to form the major (+)-trans-anti-[BP]-N(2)-dG adduct. Bulky adducts such as (+)-trans-anti-[BP]-N(2)-dG primarily block DNA replication, but are occasionally bypassed and cause mutations if paired with an incorrect base. In vitro standing-start primer-extension assays show that the preferential insertion of A opposite (+)-trans-anti-[BP]-N(2)-dG is independent of the sequence context, but the primer is extended preferentially when dT is positioned opposite the damaged base in a 5'-CG*T-3' sequence context. Regardless of the base positioned opposite (+)-trans-anti-[BP]-N(2)-dG, extension of the primer past the lesion site poses the greatest block to polymerase progression. In order to gain insight into primer-extension of each base opposite (+)-trans-anti-[BP]-N(2)-dG, we carried out molecular modeling and 1.25 ns unrestrained molecular dynamics simulations of the adduct in the +1 position of the template within the replicative pol I family T7 DNA polymerase. Each of the four bases was modeled at the 3' terminus of the primer, incorporated opposite the adduct, and the next-to-be replicated base was in the active site with its Watson-Crick partner as the incoming nucleotide. As in our studies of nucleotide incorporation, (+)-trans-anti-[BP]-N(2)-dG was modeled in the syn conformation in the +1 position, with the BP moiety on the open major groove side of the primer-template duplex region, leaving critical protein-DNA interactions intact. The present work revealed that the efficiency of primer-extension past this bulky adduct opposite each of the four bases in the 5'-CG*T-3' sequence can be rationalized by the stability of interactions between the polymerase protein, primer-template DNA and incoming nucleotide. However, the relative stabilization of each nucleotide opposite (+)-trans-anti-[BP]-N(2)-dG in the +1 position (T > G > A > or = C) differed from that when the adduct and partner were the nascent base-pair (A > T > or = G > C). In addition, extension past (+)-trans-anti-[BP]-N(2)-dG may pose a greater block to a high fidelity DNA polymerase than does nucleotide incorporation opposite the adduct because the presence of the modified base-pair in the +1 position is more disruptive to the polymerase-DNA interactions than it is within the active site itself. The dN:(+)-trans-anti-[BP]-N(2)-dG base-pair is strained to shield the bulky aromatic BP moiety from contact with the solvent in the +1 position, causing disruption of protein-DNA interactions that would likely result in decreased extension of the base-pair. These studies reveal in molecular detail the kinds of specific structural interactions that determine the function of a processive DNA polymerase when challenged by a bulky DNA adduct.

Bacteriophage T7↗

A combination of docking/dynamics simulations and pharmacophoric modeling to discover new dual c-Src/Abl kinase inhibitors.

A computational protocol was applied to identify molecular scaffolds untested toward the c-Src tyrosine kinase. A combination of docking and dynamics calculations allowed us to build three-dimensional models of the complexes between Src and several of its known inhibitors. Interactions most contributing to activity of the inhibitors, in terms of hydrogen bonds and hydrophobic contacts, were codified into pharmacophoric models that were in turn applied to perform a search of commercially available compounds within the Asinex database. As a result, we identified 1,3,4-thiadiazoles and pyrazolydine-3,5-diones showing inhibitory activity in the submicromolar range in a cell-free assay toward Src. Moreover, since several of the compounds used to generate pharmacophores were also known as Abl inhibitors, we tested the identified hits toward Abl tyrosine kinase, finding activity in the submicromolar range. Such biological data suggested that the computational protocol is an efficient tool for identifying new hits toward both Src and Abl.

CSK Tyrosine-Protein Kinase↗

Virtual screening of novel noncovalent inhibitors for SARS-CoV 3C-like proteinase.

The SARS coronavirus 3C-like proteinase is considered as a potential drug design target for the treatment of severe acute respiratory syndrome (SARS). Owing to the lack of available drugs for the treatment of SARS, the discovery of inhibitors for SARS coronavirus 3C-like proteinase that can potentially be optimized as drugs appears to be highly desirable. We have built a "flexible" three-dimensional model for SARS 3C-like proteinase by homology modeling and multicanonical molecular dynamics method and used the model for virtual screening of chemical databases. After Dock procedures, strategies including pharmocophore model, consensus scoring, and "drug-like" filters were applied in order to accelerate the process and improve the success rate of virtual docking screening hit lists. Forty compounds were purchased and tested by HPLC and colorimetric assay against SARS 3C-like proteinase. Three of them including calmidazolium, a well-known antagonist of calmodulin, were found to inhibit the enzyme with an apparent K(i) from 61 to 178 microM. These active compounds and their binding modes provide useful information for understanding the binding sites and for further selective drug design against SARS and other coronavirus.

Binding Sites↗

An altered peptide ligand for naïve cytotoxic T lymphocyte epitope of TRP-2(180-188) enhanced immunogenicity.

Tyrosinase-related protein-2 (TRP-2) is a non-mutated melanocyte differentiation antigen. The TRP-2-recognizing CD8(+) T cells can evoke immune responses to melanoma in both humans and mice. Developing epitopes with amino acid replacements in their sequences might improve the low immunogenicity against this 'self' tumor antigen. We designed altered peptide ligands (APLs) of TRP-2((180-188)) (SVYDFFVWL) with preferred primary and auxiliary HLA-A*0201 molecule anchor residue replacement. These APLs were screened for MHC-affinity by affinity prediction plots and molecular dynamics simulation, and analyzed in vitro for stability and binding-affinity to molecular HLA-A*0201. We also investigated the CTLs activities induced by TRP-2 wild-type epitope and the APLs both in vitro in human PBMCs and HLA-A2.1/K(b) transgenic mice. The results indicate that TRP-2 2M analog simultaneously had stronger binding-affinity and a lower dissociation rate to HLA-A*0201, than wild-type peptide. In addition, the analog 2M was superior to other APLs and wild-type epitope in terms of immunological efficacy ex vivo as measured by the ELISPOT assays of IFN-gamma and granzyme B. These results demonstrate that TRP-2 2M is an agonist epitope that can induce anti-tumor immunity superior to its wild-type epitope, and has potential application in peptide-mediated immunotherapy.

Amino Acid Substitution↗

Discovery of NAT-6-321056 as a novel modulator of VEGFR2 signaling to suppress tumor angiogenesis.

Vascular endothelial growth factor receptor 2 (VEGFR2) is a master regulator of angiogenesis and cancer progression. However, current VEGFR2 modulators face significant challenges, including off-target toxicity and acquired resistance, underscoring the urgent need for novel therapeutic agents with improved efficacy and safety profiles. Here, we reported that virtual screening of 39,442 natural products from the ZINC natural products-derived library, coupled with molecular docking and molecular dynamics (MD) simulations to evaluate the binding stability of candidate compounds, identified NAT-6-321056 as a highly promising modulator of VEGFR2 signaling. Biological evaluations demonstrated that NAT-6-321056 exerted potent inhibition on the growth of a broad spectrum of cancer cells, including both solid tumors and hematological malignancies. In EA.hy 926 endothelial cells and SK-N-DZ neuroblast cells, the compound significantly suppressed proliferation, migration, and invasion. Microscale thermophoresis (MST) confirmed direct binding of NAT-6-321056 to VEGFR2 with favorable affinity. Kinase profiling against a panel of 33 kinases indicated that NAT-6-321056 exhibited a multi-kinase modulation profile. Mechanistic studies revealed that NAT-6-321056 suppressed the expression of hypoxia-inducible factor 1-alpha (HIF-1α) and was associated with reduced VEGFR2 phosphorylation and attenuation of the downstream ERK/JNK/AKT signaling pathways. Moreover, NAT-6-321056 exhibited robust in vivo anti-angiogenic effects in both the chick chorioallantoic membrane (CAM) assay and transgenic zebrafish vascular fluorescence imaging models. Computational absorption, distribution, metabolism, excretion, and toxicity (ADMET) prediction suggested acceptable drug-like properties. Collectively, these findings demonstrated that NAT-6-321056 is a promising modulator of VEGFR2 signaling with potent anti-angiogenic activity and represents a viable candidate for cancer therapy.

Vascular Endothelial Growth Factor Receptor-2↗

Prediction of HIV-1 protease inhibitor resistance using a protein-inhibitor flexible docking approach.

Emergence of drug resistance remains one of the most challenging issues in the treatment of HIV-1 infection. Here we focus on resistance to HIV-1 protease inhibitors (PIs) at a molecular level, which can be analysed genotypically or phenotypically. Genotypic assays are based on the analysis of mutations associated with reduced drug susceptibility, but are problematic because of the numerous mutations and mutational patterns that confer drug resistance. Phenotypic resistance or susceptibility can be experimentally evaluated by measuring the amount of free drug bound to HIV-1 protease molecules, but this procedure is expensive and time-consuming. To overcome these problems, we have developed a docking protocol that takes protein-inhibitor flexibility into account to predict phenotypic drug resistance. For six FDA-approved Pls and a total of 1792 HIV-1 protease sequence mutants, we used a combination of inhibitor flexible docking and molecular dynamics (MD) simulations to calculate protein-inhibitor binding energies. Prediction results were expressed as fold changes of the calculated inhibitory constant (Ki), and the samples predicted to have fold-increase in calculated Ki above the fixed cut-off were defined as drug resistant. Our combined docking and MD protocol achieved accuracies ranging from 72-83% in predicting resistance/susceptibility for five of the six drugs evaluated. Evaluating the method only on samples where our predictions concurred with established knowledge-based methods resulted in increased accuracies of 83-94% for the six drugs. The results suggest that a physics-based approach, which is readily applicable to any novel PI and/or mutant, can be used judiciously with knowledge-based approaches that require experimental training data to devise accurate models of HIV-1 Pl resistance prediction.

Binding Sites↗

Vitamin D receptor agonists specifically modulate the volume of the ligand-binding pocket.

Existing crystal structure data has indicated that 1alpha,25-dihydroxyvitamin D(3) (1alpha,25(OH)(2) D(3)) and its analogues bind the ligand-binding pocket (LBP) of the human vitamin D receptor in a very similar fashion. Because docking of a ligand into the LBP is a more flexible process than crystallography can monitor, we analyzed 1alpha,25(OH)(2)D(3), its 20-epi derivative MC1288, the two side-chain analogues Gemini and Ro43-83582 (a hexafluoro-derivative) by molecular dynamics simulations in a complex with the vitamin D receptor ligand-binding domain and a co-activator peptide. Superimposition of the structures showed that the side chain of MC1288, the first side chain of the conformation II of Gemini, the second side chain of Ro43-83582 in conformation I and the first side chain of Ro43-83582 in conformation II take the same agonistic position as the side chain of 1alpha,25(OH)(2)D(3). Compared with the LBP of the natural hormone MC1288 reduced the volume by 17%, and Gemini expanded it by 19%. The shrinking of the LBP of MC1288 and its expansion to accommodate the second side chain of Gemini or Ro43-83582 is the combined result of minor movements of more than 30 residues and major movements of a few critical amino acids. The agonist-selective recognition of anchoring OH groups by the conformational flexible residues Ala-303, Leu-309, and His-397 was confirmed by in vitro assays. In summary, variations in the volume of agonists lead to adaptations in the volume of the LBP and alternative contacts of anchoring OH-groups.

Amino Acids↗

On the conformation of the COOH-terminal domain of the large mechanosensitive channel MscL.

COOH-terminal (S3) domains are conserved within the MscL family of bacterial mechanosensitive channels, but their function remains unclear. The X-ray structure of MscL from Mycobacterium tuberculosis (TbMscL) revealed cytoplasmic domains forming a pentameric bundle (Chang, G., R.H. Spencer, A.T. Lee, M.T. Barclay, and D.C. Rees. 1998. SCIENCE: 282:2220-2226). The helices, however, have an unusual orientation in which hydrophobic sidechains face outside while charged residues face inside, possibly due to specific crystallization conditions. Based on the structure of pentameric cartilage protein, we modeled the COOH-terminal region of E. coli MscL to better satisfy the hydrophobicity criteria, with sidechains of conserved aliphatic residues all inside the bundle. Molecular dynamic simulations predicted higher stability for this conformation compared with one modeled after the crystal structure of TbMscL, and suggested distances for disulfide trapping experiments. The single cysteine mutants L121C and I125C formed dimers under ambient conditions and more so in the presence of an oxidant. The double-cysteine mutants, L121C/L122C and L128C/L129C, often cross-link into tetrameric and pentameric structures, consistent with the new model. Patch-clamp examination of these double mutants under moderately oxidizing or reducing conditions indicated that the bundle cross-linking neither prevents the channel from opening nor changes thermodynamic parameters of gating. Destabilization of the bundle by replacing conservative leucines with small polar residues, or complete removal of COOH-terminal domain (Delta110-136 mutation), increased the occupancy of subconducting states but did not change gating parameters substantially. The Delta110-136 truncation mutant was functional in in vivo osmotic shock assays; however, the amount of ATP released into the shock medium was considerably larger than in controls. The data strongly suggest that in contrast to previous gating models (Sukharev, S., M. Betanzos, C.S. Chiang, and H.R. Guy. 2001a. NATURE: 409:720-724.), S3 domains are stably associated in both closed and open conformations. The bundle-like assembly of cytoplasmic helices provides stability to the open conformation, and may function as a size-exclusion filter at the cytoplasmic entrance to the MscL pore, preventing loss of essential metabolites.

Amino Acid Sequence↗

Structure of the 1,4-bis(2'-deoxyadenosin-N6-yl)-2R,3R-butanediol cross-link arising from alkylation of the human N-ras codon 61 by butadiene diepoxide.

The solution structure of the 1,4-bis(2'-deoxyadenosin-N(6)-yl)-2R,3R-butanediol cross-link arising from N(6)-dA alkylation of nearest-neighbor adenines by butadiene diepoxide (BDO(2)) was determined in the oligodeoxynucleotide 5'-d(CGGACXYGAAG)-3'.5'-d(CTTCTTGTCCG)-3'. This oligodeoxynucleotide contained codon 61 (underlined) of the human N-ras protooncogene. The cross-link was accommodated in the major groove of duplex DNA. At the 5'-side of the cross-link there was a break in Watson-Crick base pairing at base pair X(6).T(17), whereas at the 3'-side of the cross-link at base pair Y(7).T(16), base pairing was intact. Molecular dynamics calculations carried out using a simulated annealing protocol, and restrained by a combination of 338 interproton distance restraints obtained from (1)H NOESY data and 151 torsion angle restraints obtained from (1)H and (31)P COSY data, yielded ensembles of structures with good convergence. Helicoidal analysis indicated an increase in base pair opening at base pair X(6).T(17), accompanied by a shift in the phosphodiester backbone torsion angle beta P5'-O5'-C5'-C4' at nucleotide X(6). The rMD calculations predicted that the DNA helix was not significantly bent by the presence of the four-carbon cross-link. This was corroborated by gel mobility assays of multimers containing nonhydroxylated four-carbon N(6),N(6)-dA cross-links, which did not predict DNA bending. The rMD calculations suggested the presence of hydrogen bonding between the hydroxyl group located on the beta-carbon of the four-carbon cross-link and T(17) O(4), which perhaps stabilized the base pair opening at X(6).T(17) and protected the T(17) imino proton from solvent exchange. The opening of base pair X(6).T(17) altered base stacking patterns at the cross-link site and induced slight unwinding of the DNA duplex. The structural data are interpreted in terms of biochemical data suggesting that this cross-link is bypassed by a variety of DNA polymerases, yet is significantly mutagenic [Kanuri, M., Nechev, L. V., Tamura, P. J., Harris, C. M., Harris, T. M., and Lloyd, R. S. (2002) Chem. Res. Toxicol. 15, 1572-1580].

Alkylating Agents↗

A dynamical model of kinesin-microtubule motility assays.

A two-dimensional stochastic model for the dynamics of microtubules in gliding-assay experiments is presented here, which includes the viscous drag acting on the moving fiber and the interaction with the kinesins. For this purpose, we model kinesin as a spring, and explicitly use parameter values to characterize the model from experimental data. We numerically compute the mean attachment lifetimes of all motors, the total force exerted on the microtubules at all times, the effects of a distribution in the motor speeds, and also the mean velocity of a microtubule in a gliding assay. We find quantitative agreement with the results of J. Howard, A. J. Hudspeth, and R. D. Vale, Nature. 342:154-158. We perform additional numerical analysis of the individual motors, and show how cancellation of the forces exerted by the many motors creates a resultant longitudinal force much smaller than the maximum force that could be exerted by a single motor. We also examine the effects of inhomogeneities in the motor-speeds. Finally, we present a simple theoretical model for microtubules dynamics in gliding assays. We show that the model can be analytically solved in the limit of few motors attached to the microtubule and in the opposite limit of high motor density. We find that the speed of the microtubule goes like the mean speed of the motors in good quantitative agreement with the experimental and numerical results.

Biomechanical Phenomena↗

Novel carbohydrate binding site recognizing blood group A and B determinants in a hybrid of cholera toxin and Escherichia coli heat-labile enterotoxin B-subunits.

The B-subunits of cholera toxin (CTB) and Escherichia coli heat-labile enterotoxin (LTB) are structurally and functionally related. However, the carbohydrate binding specificities of the two proteins differ. While both CTB and LTB bind to the GM1 ganglioside, LTB also binds to N-acetyllactosamine-terminated glycoconjugates. The structural basis of the differences in carbohydrate recognition has been investigated by a systematic exchange of amino acids between LTB and CTB. Thereby, a CTB/LTB hybrid with a gain-of-function mutation resulting in recognition of blood group A and B determinants was obtained. Glycosphingolipid binding assays showed a specific binding of this hybrid B-subunit, but not CTB or LTB, to slowly migrating non-acid glycosphingolipids of human and animal small intestinal epithelium. A binding-active glycosphingolipid isolated from cat intestinal epithelium was characterized by mass spectrometry and proton NMR as GalNAcalpha3(Fucalpha2)Galbeta4(Fucalpha3)Glc NAcbeta3Galbeta4Glc NAcbeta3Galbeta4Glcbeta1Cer. Comparison with reference glycosphingolipids showed that the minimum binding epitope recognized by the CTB/LTB hybrid was Galalpha3(Fucalpha2)Galbeta4(Fucalpha3)GlcNAc beta. The blood group A and B determinants bind to a novel carbohydrate binding site located at the top of the B-subunit interfaces, distinct from the GM1 binding site, as found by docking and molecular dynamics simulations.

ABO Blood-Group System↗

Active-site arginines differentially control Cas12a DNA cleavage and specificity.

Cas12a is a CRISPR-Cas nuclease with biochemical features that make it useful for genome editing and nucleic acid diagnostics. However, its off-target and non-specific trans and CRISPR RNA-independent DNA cleavages can reduce the accuracy and limit applications requiring high fidelity. Here, we analyzed the role of two conserved arginine residues, R918 and R921, found in the RuvC active site pocket of Francisella novicida Cas12a. Through amino acid substitutions, biochemical assays, kinetic analysis, and computational study, we establish that a positive charge at 921 is required for CRISPR RNA-dependent DNA cleavage (cis cleavage), whereas R918 primarily enhances cleavage efficiency. Replacing R918 with lysine or alanine eliminates trans activity while retaining cis cleavage, whereas replacing R921 with lysine eliminates trans activity and replacing with alanine abolishes cis and trans cleavages. Furthermore, these changes significantly decrease RNA-independent cleavage and improve mismatch discrimination during cis cleavage, especially at PAM-distal sites. Structural analysis shows that R918 assists in the conversion of the lid covering the RuvC active site to an alpha helical form, while R921 stabilizes the DNA in the active site. Molecular dynamics simulations reveal that while R921 is critical in supporting the positioning of scissile phosphate, R918 is essential in maintaining catalytic-site organization through lid's conformational change as well as in positioning DNA through its role in stabilizing the active site framework. Together, our results highlight the importance of R918 and R921 in Cas12a's activity and the potential of modifying active pocket residues to reduce unwanted DNA cleavage while increasing on-target specificity.

CRISPR-Cas↗

Monosaccharide-linked inhibitors of O(6)-methylguanine-DNA methyltransferase (MGMT): synthesis, molecular modeling, and structure-activity relationships.

A series of potential inhibitors of the human DNA repair protein O(6)-methylguanine-DNA methyltransferase (MGMT) were synthesized, characterized in detail by NMR, and tested for their ability to deplete MGMT activity in vitro. The new compounds, omega-[O(6)-R-guan-9-yl]-(CH(2))(n)-beta-d-glucosides with R = benzyl or 4-bromothenyl and omega = n = 2, 4,. 12, were compared with the established inhibitors O(6)-benzylguanine (O(6)-BG), 8-aza-O(6)-benzylguanine (8-aza-BG), and O(6)-(4-bromothenyl)guanine (4-BTG), which exhibit in an in vitro assay IC(50) values of 0.62, 0.038, and 0.009 microM, respectively. Potential advantages of the glucosides are improved water solubility and selective uptake in tumor cells. The 4-BTG glucosides with n = 2, 4, 6 show moderate inhibition with an IC(50) of ca. 0.5 microM, while glucosides derived from BG and 8-aza-BG showed significantly poorer inhibition compared to the parent compounds. The 4-BTG glucosides with n = 8, 10, 12 were effective inhibitors with IC(50) values of ca. 0.03 microM. To understand this behavior, extensive molecular modeling studies were performed using the published crystal structure of MGMT (PDB entry: ). The inhibitor molecules were docked into the BG binding pocket, and molecular dynamics simulations with explicit water molecules were carried out. Stabilization energies for the interactions of specific regions of the inhibitor and individual amino acid residues were calculated. The alkyl spacer is located in a cleft along helix 6 of MGMT. With increasing spacer length there is increasing interaction with several amino acid residues which play an important role in the proposed nucleotide flipping mechanism required for DNA repair.

Cell-Free System↗