Search

Participating in the Janus kinase-signal transducer and activator of transcription (JAK-STAT) pathway, TYK2 emerges as a promising therapy target in controlling various autoimmune diseases, including psoriasis and multiple sclerosis. Deucravacitinib (DEU) is a novel oral TYK2-specific inhibitor approved in 2022 that is clinically effective in moderate to severe psoriasis trials. Upon the AlphaFold2 predicted TYK2 pseudokinase domain (JH2) and kinase domain (JH1), we explored the details of the underlined allosteric inhibition mechanism on TYK2 JH2-JH1 with the aid of molecular dynamics simulation. Our results suggest that the allosteric inhibition of DEU on TYK2 is accomplished by affecting the JH2-JH1 interface and hampering the state transition and ATP binding in JH1. Particularly, DEU binding stabilized the autoinhibitory interface between JH2 and JH1 while disrupting the formation of the activation interface. As a result, the negative regulation of JH2 on JH1 was greatly enhanced. These findings offer additional details on the pseudokinase-dependent autoinhibition of the JAK kinase domain and provide theoretical support for the JH2-targeted drug discovery in JAK members., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Urinary tract infections (UTIs), primarily caused by uropathogenic Escherichia coli (UPEC), have high morbidity and recurrence rates. Resistance to levofloxacin hydrochloride (LEV), a commonly used treatment for UTIs, is increasingly problematic, exacerbated by biofilm formation mediated by interactions between cyclic di-GMP (c-di-GMP or CDG) and YcgR. In this study, we identified three caffeoylquinic acid compounds from Pyrrosia lingua -chlorogenic acid (CGA), sibiricose A5 (Si-A5), and 3- O -caffeoylquinic acid methyl ester (CAM)-that target YcgR through molecular docking. Biological assays revealed that combining these compounds with levofloxacin hydrochloride significantly enhanced antibacterial activity against standard UPEC strains in a concentration-dependent manner and clinically isolated UPEC strains. Notably, chlorogenic acid and sibiricose A5, when used with levofloxacin hydrochloride, enhanced intracellular c-di-GMP levels and swimming motility, significantly reduced YcgR gene expression, and effectively inhibited biofilm formation of UPEC at multiple time points. Additionally, molecular dynamics simulations elucidated the strong binding of these compounds to YcgR, underscoring the critical roles of residues, such as Arg118 and Asp145. This research serves as a foundation for tackling antibiotic resistance and developing innovative therapeutics for UTIs.

Parathyroid hormone (PTH) type 1 receptor (PTH1R), as a typical class B1 G protein-coupled receptor (GPCR), is responsible for regulating bone turnover and maintaining calcium homeostasis, and its dysregulation has been implicated in the development of several diseases. The extracellular domain (ECD) of PTH1R is crucial for the recognition and binding of ligands, and the receptor may exhibit an autoinhibited state with the closure of the ECD in the absence of ligands. However, the correlation between ECD conformations and PTH1R activation remains unclear. Thus, this study combines enhanced sampling molecular dynamics (MD) simulations and Markov state models (MSMs) to reveal the possible relevance between the ECD conformations and the activation of PTH1R. First, 22 intermediate structures are generated from the autoinhibited state to the active state and conducted for 10 independent 200 ns simulations each. Then, the MSM is constructed based on the cumulative 44 μs simulations with six identified microstates. Finally, the potential interplay between ECD conformational changes and PTH1R activation as well as cryptic allosteric pockets in the intermediate states during receptor activation is revealed. Overall, our findings reveal that the activation of PTH1R has a specific correlation with ECD conformational changes and provide essential insights for GPCR biology and developing novel allosteric modulators targeting cryptic sites.

