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

Inhibitor discovery targeting the intermediate structure of beta-amyloid peptide on the conformational transition pathway: implications in the aggregation mechanism of beta-amyloid peptide.

Abeta peptides cleaved from the amyloid precursor protein are the main components of senile plaques in Alzheimer's disease. Abeta peptides adopt a conformation mixture of random coil, beta-sheet, and alpha-helix in solution, which makes it difficult to design inhibitors based on the 3D structures of Abeta peptides. By targeting the C-terminal beta-sheet region of an Abeta intermediate structure extracted from molecular dynamics simulations of Abeta conformational transition, a new inhibitor that abolishes Abeta fibrillation was discovered using virtual screening in conjunction with thioflavin T fluorescence assay and atomic force microscopy determination. Circular dichroism spectroscopy demonstrated that the binding of the inhibitor increased the beta-sheet content of Abeta peptides either by stabilizing the C-terminal beta-sheet conformation or by inducing the intermolecular beta-sheet formation. It was proposed that the inhibitor prevented fibrillation by blocking interstrand hydrogen bond formation of the pleated beta-sheet structure commonly found in amyloid fibrils. The study not only provided a strategy for inhibitor design based on the flexible structures of amyloid peptides but also revealed some clues to understanding the molecular events involved in Abeta aggregation.

Amyloid↗

Trisubstituted acridines as G-quadruplex telomere targeting agents. Effects of extensions of the 3,6- and 9-side chains on quadruplex binding, telomerase activity, and cell proliferation.

The synthesis is reported of a group of 3,6,9-trisubstituted acridine compounds as telomeric quadruplex-stabilizing ligands with systematic variations at the 3-, 6-, and 9-positions. A new microwave-assisted methodology has been developed for trisubstituted acridine synthesis. Structure-activity relationships are reported using surface plasmon resonance and a fluorescence melting assay to examine quadruplex binding, together with a telomerase inhibition assay. These reveal relationships between G-quadruplex stabilization and telomerase inhibition and optimal 3,6- and 9-substituent side-chain lengths for maximal activity. Qualitative molecular modeling using molecular dynamics simulations has been undertaken on four quadruplex-DNA complexes. Long-term exposure of MCF7 cancer cells to a subset of the most active compounds, at doses lower than the IC(50) values, showed that one compound produced a marked decrease in population growth, accompanied by senescence, which is consistent with telomere targeting by this agent.

Acridines↗

Inhibition of heme detoxification processes underlies the antimalarial activity of terpene isonitrile compounds from marine sponges.

A series of terpene isonitriles, isolated from marine sponges, have previously been shown to exhibit antimalarial activities. Molecular modeling studies employing 3D-QSAR with receptor modeling methodologies performed with these isonitriles showed that the modeled molecules could be used to generate a pharmacophore hypothesis consistent with the experimentally derived biological activities. It was also shown that one of the modeled compounds, diisocyanoadociane (4), as well as axisonitrile-3 (2), both of which have potent antimalarial activity, interacts with heme (FP) by forming a coordination complex with the FP iron. Furthermore, these compounds were shown to inhibit sequestration of FP into beta-hematin and to prevent both the peroxidative and glutathione-mediated destruction of FP under conditions designed to mimic the environment within the malaria parasite. By contrast, two of the modeled diterpene isonitriles, 7-isocyanoamphilecta-11(20),15-diene (12) and 7-isocyano-15-isothiocyanatoamphilecta-11(20)-ene (13), that displayed little antimalarial activity also showed little inhibitory activity in these FP detoxification assays. These studies suggest that the active isonitrile compounds, like the quinoline antimalarials, exert their antiplasmodial activity by preventing FP detoxification. Molecular dynamics simulations performed with diisocyanoadociane (4) and axisonitrile-3 (2) allowed their different binding to FP to be distinguished.

Animals↗

Amino acid modification in the HIV-1 Tat basic domain: insights from molecular dynamics and in vivo functional studies.

