SEARCH · PubMed Health

Results for “Molecular docking and molecular dynamics simulation”

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

PubMed2023

Discussion on the mechanism of Lingguizhugan Decoction in treating hypertension based on network pharmacology and molecular simulation technology.

To explore the mechanism of Lingguizhugan Decoction in treating hypertension based on network pharmacology and molecular simulation. The active ingredients and potential targets were screened by the Systematic Pharmacological Analysis Platform of Traditional Chinese Medicine (TCMSP). Hypertension-related targets were obtained from OMIM and GeneCards databases. Common targets between drug and hypertension were screened in the Venny platform. A protein-protein interaction (PPI) network was constructed in the STRING database using intersection targets. Key targets in PPI network were analyzed by Cytoscape. R language program was used for Gene Ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Finally, the binding abilities of the main active ingredients to critical targets were verified by molecular simulation. Naringenin, quercetin, kaempferol, and β-sitosterol in Lingguizhugan Decoction, and potential targets such as STAT3, AKT1, TNF, IL6, JUN, PTGS2, MMP9, CASP3, TP53, and MAPK3, were screened out. KEGG Enrichment analysis revealed that the common targets of Lingguizhugan Decoction and hypertension are mainly involved in the lipid and atherosclerosis signaling pathway, AGE-RAGE signaling pathway in diabetic complications, fluid shear stress and atherosclerosis, and IL17 signaling pathway. The molecular simulation results showed that naringenin-MAPK3, quercetin-MMP9, quercetin-PTGS2, and quercetin-TP53 were the top four in the docking scores. Naringenin-MAPK3 and quercetin-MMP9 were stable, with binding free energies of -27.97 ± 1.41 kcal/mol and -21.15 ± 3.17 kcal/mol, respectively. The possible mechanism of Lingguizhugan Decoction in treating hypertension is characterized of multi-component, multi-target, and multi-pathway.Communicated by Ramaswamy H. Sarma.

Network Pharmacology↗

PubMed2006

N-methylpurine DNA glycosylase and 8-oxoguanine dna glycosylase metabolize the antiviral nucleoside 2-bromo-5,6-dichloro-1-(beta-D-ribofuranosyl)benzimidazole.

The rapid in vivo degradation of the potent human cytomegalovirus inhibitor 2-bromo-5,6-dichloro-1-(beta-D-ribofuranosyl)benzimidazole (BDCRB) compared with a structural L-analog, maribavir (5,6-dichloro-2-(isopropylamino)-1-beta-L-ribofuranosyl-1H-benzimidazole), has been attributed to selective glycosidic bond cleavage. An enzyme responsible for this selective BDCRB degradation, however, has not been identified. Here, we report the identification of two enzymes, 8-oxoguanine DNA glycosylase (OGG1) and N-methylpurine DNA glycosylase (MPG), that catalyze N-glycosidic bond cleavage of BDCRB and its 2-chloro homolog, 2,5,6-trichloro-1-(beta-D-ribofuranosyl)benzimidazole, but not maribavir. To our knowledge, this is the first demonstration that free nucleosides are substrates of OGG1 and MPG. To understand how these enzymes might process BDCRB, docking and molecular dynamics simulations were performed with the native human OGG1 crystal coordinates. These studies showed that OGG1 was not able to bind a negative control, guanosine, yet BDCRB and maribavir were stabilized through interactions with various binding site residues, including Phe319, His270, Ser320, and Asn149. Only BDCRB, however, achieved orientations whereby its anomeric carbon, C1', could undergo nucleophilic attack by the putative catalytic residue, Lys249. Thus, in silico observations were in perfect agreement with experimental observations. These findings implicate DNA glycosylases in drug metabolism.

Antiviral Agents↗

PubMed2005

Orientational dynamics and dye-DNA interactions in a dye-labeled DNA aptamer.

