SARS-CoV-2 spike fusion peptide trans interaction with phosphatidylserine lipid triggers membrane fusion for viral entry.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike is the fusion machine for host cell entry. Still, the mechanism by which spike protein interacts with the target lipid membrane to facilitate membrane fusion during entry is not fully understood. Here, using steady-state membrane fusion and single-molecule fluorescence resonance energy transfer imaging of spike trimers on the surface of SARS-CoV-2 pseudovirion, we directly show that spike protein interacts with phosphatidylserine (PS) lipid in the target membrane for mediating fusion. We observed that the fusion peptide of the spike S2 domain interacts with the PS lipid of the target membrane. Low pH and Ca2+ trigger the spike conformational change and bring fusion peptide in close proximity to the PS lipid of the membrane. The binding of the spike with PS lipid of its viral membrane (cis interaction) impedes the fusion activation. PS on the target membrane promotes spike binding via trans interaction, prevents the cis interaction, and accelerates fusion. Sequestering or absence of PS lipid abrogates the spike-mediated fusion process and restricts SARS-CoV-2 infectivity. We found that PS-dependent interaction for fusion is conserved across all the SARS-CoV-2 spike variants of concern (D614G, Alpha, Beta, Delta, and Omicron). Our study suggests that PS lipid is indispensable for SARS-CoV-2 spike-mediated virus and target membrane fusion for entry, and restricting PS interaction with spike inhibits the SARS-CoV-2 spike-mediated entry. Therefore, PS is an important cofactor and acts as a molecular beacon in the target membrane for SARS-CoV-2 entry. IMPORTANCE: The role of lipids in the host cell target membrane for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) entry is not clear. We do not know whether SARS-CoV-2 spike protein has any specificity in terms of lipid for membrane fusion reaction. Here, using in vitro reconstitution of membrane fusion assay and single-molecule fluorescence resonance energy transfer imaging of SARS-CoV-2 spike trimers on the surface of the virion, we have demonstrated that phosphatidylserine (PS) lipid plays a key role in SARS-CoV-2 spike-mediated membrane fusion reaction for entry. Membrane-externalized PS lipid strongly promotes spike-mediated membrane fusion and COVID-19 infection. Blocking externalized PS lipid with PS-binding protein or in the absence of PS, SARS-CoV-2 spike-mediated fusion is strongly inhibited. Therefore, PS is an important target for restricting viral entry and intervening spike, and PS interaction presents new targets for COVID-19 interventions.

Bifidobacterium bifidum SAM-VI Riboswitch Conformation Change Requires Peripheral Helix Formation.

The Bifidobacterium bifidum SAM-VI riboswitch undergoes dynamic conformational changes that modulate downstream gene expression. Traditional structural methods such as crystallography capture the bound conformation at high resolution, and additional efforts would reveal details from the dynamic transition. Here, we revealed a transcription-dependent conformation model for Bifidobacterium bifidum SAM-VI riboswitch. In this study, we combine small-angle X-ray scattering, chemical probing, and isothermal titration calorimetry to unveil the ligand-binding properties and conformational changes of the Bifidobacterium bifidum SAM-VI riboswitch and its variants. Our results suggest that the SAM-VI riboswitch contains a pre-organized ligand-binding pocket and stabilizes into the bound conformation upon binding to SAM. Whether the P1 stem formed and variations in length critically influence the conformational dynamics of the SAM-VI riboswitch. Our study provides the basis for artificially engineering the riboswitch by manipulating its peripheral sequences without modifying the SAM-binding core.

Conformational dynamics underlying atypical chemokine receptor 3 activation.