Tat is an essential protein of the human immunodeficiency virus type 1 (HIV-1). It activates transcription by specifically binding a stem-loop element in the viral long terminal repeat through its highly basic arginine-rich domain. Conserved lysine residues at positions 50 and 51 inside this domain have been recently reported to be the targets of post-translational modification by acetylation, and mutation of these residues has pointed out its relevance to protein function. In an attempt to shed light on the molecular basis of the functional differences found for Tat mutants we have performed a series of molecular dynamics simulations on wt Tat, Lys50/51 --> Arg50/51, Lys50/51 --> Ala50/51 and acetylated Lys50 from HIV-1 strain Z2. Theoretical results are compared with a homogeneous set of in vivo transactivation assays on the corresponding Tat mutants from the strain B2, which exhibits high structural similarity with Tat from HIV-1 strain Z2. Remarkable correlation is found between the degree of structure conservation and the transactivation capabilities of Tat mutants.

Acetylation↗

Consensus chemistry and beta-turn conformation of the active core of the insect kinin neuropeptide family.

BACKGROUND: Neuropeptides are examples of small, flexible molecules that bind to receptors and induce signal transduction, thereby eliciting biological activity. The multifunctional insect kinin neuropeptides retain full activity when reduced to only their carboxy-terminal pentapeptide (Phe1-X2-X3-Trp4-Gly5-NH2), thereby allowing extensive structure-function studies and conformational analysis. RESULTS: A combined experimental and theoretical analysis of the insect kinin carboxy-terminal pentapeptide was used to probe the role of each residue, define the bioactive conformation, and design a constrained bioactive analog. Coupling receptor-binding data with two biological activity assays allowed receptor binding and signal transduction to be differentiated. A preferred beta-turn conformation, found for residues 1-4 by molecular dynamics simulations, was tested by designing a conformationally restricted cyclic hexapeptide. This cyclic analog showed a preference for the beta-turn conformation, as shown by a conformational search and nuclear magnetic resonance spectroscopy, and it showed stronger receptor binding but decreased activity relative to highly active linear analogs. CONCLUSIONS: Each residue of the insect kinin carboxy-terminal pentapeptide has a distinct role in conformational preference, specific receptor interactions or signal transduction. The beta-turn preference of residues Phe1-X2-X3-Trp4 implicates this as the bioactive conformation. The amidated carboxyl terminus, required for activity in many neuropeptide families, may be generally important for signal transduction and its inclusion may therefore be essential for agonist design.

Animals↗

A Computational Workflow for Prioritizing Microbial Metabolite-Associated Host Genes in Constipation-Predominant Irritable Bowel Syndrome.

No standardized computational pipeline exists for systematically prioritizing microbial metabolite-associated host genes and protein-ligand complexes from publicly available chemical, genomic, and structural databases. This article describes an eight-stage workflow that accepts a user-defined set of gut microbiota-derived metabolites and produces a ranked shortlist of candidate metabolite-associated host genes, enriched biological pathways, and structurally prioritized protein-ligand complexes for experimental follow-up. The pipeline integrates (i) chemoinformatic metabolite profiling; (ii) multi-database candidate target prediction using protein-chemical interaction and ligand-based target-prediction tool and a molecular docking program; (iii) differential gene expression analysis of publicly available transcriptomic data; (iv) target-differentially expressed gene overlap; (v) protein-protein interaction network construction and pathway enrichment; (vi) molecular docking with a molecular docking program; (vii) 200 ns molecular dynamics simulation using a molecular dynamics engine with a protein force field used for molecular dynamics simulations; and (viii) MM-PBSA binding free-energy estimation. As a worked example, nine gut microbiota-derived or microbiota-modified metabolites representing short-chain fatty acids, bile acids, tryptophan-derived metabolites, and urolithin A were processed using the public IBS-C rectal mucosal transcriptomic dataset GSE36701. The workflow ranked 17 unique predicted metabolite-associated genes that were differentially expressed in this dataset. Docking, molecular dynamics simulation, and MM-PBSA analyses structurally prioritized five metabolite-protein complexes: lithocholic acid-VDR, lithocholic acid-NR1H4/FXR, ursodeoxycholic acid-NR1H4/FXR, tryptamine-HTR2A (simulated in an explicit 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer), and urolithin A-CASP3. The protocol is designed to be adaptable to other metabolite sets, disease transcriptomic datasets, and target classes; all outputs are hypothesis-generating computational predictions that require independent transcriptomic replication, protein-level validation, and functional ligand-response assays before causal or therapeutic conclusions can be drawn.