We report the picosecond and nanosecond timescale rotational dynamics of a dye-labeled DNA oligonucleotide or "aptamer" designed to bind specifically to immunoglobulin E. Rotational dynamics in combination with fluorescence lifetime measurements provide information about dye-DNA interactions. Comparison of Texas Red (TR), fluorescein, and tetramethylrhodamine (TAMRA)-labeled aptamers reveals surprising differences with significant implications for biophysical studies employing such conjugates. Time-resolved anisotropy studies demonstrate that the TR- and TAMRA-aptamer anisotropy decays are dominated by the overall rotation of the aptamer, whereas the fluorescein-aptamer anisotropy decay displays a subnanosecond rotational correlation time much shorter than that expected for the overall rotation of the aptamer. Docking and molecular dynamics simulations suggest that the low mobility of TR is a result of binding in the groove of the DNA helix. Additionally, associated anisotropy analysis of the TAMRA-aptamer reveals both quenched and unquenched states that experience significant coupling to the DNA motion. Therefore, quenching of TAMRA by guanosine must depend on the configuration of the dye bound to the DNA. The strong coupling of TR to the rotational dynamics of the DNA aptamer, together with the absence of quenching of its fluorescence by DNA, makes it a good probe of DNA orientational dynamics. The understanding of the nature of dye-DNA interactions provides the basis for the development of bioconjugates optimized for specific biophysical measurements and is important for the sensitivity of anisotropy-based DNA-protein interaction studies employing such conjugates.

Anisotropy↗

PubMed2026

Genome-Wide Analysis of Triticum aestivum Root Meristem Growth Factor (RGF) Gene Family Highlights TaRGF5 as a Putative Component of Root-Associated Signaling.

Wheat (Triticum aestivum), a key global crop, faces rising drought stress that limits root growth and water uptake. Root meristem growth factors (RGFs) are small peptides that regulate root stem cell maintenance, meristem activity, and lateral root formation in model plants, yet the RGF gene family remains unexplored in wheat. Here, we performed a comprehensive genome-wide analysis of the TaRGF gene family, identifying 15 genes distributed across the A, B, and D subgenomes and classified into five homeologous groups (TaRGF1-TaRGF5), predominantly located on chromosomes 2 and 6. All TaRGFs contained a characteristic RGF motif, with dibasic cleavage sites and Asp-Tyr motifs indicating conserved maturation mechanisms. Based on the phylogenetic analysis, the TaRGF5 homeologs showed the highest similarity to Arabidopsis thaliana RGF5. Tested RNA-seq data revealed predominantly root-enriched expression for all TaRGF genes, with TaRGF5 exhibiting the most root-preferential and downregulation under drought stress. Quantitative real-time PCR (qRT-PCR) confirmed that drought stress suppressed the expression of TaRGF5A, TaRGF5B, and TaRGF5D in roots of wheat cultivar Sids-13 across all time points, unlike the higher accumulation seen in controls. Promoter analysis predicted a unique BES1 transcription factor binding site exclusively in TaRGF5B, linking brassinosteroid signaling to peptide-mediated root regulation. Structural modeling and molecular docking predicted an interaction between wheat TaRGF5 homeologs and root growth factor-insensitive receptor kinase (TaRGI3), characterized by conserved sulfotyrosine-mediated binding and favorable interaction energetics. Based on this characterization of the wheat RGF gene family, particularly the potential role of TaRGF5 in root development and drought-adaptation signaling, we propose targeting this gene for functional analysis to improve wheat resilience under water-limited conditions.

Triticum↗

PubMed2026

CCT2 defines a highly cisplatin-resistant and poor-prognosis subtype of lung adenocarcinoma.