Atypical Chemokine Receptor 3 (ACKR3) belongs to the G protein-coupled receptor family but it does not signal through G proteins. The structural properties that govern the functional selectivity and the conformational dynamics of ACKR3 activation are poorly understood. Here, we combined hydrogen/deuterium exchange mass spectrometry, site-directed mutagenesis, and molecular dynamics simulations to examine the binding mode and mechanism of action of ACKR3 ligands of different efficacies. Our results show that activation or inhibition of ACKR3 is governed by intracellular conformational changes of helix 6, intracellular loop 2, and helix 7, while the DRY motif becomes protected during both processes. Moreover, we identified the binding sites and the allosteric modulation of ACKR3 upon ß-arrestin 1 binding. In summary, this study highlights the structure-function relationship of small ligands, the binding mode of ß-arrestin 1, the activation dynamics, and the atypical dynamic features in ACKR3 that may contribute to its inability to activate G proteins.

Unleashing the Power of Evolution in Xylanase Engineering: Investigating the Role of Distal Mutation Regulation.

The drive to enhance enzyme performance in industrial applications frequently clashes with the practical limitations of exhaustive experimental screening, underscoring the urgency for more refined and strategic methodologies in enzyme engineering. In this study, xylanase Xyl-1 was used as the model, coupling evolutionary insights with energy functions to obtain theoretical potential mutants, which were subsequently validated experimentally. We observed that mutations in the nonloop region primarily aimed at enhancing stability and also encountered selective pressure for activity. Notably, mutations in this region simultaneously boosted the Xyl-1 stability and activity, achieving a 65% success rate. Using a greedy strategy, mutant M4 was developed, achieving a 12 °C higher melting temperature and doubled activity. By integration of spectroscopy, crystallography, and quantum mechanics/molecular mechanics molecular dynamics, the mechanism behind the enhanced thermal stability of M4 was elucidated. It was determined that the activity differences between M4 and the wild type were primarily driven by dynamic factors influenced by distal mutations. In conclusion, the study emphasizes the pivotal role of evolution-based approaches in augmenting the stability and activity of the enzymes. It sheds light on the unique adaptive mechanisms employed by various structural regions of proteins and expands our understanding of the intricate relationship between distant mutations and enzyme dynamics.

Incorporating Sequence-Dependent DNA Shape and Dynamics into Transcriptome Data Analysis.

Differentially expressed genes in a cellular context may be co-regulated by the same transcription factor. However, in the absence of a concurrent transcription factor binding data, such interactions are difficult to detect, especially at the single cell expression level. Motif enrichments in such genes can be used to gain insight into differential expressions caused by the shared upstream TFs. However, it is now established that many genes are co-regulated by the same TF due to a shared DNA shape or sequence-dependent conformational dynamics instead of sequence motif. In this work, we demonstrate how, starting from a gene expression data, such DNA shape and dynamics signatures can be potentially detected using publicly available tools, including DynaSeq, developed in our group for predicting the sequence-dependent components of these DNA shape features.

Conformational Dynamics of Influenza A Virus Ribonucleoprotein Complexes during RNA Synthesis.

Viral ribonucleoproteins (vRNPs) are the cornerstones of viral proliferation, as they form the macromolecular complexes that are responsible for the transcription and replication of most single-stranded RNA viruses. The influenza A virus (IAV) polymerase catalyzes RNA synthesis within the context of vRNPs where genomic viral RNA (vRNA) is packaged by the viral nucleoprotein (NP). We used high-speed atomic force microscopy and electron microscopy to study the conformational dynamics of individual IAV recombinant RNPs (rRNPs) during RNA synthesis. The rRNPs present an annular organization that allows for the real-time tracking of conformational changes in the NP-vRNA template caused by the advancing polymerase. We demonstrate that the rRNPs undergo a well-defined conformational cycle during RNA synthesis, which can be interpreted in light of previous transcription models. We also present initial estimations of the average RNA synthesis rate in the rRNP and its dependence on the nucleotide concentration and stability of the nascent RNA secondary structures. Furthermore, we provide evidence that rRNPs can perform consecutive cycles of RNA synthesis, accounting for their ability to recycle and generate multiple copies of RNA.

Hidden Structural States of Proteins Revealed by Conformer Selection with AlphaFold-NMR.