Irritable Bowel Syndrome↗

Sanguinarine as a multi-target therapeutic candidate for laryngeal cancer: insights from network pharmacology, molecular dynamics and in vitro validation.

OBJECTIVE: To identify the core targets and elucidate the potential molecular mechanisms of sanguinarine (SA) against laryngeal squamous cell carcinoma (LSCC), and to validate its antitumor effects in vitro. METHODS: Potential targets of SA were predicted using SwissTargetPrediction, TargetNet, and SuperPred and intersected with LSCC-related targets obtained from the GeneCards, OMIM, and DISEASES databases. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed. A protein-protein interaction (PPI) network was constructed using the STRING database (combined score > 0.900), and topological parameters including degree centrality (DC), betweenness centrality (BC), closeness centrality (CC), eigenvector centrality (EC), and local average connectivity (LAC) were calculated in Cytoscape to identify core genes based on median thresholds. Molecular docking and 100-ns molecular dynamics (MD) simulations were conducted for epidermal growth factor receptor (EGFR), Phosphatidylinositide-3-kinase catalytic subunit alpha (PIK3CA), phosphatidylinositol-4,5-biphosphate 3-kinase catalytic subunit β (PIK3CB), phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit delta (PIK3CD), and Non-Receptor Tyrosine Kinase (SRC). The effects of SA on LSCC were evaluated using CCK-8, colony formation, Transwell migration, and wound-healing assays in TU177 cells and TU212. RESULTS: A total of 213 common targets were identified, which were significantly enriched in PI3K-Akt signaling, EGFR tyrosine kinase inhibitor resistance, and adhesion- and migration-related pathways. The PPI network comprised 259 nodes and 259 edges, from which five core genes-PIK3CA, PIK3CB, PIK3CD, EGFR, and SRC-were identified. Molecular docking revealed strong binding affinities between SA and the PI3K family proteins (- 9.79 to - 10.96 kcal/mol), as well as EGFR (- 8.58 kcal/mol) and SRC (- 6.77 kcal/mol). MD simulations indicated greater stability of SA complexes with EGFR and PI3K family members compared with SRC. In vitro assays demonstrated that SA significantly inhibited TU177 cell and TU212 cell proliferation, colony formation, and migration. CONCLUSION: SA may exert anti-laryngeal cancer effects through synergistic multi-target inhibition centered on the EGFR/SRC/PI3K signaling axis, highlighting its potential as a promising therapeutic candidate for LSCC.

Humans↗

Identification of a novel glucocorticoid receptor mutation in budesonide-resistant human bronchial epithelial cells.

We developed a molecular genetic model to investigate glucocorticoid receptor (GR) signaling in human bronchial epithelial cells in response to the therapeutic steroid budesonide. Based on a genetic selection scheme using the human Chago K1 cell line and integrated copies of a glucocorticoid-responsive herpes simplex virus thymidine kinase gene and a green fluorescent protein gene, we isolated five Chago K1 variants that grew in media containing budesonide and ganciclovir. Three spontaneous budesonide-resistant subclones were found to express low levels of GR, whereas two mutants isolated from ethylmethane sulfonate-treated cultures contained normal levels of GR protein. Analysis of the GR coding sequence in the budesonide-resistant subclone Ch-BdE5 identified a novel Val to Met mutation at amino acid position 575 (GRV575M) which caused an 80% decrease in transcriptional regulatory functions with only a minimal effect on ligand binding activity. Homology modeling of the GR structure in this region of the hormone binding domain and molecular dynamic simulations suggested that the GRV575M mutation would have a decreased affinity for the LXXLL motif of p160 coactivators. To test this prediction, we performed transactivation and glutathione-S-transferase pull-down assays using the p160 coactivator glucocorticoid interacting protein 1 (GRIP1)/transcriptional intermediary factor 2 and found that GRV575M transcriptional activity was not enhanced by GRIP1 in transfected cells nor was it able to bind GRIP1 in vitro. Identification of the novel GRV575M variant in human bronchial epithelial cells using a molecular genetic selection scheme suggests that functional assays performed in relevant cell types could identify subtle defects in GR signaling that contribute to reduced steroid sensitivities in vivo.

