Conformation Of Residues Research Articles

Structural Basis of the Binding Mode of the Antineoplastic Compound Motixafortide (BL-8040) in the CXCR4 Chemokine Receptor

Modulation of the CXCL12-CXCR4 signaling axis is of the utmost importance due to its central involvement in several pathological disorders, including inflammatory diseases and cancer. Among the different currently available drugs that inhibit CXCR4 activation, motixafortide-a best-in-class antagonist of this GPCR receptor-has exhibited promising results in preclinical studies of pancreatic, breast, and lung cancers. However, detailed information on the interaction mechanism of motixafortide is still lacking. Here, we characterize the motixafortide/CXCR4 and CXCL12/CXCR4 protein complexes by using computational techniques including unbiased all-atom molecular dynamics simulations. Our microsecond-long simulations of the protein systems indicate that the agonist triggers changes associated with active-like GPCR conformations, while the antagonist favors inactive conformations of CXCR4. Detailed ligand-protein analysis indicates the importance of motixafortide's six cationic residues, all of which established charge-charge interactions with acidic CXCR4 residues. Furthermore, two synthetic bulky chemical moieties of motixafortide work in tandem to restrict the conformations of important residues associated with CXCR4 activation. Our results not only elucidate the molecular mechanism by which motixafortide interacts with the CXCR4 receptor and stabilizes its inactive states, but also provide essential information to rationally design CXCR4 inhibitors that preserve the outstanding pharmacological features of motixafortide.

Genetic encoding of 3-nitro-tyrosine reveals the impacts of 14-3-3 nitration on client binding and dephosphorylation.

14-3-3 proteins are central hub regulators of hundreds of phosphorylated "client" proteins. They are subject to over 60 post-translational modifications (PTMs), yet little is known how these PTMs alter 14-3-3 function and its ability to regulate downstream signaling pathways. An often neglected, but well-documented 14-3-3 PTM found under physiological and immune-stimulatory conditions is the conversion of tyrosine to 3-nitro-tyrosine at several Tyr sites, two of which are located at sites considered important for 14-3-3 function: Y130 (β-isoform numbering) is located in the primary phospho-client peptide-binding groove, while Y213 is found on a secondary binding site that engages with clients for full 14-3-3/client complex formation and client regulation. By genetically encoding 3-nitro-tyrosine, we sought to understand if nitration at Y130 and Y213 effectively modulated 14-3-3 structure, function, and client complexation. The 1.5Å resolution crystal structure of 14-3-3 nitrated at Y130 showed the nitro group altered the conformation of key residues in the primary binding site, while functional studies confirmed client proteins failed to bind this variant of 14-3-3. But, in contrast to other client-binding deficient variants, it did not localize to the nucleus. The 1.9Å resolution structure of 14-3-3 nitrated at Y213 revealed unusual flexibility of its C-terminal α-helix resulting in domain swapping, suggesting additional structural plasticity though its relevance is not clear as this nitrated form retained its ability to bind clients. Collectively, our data suggest that nitration of 14-3-3 will alter downstream signaling systems, and if uncontrolled could result in global dysregulation of the 14-3-3 interactome.

Self-inhibited State of Venezuelan Equine Encephalitis Virus (VEEV) nsP2 Cysteine Protease: A Crystallographic and Molecular Dynamics Analysis

The Venezuelan equine encephalitis virus (VEEV) belongs to the Togaviridae family and is pathogenic to both humans and equines. The VEEV non-structural protein 2 (nsP2) is a cysteine protease (nsP2pro) that processes the polyprotein and thus it is a drug target for inhibitor discovery. The atomic structure of the VEEV nsP2 catalytic domain was previously characterized by both X-ray crystallography and computational studies. A modified nsP2pro harboring a N475A mutation in the N terminus was observed to exhibit an unexpected conformation: the N-terminal residues bind to the active site, mimicking binding of a substrate. The large conformational change of the N terminus was assumed to be induced by the N475A mutation, as N475 has an important role in stabilization of the N terminus and the active site. This conformation was first observed in the N475A mutant, but we also found it while determining a crystal structure of the catalytically active nsP2pro containing the wild-type N475 active site residue and K741A/K767A surface entropy reduction mutations. This suggests that the N475A mutation is not a prerequisite for self-inhibition. Here, we describe a high resolution (1.46 Å) crystal structure of a truncated nsP2pro (residues 463–785, K741A/K767A) and analyze the structure further by molecular dynamics to study the active and self-inhibited conformations of nsP2pro and its N475A mutant. A comparison of the different conformations of the N-terminal residues sheds a light on the interactions that play an important role in the stabilization of the enzyme.