The microtubule-associated protein tau (MAPT) has a critical role in the development and preservation of the nervous system. However, tau's dysfunction and accumulation in the human brain can lead to several neurodegenerative diseases, such as Alzheimer's disease, Down's syndrome, and frontotemporal dementia. The microtubule binding (MTB) domain plays a significant, important role in determining the tau's pathophysiology, as the core of paired helical filaments PHF6* ( 275 VQIINK 280 ) and PHF6 ( 306 VQIVYK 311 ) of R2 and R3 repeat units, respectively, are formed in this region, which promotes tau aggregation. Post-translational modifications, and in particular lysine acetylation at K280 of PHF6* and K311 of PHF6, have been previously established to promote tau misfolding and aggregation. However, the exact aggregation mechanism is not known. In this study, we established an atomic-level nucleation-extension mechanism of the separated aggregation of acetylated PHF6* and PHF6 hexapeptides, respectively, of tau. We show that the acetylation of the lysine residues promotes the formation of β-sheet enriched high-ordered oligomers. The Markov state model analysis of ac-PHF6* and ac-PHF6 aggregation revealed the formation of an antiparallel dimer nucleus which could be extended from both sides in a parallel manner to form mixed-oriented and high-ordered oligomers. Our study describes the detailed mechanism for acetylation-driven tau aggregation, which provides valuable insights into the effect of post-translation modification in altering the pathophysiology of tau hexapeptides.

The interactions of complement receptor 2 (CR2) and the degradation fragment C3d of complement component C3 mediate the innate and adaptive immune systems. Due to the importance of C3d-CR2 interaction in the design of vaccines, many studies have indicated the interactions are pH-dependent. Moreover, C3d-CR2 interactions at pH 5.0 are unknown. To investigate the molecular mechanism of pH-regulating C3d-CR2 interaction, molecular dynamics simulations for C3d-CR2 complex in different pH are performed. Our results revealed that the protonation of His9 in C3d at pH 6.0 slightly weakens C3d-CR2 association as reducing pH from 7.4 to 6.0, initiated from a key hydrogen bond formed between Gly270 and His9 in C3d at pH 6.0. When reducing pH from 6.0 to 5.0, the protonation of His33 in C3d weakens C3d-SCR1 association by changing the hydrogen-bond network of Asp36, Glu37, and Glu39 in C3d with Arg13 in CR2. In addition, the protonation of His90 significantly enhances C3d-SCR2 association. This is because the enhanced hydrogen-bond interactions of His90 with Glu63 and Ser69 of the linker change the conformations of the linker, Cys112-Asn116 and Pro87-Gly91 regions. This study uncovers the molecular mechanism of the mediation of pH on C3d-CR2 interaction, which is valuable for vaccine design., (© 2018 John Wiley & Sons A/S.)

EGCG possesses the ability of disaggregating the existing amyloid fibrils which were associated with many age-related degenerative diseases. However, the molecular mechanism of EGCG to disaggregate these fibrils is poorly known. In this work, to study the influence of EGCG on the full-length human islet amyloid polypeptide 1-37 (hIAPP1-37) oligomers, molecular dynamics simulations of hIAPP1-37 pentamer and decamer with EGCG were performed, respectively. The obtained results indicate that EGCG indeed destabilized the hIAPP1-37 oligomers. The nematic order parameter and secondary structure calculations coupled with the free-energy landscape indicate that EGCG broke the initial ordered pattern of two polymers, greatly reduced their β-sheet content and enlarged their conformational space. On this basis, three possible target sites were identified with the binding capacity order of S1>S2>S3. After a deeper analysis of each site, we found that S1 was the most possible site on which residues B-Ile26/Ala25, A-Phe23, B/C-Leu27 and E-Tyr37 played an important role for their binding. The proposal of this molecular mechanism can not only provide a prospective interaction figure between EGCG and β-sheet-rich fibrils of hIAPP1-37, but also is useful for further discovering other potential inhibitors.