Amino Acid Sequence↗

Crystal structure of an Anti-meningococcal subtype P1.4 PorA antibody provides basis for peptide-vaccine design.

In various western countries, subtype P1.4 of Neisseria meningitidis serogroup B causes the greatest incidence of meningococcal disease. To investigate the molecular recognition of this subtype, we crystallised a peptide (P1HVVVNNKVATH(P11)), corresponding to the subtype P1.4 epitope sequence of outer membrane protein PorA, in complex with a Fab fragment of the bactericidal antibody MN20B9.34 directed against this epitope. Structure determination at 1.95 A resolution revealed a unique complex of one P1.4 antigen peptide bound to two identical Fab fragments. One Fab recognises the putative epitope residues in a 2:2 type I beta-turn at residues P5NNKV(P8), whereas the other Fab binds the C-terminal residues of the peptide that we consider a crystallisation artefact. Interestingly, recognition of the P1.4 epitope peptide is mediated almost exclusively through the complementarity-determining regions of the heavy chain. We exploited the observed turn conformation for designing conformationally restricted cyclic peptides for use as a peptide vaccine. The conformational stability of the two peptide designs was assessed by molecular dynamics simulations. Unlike the linear peptide, both cyclic peptides, conjugated to tetanus toxoid as a carrier protein, elicited antibody responses in mice that recognised meningococci of subtype P1.7-2,4. Serum bactericidal assays showed that some, but not all, of the sera induced with the cyclic peptide conjugates could activate the complement system with titres that were very high compared to the titres induced by complete PorA protein in its native conformation administered in outer membrane vesicles.

Animals↗

Altered acyl chain length specificity of Rhizopus delemar lipase through mutagenesis and molecular modeling.

The acyl binding site of Rhizopus delemar prolipase and mature lipase was altered through site-directed mutagenesis to improve lipase specificity for short- or medium-chain length fatty acids. Computer-generated structural models of R. delemar lipase were used in mutant protein design and in the interpretation of the catalytic properties of the resulting recombinant enzymes. Molecular dynamics simulations of the double mutant, val209trp + phe112trp, predicted that the introduction of trp112 and trp209 in the acyl binding groove would sterically hinder the docking of fatty acids longer than butyric acid. Assayed against a mixture of triacylglycerol substrates, the val209trp + phe112trp mature lipase mutant showed an 80-fold increase in the hydrolysis of tributyrin relative to the hydrolysis of tricaprylin while no triolein hydrolysis was detected. By comparison, the val94Trp mutant, predicted to pose steric or geometric constraints for docking fatty acids longer than caprylic acid in the acyl binding groove, resulted in a modest 1.4-fold increase in tricaprylin hydrolysis relative to the hydrolysis of tributyrin. Molecular models of the double mutant phe95asp + phe214arg indicated the creation of a salt bridge between asp95 and arg214 across the distal end of the acyl binding groove. When challenged with a mixture of triacylglycerols, the phe95asp + phe214arg substitutions resulted in an enzyme with 3-fold enhanced relative activity for tricaprylin compared to triolein, suggesting that structural determinants for medium-chain length specificity may reside in the distal end of the acyl binding groove. Attempts to introduce a salt bridge within 8 A of the active site by the double mutation leu146lys + ser115asp destroyed catalytic activity entirely. Similarly, the substitution of polar Gln at the rim of the acyl binding groove for phe112 largely eliminated catalytic activity of the lipase.

Caprylates↗

Comparative conformational analysis and in vitro pharmacological evaluation of three cyclic hexapeptide NK-2 antagonists.