In silico scanning of structural and functional deleterious nsSNPs in Arabidopsis thaliana’s SOG1 protein, using molecular dynamic simulation approaches

Suppressor of gamma response 1 (SOG1) is a member of the NAC domain family transcription factors of the DNA damage response (DDR) signaling in the plant’s genome. SOG1 is directly involved in transcriptional response to DNA damage, cell cycle checkpoints and ATR or ATM-mediated activation of the DNA damage responses and repair functioning in programmed cell death and regulation of end reduplication. Different mutations in the SOG1 protein lead to severe diseases and, ultimately, cell death. Single nucleotide polymorphisms (SNPs) are an important type of genetic alteration that cause different diseases or programmed cell death. The current study applied different computational approaches to Arabidopsis thaliana L. SOG1 protein to identify the potential deleterious nsSNPs and monitor their impact on the structure, function and protein stability. Various bioinformatics tools were applied to analyze the retrieved 34 nsSNPs and interestingly extracted four deleterious nsSNPs, that is, ensvath13968004 (Q166L), tmp18998388 (P159L), ensvath01103049 (K199N) and tmp18998295 (Y190F). For example, homology modeling, conservation and conformational analysis of the mutant’s models were considered to scrutinize the deviations of these variants from the native SOG1 structure. All atoms molecular dynamic simulation confirmed the significance of these mutations on the protein stability, residual and structural conformation, compactness, surface conformation, dominant motion, Gibbs free energy distribution and dynamic effects. Similarly, protein–protein interaction revealed that SOG1 operates as a hub-linking cluster of various proteins, and any changes in the SOG1 might result in the disassociation of several signal transduction cascades. Communicated by Ramaswamy H. Sarma

An in silico comparative study of curcumin and 2-deoxyuridine nucleoside derivatives: Reveals the role of angiogenin in ER stress-induced apoptosis signaling.

Angiogenin (ANG) protein plays a crucial role in angiogenesis, neovascularization, and cancer metastasis in NSCLC (non-small cell lung cancer) via non-coding tiRNA. It protects the cell under ER (endoplasmic reticulum) stress-induced apoptosis through the translational reprogramming process. Although B82 (Curcumin derivatives) induces ER stress-induced apoptosis, its mechanism of action was not studied. Therefore, it was hypothesized that the ribonucleolytic activity of ANG may be regulated by B82, resulting in modulated ER stress signaling for apoptosis. Hence, we designed and proposed a synthesis scheme for RNA-based anti-angiogenic derivatives of 2-deoxyuridine nucleoside forming peptide bond with amino acids like serine (Ser-3) and para-hydroxy-phenyl glycine (Normtyr-1) and compared B82 with them to know the binding affinity with ANG, anti-angiogenic potential, and its probable mechanism of anti-RNase activity through MD simulation study. Therefore, using Gromos96 43a1 and 43a2 force fields, MD simulation was performed to investigate binding affinity, ligand-induced molecular surface area change, conformational change, and dynamics of catalytic site residues to predict ligand binding to ANG in this study. The obtained binding free energy (∆Gbind ) result showed the total average ∆Gbind as -113.480 ± 1.682 (Normtyr-1) > -53.038 ± 33.069 (B82) > -27.909 ± 16.438 (Ser-3)kJ/mole specify role of B82 in regulating ER stress signaling induced apoptosis through ANG ribonucleolytic activity inhibition, suitability of 43a2 force fields and methodology in ligand screening. It shows the crucial role of Leu115 and His13 residue involvement in total ∆Gbind contribution. Hence, based on the MD result, novel conformation of catalytic residues, and ∆Gbind , a promising combination candidate could be proposed for metastatic NSCLC therapy.

A conformational equilibrium in the nitrogenase MoFe protein with an α-V70I amino acid substitution illuminates the mechanism of H2 formation.