Recent advances in molecular modeling using deep learning can revolutionize our understanding of dynamic protein structures. NMR is particularly well-suited for determining dynamic features of biomolecular structures. The conventional process for determining biomolecular structures from experimental NMR data involves its representation as conformation-dependent restraints, followed by generation of structural models guided by these spatial restraints. Here we describe an alternative approach: generating a distribution of realistic protein conformational models using artificial intelligence-(AI-) based methods and then selecting the sets of conformers that best explain the experimental data. We applied this conformational selection approach to redetermine the solution NMR structure of the enzyme Gaussia luciferase. First, we generated a diverse set of conformer models using AlphaFold2 (AF2) with an enhanced sampling protocol. The models that best-fit NOESY and chemical shift data were then selected with a Bayesian scoring metric. The resulting models include features of both the published NMR structure and the standard AF2 model generated without enhanced sampling. This "AlphaFold-NMR" protocol also generated an alternative "open" conformational state that fits nearly as well to the overall NMR data but accounts for some NOESY data that is not consistent with first "closed" conformational state; while other NOESY data consistent with this second state are not consistent with the first conformational state. The structure of this "open" structural state differs from that of the "closed" state primarily by the position of a thumb-shaped loop between α-helices H5 and H6, revealing a cryptic surface pocket. These alternative conformational states of Gluc are supported by "double recall" analysis of NOESY data and AF2 models. Additional structural states are also indicated by backbone chemical shift data indicating partially-disordered conformations for the C-terminal segment. Considered as a multistate ensemble, these multiple states of Gluc together fit the NOESY and chemical shift data better than the "restraint-based" NMR structure and provide novel insights into its structure-dynamic-function relationships. This study demonstrates the potential of AI-based modeling with enhanced sampling to generate conformational ensembles followed by conformer selection with experimental data as an alternative to conventional restraint satisfaction protocols for protein NMR structure determination.

Ultrafast conformation-dependent charge transfer in N, N, N', N'-tetramethyl-1,3-propanediamine: Effect of flexible carbon skeleton on electron lone pair interactions.