Cisplatin-based chemotherapy is a standard treatment for lung adenocarcinoma (LUAD), yet acquired cisplatin resistance remains a marked cause of treatment failure. The molecular mechanisms driving cisplatin resistance in LUAD have not been fully elucidated. The present study integrated bulk transcriptomic data, genomic mutation profiles and single-cell RNA sequencing data to systematically investigate cisplatin resistance in LUAD. Resistance-associated genes were identified through differential expression, survival analysis and database integration. Unsupervised clustering was used to define cisplatin resistance-associated subtypes. Functional characteristics were explored using pathway enrichment, immune infiltration, tumor mutation burden and weighted gene co-expression network analysis. A machine learning framework incorporating 101 algorithms was applied to identify key genes and construct a prognostic model. Single-cell analyses and in vitro experiments were performed to validate the biological role of the core gene. Molecular docking and molecular dynamics simulations were conducted to identify potential therapeutic compounds. A total of two molecular subtypes with distinct cisplatin resistance levels and prognostic outcomes were identified. The high-resistance subtype exhibited enhanced cell cycle activity, DNA repair signaling and immune heterogeneity. Machine learning analysis revealed a five-gene signature, with chaperonin-containing TCP1 subunit 2 (CCT2) emerging as a key regulator of cisplatin resistance. Single-cell analyses showed that CCT2 was predominantly enriched in resistant epithelial cell subpopulations. Functional experiments demonstrated that CCT2 knockdown significantly inhibited cell proliferation and enhanced cisplatin sensitivity in LUAD cell lines. A number of candidate compounds targeting CCT2 exhibited stable binding in silico. The present findings identified CCT2 as a key mediator of cisplatin resistance in LUAD and provided potential therapeutic strategies to overcome chemotherapy resistance.

chaperonin-containing TCP-1 subunit 2↗

PubMed2004

Loop flexibility and solvent dynamics as determinants for the selective inhibition of cyclin-dependent kinase 4: comparative molecular dynamics simulation studies of CDK2 and CDK4.

The design and discovery of selective cyclin-dependent kinase 4 (CDK4) inhibitors have been actively pursued in order to develop therapeutic cancer treatments. By means of a consecutive computational protocol involving homology modeling, docking experiments, and molecular dynamics simulations, we examine the characteristic structural and dynamic properties that distinguish CDK4 from CDK2 in its complexation with selective inhibitors. The results for all three CDK4-selective inhibitors under investigation show that the large-amplitude motion of a disordered loop of CDK4 is damped out in the presence of the inhibitors whereas their binding in the CDK2 active site has little effect on the loop flexibility. It is also found that the binding preference of CDK4- selective inhibitors for CDK4 over CDK2 stems from the reduced solvent accessibility in the active site of the former due to the formation of a stable hydrogen-bond triad by the Asp99, Arg101, and Thr102 side chains at the top of the active-site gorge. Besides the differences in loop flexibility and solvent accessibility, the dynamic stabilities of the hydrogen bonds between the inhibitors and the side chain of the lysine residue at the bottom of the active site also correlate well with the relative binding affinities of the inhibitors for the two CDKs. These results highlight the usefulness of this computational approach in evaluating the selectivity of a CDK inhibitor, and demonstrate the necessity of considering protein flexibility and solvent effects in designing new selective CDK4-selective inhibitors.

Amino Acid Sequence↗

PubMed2025

Integrated network pharmacology, molecular docking and experimental validation to investigate the mechanism of tannic acid in nasopharyngeal cancer.

Tannic acid (TA) is the primary bioactive component in the gallnut (Galla chinensis) and has exhibited the anticancer effects. However, the mechanism of its anti-cancer activity in nasopharyngeal carcinoma (NPC) remains unclear. This research aims to explore the underlying mechanism of TA in the treatment of nasopharyngeal cancer using network pharmacology, molecular docking and experimental validation. Firstly, the targets of TA and NPC were predicted and collected through databases, and the intersection targets were identified. Subsequently, protein-protein interaction (PPI) network analysis, Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes Genomes (KEGG) pathway enrichment analysis, molecular docking and molecular dynamics (MD) simulation were conducted to uncover the potential mechanisms of TA in treatment of NPC. Finally, in vitro experiments were utilized to verify the mechanism of TA with anticancer activity in NPC. The results of network pharmacology revealed 42 intersection targets between NPC-related targets and TA-related targets. The phosphoinositide 3-kinase (PI3K)/protein kinase B (AKT) signaling was identified as the main target pathway of TA against NPC. Additionally, molecular docking and MD simulation confirmed the closely binding affinities of TA with AKT1. Furthermore, the results of in vitro experiments demonstrated that TA exerts anticancer activity against NPC by targeting the PI3K/AKT signaling pathway, leading to the suppression of cell proliferation. TA is a promising therapeutic candidate for NPC through PI3K/AKT signaling pathway. These results provide insights into the clinical application of TA, particularly when considered in combination with other therapeutic modalities.