Study of α-V70I-substituted nitrogenase MoFe protein identified Fe6 of FeMo-cofactor (Fe7S9MoC-homocitrate) as a critical N2 binding/reduction site. Freeze-trapping this enzyme during Ar turnover captured the key catalytic intermediate in high occupancy, denoted E4(4H), which has accumulated 4[e-/H+] as two bridging hydrides, Fe2-H-Fe6 and Fe3-H-Fe7, and protons bound to two sulfurs. E4(4H) is poised to bind/reduce N2 as driven by mechanistically-coupled H2 reductive-elimination of the hydrides. This process must compete with ongoing hydride protonation (HP), which releases H2 as the enzyme relaxes to state E2(2H), containing 2[e-/H+] as a hydride and sulfur-bound proton; accumulation of E4(4H) in α-V70I is enhanced by HP suppression. EPR and 95Mo ENDOR spectroscopies now show that resting-state α-V70I enzyme exists in two conformational states, both in solution and as crystallized, one with wild type (WT)-like FeMo-co and one with perturbed FeMo-co. These reflect two conformations of the Ile residue, as visualized in a reanalysis of the X-ray diffraction data of α-V70I and confirmed by computations. EPR measurements show delivery of 2[e-/H+] to the E0 state of the WT MoFe protein and to both α-V70I conformations generating E2(2H) that contains the Fe3-H-Fe7 bridging hydride; accumulation of another 2[e-/H+] generates E4(4H) with Fe2-H-Fe6 as the second hydride. E4(4H) in WT enzyme and a minority α-V70I E4(4H) conformation as visualized by QM/MM computations relax to resting-state through two HP steps that reverse the formation process: HP of Fe2-H-Fe6 followed by slower HP of Fe3-H-Fe7, which leads to transient accumulation of E2(2H) containing Fe3-H-Fe7. In the dominant α-V70I E4(4H) conformation, HP of Fe2-H-Fe6 is passively suppressed by the positioning of the Ile sidechain; slow HP of Fe3-H-Fe7 occurs first and the resulting E2(2H) contains Fe2-H-Fe6. It is this HP suppression in E4(4H) that enables α-V70I MoFe to accumulate E4(4H) in high occupancy. In addition, HP suppression in α-V70I E4(4H) kinetically unmasks hydride reductive-elimination without N2-binding, a process that is precluded in WT enzyme.

Molecular Insights into the Calcium Binding in Troponin C through a Molecular Dynamics Study

Calcium-binding proteins play critical roles in various biological processes such as signal transduction, cell growth, and transcription factor regulation. Ion binding and target binding of Ca2+-binding proteins are highly related. Therefore, understanding the ion binding mechanism will benefit the relevant inhibitor design toward the Ca2+-binding proteins. The EF-hand is the typical ion binding motif in Ca2+-binding proteins. Previous studies indicate that the ion binding affinity of the EF-hand increases with the peptide length, but this mechanism has not been fully understood. Herein, using molecular dynamics simulations, thermodynamic integration calculations, and molecular mechanics Poisson-Boltzmann surface area analysis, we systematically investigated four Ca2+-binding peptides containing the EF-hand loop in site III of rabbit skeletal troponin C. These four peptides have 13, 21, 26, and 34 residues. Our simulations reproduced the observed trend that the ion binding affinity increases with the peptide length. Our results implied that the E-helix motif preceding the EF-hand loop, likely the Phe99 residue in particular, plays a significant role in this regulation. The E-helix has a significant impact on the backbone and side-chain conformations of the Asp103 residue, rigidifying important hydrogen bonds in the EF-hand and decreasing the solvent exposure of the Ca2+ ion, hence leading to more favorable Ca2+ binding in longer peptides. The present study provides molecular insights into the ion binding in the EF-hand and establishes an important step toward elucidating the responses of Ca2+-binding proteins toward the ion and target availability.

Structural basis of transcription regulation by CNC family transcription factor, Nrf2.

Several basic leucine zipper (bZIP) transcription factors have accessory motifs in their DNA-binding domains, such as the CNC motif of CNC family or the EHR motif of small Maf (sMaf) proteins. CNC family proteins heterodimerize with sMaf proteins to recognize CNC-sMaf binding DNA elements (CsMBEs) in competition with sMaf homodimers, but the functional role of the CNC motif remains elusive. In this study, we report the crystal structures of Nrf2/NFE2L2, a CNC family protein regulating anti-stress transcriptional responses, in a complex with MafG and CsMBE. The CNC motif restricts the conformations of crucial Arg residues in the basic region, which form extensive contact with the DNA backbone phosphates. Accordingly, the Nrf2-MafG heterodimer has approximately a 200-fold stronger affinity for CsMBE than canonical bZIP proteins, such as AP-1 proteins. The high DNA affinity of the CNC-sMaf heterodimer may allow it to compete with the sMaf homodimer on target genes without being perturbed by other low-affinity bZIP proteins with similar sequence specificity.