The conformational analysis of three cyclic hexapeptides is presented. Cyclo-(-Gln6-Trp7-Phe8-Gly9-Leu10-D-Met11-) (1) and cyclo-(-Gln6-Trp7-Phe8-Gly9-Leu10-Met11-) (2) are NK-2 antagonists in the hamster trachea assay, whereas cyclo-(-Gln6-Trp7-Phe8-(R)-Gly9-[ANC-2]Leu10-Met11+ ++-) (3), where Gly9[ANC-2]Leu10 represents (2S)-2-((3R)-3-amino-2-oxo-1-pyrrolidinyl)-4-methylpentanoyl, is inactive as agonist and antagonist in this assay. In DMSO, the NMR results cannot be interpreted as being consistent with a single conformation. However, the combined interpretation of results from NMR spectroscopy, restrained molecular dynamics simulations with application of proton-proton distance information from ROESY spectra, and pharmacological results leads to a reduced number of conformational domains for each peptide, which can be compared with each other and may be classified as responsible for their biological activity. Trying to match the conformational domains approximately with regular beta- and gamma-turns, we find a gamma n-turn at the position of the methionine occurring in all peptides. For the active peptides 1 and 2 we arrive at an inverse gamma i-turn at Phe8, and beta I'- or beta II-turns with Gly9 and Leu10 at the corner positions, these beta-turns having a similar topology with respect to the linking peptide unit. Other conformational domains common to only 1 and 2 support their classification as responsible for the biological activity.

Amino Acid Sequence↗

Design and characterization of a multisite fluorescence energy-transfer system for protein folding studies: a steady-state and time-resolved study of yeast phosphoglycerate kinase.

A multisite distance-based fluorescence resonance energy-transfer assay system was developed for the study of protein folding reactions. Single- and double-cysteine substitution mutagenesis was utilized to place sulfhydryl residues throughout the tertiary structure of the bidomain enzyme yeast phosphoglycerate kinase (PGK). These reactive cysteines were covalently modified with extrinsic donor [5-[[2-(2-iodoacetamido)ethyl]amino]-1-naphthalenesulfonic acid] and acceptor (5-iodoacetamidofluorescein) fluorescent labels. A detailed experimental strategy was followed, which revealed that, when these relatively large extrinsic fluorescent labels are covalently attached to properly selected solvent-exposed residues, they do not affect the intrinsic stability of the protein. The PGK crystal structure was combined with molecular dynamics simulations of the dyes built into the protein and time-resolved anisotropy experiments, in order to estimate a more realistic orientation factor, *, for each donor/acceptor pair. Time-resolved and steady-state fluorescence energy-transfer experiments revealed that this distance assay, spanning six different donor-acceptor distances, is linear and accurate (to within 10-20%) over the range of 30-70 A. This distance assay system for PGK allows for the measurement of long-range changes in intra- and interdomain spatial organization during protein folding reactions. The approach which we have developed can be applied to any protein system in which unique one- and two-site cysteine residues can be engineered into a protein. In the following paper [Lillo, M. P., et al. (1997) Biochemistry 36, 11273-11281], these multisite energy-transfer pairs are utilized for stopped-flow unfolding studies.

Crystallography, X-Ray↗

Insight into catalytically relevant correlated motions in human purine nucleoside phosphorylase.