Molecular Docking Simulation↗

PubMed2006

CCR5 interactions with the variable 3 loop of gp120.

The G-protein coupled receptor CCR5 functions pathologically as the primary co-receptor for macrophage tropic (R5) strains of HIV-1. The interactions responsible for co-receptor activity are unknown. Molecular-dynamics simulations of the extracellular and adjacent transmembrane domains of CCR5 were performed with explicit solvation utilizing a rhodopsin-based homology model. The functional unit of co-receptor binding was constructed via docking and molecular-dynamics simulation of CCR5 and the variable 3 loop of gp120, which is a dominant determinant of co-receptor utilization. The variable 3 loop was demonstrated to interact primarily with the amino terminus and the second extracellular loop of CCR5, providing novel structural information regarding the co-receptor-binding site. Alanine mutants that alter chemokine binding and co-receptor activity were examined. Molecular-dynamics simulations with and without the variable 3 loop of gp120 were able to rationalize the activities of these mutants successfully, providing support for the proposed model. Based on these results, the global complex of CCR5, gp120 including the V3 loop and CD4, was investigated. The utilization of computational analysis, in combination with molecular biological data, provides a powerful approach for understanding the use of CCR5 as a co-receptor by HIV-1.

Alanine↗

PubMed2004

A structural and dynamic model for the interaction of interleukin-8 and glycosaminoglycans: support from isothermal fluorescence titrations.

Binding of interleukin-8 (IL-8) to glycosaminoglycans (GAGs) on the surface of endothelial cells is crucial for the recruitment of neutrophils to an inflammatory site. Deriving structural knowledge about this interaction from in silico docking experiments has proved difficult because of the high flexibility and the size of GAGs. Therefore, we developed a docking method that takes into account ligand and protein flexibility by running approximately 15,000 molecular dynamics simulations of the docking event with different initial orientations of the binding partners. The method was shown to successfully reproduce the residues of basic fibroblast growth factor involved in GAG binding. Docking of a heparin hexasaccharide to IL-8 gave an interaction interface involving the basic residues His18, Lys20, Arg60, Lys64, Lys67, and Arg68. By subjecting IL-8 single-site mutants, in which these amino acids were replaced by alanine, to isothermal fluorescence titrations, the affinities for heparin were determined to be wtIL-8 > IL-8(H18A) >> IL-8(R68A) > IL-8(K67A) >> IL-8(K20A) > IL-8(R60A) >> IL-8(K64A). A comparison with the binding energies calculated from the model revealed high values for wtIL-8 and the H18A mutant and significantly lower but similar energies for the remaining mutants. Connecting the two fully sulfated hexasaccharides bound to each of the two IL-8 monomers in the dimeric chemokine by an N-acetylated dodecasaccharide gave a complex structure in which the GAG molecule aligned in a parallel fashion to the N-terminal alpha-helices of IL-8 like a horseshoe. A 5-ns molecular dynamics simulation of this complex confirmed its structural stability and revealed a reorientation in both binding sites where a disaccharide became the central binding unit. Isothermal fluorescence titration experiments using differently sulfated heparin disaccharides confirmed that a single disaccharide can indeed bind IL-8 with high affinity.

Amino Acid Substitution↗

PubMed2005

Dynamic mechanism for the autophosphorylation of CheA histidine kinase: molecular dynamics simulations.