Insights from the Structure of an Active Form of Bacillus thuringiensis Cry5B.

The crystal protein Cry5B, a pore-forming protein produced by the soil bacterium Bacillus thuringiensis, has been demonstrated to have excellent anthelmintic activity. While a previous structure of the three-domain core region of Cry5B(112-698) had been reported, this structure lacked a key N-terminal extension critical to function. Here we report the structure of Cry5B(27-698) containing this N-terminal extension. This new structure adopts a distinct quaternary structure compared to the previous Cry5B(112-698) structure, and also exhibits a change in the conformation of residues 112-140 involved in linking the N-terminal extension to the three-domain core by forming a random coil and an extended α-helix. A role for the N-terminal extension is suggested based on a computational model of the tetramer with the conformation of residues 112-140 in its alternate α-helix conformation. Finally, based on the Cry5B(27-698) structure, site-directed mutagenesis studies were performed on Tyr495, which revealed that having an aromatic group or bulky group at this residue 495 is important for Cry5B toxicity.

Characterization of Deubiquitinase Catalytic State Using a Structure-Based Approach.

The activity of deubiquitinases (DUBs) is tightly regulated in eukaryotes via various mechanisms. One of the regulatory strategies is substrate-induced catalytic triad rearrangement, where ubiquitin-binding helps the DUB adopt an active conformation for catalysis. The crystal structure of the apo form of such a DUB, when not bound to ubiquitin, reveals an inactive conformation of the catalytic residues, necessitating the structure of the ubiquitin-bound form to visualize the active state of the DUB. Comparing the apo and ubiquitin-bound structures reveals conformational changes leading to catalytic activation. To capture the deubiquitinase in its ubiquitin-bound form, a series of activity-based ubiquitin probes (Ub-ABPs) harboring C-terminal electrophiles were designed to react with the catalytic nucleophile of cysteine protease DUBs. The resulting covalently linked DUB-ubiquitin complex is amendable for structural studies to probe the DUB-ubiquitin interface and the potential conformational change of the DUB. Here, we present a detailed protocol for the generation and purification of ubiquitin carboxy-terminal hydrolase L1 (UCHL1) in complex with a Ub-ABP, ubiquitin-vinyl methyl ester (UbVME), and the subsequent structural analysis to characterize the catalytic state of the DUB.

Molecular Dynamics Investigation on the Effects of Protonation and Lysyl Hydroxylation on Sulfilimine Cross-links in Collagen IV.

Collagen IV networks are an essential component of basement membranes that are important for their structural integrity and thus that of an organism's tissues. Improper functioning of these networks has been associated with several diseases. Cross-links, such as sulfilimine bonds interconnecting NC1 domains, are critical for forming and mechanically stabilizing these collagen IV networks. More specifically, the sulfilimine cross-links form between methionine (Met93) and lysine/hydroxylsine (Lys211/Hyl211) residues of NC1 domains. Therefore, the dynamic nature of the sulfilimine bond in collagen IV is crucial for network formation. To understand the dynamic nature of a neutral and protonated sulfilimine bond in collagen IV, we performed molecular dynamics (MD) simulations on four sulfilimine cross-linked systems (i.e., Met93S-NLys211, Met93S-NHLys211 +, Met93S-NHyl211, and Met93S-NHHyl211 +) of collagen IV. The MD results showed that the neutral Met93S-NLys211 system has the smallest protein backbone and showed the cross-linked residues' RMSD value. The conformational change analyses showed that the conformations of the sulfilimine cross-linked residues take on a U-shape for the Met93S-NHyl211 and Met93S-HNHyl211 + systems, whereas the conformations of the sulfilimine cross-linked residues are more open for the Met93S-NLys211, and Met93S-NHLys211 + systems. Protonation is a crucial biochemical process to stabilize the protein structure or the biological cross-links. Furthermore, the protonation of the sulfilimine bond could potentially influence hydrogen bond interaction with near amino acid residues, and according to water distribution analyses, the sulfilimine bond can potentially exist in one or more protonation states.

Comparison of allosteric signaling in DnaK and BiP using mutual information between simulated residue conformations.