Prion diseases are marked by cerebral accumulation of the abnormal isoform of the prion protein. A fragment of prion protein composed of residues 106-126 (PrP106-126) exhibits similar properties to full length prion and plays a key role in the conformational conversion from cellular prion to its pathogenic pattern. Soluble oligomers of PrP106-126 have been proposed to be responsible for neurotoxicity. However, the monomeric conformational space and initial oligomerization of PrP106-126 are still obscure, which are very important for understanding the conformational conversion of PrP106-126. In this study, replica exchange molecular dynamics simulations were performed to investigate monomeric and dimeric states of PrP106-126 in implicit solvent. The structural diversity of PrP106-126 was observed and this peptide did not acquire stable structure. The dimeric PrP106-126 also displayed structural diversity and hydrophobic interaction drove the dimerization. To further study initial oligomerization of PrP106-126, 1 µs conventional molecular dynamics simulations of trimer and tetramer formation were carried out in implicit solvent. We have observed the spontaneous formation of several basic oligomers and stable oligomers with high β-sheet contents were sampled in the simulations of trimer and tetramer formation. The β-hairpin formed in hydrophobic tail of PrP106-126 with residues 118-120 in turn may stabilize these oligomers and seed the formation oligomers. This study can provide insight into the detailed information about the structure of PrP106-126 and the dynamics of aggregation of monomeric PrP106-126 into oligomers in atomic level.

Prion diseases are associated with misfolding and aggregation of prion protein (PrP). Cellular prion protein contains a disulfide bond linking Cys residues at positions 179 and 214. It has been proposed that this disulfide bond plays an important role in the conversion between cellular (PrP(C)) and the scrapie form of prion protein (PrP(Sc)). To probe the role of this disulfide bond in the stability and folding of prion protein, we employed molecular dynamics simulations to study the reduced prion protein and a variant of PrP in which the two cysteines were replaced by alanines residues. The simulations highlighted the changes that occurred upon breakage of the disulfide bond. Breakage of the disulfide bond resulted in a shift of H1, elongation of the native β-sheet and perturbation of the hydrophobic core of huPrP. The changes are similar to the conformational transitions of prion protein in low pH, in denaturing conditions or with pathogenic mutations, which indicate that rupture of the disulfide bond may lead to the misfolding of prion protein.

Background: The formation of amyloid fibrils is associated with many age-related degenerative diseases. Nevertheless, the molecular mechanism that directs the nucleation of these fibrils is not fully understood., Methods: Here, we performed MD simulations for the NFGAILS motif of hIAPP associated with the type II diabetes to estimate the stabilities of hIAPP22-28 protofibrils with different sizes: from 2 to 16 chains. In addition, to study the initial self-assembly stage, 4 and 8 IAPP22-28 chains in explicit solvent were also simulated., Results: Our results indicate that the ordered protofibrils with no more than 16 hIAPP22-28 chains will be structurally stable in two layers, while one-layer or three-layer models are not stable as expected. Furthermore, the oligomerization simulations show that the initial coil structures of peptides can quickly aggregate and convert to partially ordered β-sheet-rich oligomers., Conclusions: Based on the obtained results, we found that the stability of an IAPP22-28 oligomer was not only related with its size but also with its morphology. The driving forces to form and stabilize an oligomer are the hydrophobic effects and backbone H-bond interaction. Our simulations also indicate that IAPP22-28 peptides tend to form an antiparallel strand orientation within the sheet., General Significance: Our finding can not only enhance the understanding about potential mechanisms of hIAPP nuclei formation and the extensive structural polymorphisms of oligomers, but also provide valuable information to develop potential β-sheet formation inhibitors against type II diabetes., (© 2013.)