The catalytic site of the homotrimeric enzyme human purine nucleoside phosphorylase enzyme (hPNP) features residue F200 and the 241-265 loop directly skirting the purine base and a residue belonging to the adjacent monomer, F159, immediately conterminous to the ribosyl moiety. Crystallographic B-factors of apo human purine nucleoside phosphorylase, and hPNP complexed with substrate/transition state (TS) analogues, show that residue E250 is the centroid of a highly mobile loop region. Furthermore, superimposition of apo hPNP and hPNP complexed with TS analogue Immucillin-H shows a tightening of the active site, caused by the ligand-dependent 241-265 loop rearrangement taking place upon substrate/inhibitor binding, suggesting a putative dynamic role of the loop in binding/catalysis. However, crystallographic structures reveal only average atomic positions, and more detailed information is needed to discern the dynamic behavior of hPNP. The Essential Dynamics (ED) method is used here to investigate the existence of correlated motions in hPNP and consequently proposes mutagenesis assays to estimate the relative importance of these motions in the phosphorolytic efficiency of the reaction catalyzed by hPNP. We compare the concerted motions obtained from multiple molecular dynamics simulations of apo and Michaelis complex of hPNP both in vacuo and in solution. The results of the principal component analysis for the apo hPNP indicate the existence of strong correlations predominantly in the vicinity of residue F159. However, for the Michaelis complex, concerted motions are seen mostly around both active site residue F200 and loop residue E250. Additionally, for a simulation depicting the relaxation of tight complexed hPNP with a TS analogue, toward its relaxed apo form (after removal of the TS analog), a combination of the apo hPNP and Michaelis complex motions is found, with prominent concerted modes centered around neighboring residues F159, F200, and E250. Finally, we probed the extent to which these concerted motions bear an intrinsic catalytic role by performing experimental site-directed mutagenesis on some residues, followed by kinetic analysis. The F159G and F200G mutants displayed a strong increase in K(M) and modest decrease in k(cat), suggesting that these concerted motions may provide dynamical roles in substrate binding and/or catalysis. However, further structural data for the hPNP mutants are needed to confirm our hypothesis.

Binding Sites↗

A two-factor authentication mechanism licenses pilins for pilus assembly in gram-positive bacteria.

Gram-positive bacteria display virulence-associated pili that facilitate adhesion and biofilm formation. These pili are covalently polymerized by class C sortase enzymes, which selectively recognize their cognate pilin substrates amid numerous cell wall sorting signal (CWSS)-bearing proteins. The molecular basis for this stringent substrate specificity has remained unclear. Here, we develop a rapid, quantitative fluorescence-activated cell sorting assay to monitor pilus assembly in Corynebacterium diphtheriae, enabling high-throughput analysis of SpaA pilin and SrtA sortase variants. Using this platform, together with molecular modeling and dynamics simulations, we show that SrtA engages nearly the entire SpaA CWSS to form a membrane-embedded complex that incorporates not only the LPXTG motif but also its connector and transmembrane helix elements. Formation of this interface displaces an inhibitory active-site lid and activates the enzyme to load the pilin substrate. Systematic CWSS swapping experiments and deep mutational scanning further support this model, demonstrating that noncognate pilins are excluded because they fail to form the required interface. Conversely, SrtA variants with an artificially unlatched lid bypass the need for this interface, indicating that membrane-driven complex formation is important for substrate licensing. Together, these findings define a "two-factor authentication" mechanism for pilus assembly in gram-positive bacteria: class C sortases first verify pilin identity by forming a membrane-embedded interface that activates the enzyme, then they recognize the LPXTG motif to initiate loading and crosslinking. This work provides a unified molecular framework for selective pilin incorporation in gram-positive bacteria and identifies potential vulnerabilities in the licensing machinery that may be exploited therapeutically.

Fimbriae, Bacterial↗

Insight into the inhibition of human choline kinase: homology modeling and molecular dynamics simulations.

A homology model of human choline kinase (CK-alpha) based on the X-ray crystallographic structure of C. elegans choline kinase (CKA-2) is presented. Molecular dynamics simulations performed on CK-alpha confirm the quality of the model, and also support the putative ATP and choline binding sites. A good correlation between the MD results and reported CKA-2 mutagenesis assays has been found for the main residues involved in catalytic activity. Preliminary docking studies performed on the CK-alpha model indicate that inhibitors can bind to the binding sites of both substrates (ATP and choline). A possible reason for inhibition of choline kinase by Ca(2+) ion is also proposed.

Amino Acid Sequence↗

Selection of peptides that bind to the HLA-A2.1 molecule by molecular modelling.