The heat shock protein 70 kDa (Hsp70) chaperone system serves as a critical component of protein quality control across a wide range of prokaryotic and eukaryotic organisms. Divergent evolution and specialization to particular organelles have produced numerous Hsp70 variants which share similarities in structure and general function, but differ substantially in regulatory aspects, including conformational dynamics and activity modulation by cochaperones. The human Hsp70 variant BiP (also known as GRP78 or HSPA5) is of therapeutic interest in the context of cancer, neurodegenerative diseases, and viral infection, including for treatment of the pandemic virus SARS-CoV-2. Due to the complex conformational rearrangements and high sequential variance within the Hsp70 protein family, it is in many cases poorly understood which amino acid mutations are responsible for biochemical differences between protein variants. In this study, we predicted residues associated with conformational regulation of human BiP and Escherichia coli DnaK. Based on protein structure networks obtained from molecular dynamics simulations, we analyzed the shared information between interaction timelines to highlight residue positions with strong conformational coupling to their environment. Our predictions, which focus on the binding processes of the chaperone's substrate and cochaperones, indicate residues filling potential signaling roles specific to either DnaK or BiP. By combining predictions of individual residues into conformationally coupled chains connecting ligand binding sites, we predict a BiP specific secondary signaling pathway associated with substrate binding. Our study sheds light on mechanistic differences in signaling and regulation between Hsp70 variants, which provide insights relevant to therapeutic applications of these proteins.

Structural perturbations of substrate binding and oxidation state changes in a lytic polysaccharide monooxygenase

LPMOs are enzymes which catalyse the oxidation of a C-H bond within polysaccharides, leading to their oxidative cleavage. To achieve this, LPMOs employ highly reactive oxidising intermediates, the generation of which is likely coupled to substrate binding to the enzyme. The nature of this coupling is unknown. Here we report a statistical comparison for four three-dimensional structures of an AA9 LPMO crystallised in the same space group but in different oxidation and substrate-binding states, to determine which significant structural perturbations occur at the enzyme upon either oxidation state change or the binding of substrate. In a novel step, we determine the global random error associated with the positional coordinates of atoms using the method of moments to ascertain the statistical estimators of Gaussian distributions of pairwise RMS differences between individual atoms in different structures. The results show that a change in the oxidation state of the copper leads to no significant structural changes, and that binding of the substrate leads to a single change in the conformation of a tryptophan residue. This tryptophan has previously been identified as part of a charge transfer pathway between the active site and the external surface of the protein, and the structural change identified herein may be part of the substrate-enzyme coupling mechanism.Graphical abstract

Simultaneously improving the specific activity and thermostability of α-amylase BLA by rational design.

Higher activity and alkaline α-amylases are desired for textile desizing and detergent additive. Here, rational design was used to improve the specific activity and thermostability of the α-amylase BLA from Bacillus licheniformis. Seventeen mutants of BLA were designed based on sequence consensus analysis and folding free energy calculation, and characterized by measuring their respective activity and thermostability at pH 8.5. Among them, mutant Q360C exhibited nearly threefold improved activity than that of wild-type and retained a higher residual activity (75% vs 59% for wild-type) after preincubation at 70℃ for 30min. The modeled structures and molecular dynamics simulations analysis demonstrated that the enhanced hydrophobic interaction near residue 360 and reduced disturbance to the conformation of catalytic residues are the possible reasons for the improved thermostability and activity of Q360C. The results suggest that 360th of BLA may act as a hotspot for engineering other enzymes in the GH13 superfamily.

Interdisciplinary in silico studies to understand in-depth molecular level mechanism of drug resistance involving NS3-4A protease of HCV

Hepatitis C virus (HCV) causes hepatitis, a life-threatening disease responsible for liver cirrhosis. Urgent measures have been taken to develop therapeutics against this deadly pathogen. NS3/4A protease is an extremely important target. A series of inhibitors have been developed against this viral protease including Faldaprevir. Unfortunately, the error-prone viral RNA polymerase causes the emergence of resistance, thereby causing reduced effectiveness of those peptidomimetic inhibitors. Among the drug resistant variants, three single amino acid residues (R155, A156 and D168) are notable for their presence in clinical isolates and also their effectivity against most of the known inhibitors in clinical development. Therefore, it is crucial to understand the mechanistic role of those drug resistant variants while designing potent novel inhibitors. In this communication, we have deeply analyzed through using in silico studies to understand the molecular mechanism of alteration of inhibitor binding between wild type and its R155K, A156V and D168V variants. Principal component analysis was carried to identify the backbone fluctuations of important residues in HCV NS3/4A responsible for the inhibitor binding and maintaining drug resistance. Free energy landscape as a function of the principal components has been used to identify the stability and conformation of the key residues that regulate inhibitor binding and their impact in developing drug resistance. Our findings are consistent with the trend of experimental results. The observations are also true in case of other Faldaprevir-like peptidomimetic inhibitors. Understanding this binding mechanism would be significant for the development of novel inhibitors with less susceptibility towards drug resistance. Communicated by Ramaswamy H. Sarma

