Microsecond Molecular Dynamics Research Articles

Elucidation of the Decomplexation and Reverse Migration Mechanism of Uranyl Ions from the Oil Phase to the Aqueous Phase Using Microsecond (0.2 μs) Molecular Dynamics Simulations.

Interphase ion transfer is ubiquitous in chemistry, physics, biology, and various engineering sciences. Ion transfer from the aqueous phase to the oil phase or vice versa is a complex chemical phenomenon, and its fundamental understanding is crucial for efficient and economical mass transfer. This ion transfer is much more complex for radionuclide metal ions. Therefore, an attempt has been made to elucidate the decomplexation and mechanism of reverse migration of uranyl ions from the loaded oil phase to the aqueous phase using large time-scale molecular dynamics simulations (microseconds) and density functional theory (DFT). The strength of the metal-ligand complex is the key factor for stripping among other mass transfer parameters. Stronger the metal-ligand complex, the lower the tendency for reverse migration into the aqueous phase, which was demonstrated by three different ligands for reverse migration using water as the stripping agent. The interaction energy of the metal-ligand complex has been calculated using DFT. The reverse migration is validated by the distance between the centers of mass. Higher interface thickness and lower interfacial tension belong to N,N,N',N'-tetra-octyldiglycolamide (TODGA), which are favorable for mass transfer across the liquid-liquid interface. The computed distribution coefficient (KD) is favorable for tributyl phosphate (TBP) and tri-iso-amyl phosphate (TiAP), whereas TODGA shows a higher KD, indicating a less favorable value for reverse migration. The distance between the uranyl ion and the TODGA molecule confirms that the entrapment of uranyl ions in the TODGA phase might be attributed to aggregation. Further, the aggregation tendency of TODGA molecules reduces the back extraction and the recovery of uranyl ions into the aqueous phase. The present MD results might be important for predicting an efficient back-extracting agent for the recovery of the metal ions from the oil phase.

Cobaltabis(dicarbollide) Interaction with DNA Resolved at the Atomic Scale.

Boron neutron capture therapy represents a promising avenue for cancer treatment that requires nontoxic drugs with a high boron content efficiently distributed into cancerous cells. The metallacarborane o-cobaltabis(dicarbollide) ([COSAN]-) fulfills these requirements and constitutes an attractive candidate. Nevertheless, the interaction of this promising drug with nucleic acids, the assumed target of the biological damage, is poorly understood since contradictory results are reported in the literature. This work establishes the DNA/[COSAN]- interaction strength, mechanism, and time scale at the atomistic level by using a combination of microsecond-molecular dynamics and hybrid quantum mechanics/molecular mechanics simulations and by quantifying the absolute binding free energy. Results show that the DNA/[COSAN]- interaction is highly dependent on the ionic strength of the medium. A relatively weak DNA major groove binding (ΔGbind= -2.49 kcal/mol) driven mostly by dihydrogen B-H···H-N bonding is observed in the simulations only at a high NaCl concentration, whereas DNA intercalation mode is deemed highly unlikely.

Building predictive Markov models of ion channel permeation from molecular dynamics simulations

Molecular dynamics (MD) simulation of biological processes has always been a challenging task due to the long timescales of the processes involved and the large amount of output data to handle. Markov state models (MSMs) have been introduced as a powerful tool in this area of research, as they provide a mechanistically comprehensible synthesis of the large amount of MD data and, at the same time, can be used to rapidly estimate experimental properties of biological processes. Herein, we propose a method for building MSMs of ion channel permeation from MD trajectories, which directly evaluates the current flowing through the channel from the model’s transition matrix (T), which is crucial for comparing simulations and experimental data. This is achieved by including in the model a flux matrix that summarizes information on the charge moving across the channel between each pair of states of the MSM and can be used in conjunction with T to predict the ion current. A procedure to drastically reduce the number of states in the MSM while preserving the estimated ion current is also proposed. Applying the method to the KcsA channel returned an MSM with five states with significant equilibrium occupancy, capable of accurately reproducing the single-channel ion current from microsecond MD trajectories.

Parsing Dynamics of Protein Backbone NH and Side-Chain Methyl Groups using Molecular Dynamics Simulations.