Cytotoxic T lymphocytes recognize antigenic peptides in association with major histocompatibility complex class I proteins. Although a large set of class I binding peptides has been described, it is not yet easy to search for potentially antigenic peptides without synthesis of a panel of peptides, and subsequent binding assays. In order to predict HLA-A2.1-restricted antigenic epitopes, a computer model of the HLA-A2.1 molecule was established using X-ray crystallography data. In this model nonameric peptide sequences were aligned. In a molecular dynamics (MD) simulation with two sets of peptides known to be presented by HLA-A2.1, it was important to know the anchor amino acid residue preference and the distance between the anchor residues. We show here that the peptides bound to the HLA-A2.1 model structure possess a side chain of C-terminal anchor residue oriented into the binding groove with different distances between the two anchor residues from 15 to 21A. We also synthesized a set of nonamer peptides containing amino acid sequences of Hepatitis B virus protein that were selected on the basis of previously described HLA-A2.1 specific motifs. When results obtained from the MD simulation were compared with functional binding assays using the TAP-deficient cell line T2, it was evident that the MD simulation method improves prediction of the HLA-A2.1 binding epitope sequence. These results suggest that this approach can provide a way to predict peptide epitopes and search for antigenic regions in sequences in a variety of antigens without screening a large number of synthetic peptides.

Amino Acid Sequence↗

Biochemical and molecular characterization of the novel BRAF(V599Ins) mutation detected in a classic papillary thyroid carcinoma.

Activating mutations of the BRAF gene are the most common genetic alterations in papillary thyroid carcinomas (PTCs) and the T1799A transversion, resulting in BRAFV600E, appeared virtually unique in this cancer type. Here, we report on the identification in a classic PTC of a novel BRAF mutation, namely a 1795GTT insertion, resulting in BRAFV599Ins, and describe its biochemical and molecular characterization. Kinase assays carried out on BRAFV599Ins and BRAFV600E revealed a three- to five-fold increase in the enzymatic activity of both mutants with respect to BRAFWT. Similarly, evaluation of BRAF-induced phosphorylation of MEK, MAPK and RSK revealed a significant MAPK cascade activation in cells expressing BRAFV599Ins or BRAFV600E, but not in cells expressing BRAFWT. Molecular dynamic simulations showed a destabilization of the inactive conformation of the enzyme in both BRAFV599Ins and BRAFV600E mutants, but not in BRAFWT. The analysis of the interaction energies inside the catalytic site allowed to demonstrate the presence of repulsive electrostatic forces acting on the activation loop and moving from inward to outward of the mutant enzymes. Finally, focus assays in NIH-3T3 cells confirmed a high transformation rate in the cells transfected either with BRAFV599Ins or BRAFV600E. In conclusion, this study demonstrated that BRAFV599Ins, as BRAFV600E, is a 'gain of function' mutation, characterized by a constitutive catalytic activation, which accounts for its causative role in the studied PTC.

Aged↗

Lysine 183 and glutamic acid 157 of the TSH receptor: two interacting residues with a key role in determining specificity toward TSH and human CG.

A naturally occurring mutation in the ectodomain of the TSH receptor (TSHr), K183R, has been described recently in a familial case of gestational hyperthyroidism. Hyperthyroidism was explained by the widening of the specificity of the mutant receptor toward human CG (hCG). In the present study, we attempted to understand in molecular terms the structure-phenotype relationships of this mutant in light of the available structural model of TSHr ectodomain established on the template of the atomic structure of the porcine ribonuclease inhibitor. To this aim, we studied by site-directed mutagenesis and functional assays in transfected COS cells the effects of substituting amino acids with different physicochemical properties for lysine 183. Unexpectedly, all TSHr mutants displayed widening of their specificity toward hCG. Molecular dynamics simulations suggested that the gain of function would be secondary to the release of a nearby glutamate residue (E157) from a salt bridge with K183. This hypothesis was supported by further site-directed mutagenesis experiments showing that the presence of an acidic residue in position 157, or in its vicinity, was required to observe the increase in sensitivity to hCG (an acidic residue in position 183 can partially fulfill the role of a free acidic residue in position 157 when tested on the background of a E157A mutant). Our results suggest also that additional natural mutations (especially K183M, N, or Q) in position 183 of TSHr are expected to be found in gestational hyperthyroidism.

Amino Acid Substitution↗