Chemical Conformation of the Essential Glutamate Site of the c-Ring within Thermophilic Bacillus FoF1-ATP Synthase Determined by Solid-State NMR Based on its Isolated c-Ring Structure.

Proton translocation through the membrane-embedded Fo component of F-type ATP synthase (FoF1) is facilitated by the rotation of the Fo c-subunit ring (c-ring), carrying protons at essential acidic amino acid residues. Cryo-electron microscopy (Cryo-EM) structures of FoF1 suggest a unique proton translocation mechanism. To elucidate it based on the chemical conformation of the essential acidic residues of the c-ring in FoF1, we determined the structure of the isolated thermophilic Bacillus Fo (tFo) c-ring, consisting of 10 subunits, in membranes by solid-state NMR. This structure contains a distinct proton-locking conformation, wherein Asn23 (cN23) CγO and Glu56 (cE56) CδOH form a hydrogen bond in a closed form. We introduced stereo-array-isotope-labeled (SAIL) Glu and Asn into the tFoc-ring to clarify the chemical conformation of these residues in tFoF1-ATP synthase (tFoF1). Two well-separated 13C signals could be detected for cN23 and cE56 in a 505 kDa membrane protein complex, respectively, thereby suggesting the presence of two distinct chemical conformations. Based on the signal intensity and structure of the tFoc-ring and tFoF1, six pairs of cN23 and cE56 surrounded by membrane lipids take the closed form, whereas the other four in the a-c interface employ the deprotonated open form at a proportion of 87%. This indicates that the a-c interface is highly hydrophilic. The pKa values of the four cE56 residues in the a-c interface were estimated from the cN23 signal intensity in the open and closed forms and distribution of polar residues around each cE56. The results favor a rotation of the c-ring for ATP synthesis.

Molecular mimicry between tumor associated antigens and microbiota-derived epitopes

BackgroundThe gut microbiota profile is unique for each individual and are composed by different bacteria species according to individual birth-to-infant transitions. In the last years, the local and systemic effects of microbiota on cancer onset, progression and response to treatments, such as immunotherapies, has been extensively described. Here we offer a new perspective, proposing a role for the microbiota based on the molecular mimicry of tumor associated antigens by microbiome-associated antigens.MethodsIn the present study we looked for homology between published TAAs and non-self microbiota-derived epitopes. Blast search for sequence homology was combined with extensive bioinformatics analyses.ResultsSeveral evidences for homology between TAAs and microbiota-derived antigens have been found. Strikingly, three cases of 100% homology between the paired sequences has been identified. The predicted average affinity to HLA molecules of microbiota-derived antigens is very high (< 100 nM). The structural conformation of the microbiota-derived epitopes is, in general, highly similar to the corresponding TAA. In some cases, it is identical and contact areas with both HLA and TCR chains are indistinguishable. Moreover, the spatial conformation of TCR-facing residues can be identical in paired TAA and microbiota-derived epitopes, with exactly the same values of planar as well as dihedral angles.ConclusionsThe data reported in the present study show for the first time the high homology in the linear sequence as well as in structure and conformation between TAAs and peptides derived from microbiota species of the Firmicutes and the Bacteroidetes phyla, which together account for 90% of gut microbiota. Cross-reacting CD8+ T cell responses are very likely induced. Therefore, the anti-microbiota T cell memory may turn out to be an anti-cancer T cell memory, able to control the growth of a cancer developed during the lifetime if the expressed TAA is similar to the microbiota epitope. This may ultimately represent a relevant selective advantage for cancer patients and may lead to a novel preventive anti-cancer vaccine strategy.

Folding molecular dynamics simulation of T‐peptide, a HIV viral entry inhibitor: Structure, dynamics, and comparison with the experimental data