Experimental NMR spectroscopy and theoretical molecular dynamics (MD) simulations provide complementary insights into protein conformational dynamics and hence into biological function. The present work describes an extensive set of backbone NH and side-chain methyl group generalized order parameters for the Escherichia coli ribonuclease HI (RNH) enzyme derived from 2-μs microsecond MD simulations using the OPLS4 and AMBER-FF19SB force fields. The simulated generalized order parameters are compared with values derived from NMR 15N and 13CH2D spin relaxation measurements. The squares of the generalized order parameters, S2 for the N-H bond vector and Saxis2 for the methyl group symmetry axis, characterize the equilibrium distribution of vector orientations in a molecular frame of reference. Optimal agreement between simulated and experimental results was obtained by averaging S2 or Saxis2 calculated by dividing the simulated trajectories into 50 ns blocks (∼five times the rotational diffusion correlation time for RNH). With this procedure, the median absolute deviations (MAD) between experimental and simulated values of S2 and Saxis2 are 0.030 (NH) and 0.061 (CH3) for OPLS4 and 0.041 (NH) and 0.078 (CH3) for AMBER-FF19SB. The MAD between OPLS4 and AMBER-FF19SB are 0.021 (NH) and 0.072 (CH3). The generalized order parameters for the methyl group symmetry axis can be decomposed into contributions from backbone fluctuations, between-rotamer dihedral angle transitions, and within-rotamer dihedral angle fluctuations. Analysis of the simulation trajectories shows that (i) backbone and side chain conformational fluctuations exhibit little correlation and that (ii) fluctuations within rotamers are limited and highly uniform with values that depend on the number of dihedral angles considered. Low values of Saxis2, indicative of enhanced side-chain flexibility, result from between-rotamer transitions that can be enhanced by increased local backbone flexibility.

AlphaFold2 Predictions of Conformational Ensembles and Atomistic Simulations of the SARS-CoV-2 Spike XBB Lineages Reveal Epistatic Couplings between Convergent Mutational Hotspots that Control ACE2 Affinity.

In this study, we combined AlphaFold-based atomistic structural modeling, microsecond molecular simulations, mutational profiling, and network analysis to characterize binding mechanisms of the SARS-CoV-2 spike protein with the host receptor ACE2 for a series of Omicron XBB variants including XBB.1.5, XBB.1.5+L455F, XBB.1.5+F456L, and XBB.1.5+L455F+F456L. AlphaFold-based structural and dynamic modeling of SARS-CoV-2 Spike XBB lineages can accurately predict the experimental structures and characterize conformational ensembles of the spike protein complexes with the ACE2. Microsecond molecular dynamics simulations identified important differences in the conformational landscapes and equilibrium ensembles of the XBB variants, suggesting that combining AlphaFold predictions of multiple conformations with molecular dynamics simulations can provide a complementary approach for the characterization of functional protein states and binding mechanisms. Using the ensemble-based mutational profiling of protein residues and physics-based rigorous calculations of binding affinities, we identified binding energy hotspots and characterized the molecular basis underlying epistatic couplings between convergent mutational hotspots. Consistent with the experiments, the results revealed the mediating role of the Q493 hotspot in the synchronization of epistatic couplings between L455F and F456L mutations, providing a quantitative insight into the energetic determinants underlying binding differences between XBB lineages. We also proposed a network-based perturbation approach for mutational profiling of allosteric communications and uncovered the important relationships between allosteric centers mediating long-range communication and binding hotspots of epistatic couplings. The results of this study support a mechanism in which the binding mechanisms of the XBB variants may be determined by epistatic effects between convergent evolutionary hotspots that control ACE2 binding.

Unveiling the Mechanisms of EGCG-p53 Interactions through Molecular Dynamics Simulations.