Recent advances in nanotechnologies have led to wide use of nanomaterials in biomedical field. However, nanoparticles are found to interfere with protein misfolding and aggregation associated with many human diseases. It is still a controversial issue whether nanoparticles inhibit or promote protein aggregation. In this study, we used molecular dynamics simulations to explore the effects of three kinds of carbon nanomaterials including graphene, carbon nanotube and C₆₀ on the aggregation behavior of islet amyloid polypeptide fragment 22-28 (IAPP₂₂₋₂₈). The diverse behaviors of IAPP₂₂₋₂₈ peptides on the surfaces of carbon nanomaterials were studied. The results suggest these nanomaterials can prevent β-sheet formation in differing degrees and further affect the aggregation of IAPP₂₂₋₂₈. The π-π stacking and hydrophobic interactions are different in the interactions between peptides and different nanoparticles. The subtle differences in the interaction are due to the difference in surface curvature and area. The results demonstrate the adsorption interaction has competitive advantages over the interactions between peptides. Therefore, the fibrillation of IAPP₂₂₋₂₈ may be inhibited at its early stage by graphene or SWCNT. Our study can not only enhance the understanding about potential effects of nanomaterials to amyloid formation, but also provide valuable information to develop potential β-sheet formation inhibitors against type II diabetes.

The central event in the pathogenesis of prion protein (PrP) is a profound conformational change from its α-helical (PrP(C)) to its β-sheet-rich isoform (PrP(Sc)). Many single amino acid mutations of PrP are associated with familial prion diseases, such as D202N, E211Q, and Q217R mutations located at the third native α-helix of human PrP. In order to explore the underlying structural and dynamic effects of these mutations, we performed all-atom molecular dynamics (MD) simulations for the wild-type (WT) PrP and its mutants. The obtained results indicate that these amino acid substitutions have subtle effects on the protein structures, but show large changes of the overall electrostatic potential distributions. We can infer that the changes of PrP electrostatic surface due to the studied mutations may influence the intermolecular interactions during the aggregation process. In addition, the mutations also affect the thermodynamic stabilities of PrP., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Background: Prion diseases are associated with a conformational switch for PrP from PrP(C) to PrP(Sc). Many genetic mutations are linked with prion diseases, such as mutations T188K/R/A with fCJD., Scope of Review: MD simulations for the WT PrP and its mutants were performed to explore the underlying dynamic effects of T188 mutations on human PrP. Although the globular domains are fairly conserved, the three mutations have diverse effects on the dynamics properties of PrP, including the shift of H1, the elongation of native β-sheet and the conversion of S2-H2 loop to a 3(10) helix., Major Conclusions: Our present study indicates that the three mutants for PrP may undergo different pathogenic mechanisms and the realistic atomistic simulations can provide insights into the effects of disease-associated mutations on PrP dynamics and stability, which can enhance our understanding of how mutations induce the conversion from PrP(C) to PrP(Sc). General significance Our present study helps to understand the effects of T188K/R/A mutations on human PrP: despite the three pathogenic mutations almost do not alter the native structure of PrP, but perturb its stability. This instability may further modulate the oligomerization pathways and determine the features of the PrP(Sc) assemblies., (Copyright © 2011 Elsevier B.V. All rights reserved.)

Biological membranes (biomembranes) are one of the most complicated structures that allow life to exist. Investigating their structure, dynamics, and function is crucial for advancing our knowledge of cellular mechanisms and developing novel therapeutic strategies. However, experimental investigation of many biomembrane phenomena is challenging due to their compositional and structural complexity, as well as the inherently multi‐scalar features. Computational approaches, particularly molecular dynamics (MD) simulations, have emerged as powerful tools for addressing the atomic details of biomembrane systems, driving breakthroughs in our understanding of biomembranes and their roles in cellular function. This review presents an overview of the latest advancements in related computational approaches, from force fields and model construction to MD simulations and trajectory analysis. We also discussed current hot research topics and challenges. Finally, we outline future directions, emphasizing the integration of force field development, enhanced sampling techniques, and data‐driven approaches to accelerate the growth of this field in the years to come. We aim to equip readers with an understanding of the promise and limitations of emerging computational technologies in biomembrane systems and offer valuable recommendations for future research endeavors. This article is categorized under:Structure and Mechanism > Computational Biochemistry and BiophysicsMolecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods [ABSTRACT FROM AUTHOR]