Peptide T is a synthetic octapeptide fragment, which corresponds to the region 185-192 of the gp120 HIV coat protein and functions as a viral entry inhibitor. In this work, a folding molecular dynamics simulation of peptide T in a membrane-mimicking (DMSO) solution was performed with the aim of characterizing the peptide's structural and dynamical properties. We show that peptide T is highly flexible and dynamic. The main structural characteristics observed were rapidly interconverting short helical stretches and turns, with a notable preference for the formation of β-turns. The simulation also indicated that the C-terminal part appears to be more stable than the rest of the peptide, with the most preferred conformation for residues 5-8 being a β-turn. In order to validate the accuracy of the simulations, we compared our results with the experimental NMR data obtained for the T-peptide in the same solvent. In agreement with the simulation, the NMR data indicated the presence of a preferred structure in solution that was consistent with a β-turn comprising the four C-terminal residues. An additional comparison between the experimental and simulation-derived chemical shifts also showed a reasonable agreement between experiment and simulation, further validating the simulation-derived structural characterization of the T-peptide. We conclude that peptide folding simulations produce physically relevant results even when performed with organic solvents that were not part of the force field parameterization procedure.

Structural and Biophysical Insights into SPINK1 Bound to Human Cationic Trypsin.

(1) The serine protease inhibitor Kazal type 1 (SPINK1) inhibits trypsin activity in zymogen granules of pancreatic acinar cells. Several mutations in the SPINK1 gene are associated with acute recurrent pancreatitis (ARP) and chronic pancreatitis (CP). The most common variant is SPINK1 p.N34S. Although this mutation was identified two decades ago, the mechanism of action has remained elusive. (2) SPINK1 and human cationic trypsin (TRY1) were expressed in E. coli, and inhibitory activities were determined. Crystals of SPINK1–TRY1 complexes were grown by using the hanging-drop method, and phases were solved by molecular replacement. (3) Both SPINK1 variants show similar inhibitory behavior toward TRY1. The crystal structures are almost identical, with minor differences in the mutated loop. Both complexes show an unexpected rotamer conformation of the His63 residue in TRY1, which is a member of the catalytic triad. (4) The SPINK1 p.N34S mutation does not affect the inhibitory behavior or the overall structure of the protein. Therefore, the pathophysiological mechanism of action of the p.N34S variant cannot be explained mechanistically or structurally at the protein level. The observed histidine conformation is part of a mechanism for SPINK1 that can explain the exceptional proteolytic stability of this inhibitor.

Molecular profiling of ginsenoside metabolites to identify estrogen receptor alpha activity.

20(S)-Protopanaxadiol (PPD) and 20(S)-Protopanaxatriol (PPT) are major metabolites of ginseng in humans and are considered to have estrogenic activity in cellular bioassays. In this study, we conducted in silico analyses to determine whether PPD and PPT interact with estrogen receptor alpha (ERα) and compared them with ERα agonists, partial agonists, and antagonists to identify their ERα activity. The transcriptome profile of 17β-estradiol (E2), PPD, and PPT in MCF-7 cells expressing ERα was further compared to understand the ERα activity of ginsenoside metabolites. The results showed that PPD and PPT interacted with the 1ERE, 1GWR, and 3UUD ERα proteins in the E2 interaction model, the 3ERD protein in the diethylstilbestrol (DES) interaction model, and the 1X7R protein in the genistein (GEN) interaction model. Conversely, neither the 4PP6 protein of the interaction model with the antagonist resveratrol (RES) nor the 1ERR protein of the interaction model with the antagonist raloxifene (RAL) showed the conformation of amino acid residues. When E2, PPD, and PPT were exposed to MCF-7 cells, cell proliferation and gene expression were observed. The transcriptomic profiles of E2, PPD, and PPT were compared using a knowledge-based pathway. PPD-induced transcription profiling was similar to that of E2, and the neural transmission pathway was detected in both compounds. In contrast, PPT-induced transcription profiling displayed characteristics of gene expression associated with systemic lupus erythematosus. These results suggest that ginsenoside metabolites have ERα agonist activity and exhibit neuroprotective effects and anti-inflammatory actions. However, a meta-analysis using public microarray data showed that the mother compounds GRb1 and GRg1 of PPD and PPT showed metabolic functions in insulin signaling pathways, condensed DNA repair and cell cycle pathways, and immune response and synaptogenesis. These results suggest that the ginsenoside metabolites have potent ERα agonist activity; however, their gene expression profiles may differ from those of E2.