Flexible three-carbon skeleton makes N, N, N', N'-tetramethyl-1,3-propanediamine (TMPDA) an important diamine system to investigate the conformation-dependent electron lone pair interactions and charge delocalization. The charge transfer process linked to structural motions of the three-carbon skeleton has been monitored in real time by the Rydberg electron binding energy (BE) spectra of TMPDA coupled with quantum chemical calculations. Optical excitation to the 3p state with a 200 nm pump pulse initially generated a localized charge on one of the two nitrogen atoms that may partially transfer to the other one. Rapid internal conversion (IC) from the 3p to 3s state occurred within 430 fs, resulting in an initial charge delocalized 3s_h/3s_l population ratio of 23.6 %/76.4 %. A final 3s_h/3s_l (51.9 %/48.1 %) equilibrium proceeded within about 2.64 ps. The 3s_h (TTTT+, GG'TG+ and G'GG'G+) and 3s_l (GG'GG'+ and GG'G'G+) (see text for structure definitions) are identified as the extended and folded conformers, respectively. Two types of electron lone pair interactions, i.e., through-space interaction (TSI) and through-bond interaction (TBI), are found to coexist in TMPDA to drive charge transfer. The GG'GG'+ and GG'G'G+ structures exhibit TSI, while the TTTT+ structure shows TBI. The GG'TG+ and G'GG'G+ structures exhibit both TSI and TBI. Flexible three-carbon skeleton provide more opportunities for the two N-electron lone pairs to overlap in space (i.e., TSI), making TMPDA to be favorable for the most stably folded conformation.

Ligand binding initiates single-molecule integrin conformational activation.

Integrins link the extracellular environment to the actin cytoskeleton in cell migration and adhesiveness. Rapid coordination between events outside and inside the cell is essential. Single-molecule fluorescence dynamics show that ligand binding to the bent-closed integrin conformation, which predominates on cell surfaces, is followed within milliseconds by two concerted changes, leg extension and headpiece opening, to give the high-affinity integrin conformation. The extended-closed integrin conformation is not an intermediate but can be directly accessed from the extended-open conformation and provides a pathway for ligand dissociation. In contrast to ligand, talin, which links the integrin ß-subunit cytoplasmic domain to the actin cytoskeleton, modestly stabilizes but does not induce extension or opening. Integrin activation is thus initiated by outside-in signaling and followed by inside-out signaling. Our results further imply that talin binding is insufficient for inside-out integrin activation and that tensile force transmission through the ligand-integrin-talin-actin cytoskeleton complex is required.

Structural Basis of Conformational Dynamics in the PROTAC-Induced Protein Degradation.

Pronounced conformational dynamics is unveiled upon analyzing multiple crystal structures of the same proteins recruited to the same E3 ligases by PROTACs, and yet, is largely permissive for targeted protein degradation due to the intrinsic mobility of E3 assemblies creating a large ubiquitylation zone. Mathematical modelling of ternary dynamics on ubiquitylation probability confirms the experimental finding that ternary complex rigidification need not correlate with enhanced protein degradation. Salt bridges are found to prevail in the PROTAC-induced ternary complexes, and may contribute to a positive cooperativity and prolonged half-life. The analysis highlights the importance of presenting lysines close to the active site of the E2 enzyme while constraining ternary dynamics in PROTAC design to achieve high degradation efficiency.

The effect of chrysin binding on the conformational dynamics and unfolding pathway of human serum albumin.

Studies on the interactions between ligands and proteins provide insights into how a possible medication alters the structures and activities of the target or carrier proteins. The natural flavonoid aglycone Chrysin (CHR) has demonstrated anti-inflammatory, antioxidant, antiapoptotic, neuroprotective, and antineoplastic effects, both in vitro and in vivo. In this work, we investigated the impact of CHR binding on the as-yet-unexplored conformation, dynamics, and unfolding mechanism of human serum albumin (HSA). We determined CHR binding to HSA domain-II with the association constant (Ka) of 2.70 ± 0.21 × 105 M-1. The urea-induced sequential unfolding mechanism of HSA was used to elucidate the debatable binding location of CHR. CHR binding induced both secondary and tertiary structural alterations in the protein as studied by far-UV circular dichroism and intrinsic fluorescence spectroscopy. Red edge excitation shift (REES) indicated a decrease in conformational dynamics of the protein on the complex formation. This suggested an ordered compact and spatial arrangement of the CHR-boundmolecule. The binding of CHR was found to significantly modulate the urea-induced unfolding pathway of HSA. Urea-induced unfolding pathway of HSA became a two-state process (N-U) from a three-state process (N-I-U). The interaction of CHR is found to increase the thermal stability of the protein by ∼4 °C. This study focuses on the fundamental sciences and demonstrates how prospective medication compounds can alter the dynamics and stability of protein structure.

Illuminating ligand-induced dynamics in nuclear receptors through MD simulations.

Nuclear receptors (NRs) regulate gene expression in critical physiological processes, with their functionality finely tuned by ligand-induced conformational changes. While NRs may sometimes undergo significant conformational motions in response to ligand-binding, these effects are more commonly subtle and challenging to study by traditional structural or biophysical methods. Molecular dynamics (MD) simulations are a powerful tool to bridge the gap between static protein-ligand structures and dynamical changes that govern NR function. Here, we summarize a handful of recent studies that apply MD simulations to study NRs. We present diverse methodologies for analyzing simulation data with a detailed examination of the information each method can yield. By delving into the strengths, limitations and unique contributions of these tools, this review provides guidance for extracting meaningful data from MD simulations to advance the goal of understanding the intricate mechanisms by which ligands orchestrate a range of functional outcomes in NRs.

Characterizing the Conformational Dynamics of Human SUMO2: Insights into its Interaction with Metal Ions and SIMs.

SUMO (Small Ubiquitin-like Modifiers) proteins are involved in a crucial post-translational modification commonly termed as SUMOylation. In this work, we have investigated the native-state conformational flexibility of human SUMO2 and its interaction with Cu2+ and Zn2+ ions using 15N-1H based 2D NMR spectroscopy. After SUMO1, SUMO2 is the most studied SUMO isoform in humans which shares 45 % and ~80 % similarity with SUMO1 in terms of sequence and structure, respectively. In this manuscript, we demonstrate that compared to SUMO1, several amino acids around the α1-helix region of SUMO2 access energetically similar near-native conformations. These conformations could play a crucial role in SUMO2's non-covalent interactions with SUMO interaction motifs (SIMs) on other proteins. The C-terminal of SUMO2 was found to bind strongly with Cu2+ ions resulting in a trimeric structure as observed by gel electrophoresis. This interaction seems to interfere in its non-covalent interaction with a V/I-x-V/I-V/I based SIM in Daxx protein.

Four Parallel Pathways in T4 Ligase-Catalyzed Repair of Nicked DNA with Diverse Bending Angles.

The structural diversity of biological macromolecules in different environments contributes complexity to enzymological processes vital for cellular functions. Fluorescence resonance energy transfer and electron microscopy are used to investigate the enzymatic reaction of T4 DNA ligase catalyzing the ligation of nicked DNA. The data show that both the ligase-AMP complex and the ligase-AMP-DNA complex can have four conformations. This finding suggests the parallel occurrence of four ligation reaction pathways, each characterized by specific conformations of the ligase-AMP complex that persist in the ligase-AMP-DNA complex. Notably, these complexes have DNA bending angles of ≈0°, 20°, 60°, or 100°. The mechanism of parallel reactions challenges the conventional notion of simple sequential reaction steps occurring among multiple conformations. The results provide insights into the dynamic conformational changes and the versatile attributes of T4 DNA ligase and suggest that the parallel multiple reaction pathways may correspond to diverse T4 DNA ligase functions. This mechanism may potentially have evolved as an adaptive strategy across evolutionary history to navigate complex environments.

Peptide-based allosteric inhibitor targets TNFR1 conformationally active region and disables receptor-ligand signaling complex.

Tumor necrosis factor (TNF) receptor 1 (TNFR1) plays a pivotal role in mediating TNF induced downstream signaling and regulating inflammatory response. Recent studies have suggested that TNFR1 activation involves conformational rearrangements of preligand assembled receptor dimers and targeting receptor conformational dynamics is a viable strategy to modulate TNFR1 signaling. Here, we used a combination of biophysical, biochemical, and cellular assays, as well as molecular dynamics simulation to show that an anti-inflammatory peptide (FKCRRWQWRMKK), which we termed FKC, inhibits TNFR1 activation allosterically by altering the conformational states of the receptor dimer without blocking receptor-ligand interaction or disrupting receptor dimerization. We also demonstrated the efficacy of FKC by showing that the peptide inhibits TNFR1 signaling in HEK293 cells and attenuates inflammation in mice with intraperitoneal TNF injection. Mechanistically, we found that FKC binds to TNFR1 cysteine-rich domains (CRD2/3) and perturbs the conformational dynamics required for receptor activation. Importantly, FKC increases the frequency in the opening of both CRD2/3 and CRD4 in the receptor dimer, as well as induces a conformational opening in the cytosolic regions of the receptor. This results in an inhibitory conformational state that impedes the recruitment of downstream signaling molecules. Together, these data provide evidence on the feasibility of targeting TNFR1 conformationally active region and open new avenues for receptor-specific inhibition of TNFR1 signaling.

Structural biology in the age of X-ray free-electron lasers and exascale computing.

Serial femtosecond X-ray crystallography has emerged as a powerful method for investigating biomolecular structure and dynamics. With the new generation of X-ray free-electron lasers, which generate ultrabright X-ray pulses at megahertz repetition rates, we can now rapidly probe ultrafast conformational changes and charge movement in biomolecules. Over the last year, another innovation has been the deployment of Frontier, the world's first exascale supercomputer. Synergizing extremely high repetition rate X-ray light sources and exascale computing has the potential to accelerate discovery in biomolecular sciences. Here we outline our perspective on each of these remarkable innovations individually, and the opportunities and challenges in yoking them within an integrated research infrastructure.

Mispacking of the F87 sidechain drives aggregation-promoting conformational fluctuations in the subunit interfaces of the transthyretin tetramer.

Aberrant formation and deposition of human transthyretin (TTR) aggregates causes transthyretin amyloidosis. To initialize aggregation, transthyretin tetramers must first dissociate into monomers that partially unfold to promote entry into the aggregation pathway. The native TTR tetramer (T) is stabilized by docking of the F87 sidechain into an interfacial cavity enclosed by several hydrophobic residues including A120. We have previously shown that an alternative tetramer (T*) with mispacked F87 sidechains is more prone to dissociation and aggregation than the native T state. However, the molecular basis for the reduced stability in T* remains unclear. Here we report characterization of the A120L mutant, where steric hindrance is introduced into the F87 binding site. The X-ray structure of A120L shows that the F87 sidechain is displaced from its docking site across the subunit interface. In A120S, a naturally occurring pathogenic mutant that is less aggregation-prone than A120L, the F87 sidechain is correctly docked, as in the native TTR tetramer. Nevertheless, 19F-NMR aggregation assays show an elevated population of a monomeric aggregation intermediate in A120S relative to a control containing the native A120, due to accelerated tetramer dissociation and slowed monomer tetramerization. The mispacking of the F87 sidechain is associated with enhanced exchange dynamics for interfacial residues. At 298 K, the T* populations of various naturally occurring mutants fall between 4-7% (ΔG ~ 1.5-1.9 kcal/mol), consistent with the free energy change expected for undocking and solvent exposure of one of the four F87 sidechains in the tetramer (ΔG ~ 1.6 kcal/mol). Our data provide a molecular-level picture of the likely universal F87 sidechain mispacking in tetrameric TTR that promotes interfacial conformational dynamics and increases aggregation propensity.

Mutations influence the conformational dynamics of the GDP/KRAS complex.

Mutations near allosteric sites can have a significant impact on the function of KRAS. Three specific mutations, K104Q, G12D/K104Q, and G12D/G75A, which are located near allosteric positions, were selected to investigate the molecular mechanisms behind mutation-induced influences on the activity of KRAS. Gaussian accelerated molecular dynamics (GaMD) simulations followed by the principal component analysis (PCA) were performed to improve the sampling of conformational states. The results revealed that these mutations significantly alter the structural flexibility, correlated motions, and dynamic behavior of the switch regions that are essential for KRAS binding to effectors or regulators. Furthermore, the mutations have a significant impact on the hydrogen bonding interactions between GDP and the switch regions, as well as on the electrostatic interactions of magnesium ions (Mg2+) with these regions. Our results verified that these mutations strongly influence the binding of KRAS to its effectors or regulators and allosterically regulate the activity. We believe that this work can provide valuable theoretical insights into a deeper understanding of KRAS function.Communicated by Ramaswamy H. Sarma.

Deciphering the molecular choreography of Janus kinase 2 inhibition via Gaussian accelerated molecular dynamics simulations: a dynamic odyssey.

The Janus kinases (JAK) are crucial targets in drug development for several diseases. However, accounting for the impact of possible structural rearrangements on the binding of different kinase inhibitors is complicated by the extensive conformational variability of their catalytic kinase domain (KD). The dynamic KD contains mainly four prominent mobile structural motifs: the phosphate-binding loop (P-loop), the αC-helix within the N-lobe, the Asp-Phe-Gly (DFG) motif, and the activation loop (A-loop) within the C-lobe. These distinct structural orientations imply a complex signal transmission path for regulating the A-loop's flexibility and conformational preference for optimal JAK function. Nevertheless, the precise dynamical features of the JAK induced by different types of inhibitors still remain elusive. We performed comparative, microsecond-long, Gaussian accelerated molecular dynamics simulations in triplicate of three phosphorylated JAK2 systems: the KD alone, type-I ATP-competitive inhibitor (CI) bound KD in the catalytically active DFG-in conformation, and the type-II inhibitor (AI) bound KD in the catalytically inactive DFG-out conformation. Our results indicate significant conformational variations observed in the A-loop and αC helix motions upon inhibitor binding. Our studies also reveal that the DFG-out inactive conformation is characterized by the closed A-loop rearrangement, open catalytic cleft of N and C-lobe, the outward movement of the αC helix, and open P-loop states. Moreover, the outward positioning of the αC helix impacts the hallmark salt bridge formation between Lys882 and Glu898 in an inactive conformation. Finally, we compared their ligand binding poses and free energy by the MM/PBSA approach. The free energy calculations suggested that the AI's binding affinity is higher than CI against JAK2 due to an increased favorable contribution from the total non-polar interactions and the involvement of the αC helix. Overall, our study provides the structural and energetic insights crucial for developing more promising type I/II JAK2 inhibitors for treating JAK-related diseases.

Studying Conformational Properties of Transmembrane Domain of KCNE3 in a Lipid Bilayer Membrane Using Molecular Dynamics Simulations.

KCNE3 is a single-pass integral membrane protein that regulates numerous voltage-gated potassium channel functions such as KCNQ1. Previous solution NMR studies suggested a moderate degree of curved α-helical structure in the transmembrane domain (TMD) of KCNE3 in lyso-myristoylphosphatidylcholine (LMPC) micelles and isotropic bicelles with the residues T71, S74 and G78 situated along the concave face of the curved helix. During the interaction of KCNE3 and KCNQ1, KCNE3 pushes its transmembrane domain against KCNQ1 to lock the voltage sensor in its depolarized conformation. A cryo-EM study of KCNE3 complexed with KCNQ1 in nanodiscs suggested a deviation of the KCNE3 structure from its independent structure in isotropic bicelles. Despite the biological significance of KCNE3 TMD, the conformational properties of KCNE3 are poorly understood. Here, all atom molecular dynamics (MD) simulations were utilized to investigate the conformational dynamics of the transmembrane domain of KCNE3 in a lipid bilayer containing a mixture of POPC and POPG lipids (3:1). Further, the effect of the interaction impairing mutations (V72A, I76A and F68A) on the conformational properties of the KCNE3 TMD in lipid bilayers was investigated. Our MD simulation results suggest that the KCNE3 TMD adopts a nearly linear α helical structural conformation in POPC-POPG lipid bilayers. Additionally, the results showed no significant change in the nearly linear α-helical conformation of KCNE3 TMD in the presence of interaction impairing mutations within the sampled time frame. The KCNE3 TMD is more stable with lower flexibility in comparison to the N-terminal and C-terminal of KCNE3 in lipid bilayers. The overall conformational flexibility of KCNE3 also varies in the presence of the interaction-impairing mutations. The MD simulation data further suggest that the membrane bilayer width is similar for wild-type KCNE3 and KCNE3 containing mutations. The Z-distance measurement data revealed that the TMD residue site A69 is close to the lipid bilayer center, and residue sites S57 and S82 are close to the surfaces of the lipid bilayer membrane for wild-type KCNE3 and KCNE3 containing interaction-impairing mutations. These results agree with earlier KCNE3 biophysical studies. The results of these MD simulations will provide complementary data to the experimental outcomes of KCNE3 to help understand its conformational dynamic properties in a more native lipid bilayer environment.