Phthalate esters (PAEs) are a ubiquitous kind of environmental endocrine that disrupt chemicals, causing environmental and health issues. EstJ6 is an effective phthalate-degrading hydrolase, and its mutant with a combination of three non-conservative distal mutations has an improved activity against PAEs with unknown molecular mechanisms. Herein, we attempt to fill the significant gap between distal mutations and the activity of this enzyme using computational approaches. We found that mutations resulted in a redistribution of the enzyme's preexisting conformational states and dynamic changes of key functional regions, especially the lid over the active site. The outward motion of the lid upon the mutations made it easier for substrates or products to enter or exit. Additionally, a stronger substrate binding affinity and conformational rearrangements of catalytic reaction-associated residues in the mutant, accompanied by the strengthened communication within the protein, could synergistically contribute to the elevated catalytic efficiency. Finally, an attempt was made to improve the thermostability of EstJ6 upon introducing a distal disulfide bond between residues A23 and A29, and the simulation results were as expected. Together, our work explored the allosteric effects caused by distal mutations, which could provide insights into the rational design of esterases for industrial applications in the future. [ABSTRACT FROM AUTHOR]

The parathyroid hormone type 1 receptor (PTH1R) acts as a canonical class B G protein-coupled receptor, regulating crucial functions including calcium homeostasis and bone formation. The identification and development of PTH1R non-peptide allosteric modulators have obtained widespread attention. It has been found that a negative allosteric modulator (NAM) could inhibit the activation of PTH1R, but the implied mechanism remains unclear. Herein, extensive molecular dynamics simulations together with multiple analytical approaches are utilized to unravel the mechanism of PTH1R allosteric inhibition. The results suggest that the binding of NAM destabilizes the structure of the PTH1R–PTH–spep/qpep (the C terminus of Gs/Gq proteins) complexes. Moreover, the presence of NAM weakens the binding of PTH/peps (spep and qpep) and PTH1R. The intra- and inter-molecular couplings are also weakened in PTH1R upon NAM binding. Interestingly, compared with our previous study of the positive allosteric effects induced by extracellular Ca2+, the enhanced correlation between the PTH and G-protein binding sites is significantly reduced by the replacement of this negative allosteric regulator. Our findings might contribute to the development of new therapeutic agents for diseases caused by the abnormal activation of PTH1R. [ABSTRACT FROM AUTHOR]

BACKGROUND: Amitraz is one of the critical acaricides/insecticides for effective control of pest infestation of Varroa destructor mite, a devastating parasite of Apis mellifera, because of its low toxicity to honeybees. Previous assays verified that a typical G protein‐coupled receptor, β‐adrenergic‐like octopamine receptor (Octβ2R), is the unique target of amitraz, but the honeybee Octβ2R resists to amitraz. However, the underlying molecular mechanism of the enhanced sensitivity or toxicity of amitraz to mutated honeybee Octβ2RE208V/I335T/I350V is not fully understood. Here, molecular dynamics simulations are employed to explore the implied mechanism of the enhanced sensitivity to amitraz in mutant honeybee Octβ2R. RESULTS: We found that amitraz binding stabilized the structure of Octβ2R, particularly the intracellular loop 3 associated with the Octβ2R signaling. Then, it was further demonstrated that both mutations and ligand binding resulted in a more rigid and compact amitraz binding site, as well as the outward movement of the transmembrane helix 6, which was a prerequisite for G protein coupling and activation. Moreover, mutations were found to promote the binding between Octβ2R and amitraz. Finally, community analysis illuminated that mutations and amitraz strengthened the residue–residue communication within the transmembrane domain, which might facilitate the allosteric signal propagation and activation of Octβ2R. CONCLUSION: Our results unveiled structural determinants of improved sensitivity in the Octβ2R‐amitraz complex and may contribute to further structure‐based drug design for safer and less toxic selective insecticides. © 2022 Society of Chemical Industry. [ABSTRACT FROM AUTHOR]