Green tea consumption is associated with protective and preventive effects against various types of cancer. These effects are attributed to polyphenols, particularly epigallocatechin-3-gallate (EGCG). EGCG acts by directly inhibiting tumor suppressor protein p53. The binding mechanism by which EGCG inhibits p53 activity is associated with residues Trp23-Lys24 and Pro47-Thr55 within the p53 N-terminal domain (NTD). However, the structural and thermodynamic aspects of the interaction between EGCG and p53 are poorly understood. Therefore, based on crystallographic data, we combine docking, molecular dynamics (MD) simulations, and molecular mechanics generalized Born surface area approaches to explore the intricacies of the EGCG-p53 binding mechanism. A triplicate microsecond MD simulation for each system is initially performed to capture diverse p53 NTD conformations. From the start, the most populated cluster of the second run (R2-1) stands out due to a unique opening between Trp23 and Trp53. During MD simulations, this conformation allows EGCG to sustain a high level of stability and affinity while interacting with both regions of interest and deepening the binding pocket. Structural analysis emphasizes the significance of pyrogallol motifs in EGCG binding. Therefore, the conformational shift in this gap is pivotal, enabling EGCG to impede p53 interactions and manifest its anticancer properties. These findings enhance the present comprehension of the anticancer properties of green tea polyphenols and pave the way for future therapeutic developments.

Dissecting the effect of ALS mutation S375G on the conformational properties and aggregation dynamics of TDP-43370-375 fragment

The aggregation of transactive response deoxyribonucleic acid (DNA) binding protein of 43 kDa (TDP-43) into ubiquitin-positive inclusions is closely associated with amyotrophic lateral sclerosis (ALS), frontotemporal lobar degeneration, and chronic traumatic encephalopathy. The 370–375 fragment of TDP-43 (370GNNSYS375, TDP-43370-375), the amyloidogenic hexapeptides, can be prone to forming pathogenic amyloid fibrils with the characteristic of steric zippers. Previous experiments reported the ALS-associated mutation, serine 375 substituted by glycine (S375G) is linked to early onset disease and protein aggregation of TDP-43. Based on this, it is necessary to explore the underlying molecular mechanisms. By utilizing all-atom molecular dynamics (MD) simulations of 102 μs in total, we investigated the impact of S375G mutation on the conformational ensembles and oligomerization dynamics of TDP-43370-375 peptides. Our replica exchange MD simulations show that S375G mutation could promote the unstructured conformation formation and induce peptides to form a loose packed oligomer, thus inhibiting the aggregation of TDP-43370-375. Further analyses suggest that S375G mutation displays a reduction effect on the number of total hydrogen bonds and contacts among TDP-43370-375 peptides. Hydrogen bonding and polar interactions among TDP-43370-375 peptides, as well as Y374-Y374 π-π stacking interaction, are attenuated by S375G mutation. Additional microsecond MD simulations demonstrate that S375G mutation could prohibit the conformational conversion to β-structure-rich aggregates and possess an inhibitory effect on the oligomerization dynamics of TDP-43370-375. This study offers for the first time of molecular insights into the S375G mutation affecting the aggregation of TDP-43370-375 at the atomic level, and may open new avenues in the development of future site-specific mutation therapeutics.

Abstract 4651: A macrocyclic kinase inhibitor overcomes the compound mutations that cause the resistance to all approved EGFR inhibitors in EGFR-positive lung cancer

Abstract EGFR mutant cancers account for 20-30 % of non-small cell lung cancers and show marked sensitivity to EGFR tyrosine kinase inhibitors (TKIs). However, the tumor almost inevitably relapses due to acquired resistance. To date, several EGFR-TKIs have been approved, and the third-generation EGFR-TKI, osimertinib, has been widely used both as first-line therapy and for patients with the EGFR-T790M mutation after 1st or 2nd generation EGFR-TKI therapy. However, the EGFR-T790M+C797S compound mutation confers resistance to all approved EGFR-TKIs. In our previous study, we experimentally demonstrated that brigatinib with anti-EGFR antibody was effective against osimertinib-resistant EGFR-C797S and EGFR-T790M+C797S mutants, and is currently undergoing clinical studies. However, the tumor would relapse because of the acquisition of additional resistance mutations. Here, we first demonstrated the binding mode of brigatinib to the EGFR-T790M/C797S mutant by crystal structure analysis and predicted brigatinib-resistant mutations through a cell-based assay, including N-ethyl-N-nitrosourea (ENU) mutagenesis. We found that clinically reported L718 and G796 compound mutations appeared, consistent with their proximity to the binding site of brigatinib, and brigatinib-resistant quadruple mutants, such as EGFR-activating mutation/T790M/C797S/L718M, were resistant to all clinically available EGFR-TKIs. BI-4020, a fourth-generation macrocyclic EGFR inhibitor plus anti-EGFR antibody, overcomes the quadruple and major EGFR-activating mutants in vitro and in vivo, but not the minor mutants such as L747P. Microsecond molecular dynamics simulations revealed the binding mode and affinity between BI-4020 and the EGFR mutants. This study identified potential therapeutic strategies using new-generation macrocyclic EGFR inhibitor to overcome the emerging resistance mutants. Citation Format: Ryohei Katayama, Ken Uchibori, Makoto Nishio, Naoya Fujita. A macrocyclic kinase inhibitor overcomes the compound mutations that cause the resistance to all approved EGFR inhibitors in EGFR-positive lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 4651.