The two-component system (TCS) is an important signal transduction component for most bacteria. This signaling pathway is mediated by histidine kinases via autophosphorylation between P1 and P4 domains. Taking chemotaxis protein CheA as a model of TCS, the autophosphorylation mechanism of the TCS histidine kinases has been investigated in this study by using a computational approach integrated homology modeling, ligand-protein docking, protein-protein docking, and molecular dynamics (MD) simulations. Four nanosecond-scale MD simulations were performed on the free P4 domain, P4-ATP, P4-TNPATP, and P1-P4-ATP complexes, respectively. Upon its binding to the binding pocket of P4 with a folded conformation, ATP gradually extends to an open state with help from a water molecule. Meanwhile, ATP forms two hydrogen bonds with His413 and Lys494 at this state. Because of the lower energy of the folded conformations, ATP shrinks back to its folded conformations, leading to the rupture of the hydrogen bond between ATP and Lys494. Consequently, Lys494 moves away from the pocket entrance, resulting in an open of the ATP lid of P4. It is the open state of P4 that can bind tightly to P1, where the His45 of P1 occupies a favorable position for its autophosphorylation from ATP. This indicates that ATP is not only a phosphoryl group donor but also an activator for CheA phosphorylation. Accordingly, a mechanism of the autophosphorylation of CheA is proposed as that the ATP conformational switch triggers the opening of the ATP lid of P4, leading to P1 binding tightly, and subsequently autophosphorylation from ATP to P1.

Adenosine Triphosphate↗

PubMed2004

A computational study of the binding of propidium to the peripheral anionic site of human acetylcholinesterase.

Combined docking and molecular dynamics (MD) simulations were carried out in order to investigate the binding mode of propidium at the human acetylcholinesterase (HuAChE) peripheral site. Two different docking protocols followed by cluster analyses were performed, allowing the identification of five high-populated and low-energy configuration families. To dynamically explore the behavior of the ligand at the peripheral HuAChE binding site, six complexes (five low-energy and one high-energy) were submitted to 2.5 ns of MD simulations. The representative propidium/HuAChE binding modes were chosen on the basis of both the docking energy score and the dynamic stability of the complexes throughout the MD simulations. The most stable poses of propidium at HuAChE PAS were similar to those experimentally determined with the murine enzyme. We therefore suggest that the present modeling protocol might be used in the dynamic investigation of the interactions of a small-molecule inhibitor with a surface-like binding site of a protein. Finally, because of the biological relevance of the target studied here, the present results can be of interest for the rational design of molecules potentially useful in the treatment of the Alzheimer's disease.

Acetylcholinesterase↗

PubMed2005

Regiospecificity of human cytochrome P450 1A1-mediated oxidations: the role of steric effects.

Cytochrome P450 1A1 oxidizes a diverse range of substrates, including the procarcinogenic xenobiotic benzo[a]pyrene (B[a]P) and endogenous fatty acid precursors of prostaglandins, such as arachidonic acid (AA) and eicosapentaenoic acid (EA). We have investigated the extent to which enzyme-substrate interactions govern regio- and stereoselectivity of oxidation of these compounds by using docking and molecular dynamics (MD) simulations to examine the likelihood of substrate oxidation at various sites. Due to structural differences between the substrates analyzed, B[a]P and its diols (planar, rigid), and the fatty acids AA and EA (long, flexible), different docking strategies were required. B[a]P, B[a]P-7,8-diols, (+) 7S,8S- and (-) 7R,8R-diols, were docked into the active site of a homology model of P450 1A1 using an automated routine, Affinity (Accelrys, San Diego, CA). AA and EA, on the other hand, required a series of restrained MD simulations to obtain a variety of productive binding modes. All complexes were evaluated by MD-based in silico site scoring to predict product profiles based on certain geometric criteria, such as angle and distance of a given substrate atom from the ferryl oxygen. For all substrates studied, the in vitro profiles were generally reflected by the in silico scores, which suggests that steric factors play a key role in determining regiospecificity in P450 1A1-mediated oxidations. We have also shown that molecular dynamics simulations may be very useful in determination of product profiles for structurally diverse substrates of P450 enzymes.