Epidermal growth factor receptor (EGFR) is an intensively focused target for anti-tumor compounds used in non-small cell lung cancer (NSCLC) therapy. Compared to the classical activating mutations, there are still many uncommon EGFR mutations associated with poorer responses to EGFR inhibitors. A detailed understanding of the molecular basis for multiple EGFR mutants exhibiting diverse responses to inhibitors is of critical importance for related drug development. Herein, we explored the molecular determinants contributing to the distinct responses of EGFR with a single rare mutation (G719S) or combined mutations (G719S/L858R and G719S/l861Q) to Gefitinib (IRE). Our results indicated that interactions, formed within the tetrad of residues S768 (in the αC-helix), D770 (in the αC-β4 loop), Y827 (in the αE-helix), and R831 (in the catalytic loop), play an important role in the stability of αC-helix and the maintenance of K745–E762 salt bridge in the absence of IRE, which are weakened in the EGFRG719S system and enhanced in the EGFRG719S/L858R system upon IRE binding. Besides, the introduced hydrogen bonds by the co-occurring mutation partner also contribute to the stability of αC-helix. The work done for inhibitor dissociation suggests that IRE exhibits a stronger binding affinity to EGFRG719S/L858R mutant. Together, these findings provide a deeper understanding of minor mutations, which is essential for drug development targeting EGFR with less common mutations. [ABSTRACT FROM AUTHOR]

Abstract As the broad application of graphene in the biomedical field, it is urgent and important to evaluate how the graphene affects the structure and function of the proteins in our body, especially the amyloid-related proteins. Prion protein, as a typical amyloid protein, it misfolding and aggregation will lead to serious prion diseases. To explore if graphene promotes or inhibits the formation of amyloid, here, we combined the experimental and molecular dynamics (MD) simulation methods to study the influence of graphene on the globular domain of prion protein (PrP 117–231). The results from fluorescence quenching and circular dichroism spectrum showed that the addition of graphene changed the secondary structure of prion protein largely, mainly reflecting in the reduced α-helix structure and the increased coil structure, indicating graphene may strengthen the misfolding inclination of prion. To further uncover the mechanism of conformational change of prion under the induction of graphene, the all-atoms MD simulations in explicit solvent were performed. Our simulations suggest that prion protein can be quickly and tightly adsorbed onto graphene together with the weak conformational rearrangement and may reorient when approaching the surface. The Van der Waals' force drive the adsorption process. In the induction of graphene, H1 and S2-H2 loop regions of prion become unstable and prion begins to misfold partially. Our work shows that graphene can induce the misfolding of prion protein and may cause the potential risk to biosystems. Highlights • To evaluate the effect of graphene on the amyloid-related proteins is important. • The influence of graphene on prion was explored by experiments and MD simulation. • The addition of graphene changed the secondary structure of prion protein. • In the induction of graphene, H1 and S2-H2 loop regions of prion become unstable. • Graphene can induce the misfolding of prion and may cause the potential risk. [ABSTRACT FROM AUTHOR]

The coreceptor binding site of gp120 plays an important role in HIV entry into host cell. X5 and 17b are typical coreceptor binding site antibodies with the ability to broadly neutralize HIV. Thus, here, to study the neutralizing mechanism of two antibodies and identify the source of two antibodies with different neutralizing ability, we performed molecular dynamics simulations for the complexes of X5 and 17b with gp120 and CD4. The simulation results indicate X5 and 17b mainly affects CD4 and coreceptor binding sites. Specifically, for CD4 binding site (CD4bs), the binding of antibodies has different effects on CD4bs with and without CD4. However, for coreceptor binding sites, the binding of the antibodies has consistent influence on the region adjacent to loop V3 despite of the simulated systems with or without CD4. The binding of the antibodies enhances the interactions of gp120 region adjacent to loop V3 with other region of gp120, which are unfavorable for conformational rearrangements of the region adjacent to loop V3 and further binding the coreceptor. Additionally, the interactions of loop V3 and bridging sheet with X5 lead to the close motion of loop V3 in X5 bound form, which further influences the rearrangements in gp120. [ABSTRACT FROM AUTHOR]