Effect of Glucose on CH4 Hydrate Formation in Clay Nanopores and Bulk Solution: Insights from Microsecond Molecular Dynamics Simulations

Effect of Glucose on CH₄ Hydrate Formation in Clay Nanopores and Bulk Solution: Insights from Microsecond Molecular Dynamics Simulations

Local Structures in Proteins from Microsecond Molecular Dynamics Simulations: 2. The Role of Symmetry in GTPase Binding and Dimer Formation.

The Rho GTPase binding domain of Plexin-B1 (RBD) prevails in solution as dimer. Under appropriate circumstances, it binds the small GTPase Rac1 to yield the complex RBD-Rac1. Here, we study RBD dimerization and complex formation from a symmetry-based perspective using data derived from 1 μs long MD simulations. The quantities investigated are the local potentials, u(MD), prevailing at the N-H sites of the protein. These potentials are statistical in character providing an empirical description of the local structure. To establish more methodical description, a method for approximating them by explicit functions, u(simulated), was developed in the preceding article in this journal issue. These functions are combinations of analytical Wigner functions, DL,K, belonging to the D2h point group. The D2h subgroups Ag and B2u are found to dominate u(simulated); the B1u subgroup contributes in some cases. The Ag (B2u) functions have axial or rhombic symmetry. For the first time, local potentials in proteins can be quantitatively characterized in terms of their strength (rhombicity) evaluated by axial Ag (rhombic Ag and B2u) contributions. Until now, the chain-segment [β3-L3-β4] and to some extent the α2-helix have been associated with GTPase binding. Here, we find that this process causes an increase (decrease) in the potential strength of β3 and β4 (the preceding L2 loop and the remote chain-segment [(α2-helix)-(α2/β5-turn)-(β5-strand)]), suggesting effects of counterbalancing and allostery. There is evidence for the L2 loop being associated with RBD-GTPase binding. Until now only the L4 loop has been associated with RBD dimerization. The latter process is found to cause an increase (decrease) in the potential strength and rhombicity of the L4 loop (the adjacent chain-segment [(α2-helix)-(α2/β5-turn)-(β5-strand)]), suggesting counterbalancing activity. On average, the RBD dimer features stronger local potentials than RBD-Rac1. The novel information inherent in these findings is mesoscopic in character. Prospects of interest include exploring relation to atomistic force-field parameters.

Local Structures in Proteins from Microsecond Molecular Dynamics Simulations: A Symmetry-Based Perspective.

We report on a new method for the characterization of local structures in proteins based on extensive molecular dynamics (MD) simulations, here, 1 μs in length. The N-H bond of the Rho GTPase binding domain of plexin-B1 (RBD) serves as a probe and the potential, u(MD), which restricts its internal motion, as a qualifier of the local dynamic structure. u(MD) is derived from the MD trajectory as a function of the polar angles, (θ, φ), which specify the N-H orientation in the protein. u(MD) is statistical in character yielding empirical descriptions. To establish more insightful methodical descriptions, we develop a comprehensive method which approximates u(MD) by combinations of analytical Wigner functions that belong to the D2h point group. These combinations, called u(simulated), make it possible to gain a new perspective of local dynamic structures in proteins based on explicit potentials/free energy surfaces and associated probability densities, entropy, and ordering. A simpler method was developed previously using 100 ns MD simulations. In that case, the traditional "perpendicular N-H ordering" setting centered at Cα-Cα with (θ, φ) = (90, 90) and generally, featuring positive φ, prevailed. u(MD) derived from 1 μs MD simulations is considerably more complex requiring substantial model enhancement. The enhanced method applies to the well-structured sections of the RBD. It only applies partly to its loops where u(MD) extends into the negative-φ region where we detect nonperpendicular N-H ordering. This arrangement requires devising new reference structures and making substantial algorithmic changes, to be performed in future work. Here, we focus on developing the comprehensive method and using it to investigate perpendicular ordering settings. We find that secondary structures (loops) exhibit varying (virtually invariant) potentials with Ag, B2u, and B1u (Ag and B2u) D2h symmetry. Application to RBD dimerization and RBD binding to the GTPase Rac1 is described in the subsequent article. Applications to other probes, proteins, and biological functions, based on explicit local potentials, probability densities, entropy, and ordering, are possible.