Binding Sites↗

PubMed2004

Molecular docking and 3D QSAR studies on 1-amino-2-phenyl-4-(piperidin-1-yl)-butanes based on the structural modeling of human CCR5 receptor.

In the present study, we have used an approach combining protein structure modeling, molecular dynamics (MD) simulation, automated docking, and 3D QSAR analyses to investigate the detailed interactions of CCR5 with their antagonists. Homology modeling and MD simulation were used to build the 3D model of CCR5 receptor based on the high-resolution X-ray structure of bovine rhodopsin. A series of 64 CCR5 antagonists, 1-amino-2-phenyl-4-(piperidin-1-yl)-butanes, were docked into the putative binding site of the 3D model of CCR5 using the docking method, and the probable interaction model between CCR5 and the antagonists were obtained. The predicted binding affinities of the antagonists to CCR5 correlate well with the antagonist activities, and the interaction model could be used to explain many mutagenesis results. All these indicate that the 3D model of antagonist-CCR5 interaction is reliable. Based on the binding conformations and their alignment inside the binding pocket of CCR5, three-dimensional structure-activity relationship (3D QSAR) analyses were performed on these antagonists using comparative molecular field analysis (CoMFA) and comparative molecular similarity analysis (CoMSIA) methods. Both CoMFA and CoMSIA provide statistically valid models with good correlation and predictive power. The q(2)(r(cross)(2)) values are 0.568 and 0.587 for CoMFA and CoMSIA, respectively. The predictive ability of these models was validated by six compounds that were not included in the training set. Mapping these models back to the topology of the active site of CCR5 leads to a better understanding of antagonist-CCR5 interaction. These results suggest that the 3D model of CCR5 can be used in structure-based drug design and the 3D QSAR models provide clear guidelines and accurate activity predictions for novel antagonist design.

Binding Sites↗

PubMed2006

Essential structural profile of a dual functional inhibitor against cyclooxygenase-2 (COX-2) and 5-lipoxygenase (5-LOX): molecular docking and 3D-QSAR analyses on DHDMBF analogues.

It is recently proposed that compounds with equal capabilities of inhibiting COX and 5-LOX, both are key enzymes involved in the arachidonic acid (AA) cascade, are expected to be safer non-steroidal anti-inflammatory drugs (NSAIDs). To dig out helpful information in designing dual functional inhibitors against the two enzymes, homology modeling, molecular dynamics (MD) simulations, automated docking, and 3D-QSAR analyses were performed in this study on 21 COX-2/5-LOX dual inhibitors, namely, 7-tert-butyl-2,3-dihydro-3,3-dimethylbenzofuran (DHDMBF) analogues. A 3D-model of 5-LOX was built based on the high-resolution X-ray structure of rabbit reticulocyte 15-lipoxygenase. Molecular docking was then applied to locate the binding orientations and conformations of DHDMBF analogues with COX-2 and 5-LOX, respectively, leading to highly predictive CoMFA models constructed on the basis of the binding conformations with q2 values of 0.782 and 0.634 for COX-2 and 5-LOX, respectively. In addition, CoMFA field distributions were found in good agreement with the structural characteristics of the corresponding binding sites. Both the docking simulations and QSAR analyses suggest that new potent dual inhibitors should share a structural feature with a moderately bulky group at R2 position and a rather negatively charged group around the position of the carbonyl group of DHDMBFs. Therefore, the final 3D-QSAR models and the information of the inhibitor-enzyme interaction should be useful in developing new NSAIDs as anti-inflammation drugs with favorable safety profile.

Amino Acid Sequence↗

PubMed2003

An all-atom solution-equilibrated model for human extrinsic blood coagulation complex (sTF-VIIa-Xa): a protein-protein docking and molecular dynamics refinement study.