Background The interactions of complement receptor 2 (CR2) and the degradation fragment C3d of complement component C3 play important links between the innate and adaptive immune systems. Due to the importance of C3d–CR2 interaction in the design of vaccines and inhibitors, a number of studies have been performed to investigate C3d–CR2 interaction. Many studies have indicated C3d–CR2 interactions are ionic strength-dependent. Methods To investigate the molecular mechanism of C3d–CR2 interaction and the origin of effects of ionic strength, molecular dynamics simulations for C3d–CR2 complex together with the energetic and structural analysis were performed. Results Our results revealed the increased interactions between charged protein and ions weaken C3d–CR2 association, as ionic strengths increase. Moreover, ion strengths have similar effects on antigen-binding site and CR2 binding site. Meanwhile, Ala17 and Gln20 will transform between the activated and non-activated states mediated by His133 and Glu135 at different ion strengths. Conclusions Our results reveal the origins of the effects of ionic strengths on C3d–CR2 interactions are due to the changes of water, ion occupancies and distributions. General significance This study uncovers the origin of the effect of ionic strength on C3d–CR2 interaction and deepens the understanding of the molecular mechanism of their interaction, which is valuable for the design of vaccines and small molecule inhibitors. [ABSTRACT FROM AUTHOR]

Virus-specific cytotoxic T lymphocytes contribute to the control of virus infections including those caused by influenza viruses. However, during the evolution of influenza A viruses, variations in cytotoxic T lymphocytes epitopes have been observed and it will affect the recognition by virus-specific cytotoxic T lymphocytes and the human virus-specific cytotoxic T lymphocytes response in vitro. Here, to gain further insights into the molecular mechanism of the virus-specific cytotoxic T lymphocytes immunity, the class I major histocompatibility complex-encoded HLA-B*3501 protein with six different NP418-426 antigenic peptides emerging from 1918 to 2009 pandemic influenza A virus were studied by molecular dynamics simulation. Dynamical and structural properties (such as atomic fluctuations, solvent-accessible surface areas, binding free energy), based on the solvated protein-peptide complexes, were compared. Free energy calculations emphasized the important role of the secondary anchors (positions 2 and 9) in influencing the binding of MHC-I with antigenic non-apeptides. Furthermore, major interactions with peptides were gained from HLA-B*3501 residues: Tyr7, Ile66, Lys146, Trp147, and Tyr159. Detailed analysis could help to understand how different NP418-426 mutants effectively bind with the HLA-B*3501. [ABSTRACT FROM AUTHOR]

The outer membrane (OM) is the first defence for Gram-negative bacteria against their environments making it important in strain improvement for sustainable biobutanol production. While modifying the OM structure using chemical additives could enhance microbial viability, there are currently no model systems that accurately describe OM responses to butanol. Here, we experimentally determined that reducing the lipopolysaccharide (LPS) core length and charge increased Escherichia coli sensitivity to butanol. In silico models were built to describe how OM structure contributes to its ability to withstand butanol under conditions of exposure and production. Consistent with experiments, resistance to ingress of butanol into OMs correlates with both core length and charge, where a lower charge density is more conducive to butanol assimilation. The core length and branching correlate with the lateral spacing of the lipids, suggestive of a role of them in maintaining OM fluidity. In contrast to systems with short-length LPS cores, butanol intercalation into OMs with longer LPS cores increases membrane order and rigidity, which might be due to their more porous internal structure. These findings will assist the development of more butanol-tolerant biobutanol-producing bacteria, where thicker, more compact and less polar LPS-core surfaces reinforce the integrity of OMs and further improve resilience in extreme environments. Unlabelled Image • Our OM model describes the effects of LPS charge and length on E. coli sensitivity to butanol • Outer membrane acyl chain fluidity is inversely proportional to head group packing • Butanol can increase acyl chain fluidity by disrupting head group interactions • Effect of butanol intercalation on acyl chain fluidity depends on LPS core length and charge • A model is described to inform OM interventions for increasing the tolerance of Gram-negative bacteria to butanol [ABSTRACT FROM AUTHOR]