Differential behavior of conformational dynamics in active and inactive states of cannabinoid receptor 1 revealed by microsecond molecular dynamics simulation

Differential behavior of conformational dynamics in active and inactive states of cannabinoid receptor 1 revealed by microsecond molecular dynamics simulation

Complete Virion Simulated: All-Atom Model of an MS2 Bacteriophage with Native Genome.

For the first time, a complete all-atom molecular dynamics (MD) model of a virus, bacteriophage MS2, in its entirety, including a protein outer shell, native genomic RNA with necessary divalent ions, and surrounding explicit aqueous solution with ions at physiological concentration, was built. The model is based on an experimentally measured cryo-EM structure, which was substantially augmented by reconstructing missing or low-resolution parts of the measured density (where the atomistic structure cannot be fit unambiguously). The model was tested by a quarter of a microsecond MD run, and various biophysical characteristics are obtained and analyzed. The developed methodology of building the model can be used for reconstructing other large biomolecular structures when experimental data are fragmented and/or of varying resolution, while the model itself can be used for studying the biology of MS2, including the dynamics of its interaction with the host bacteria.

Structure and dynamics of whole-sequence homology model of ORF3a protein of SARS-CoV-2: An insight from microsecond molecular dynamics simulations

The ORF3a is a large accessory protein in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which plays an important role in virulence and viral replication; especially in inflammasome activation and apoptosis. However,, the existing cryo-EM structure of SARS-CoV-2 ORF3a is incomplete, . making it challenging to understand its structural and functional features. The aim of this study is to investigate the dynamic behaviors of the full-sequence homology model of ORF3a and compare it with the cryo-EM structure using microsecond molecular dynamics simulations. The previous studies indicated that the unresolved residues of the cryo-EM structure are not only involved in the pathogenesis of the SARS-CoV-2 but also exhibit a significant antigenicity. The dynamics scenario of homology model revealed higher RMSD, Rg, and SASA values with stable pattern when compared to the cryo-EM structure. Moreover, the RMSF analysis demonstrated higher fluctuations at specific positions (1-43, 97-110, 172-180, 219-243) in the model structure, whereas the cryo-EM structure displayed lower overall drift (except 1-43) in comparison to the model structure.Secondary structural features indicated that a significant unfolding in the transmembrane domains and β-strand at positions 166 to 172, affecting the stability and compactness of the cryo-EM structure , whereas the model exhibited noticeable unfolding in transmembrane domains and small-coiled regions in the N-terminal. , The results from molecular docking and steered molecular dynamics investigations showed the model structure had a greater number of non-bonding interactions, leading to enhanced stability when compared to the cryo-EM structure. Consequently, higher forces were necessary for unbinding of the baricitinib and ruxolitinib inhibitors from the model structure.. Our findings can help better understanding of the significance of unresolved residues at the molecular level. Additionally, this information can guide researchers for experimental endeavors aimed at completing the full-sequence structure of the ORF3a. Communicated by Ramaswamy H. Sarma

Structural and functional effects of the L84S mutant in the SARS-COV-2 ORF8 dimer based on microsecond molecular dynamics study