Tissue factor (TF)-bound factor (F)VIIa plays a critical role in activating FX, an event that rapidly results in blood coagulation. Despite recent advances in the structural information about soluble TF (sTF)-bound VIIa and Xa individually, the atomic details of the ternary complex are not known. As part of our long-term goal to provide a structural understanding of the extrinsic blood coagulation pathway, we built an all atom solution-equilibrated model of the human sTF-VIIa-Xa ternary complex using protein-protein docking and molecular dynamics (MD) simulations. The starting structural coordinates of sTF-VIIa and Xa were derived from dynamically equilibrated solution structures. Due to the flexible nature of the light-chain of the Xa molecule, a three-stage docking approach was employed in which SP (Arg195-Lys448)/EGF2 (Arg86-Arg139), EGF1 (Asp46-Thr85) and GLA (Ala1-Lys45) domains were docked in a sequential manner. The rigid-body docking approach of the FTDOCK method in conjunction with filtering based on biochemical knowledge from experimental site-specific mutagenesis studies provided the strategy. The best complex obtained from the docking experiments was further refined using MD simulations for 3 ns in explicit water. In addition to explaining most of the known experimental site-specific mutagenesis data pertaining to sTF-VIIa, our model also characterizes likely enzyme-binding exosites on FVIIa and Xa that may be involved in the ternary complex formation. According to the equilibrated model, the 140s loop of VIIa serves as the key recognition motif for complex formation. Stable interactions occur between the FVIIa 140s loop and the FXa -strand B2 region near the sodium-binding domain, the 160 s loop and the N-terminal activation loop regions. The helical-hydrophobic stack region that connects the GLA and EGF1 domains of VIIa and Xa appears to play a potential role in the membrane binding region of the ternary complex. The proposed model may serve as a reasonable structural basis for understanding the exosite-mediated substrate recognition of sTF-VIIa and to advance understanding of the TFPI-mediated regulatory pathway of the extrinsic blood coagulation cascade.

Binding Sites↗

PubMed2005

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↗

PubMed1999

Ligand binding by antibody IgE Lb4: assessment of binding site preferences using microcalorimetry, docking, and free energy simulations.

Antibody IgE Lb4 interacts favorably with a large number of different compounds. To improve the current understanding of the structural basis of this vast cross-reactivity, the binding of three dinitrophenyl (DNP) amino acids (DNP-alanine, DNP-glycine, and DNP-serine) is investigated in detail by means of docking and molecular dynamics free energy simulations. Experimental binding energies obtained by isothermal titration microcalorimetry are used to judge the results of the computational studies. For all three ligands, the docking procedure proposes two plausible subsites within the binding region formed by the antibody CDR loops. By subsequent molecular dynamics simulations and calculations of relative free energies of binding, one of these subsites, a tyrosine-surrounded pocket, is revealed as the preferred point of complexation. For this subsite, results consistent with experimental observations are obtained; DNP-glycine is found to bind better than DNP-serine, and this, in turn, is found to bind better than DNP-alanine. The suggested binding mode makes it possible to explain both the moderate binding affinity and the differences in binding energy among the three ligands.

Animals↗

PubMed2023

Structural and functional insights into a novel homozygous missense pathogenic variant in CUL7 identified in consanguineous Pakistani family.

3M syndrome is a rare genetic familial disorder characterized by short stature, growth retardation, facial dysmorphism, skeletal abnormalities, fleshy protruding heels, and normal intelligence, caused by mutations in the CUL7, OBSL1 and CCDC8 genes. In the present study, a novel homozygous missense variant of CUL7 (NP_001161842.1, c.4493T > C, p.L1498P) has been identified in a consanguineous Pakistani family by whole exome sequencing. In silico structural evaluation, molecular docking and simulation studies of mutant CUL7 provides substantial evidence about its crucial role in the progression of discussed ailment. The newly discovered variant significantly altered the protein's three dimensional structure, leading to abnormal interaction with binding proteins. This computational and experimental investigation provides useful information to drug developers for the synthesis of novel therapeutics against the discussed ailment.Communicated by Ramaswamy H. Sarma.

Humans↗