The L84S mutation has been observed frequently in the ORF8 protein of SARS-CoV-2, which is an accessory protein involved in various important functions such as virus propagation, pathogenesis, and evading the immune response. However, the specific effects of this mutation on the dimeric structure of ORF8 and its impacts on interactions with host components and immune responses are not well understood. In this study, we performed one microsecond molecular dynamics (MD) simulation and analyzed the dimeric behavior of the L84S and L84A mutants in comparison to the native protein. The MD simulations revealed that both mutations caused changes in the conformation of the ORF8 dimer, influenced protein folding mechanisms, and affected the overall structural stability. In particular, the 73YIDI76 motif has found to be significantly affected by the L84S mutation, leading to structural flexibility in the region connecting the C-terminal β4 and β5 strands. This flexibility might be responsible for virus immune modulation. The free energy landscape (FEL) and principle component analysis (PCA) have also supported our investigation. Overall, the L84S and L84A mutations affect the ORF8 dimeric interfaces by reducing the frequency of protein–protein interacting residues (Arg52, Lys53, Arg98, Ile104, Arg115, Val117, Asp119, Phe120, and Ile121) in the ORF8 dimer. Our findings provide detail insights for further research in designing structure-based therapeutics against the SARS-CoV-2. Communicated by Ramaswamy H. Sarma

IDH3γ functions as a redox switch regulating mitochondrial energy metabolism and contractility in the heart

Redox signaling and cardiac function are tightly linked. However, it is largely unknown which protein targets are affected by hydrogen peroxide (H2O2) in cardiomyocytes that underly impaired inotropic effects during oxidative stress. Here, we combine a chemogenetic mouse model (HyPer-DAO mice) and a redox-proteomics approach to identify redox sensitive proteins. Using the HyPer-DAO mice, we demonstrate that increased endogenous production of H2O2 in cardiomyocytes leads to a reversible impairment of cardiac contractility in vivo. Notably, we identify the γ-subunit of the TCA cycle enzyme isocitrate dehydrogenase (IDH)3 as a redox switch, linking its modification to altered mitochondrial metabolism. Using microsecond molecular dynamics simulations and experiments using cysteine-gene-edited cells reveal that IDH3γ Cys148 and 284 are critically involved in the H2O2-dependent regulation of IDH3 activity. Our findings provide an unexpected mechanism by which mitochondrial metabolism can be modulated through redox signaling processes.

Allosteric inhibition of TEM-1 β lactamase: Microsecond molecular dynamics simulations provide mechanistic insights.

β-lactam antibiotics target DD-transpeptidases, enzymes that perform the last step of bacterial cell-wall synthesis. To block the antimicrobial activity of these antibiotics, bacteria have evolved lactamases that render them inert. Among these, TEM-1, a class A lactamase, has been extensively studied. In 2004, Horn et al. described a novel allosteric TEM-1 inhibitor, FTA, that binds distant from the TEM-1 orthosteric (penicillin-binding) pocket. TEM-1 has subsequently become a model for the study of allostery. In the present work, we perform molecular dynamics simulations of FTA-bound and FTA-absent TEM-1, totaling ~3μS, that provide new insight into TEM-1 inhibition. In one of the simulations, bound FTA assumed a conformation different than that observed crystallographically. We provide evidence that the alternate pose is physiologically plausible and describe how it impacts our understanding of TEM-1 allostery.

Adenine Fine-Tunes DNA Photolyase's Repair Mechanism.

The comparative study of DNA repair by mesophilic and extremophilic photolyases helps us understand the evolution of these enzymes and their role in preserving life on our changing planet. The mechanism of repair of cyclobutane pyrimidine dimer lesions in DNA by electron transfer from the flavin adenine dinucleotide cofactor is the subject of intense interest. The role of adenine in mediating this process remains unresolved. Using microsecond molecular dynamics simulations, we find that adenine mediates the electron transfer in both mesophile and extremophile DNA photolyases through a similar mechanism. In fact, in all photolyases studied, the molecular conformations with the largest electronic couplings between the enzyme cofactor and DNA show the presence of adenine in 10-20% of the strongest-coupling tunneling pathways between the atoms of the electron donor and acceptor. Our theoretical analysis finds that adenine serves the critical role of fine-tuning rather than maximizing the donor-acceptor coupling within the range appropriate for the repair function.

SARS-CoV-2 envelope protein attain Kac mediated dynamical interaction network to adopt ‘histone mimic’ at BRD4 interface

Interface mimicry, achieved by recognition of host-pathogen interactions, is the basis by which pathogen proteins can hijack the host machinery. The envelope (E) protein of SARS-CoV-2 is reported to mimic the histones at the BRD4 surface via establishing the structural mimicry; however, the underlying mechanism of E protein mimicking the histones is still elusive. To explore the mimics at dynamic and structural residual network level an extensive docking, and MD simulations were carried out in a comparative manner between complexes of H3-, H4-, E-, and apo-BRD4. We identified that E peptide is able to attain an ‘interaction network mimicry’, as its acetylated lysine (Kac) achieves orientation and residual fingerprint similar to histones, including water-mediated interactions for both the Kac positions. We identified Y59 of E, playing an anchor role to escort lysine positioning inside the binding site. Furthermore, the binding site analysis confirms that E peptide needs a higher volume, similar to the H4-BRD4 where both the lysine’s (Kac5 and Kac8) can accommodate nicely, however, the position of Kac8 is mimicked by two additional water molecules other than four water-mediated bridging’s, strengthening the possibility that E peptide could hijack host BRD4 surface. These molecular insights seem pivotal for mechanistic understanding and BRD4-specific therapeutic intervention. KEY POINTS Molecular mimicry is reported in hijacking and then outcompeting the host counterparts so that pathogens can rewire their cellular function by overcoming the host defense mechanism. The molecular recognition process is the basis of molecular mimicry. The E peptide of SARS-CoV-2 is reported to mimic host histone at the BRD4 surface by utilizing its C-terminally placed acetylated lysine (Kac63) to mimic the N-terminally placed acetylated lysine Kac5GGKac8 histone (H4) by interaction network mimicry identified through microsecond molecular dynamics (MD) simulations and post-processing extensive analysis. There are two steps to mimic: firstly, tyrosine residues help E to anchor at the BRD4 surface to position Kac and increase the volume of the pocket. Secondary, after positioning of Kac, a common durable interaction network N140:Kac5; Kac5:W1; W1:Y97; W1:W2; W2:W3; W3:W4; W4:P82 is established between Kac5, with key residues P82, Y97, N140, and four water molecules through water mediate bridge. Furthermore, the second acetylated lysine Kac8 position and its interaction as polar contact with Kac5 were also mimicked by E peptide through interaction network P82:W5; W5:Kac63; W5:W6; W6:Kac63. The binding event at BRD4/BD1 seems an induced-fit mechanism as a bigger binding site volume was identified at H4-BRD4 on which E peptide attains its better stability than H3-BRD4. We identified the tyrosine residue Y59 of E that acts like an anchor on the BRD4 surface to position Kac inside the pocket and attain the interaction network by using aromatic residues of the BRD4 surface. Communicated by Ramaswamy H. Sarma

Fibroblast growth factor 5 (FGF5) and its missense mutant FGF5-H174 underlying trichomegaly: a molecular dynamics simulation investigation

The missense mutation Y174H of FGF5 (FGF5-H174) had been associated with trichomegaly, characterized by abnormally long and pigmented eyelashes. The amino acid tyrosine (Tyr/Y) at position 174 is conserved across many species, proposedly holding important characteristics for the functions of FGF5. Microsecond molecular dynamics simulations along with protein-protein docking and residue interacting network analysis were employed to investigate the structural dynamics and binding mode of both wild-type (FGF5-WT) and its mutated counterpart (FGF5-H174). It was found that the mutation decreased number of hydrogen bonds within the protein, sheet secondary structure, interaction of residue 174 with other residues, and number of salt-bridges. On the other hand, the mutation increased solvent accessible surface area, number of hydrogen bonds between the protein and solvent, coil secondary structure, protein C-alpha backbone root mean square deviation, protein residue root mean square fluctuations, as well as occupied conformational space. In addition, protein-protein docking integrated with molecular dynamics simulations and molecular mechanics - Poisson-Boltzmann surface area (MM/PBSA) binding energy calculation demonstrated that the mutated variant possessed stronger binding affinity towards fibroblast growth factor receptor 1 (FGFR1). However, residue interaction network analysis demonstrated that the binding mode of the FGFR1-FGF5-H174 complex was drastically different from that of the FGFR1-FGF5-WT complex. In conclusion, the missense mutation conferred more instability within itself and stronger binding affinity towards FGFR1 with distinctively altered binding mode or residue connectivity. These findings might help explain the decreased pharmacological activity of FGF5-H174 towards FGFR1, underlying trichomegaly. Communicated by Ramaswamy H. Sarma