GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation

Identification of potent BRD4-BD1 inhibitors using classical and steered molecular dynamics based free energy analysis

In the present work a combination of traditional and steered molecular dynamics based techniques were employed to identify potential inhibitors against the human BRD4 protein (BRD4- BD1); an established drug target for multiple illnesses including various malignancies. Quinoline derivatives that were synthesized in-house were tested for their potential as new BRD4-BD1 inhibitors. Initially molecular docking experiments were performed to determine the binding poses of BRD4-BD1 inhibitors. To learn more about the thermodynamics of inhibitor binding to the BRD4-BD1 active site, the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) free energy calculations were conducted afterwards. The findings of the MM-PBSA analysis were further reinforced by performing steered umbrella sampling simulations which revealed crucial details about the binding/unbinding process of the most potent quinoline derivatives at the BRD4-BD1 active site. We report a novel quinoline derivative which can be developed into a fully functional BRD4-BD1 inhibitor after experimental validation. The identified compound (4 g) shows better properties than the standard BRD4-BD1 inhibitors considered in the study. The study also highlights the crucial role of Gln78, Phe79, Trp81, Pro82, Phe83, Gln84, Gln85, Val87, Leu92, Leu94, Tyr97, Met105, Cys136, Asn140, Ile146 and Met149 in inhibitor binding. The study provides a possible lead candidate and key amino acids involved in inhibitor recognition and binding at the active site of BRD4-BD1 protein. The findings might be of significance to medicinal chemists involved in the development of potent BRD4-BD1 inhibitors.

Monitoring the interactions between POPG phospholipid bilayer and amyloid-forming protein human cystatin C. Does the bilayer influence the oligomeric state and structure of the protein?

A biological membrane is a structure characteristic for various cells and organelles present in almost all living organisms. Even though, it is one of the most common structures in organisms, where it serves crucial functions, a phospholipid bilayer may also take part in pathological processes leading to severe diseases. Research indicates that biological membranes have a profound impact on the pathological processes of oligomerization of amyloid-forming proteins. These processes are a hallmark of amyloid diseases, a group of pathological states involving, e.g., Parkinson's or Alzheimer's disease. Even though amyloidogenic diseases reap the harvest in modern societies, especially in elderly patients, the mechanisms governing the amyloid deposition are not clearly described. Therefore, the presented study focuses on the description of interactions between a model biological membrane (POPG) and one of amyloid forming proteins - human cystatin C. For the purpose of the study molecular dynamics simulations were applied to confirm interactions between the protein and POPG membrane. Next the NMR techniques were used to verify how the data obtained in solution compared to MD simulations and determine fragments of the protein responsible for interactions with POPG. Finally, circular dichroism was used to monitor the changes in secondary structure of the protein and size exclusion chromatography was used to monitor its oligomerization process. Obtained data indicates that the protein interacts with POPG submerging itself into the bilayer with the AS region. However, the presence of POPG bilayer does not significantly affect the structure or oligomerization process of human cystatin C.

A simulation-based approach to target Zika virus RNA-dependent RNA polymerase with marine compounds for antiviral development

Despite significant efforts, currently, there is no particular drug available to treat Zika virus (ZIKV) infection, highlighting the urgent need for effective therapeutic interventions. To identify putative inhibitors of the ZIKV RdRp protein's RNA binding function, the present study applied an extensive in-silico drug discovery methodology. The initial phase involved virtual screening using Lipinski's rule of five as a filter, ensuring the selection of molecules with favorable pharmacokinetic properties. This process yielded 238 compounds with promising docking scores, ranging from -6.0 to -7.48 kcal/mol, indicative of their potential binding affinity to the ZIKV RdRp. To refine the selection, these compounds underwent a re-docking process, comparing their binding energies with a reference molecule known for its inhibitory action against RdRp. Remarkably, five compounds, labeled CMNPD30598, CMNPD27464, CMNPD25971, CMNPD27444, and CMNPD16599, demonstrated superior re-docking energies compared to the reference, suggesting a stronger interaction with the RdRp allosteric site. Subsequent molecular dynamics (MD) simulations provided insights into the stability of these complexes over time, reinforcing their potential as RdRp inhibitors. Additionally, the calculation of free binding energies and principal component analysis (PCA) of the free energy landscape offered a deeper understanding of the binding dynamics and energetics. This study not only highlights the utility of marine fungi compounds in antiviral drug discovery but also showcases the power of computational tools in identifying novel therapeutics. The identified compounds represent promising candidates for further experimental validation and development as ZIKV RdRp inhibitors.Communicated by Ramaswamy H. Sarma.

Hydrogen bonding in glassy trehalose-water system: Insights from density functional theory and molecular dynamics simulations

We report a detailed density functional theory and molecular dynamics study of hydrogen bonding between trehalose and water, with a special emphasis on interactions in the amorphous solid state. For comparison, water-water interactions in water dimers and tetramers are evaluated using quantum calculations. The results show that the hydrogen bonding energy is dependent not only on the geometry (bond length and angle) but also on the local environment of the hydrogen bond. This is seen in quantum calculations of complexes in vacuum as well as in amorphous solid states with periodic boundary conditions. The temperature-induced glass transition in the trehalose-water system was studied using molecular dynamics simulations with varying cooling and heating rates. The obtained parameters of the glass transition are in good agreement with the experiments. Moreover, the dehydration of trehalose in the glassy state was investigated through a gradual dehydration with multiple small steps under isothermal conditions. From these simulations, the values of water sorption energy at different temperatures were obtained. The partial molar enthalpy of mixing of water value of -18 kJ/mol found in calorimetric experiments was accurately reproduced in these simulations. These findings are discussed in light of the hydrogen bonding data in the system. We conclude that the observed exothermic effect is due to different responses of liquid and glassy matrices to perturbations associated with the addition or removal of water molecules.

Understanding a point mutation signature D54K in the caspase activation recruitment domain of NOD1 capitulating concerted immunity via atomistic simulation

Point mutation D54K in the human N-terminal caspase recruitment domain (CARD) of nucleotide-binding oligomerization domain -1 (NOD1) abrogates an imperative downstream interaction with receptor-interacting protein kinase (RIPK2) that entails combating bacterial infections and inflammatory dysfunction. Here, we addressed the molecular details concerning conformational changes and interaction patterns (monomeric-dimeric states) of D54K by signature-based molecular dynamics simulation. Initially, the sequence analysis prioritized D54K as a pathogenic mutation, among other variants, based on a sequence signature. Since the mutation is highly conserved, we derived the distant ortholog to predict the sequence and structural similarity between native and mutant. This analysis showed the utility of 33 communal core residues associated with structural-functional preservation and variations, concurrently served to infer the cryptic hotspots Cys39, Glu53, Asp54, Glu56, Ile57, Leu74, and Lys78 determining the inter helical fold forming homodimers for putative receptor interaction. Subsequently, the atomistic simulations with free energy (MM/PB(GB)SA) calculations predicted structural alteration that takes place in the N-terminal mutant CARD where coils changed to helices (45 α3- L4-α4-L6- α683) in contrast to native (45T2-L4-α4-L6-T483). Likewise, the C-terminal helices 93T1-α7105 connected to the loops distorted compared to native 93α6-L7105 may result in conformational misfolding that promotes functional regulation and activation. These structural perturbations of D54K possibly destabilize the flexible adaptation of critical homotypic NOD1CARD-CARDRIPK2 interactions (α4Asp42-Arg488α5 and α6Phe86-Lys471α4) is consistent with earlier experimental reports. Altogether, our findings unveil the conformational plasticity of mutation-dependent immunomodulatory response and may aid in functional validation exploring clinical investigation on CARD-regulated immunotherapies to prevent systemic infection and inflammation.Communicated by Ramaswamy H. Sarma.

Molecular Simulation of SO2 Separation and Storage Using a Cryptophane-Based Porous Liquid

A theoretical molecular simulation study of the encapsulation of gaseous SO2 at different temperature conditions in a type II porous liquid is presented here. The system is composed of cage cryptophane-111 molecules that are dispersed in dichloromethane, and it is described using an atomistic modelling of molecular dynamics. Gaseous SO2 tended to almost fully occupy cryptophane-111 cavities throughout the simulation. Calculations were performed at 300 K and 283 K, and some insights into the different adsorption found in each case were obtained. Simulations with different system sizes were also studied. An experimental-like approach was also employed by inserting a SO2 bubble in the simulation box. Finally, an evaluation of the radial distribution function of cryptophane-111 and gaseous SO2 was also performed. From the results obtained, the feasibility of a renewable separation and storage method for SO2 using porous liquids is mentioned.

Optimizing the Martini 3 Force Field Reveals the Effects of the Intricate Balance between Protein-Water Interaction Strength and Salt Concentration on Biomolecular Condensate Formation

Condensation/dissolution has become a widely acknowledged biological macromolecular assembly phenomenon in subcellular compartmentalization. The MARTINI force field offers a coarse-grained protein model with a resolution that preserves molecular details with an explicit (CG) solvent. Despite its relatively higher resolution, it can still achieve condensate formation in a reasonable computing time with explicit solvent and ionic species. Therefore, it is highly desirable to tune this force field to be able to reproduce the experimentally observed properties of the condensate formation. In this work, we studied the condensate formation of the low-sequence complexity domain of fused in sarcoma protein using a MARTINI 3 force field by systematically modifying (increasing) the protein-water interaction strength and varying the salt concentration. We found that the condensate formation is sensitive both to the protein-water interaction strength and the presence of salt. While the unmodified MARTINI force field yields a complete collapse of proteins into one dense phase (i.e., no dilute phase), we reported a range of modified protein-water interaction strength that is capable of capturing the experimentally found transfer free energy between dense and dilute phases. We also found that the condensates lose their spherical shape upon the addition of salt, especially when the protein-water interactions are weak. Interchain amino acid contact map analysis showed one explanation for this observation: the protein-protein contact fraction reduces as salt is added to systems (when the protein-water interactions are weak), consistent with electrostatic screening effects. This reduction might be responsible for the condensates becoming nonspherical upon the addition of salt by reducing the need for minimizing the interfacial area. However, as the protein-water interactions become stronger to the extent that makes the transfer free energy agree well with experimentally observed transfer free energy, we found an increase in the protein-protein contact fraction upon the addition of salt, consistent with the salting-out effects. Therefore, we concluded that there is an intricate balance between screening effects and salting-out effects upon the addition of salt and this balance is highly sensitive to the strength of protein-water interactions.

Histone lactylation inhibits RARγ expression in macrophages to promote colorectal tumorigenesis through activation of TRAF6-IL-6-STAT3 signaling

Macrophages are phenotypically and functionally diverse in the tumor microenvironment (TME). However, how to remodel macrophages with a protumor phenotype and how to manipulate them for therapeutic purposes remain to be explored. Here, we show that in the TME, RARγ is downregulated in macrophages, and its expression correlates with poor prognosis in patients with colorectal cancer (CRC). In macrophages, RARγ interacts with tumor necrosis factor receptor-associated factor 6 (TRAF6), which prevents TRAF6 oligomerization and autoubiquitination, leading to inhibition of nuclear factor κB signaling. However, tumor-derived lactate fuels H3K18 lactylation to prohibit RARγ gene transcription in macrophages, consequently enhancing interleukin-6 (IL-6) levels in the TME and endowing macrophages with tumor-promoting functions via activation of signal transducer and activator of transcription 3 (STAT3) signaling in CRC cells. We identified that nordihydroguaiaretic acid (NDGA) exerts effective antitumor action by directly binding to RARγ to inhibit TRAF6-IL-6-STAT3 signaling. This study unravels lactate-driven macrophage function remodeling by inhibition of RARγ expression and highlights NDGA as a candidate compound for treating CRC.

Structural and energetic analysis of NS5 protein inhibition by small molecules in Japanese encephalitis virus using machine learning and steered molecular dynamics approach

One of the viral diseases that affect millions of people around the world, particularly in developing countries, is Japanese encephalitis (JE). In this study, the conserved protein of this virus, that is, non-structural protein 5 (NS5), was used as a target protein for this study, and a compound library of 749 antiviral molecules was screened against NS5. The current study employed machine learning-based virtual screening combined with molecular docking. Here, three hits (24360, 123519051 and 213039) had lower binding energies (< -8 kcal/mol) than the control, S-Adenosyl-L-homocysteine (SAH). All the compounds showed significant H-bond interactions with functional residues, which were also observed by the control. Molecular dynamics simulation, MM/GBSA for binding free energy analysis, principal component analysis and free energy landscape were also performed to study the stability of the complex formation. All three compounds had similar root mean square deviation trends, which were comparable to the control, SAH. Post-MD, the 123519051-receptor complex had the highest number of H-bonds (4 to 5) after the control, out of which three exhibited the highest percentage occupancy (50%, 24% and 79%). Both docking and MD, 123519051 showed an H-bond with the residue Gly111, which was also found for the control-protein complex. 123519051 showed the lowest binding free energy with ΔGbind of -89 kJ/mol. Steered molecular dynamics depicted that 123519051 had the maximum magnitude of dissociation (1436.43 kJ/mol/nm), which was more than the control, validating its stable complex formation. This study concluded that 123519051 is a binder and could inhibit the protein NS5 of JE.Communicated by Ramaswamy H. Sarma.

Solvent-Dependent Folding Behavior of a Helix-Forming Aryl-Triazole Foldamer

Aromatic foldamers make up a novel class of bioinspired molecules that display helical conformations and have functions that rely on control over their coil-helix folding preferences. While the folding has been extensively examined by experiment, it has rarely been paired with the types of atomic level insights offered by theory. We present the results of all-atom molecular dynamics (MD) simulations to examine the role of solvent polarity on driving the helical folding behavior of the aryl-triazole foldamer. The temperature-dependent enhanced sampling technique, replica-exchange MD simulations, was employed to understand the folding phenomena. The simulation results show that in a low polarity solvent (dichloromethane), the foldamer prefers to stay in the unfolded state. The unfolded state has four dipolar triazoles (5 D) in their favored all-anti geometries and favoring anti-parallel geometries. However, an increase in solvent polarity using acetonitrile resulted in solvophobic collapse, yielding the helically folded form as the predominant state. The folded helix has an all-syn geometry, with triazoles in parallel arrangements. Intermediate conformations with a mixture of syn and anti arrangements of the triazoles are of lower abundance in both the DCM and MeCN solvents. The chiral handedness of the helix observed experimentally is assigned as left-handed by correlation with computed electronic circular dichroism spectra using time-dependent density functional theory.

Addressing Suboptimal Poses in Nonequilibrium Alchemical Calculations

Alchemical transformations can be used to quantitatively estimate absolute binding free energies at a reasonable computational cost. However, most of the approaches currently in use require knowledge of the correct (crystallographic) pose. In this paper, we present a combined Hamiltonian replica exchange nonequilibrium alchemical method that allows us to reliably calculate absolute binding free energies, even when starting from suboptimal initial binding poses. Performing a preliminary Hamiltonian replica exchange enhances the sampling of slow degrees of freedom of the ligand and the target, allowing the system to populate the correct binding pose when starting from an approximate docking pose. We apply the method on 6 ligands of the first bromodomain of the BRD4 bromodomain-containing protein. For each ligand, we start nonequilibrium alchemical transformations from both the crystallographic pose and the top-scoring docked pose that are often significantly different. We show that the method produces statistically equivalent binding free energies, making it a useful tool for computational drug discovery pipelines.

Allosteric regulation in SARS-CoV-2 spike protein

Allosteric regulation is common in protein-protein interactions and is thus promising in drug design. Previous experimental and simulation work supported the presence of allosteric regulation in the SARS-CoV-2 spike protein. Here the route of allosteric regulation in SARS-CoV-2 spike protein is examined by all-atom explicit solvent molecular dynamics simulations, contrastive machine learning, and the Ohm approach. It was found that peptide binding to the polybasic cleavage sites, especially the one at the first subunit of the trimeric spike protein, activates the fluctuation of the spike protein's backbone, which eventually propagates to the receptor-binding domain on the third subunit that binds to ACE2. Remarkably, the allosteric regulation routes starting from the polybasic cleavage sites share a high fraction (39-67%) of the critical amino acids with the routes starting from the nitrogen-terminal domains, suggesting the presence of an allosteric regulation network in the spike protein. Our study paves the way for the rational design of allosteric antibody inhibitors.

Molecular Insights into the Dynamics of Amyloid Fibril Growth: Elongation and Lateral Assembly of GNNQQNY Protofibrils

The self-assembly of peptides and proteins into β-sheet rich amyloid fibrils is linked to both functional and pathological states. In this study, the growth of fibrillar structures of the short peptide GNNQQNY, a fragment from the yeast prion Sup35 protein, was examined. Molecular dynamics simulations were used to study alternative mechanisms of fibril growth, including elongation through binding of monomers as well as fibril self-assembly into larger, more mature structures. It was found that after binding, monomers diffused along preformed fibrils toward the ends, supporting the mechanism of fibril growth via elongation. Lateral assembly of protofibrils was found to occur readily, suggesting that this could be the key to transitioning from isolated fibrils to mature multilayer structures. Overall, the work provides mechanistic insights into the competitive pathways that govern amyloid fibril growth.

Computing the frequency-dependent NMR relaxation of 1H nuclei in liquid water

We present a computational framework for reliably determining the frequency-dependent intermolecular and intramolecular nuclear magnetic resonance (NMR) dipole-dipole relaxation rates of spin 1/2 nuclei from Molecular Dynamics (MD) simulations. This approach avoids the alterations caused by the well-known finite-size effects of translational diffusion. Moreover, a procedure is derived to control and correct for effects caused by fixed distance-sampling cutoffs and periodic boundary conditions. By construction, this approach is capable of accurately predicting the correct low-frequency scaling behavior of the intermolecular NMR dipole-dipole relaxation rate and thus allows for the reliable calculation of the frequency-dependent relaxation rate over many orders of magnitude. Our approach is based on the utilization of the theory of Hwang and Freed for the intermolecular dipole-dipole correlation function and its corresponding spectral density [L.-P. Hwang and J. H. Freed, J. Chem. Phys. 63, 4017-4025 (1975)] and its combination with data from MD simulations. The deviations from the Hwang and Freed theory caused by periodic boundary conditions and sampling distance cutoffs are quantified by means of random walker Monte Carlo simulations. An expression based on the Hwang and Freed theory is also suggested for correcting those effects. As a proof of principle, our approach is demonstrated by computing the frequency-dependent intermolecular and intramolecular dipolar NMR relaxation rates of 1H nuclei in liquid water at 273 and 298 K based on the simulations of the TIP4P/2005 model. Our calculations are suggesting that the intermolecular contribution to the 1H NMR relaxation rate of the TIP4P/2005 model in the extreme narrowing limit has previously been substantially underestimated.

Targeting RdRp of SARS-CoV-2 with De Novo Molecule Generation

Viruses are known for their extremely high mutation rates, allowing them to evade both the human immune system and many forms of standard medicine. Despite this, the RNA dependent RNA polymerase (RdRp) of the RNA viruses has been largely conserved, and any significant mutation of this protein is unlikely. The recent COVID-19 pandemic presents a need for therapeutics. We have designed a de novo drug design algorithm that generates strong binding ligands from scratch, based on only the structure of the target protein's receptor. In this paper, we applied our method to target SARS-CoV-2 RdRp and generated several de novo molecules. We then chose some drug molecules based on the structural similarity to some of our strongest binding de novo molecules. Subsequently, we showed, using rigorous all-atom explicit-water free energy calculations in near-microsecond time scales using state-of-the-art well-tempered metadynamics simulations, that some of our de novo generated ligands bind more strongly to RdRp than the recent FDA approved drug remdesivir in its active form, remdesivir triphosphate (RTP). We elucidated the binding mechanism for some of the top binders and compared it with RTP. We believe that this work will be useful both by presenting lead structures for RdRp inhibition and by delivering key insights into the residues of the protein potentially involved in the binding/unbinding of these small molecule drugs, leading to more targeted studies in the future.

Marine Brown Algae-Derived Compounds as Potential Inhibitors of Japanese Encephalitis Virus RNA-Dependent RNA Polymerase

The Japanese encephalitis virus (JEV) is a mosquito-borne flavivirus that primarily affects people in Asia and seriously threatens public health. Considering the rising occurrence rates and lack of targeted antiviral treatments, it is essential to comprehend and tackle obstacles related to JEV in order to lessen its influence on world health. This investigation explores compounds derived from marine brown algae (Phaeophyceae) as potential inhibitors of JEV RNA-dependent RNA polymerase (RdRp), a critical enzyme in the virus's replication cycle. Employing the computational virtual screen approach, four compounds, i.e., CMNPD16749, CMNPD2606, CMNPD27817, and CMNPD23662, with favorable binding energies ranging from -15.7 Kcal/mol to -13.9 kcal/mol were identified. Subsequently, through molecular docking analysis, the interactions responsible for the binding stability between the target protein and hit molecules compared to the reference molecule Galidesvir were studied. Further, through extensive molecular dynamic (MD) simulation studies at 200 ns, it was confirmed that each docked complex showed acceptable dynamic stability compared to the reference molecule. These findings were further validated using MM/PBSA free binding energy calculations, PCA analysis and free energy landscape construction. These computational findings suggested that the brown algae-derived compounds may act as an antiviral drug against JEV infection and lay a crucial foundation for future experimental studies against JEV.

Unveiling the State Transition Mechanisms of Ras Proteins through Enhanced Sampling and QM/MM Simulations

In cells, wild-type RasGTP complexes exist in two distinct states: active State 2 and inactive State 1. These complexes regulate their functions by transitioning between the two states. However, the mechanisms underlying this state transition have not been clearly elucidated. To address this, we conducted a detailed simulation study to characterize the energetics of the stable states involved in the state transitions of the HRasGTP complex, specifically from State 2 to State 1. This was achieved by employing multiscale quantum mechanics/molecular mechanics and enhanced sampling molecular dynamics methods. Based on the simulation results, we constructed the two-dimensional free energy landscapes that provide crucial information about the conformational changes of the HRasGTP complex from State 2 to State 1. Furthermore, we also explored the conformational changes from the intermediate state to the product state during guanosine triphosphate hydrolysis. This study on the conformational changes involved in the HRas state transitions serves as a valuable reference for understanding the corresponding events of both KRas and NRas as well.

Therapeutic Peptides Are Preferentially Solubilized in Specific Microenvironments within PEG-PLGA Polymer Nanoparticles

Polymeric nanoparticles are a highly promising drug delivery formulation. However, a lack of understanding of the molecular mechanisms that underlie their drug solubilization and controlled release capabilities has hindered the efficient clinical translation of such technologies. Polyethylene glycol-poly(lactic-co-glycolic) acid (PEG-PLGA) nanoparticles have been widely studied as cancer drug delivery vehicles. In this letter, we use unbiased coarse-grained molecular dynamics simulations to model the self-assembly of a PEG-PLGA nanoparticle and its solubulization of the anticancer peptide, EEK, with good agreement with previously reported experimental structural data. We applied unsupervised machine learning techniques to quantify the conformations that polymers adopt at various locations within the nanoparticle. We find that the local microenvironments formed by the various polymer conformations promote preferential EEK solubilization within specific regions of the NP. This demonstrates that these microenvironments are key in controlling drug storage locations within nanoparticles, supporting the rational design of nanoparticles for therapeutic applications.

Multiscale modelling of biomolecular corona formation on metallic surfaces

In the realm of food industry, the choice of non-consumable materials used plays a crucial role in ensuring consumer safety and product quality. Aluminum is widely used in food packaging and food processing applications, including dairy products. However, the interaction between aluminum and milk content requires further investigation to understand its implications. In this work, we present the results of multiscale modelling of the interaction between various surfaces, that is (100), (110), and (111), of fcc aluminum with the most abundant milk proteins and lactose. Our approach combines atomistic molecular dynamics, a coarse-grained model of protein adsorption, and kinetic Monte Carlo simulations to predict the protein corona composition in the deposited milk layer on aluminum surfaces. We consider a simplified model of milk, which is composed of the six most abundant milk proteins found in natural cow milk and lactose, which is the most abundant sugar found in dairy. Through our study, we ranked selected proteins and lactose adsorption affinities based on their corresponding interaction strength with aluminum surfaces and predicted the content of the naturally forming biomolecular corona. Our comprehensive investigation sheds light on the implications of aluminum in food processing and packaging, particularly concerning its interaction with the most abundant milk proteins and lactose. By employing a multiscale modelling approach, we simulated the interaction between metallic aluminum surfaces and the proteins and lactose, considering different crystallographic orientations. The results of our study provide valuable insights into the mechanisms of lactose and protein deposition on aluminum surfaces, which can aid in the general understanding of protein corona formation.

In Silico Screening of Chlorogenic Acids from Plant Sources against Human Translocase-I to Identify Competitive Inhibitors to Treat Diabetes

Chlorogenic acids (CHLs) are known to competitively bind to translocase-I (T1) of the glucose-6-phosphatase (G6 Pase) system, thereby inhibiting the transport of glucose-6-phosphate (G6P). This competitive binding results in a consequential reduction in blood sugar levels. In this study, steered molecular dynamics (SMD) simulation is employed to investigate the interaction between T1 and G6P, aiming to gain insights into the binding dynamics and diffusion process of G6P through T1. A database comprising 41 CHLs sourced from various plants was developed, subjected to minimization, and screened against T1 through conventional docking methods. The docked conformations were fed into a newly developed customized scoring method incorporating contact-based weights to assess the binding affinities that systematically rank and identify the most effective competitive inhibitors. Among the screened CHLs, 1-methoxy 3,5-dicaffeoylquinic acid, 3,4 dicaffeoyl quinic acid, and 3,4,5-tricaffeoylquinic acid stood out as the top three inhibitors, showcasing crucial atomic interactions with key residues within the binding pocket of T1, and these CHLs are sourced from readily available plants, diminishing reliance on coffee as the predominant CHL source. Along with the devised scoring function, which serves as a valuable tool for virtual screening and lead optimization in drug development, this study also marks a pioneering effort as it involves the modeling of the human translocase and unravels the mechanism of binding and diffusion of G6P within human T1, providing valuable insights into the structural prerequisites for successfully inhibiting the G6P system, laying the foundation for a rational approach to drug design. This research contributes to the progress of drug discovery strategies focused on the G6P system, presenting potential therapeutic avenues for addressing metabolic disorders linked to an impaired glucose metabolism.

Exploring the Role of Globular Domain Locations on an Intrinsically Disordered Region of p53: A Molecular Dynamics Investigation

The pre-tetramerization loop (PTL) of the human tumor suppressor protein p53 is an intrinsically disordered region (IDR) necessary for the tetramerization process, and its flexibility contributes to the essential conformational changes needed. Although the IDR can be accurately simulated in the traditional manner of molecular dynamics (MD) with the end-to-end distance (EEdist) unhindered, we sought to explore the effects of restraining the EEdist to the values predicted by electron microscopy (EM) and other distances. Simulating the PTL trajectory with a restrained EEdist , we found an increased agreement of nuclear magnetic resonance (NMR) chemical shifts with experiments. Additionally, we observed a plethora of secondary structures and contacts that only appear when the trajectory is restrained. Our findings expand the understanding of the tetramerization of p53 and provide insight into how mutations could make the protein impotent. In particular, our findings demonstrate the importance of restraining the EEdist in studying IDRs and how their conformations change under different conditions. Our results provide a better understanding of the PTL and the conformational dynamics of IDRs in general, which are useful for further studies regarding mutations and their effects on the activity of p53.

Characterization and transformation of TtMYB1 transcription factor from Tritipyrum to improve salt tolerance in wheat

BACKGROUND

Common wheat (Triticum aestivum L.) is a worldwide cereal crop, which is an integral part of the diets of many countries. In addition, the MYB gene of wheat plays a role in the response to salt stress.

RESULTS

"Y1805" is a Tritipyrum variety that is relatively tolerant to salt. We used transcriptome analysis to show that the "Y1805" MYB gene was both highly expressed and sensitive to salt stress. Compared with control roots, the level of MYB expression during salt stress was higher, which rapidly decreased to control levels during the recovery process. MYB gene relative expression showed the highest levels in "Y1805" roots during salt stress, with the stems and then leaves being the next highest stressed tissues. The novel MYB gene (TtMYB1) was successfully cloned from "Y1805". It showed a coding sequence length of 783 bp with 95.79% homology with Tel2E01G633100 from Thinopyrum elongatum. TtMYB1 and MYB from Th. elongatum were clustered in the same branch using phylogenetic analysis, which indicated high similarities. The TtMYB1 gene is located in the nucleus. The coleoptile method was employed when a TtMYB1 overexpression vector was used during transformation into "1718" (common wheat). Under high salt stress, TtMYB1 leaves of overexpression lines had decreased wilting, when compared with wild-type (WT) plants. During normal conditions, salt stress, and recovery, the lengths of the roots and the heights of seedlings from the overexpression lines were found to be significantly greater than roots and seedlings of WT plants. In addition, during high salt stress, the overexpression lines showed that proline and soluble sugar levels were higher than that of WT plants, but with lower malondialdehyde levels. Forty-three proteins that interacted with TtMYB1 were identified using the yeast two-hybrid assay. Protein-protein interaction analyses indicated that most were SANT domain-containing and Wd repeat region domain-containing proteins. Among these proteins, ribosomal proteins were the main node. Abiotic stress-related terms (such as "carbonate dehydratase activity", "protein targeting peroxisomes", and "glutathione peroxidase activity") were enriched in GO analysis. In KEGG analysis, "carbohydrate metabolism", "environmental information processing", "genetic information processing", "signaling and cell precursors", and "energy metabolism" pathways were enriched.

CONCLUSION

The TtMYB1 gene might enhance salt tolerance by increasing proline and soluble sugar content and antioxidase activity in transgenic wheat. It therefore has the potential to enhance high salt tolerance in plants.

Chain-Level Analysis of Reinforced Polyethylene through Stretch-Induced Crystallization

Herein, we propose a large-scale simulation approach to perform the stretch-induced crystallization of entangled polyethylene (PE) melts. Sufficiently long (1000 ns) united-atom molecular dynamics (UAMD) simulations for 16000 chains of 1000 consecutive CH2 united-atom particles under periodic boundary conditions were performed to achieve the crystallinity observed in experiments. Before the isothermal crystallization process, we applied uniaxial stretching as pre-elongation to the embedded strain memory on the entangled PE melts. We confirmed significant differences in the morphologies of crystal domains and scattering patterns for pre-elongation ratios of 400% and 800%. The obtained scattering patterns were consistent with the experimental results. Uniaxial stretching MD simulations revealed that the elastic modulus at 800% pre-elongation was stronger than that at 400% pre-elongation. From this observation, we can derive the structure-property relationship, wherein the magnitude of the pre-elongation governs the crystal domain structures and mechanical properties.

Cs+ Promoting the Diffusion of K+ and Inhibiting the Generation of Newberyite in Struvite-K Cements: Experiments and Molecular Dynamics Simulation Calculations

Struvite-K cements, also called magnesium potassium phosphate cements (MKPCs), are applicable for particular applications, especially the immobilization of radioactive Cs+ in the nuclear industry. This work focuses on how Cs+ affects the hydration mechanism of struvite-K cements because newberyite and brucite in the hydration products are deemed to be risky products that result in cracking. Experiments and molecular dynamics simulations showed that Cs+ promoted the diffusion of K+ to the surface of MgO, which greatly facilitates the formation of more K-struvite crystals, inhibiting the formation of newberyite and brucite. A total of 0.02 M Cs+ resulted in a 40.44%, 13.93%, 60.81%, and 32.18% reduction in the amount of newberyite and brucite, and the Cs immobilization rates were 99.07%, 99.84%, 99.87%, and 99.83% when the ratios of Mg/P were 1, 3, 5, and 7, respectively. This provides new evidence of stability for struvite-K cements on radioactive Cs+ immobilization. Surprisingly, another new crystal, [CsPO3·H2O]4, was found to be a dominating Cs-containing phase in Cs-immobilizing struvite-K cements, in addition to Cs-struvite.

Discovery of Quinolinone Hybrids as Dual Inhibitors of Acetylcholinesterase and Aβ Aggregation for Alzheimer's Disease Therapy

The development of multitargeted therapeutics has evolved as a promising strategy to identify efficient therapeutics for neurological disorders. We report herein new quinolinone hybrids as dual inhibitors of acetylcholinesterase (AChE) and Aβ aggregation that function as multitargeted ligands for Alzheimer's disease. The quinoline hybrids (AM1-AM16) were screened for their ability to inhibit AChE, BACE1, amyloid fibrillation, α-syn aggregation, and tau aggregation. Among the tested compounds, AM5 and AM10 inhibited AChE activity by more than 80% at single-dose screening and possessed a remarkable ability to inhibit the fibrillation of Aβ42 oligomers at 10 μM. In addition, dose-dependent screening of AM5 and AM10 was performed, giving half-maximal AChE inhibitory concentration (IC50) values of 1.29 ± 0.13 and 1.72 ± 0.18 μM, respectively. In addition, AM5 and AM10 demonstrated concentration-dependent inhibitory profiles for the aggregation of Aβ42 oligomers with estimated IC50 values of 4.93 ± 0.8 and 1.42 ± 0.3 μM, respectively. Moreover, the neuroprotective properties of the lead compounds AM5 and AM10 were determined in SH-SY5Y cells incubated with Aβ oligomers. This work would enable future research efforts aiming at the structural optimization of AM5 and AM10 to develop potent dual inhibitors of AChE and amyloid aggregation. Furthermore, the in vivo assay confirmed the antioxidant activity of compounds AM5 and AM10 through increasing GSH, CAT, and SOD activities that are responsible for scavenging the ROS and restoring its normal level. Blood investigation illustrated the protective activity of the two compounds against lead-induced neurotoxicity through retaining hematological and liver enzymes near normal levels. Finally, immunohistochemistry investigation revealed the inhibitory activity of β-amyloid (Aβ) aggregation.

In-vitro antimicrobial activity of AF-DP protein and in-silico approach of cell membrane disruption

Microbial resistance against common antibiotics has become one of the most serious threats to human health. The increasing statistics on this problem show the necessity of finding a way to deal with it. In recent years, antimicrobial peptides with unique properties and the capability of targeting a wide range of pathogens, have been considered as a potential for replacing common antibiotics. A small chitin-binding protein with anticandidal activity was isolated from Moringa oleifera seeds by Neto and colleagues in 2017, which very much resembled antimicrobial peptides. In this study, the antimicrobial protein 'AF-DP' was identified and characterized. AF-DP was heterologously expressed, purified, and characterized, and its 3D structure was predicted. Six molecular dynamic simulations were performed to investigate how the protein interacts with Gram-negative inner and outer, Gram-positive, fungal, cancerous, and normal mammalian membranes. Also, its antimicrobial and anticancer activity was assessed in vitro via minimum inhibition concentration (MIC) and MTT assays, respectively. This protein with 111 amino acids and a total net charge (of 10.5) has been predicted to be mainly composed of alpha helix and random coils. Its MIC affecting the growth of Escherichia coli, Staphylococcus aureus, and Candida albicans was 30 µg/ml, 100 µg/ml, and 100 µg/ml, respectively; AF-DP showed anticancer activity against MCF-7 breast cancer cell line. Scanning electron microscopic analysis confirmed the creation of pores and scratches on the surface of the bacterial membrane. The results of this research show that AF-DP can be a candidate for the production of new drugs as an AMP with antimicrobial activity.Communicated by Ramaswamy H. Sarma.

Transformations in crystals of DNA-functionalized nanoparticles by electrolytes

Colloidal crystals have applications in water treatments, including water purification and desalination technologies. It is, therefore, important to understand the interactions between colloids as a function of electrolyte concentration. We study the assembly of DNA-grafted gold nanoparticles immersed in concentrated electrolyte solutions. Increasing the concentration of divalent Ca2+ ions leads to the condensation of nanoparticles into face-centered-cubic (FCC) crystals at low electrolyte concentrations. As the electrolyte concentration increases, the system undergoes a phase change to body-centered-cubic (BCC) crystals. This phase change occurs as the interparticle distance decreases. Molecular dynamics analysis suggests that the interparticle interactions change from strongly repulsive to short-range attractive as the divalent-electrolyte concentration increases. A thermodynamic analysis suggests that increasing the salt concentration leads to significant dehydration of the nanoparticle environment. We conjecture that the intercolloid attractive interactions and dehydrated states favour the BCC structure. Our results gain insight into salting out of colloids such as proteins as the concentration of salt increases in the solution.

Properties of aqueous electrolyte solutions at carbon electrodes: effects of concentration and surface charge on solution structure, ion clustering and thermodynamics in the electric double layer

Surfaces are able to control physical-chemical processes in multi-component solution systems and, as such, find application in a wide range of technological devices. Understanding the structure, dynamics and thermodynamics of non-ideal solutions at surfaces, however, is particularly challenging. Here, we use Constant Chemical Potential Molecular Dynamics (CμMD) simulations to gain insight into aqueous NaCl solutions in contact with graphite surfaces at high concentrations and under the effect of applied surface charges: conditions where mean-field theories describing interfaces cannot (typically) be reliably applied. We discover an asymmetric effect of surface charge on the electric double layer structure and resulting thermodynamic properties, which can be explained by considering the affinity of the surface for cations and anions and the cooperative adsorption of ions that occurs at higher concentrations. We characterise how the sign of the surface charge affects ion densities and water structure in the double layer and how the capacitance of the interface-a function of the electric potential drop across the double layer-is largely insensitive to the bulk solution concentration. Notably, we find that negatively charged graphite surfaces induce an increase in the size and concentration of extended liquid-like ion clusters confined to the double layer. Finally, we discuss how concentration and surface charge affect the activity coefficients of ions and water at the interface, demonstrating how electric fields in this region should be explicitly considered when characterising the thermodynamics of both solute and solvent at the solid/liquid interface.

Molecular Dynamics Study on the Effect of Polyacrylamide on Electric Field Demulsification of Oil-in-Water Emulsion

The effect of the water-soluble polymer (partially hydrolyzed polyacrylamide, HPAM) in produced water on the demulsification process of the electric field was studied by molecular dynamics simulations. By comparing the coalescence process of oil droplets in the electric field environment with or without HPAM, we find that HPAM in the water phase can promote the coalescence of nearly oil droplets but hinder the deformation and migration of oil droplets. By analyzing the radial distribution function and interaction energy between molecules, we conclude that the existence of HPAM molecules can reduce the hydrophilicity of other molecules through their strong interaction with water, and sodium ions (Na+) have strong interaction with bound water in the process of breaking away from HPAM, thus leading the movement of water molecules. At the same time, the influence of HPAM molecules located between the two oil droplets on the demulsification process was also studied. The HPAM molecules and sodium ions located between the two oil droplets also affected the coalescence process of oil droplets under an electric field by interacting with water.

Unveiling the structure-emulsifying function relationship of truncated recombinant forms of the SA01-OmpA protein opens up a new vista in bioemulsifiers

The emulsifying ability of SA01-OmpA (outer membrane protein A from Acinetobacter sp. SA01) was found to be constrained by challenges like low production efficiency and high costs associated with protein recovery from E. coli inclusion bodies, as described in our previous study. The present study sought to benefit from the advantages of the targeted truncating of SA01-OmpA protein, taking into account the reduced propensity of protein expression as inclusion bodies and cytotoxicity. Here, the structure and activity relationship of two truncated recombinant forms of SA01-OmpA protein was unraveled through a hybrid approach based on experimental data and computational methodologies, representing an innovative bioemulsifier with advantageous emulsifying activity. The recombinant truncated SA01-OmpA variants were cloned and heterologously expressed in E. coli host cells and subsequently purified. The results showed increased emulsifying activity of N-terminally truncated SA01-OmpA (NT-OmpA) compared to full-length SA01-OmpA. Molecular dynamics (MD) simulations analysis demonstrated a direct correlation between the C-terminally truncated SA01-OmpA (CT-OmpA) and its expression as inclusion bodies. Analysis of the structure-activity relationship of truncated variants of SA01-OmpA revealed that, compared to the full-length protein, deletion of the β-barrel portion from the N-terminal of SA01-OmpA increased the emulsifying activity of NT-OmpA while lowering its expression as inclusion bodies. Contrary to the full-length protein, the N-terminally truncated SA01-OmpA was not as cytotoxic, according to the MTT assay, FCM analysis, and AO/EB staining. The findings of this extensive study advance our knowledge of SA01-OmpA at the molecular level as well as the design and development of efficient bioemulsifiers.IMPORTANCEPrevious research (Shahryari et al. 2021, mSystems 6: e01175-20) introduced and characterized the SA01-OmpA protein as a multifaceted protein with a variety of functions, including maintaining cellular homeostasis under oxidative stress conditions, biofilm formation, outer membrane vesicles (OMV) biogenesis, and beneficial emulsifying capacity. By truncating the SA01-OmpA protein, the current study presents a unique method for developing protein-type bioemulsifiers. The findings indicate that the N-terminally truncated SA01-OmpA (NT-OmpA) has the potential to fully replace full-length SA01-OmpA as a novel bioemulsifier with significant emulsifying activity. This study opens up a new frontier in bioemulsifiers, shedding light on a possible relationship between the structure and activity of SA01-OmpA truncated forms.

N-(Heteroaryl)thiophene sulfonamides as angiotensin AT2 receptor ligands

Two series of N-(heteroaryl)thiophene sulfonamides, encompassing either a methylene imidazole group or a tert-butylimidazolylacetyl group in the meta position of the benzene ring, have been synthesized. An AT2R selective ligand with a Ki of 42 nM was identified in the first series and in the second series, six AT2R selective ligands with significantly improved binding affinities and Ki values of <5 nM were discovered. The binding modes to AT2R were explored by docking calculations combined with molecular dynamics simulations. Although some of the high affinity ligands exhibited fair stability in human liver microsomes, comparable to that observed with C21 undergoing clinical trials, most ligands displayed a very low metabolic stability with t½ of less than 10 min in human liver microsomes. The most promising ligand, with an AT2R Ki value of 4.9 nM and with intermediate stability in human hepatocytes (t½ = 77 min) caused a concentration-dependent vasorelaxation of pre-contracted mouse aorta.

Multiscale Modeling of Macromolecular Interactions between Tau-Amylin Oligomers and Asymmetric Lipid Nanodomains That Link Alzheimer's and Diabetic Diseases

The molecular events of protein misfolding and self-aggregation of tau and amylin are associated with the progression of Alzheimer's and diabetes, respectively. Recent studies suggest that tau and amylin can form hetero-tau-amylin oligomers. Those hetero-oligomers are more neurotoxic than homo-tau oligomers. So far, the detailed interactions between the hetero-oligomers and the neuronal membrane are unknown. Using multiscale MD simulations, the lipid binding and protein folding behaviors of hetero-oligomers on asymmetric lipid nanodomains or raft membranes were examined. Our raft membranes contain phase-separated phosphatidylcholine (PC), cholesterol, and anionic phosphatidylserine (PS) or ganglioside (GM1) in one leaflet of the lipid bilayer. The hetero-oligomers bound more strongly to the PS and GM1 than other lipids via the hydrophobic and hydrophilic interactions, respectively, in the raft membranes. The hetero-tetramer disrupted the acyl chain orders of both PC and PS in the PS-containing raft membrane, but only the GM1 in the GM1-containing raft membrane as effectively as the homo-tau-tetramer. We discovered that the alpha-helical content in the heterodimer was greater than the sum of alpha-helical contents from isolated tau and amylin monomers on both raft membranes, indicative of a synergetic effect of tau-amylin interactions in surface-induced protein folding. Our results provide new molecular insights into understanding the cross-talk between Alzheimer's and diabetes.

Ten-eleven translocation-2-mediated macrophage activation promotes liver regeneration

BACKGROUND

The remarkable regenerative capacity of the liver enables recovery after radical Hepatocellular carcinoma (HCC) resection. After resection, macrophages secrete interleukin 6 and hepatocyte growth factors to promote liver regeneration. Ten-eleven translocation-2 (Tet2) DNA dioxygenase regulates pro-inflammatory factor secretion in macrophages. In this study, we explored the role of Tet2 in macrophages and its function independent of its enzymatic activity in liver regeneration.

METHODS

The model of liver regeneration after 70% partial hepatectomy (PHx) is a classic universal model for studying reparative processes in the liver. Mice were euthanized at 0, 24, and 48 h after PHx. Enzyme-linked immunosorbent assays, quantitative reverse transcription-polymerase chain reaction, western blotting, immunofluorescence analysis, and flow cytometry were performed to explore immune cell infiltration and liver regenerative capability. Molecular dynamics simulations were performed to study the interaction between Tet2 and signal transducer and activator of transcription 1 (Stat1).

RESULTS

Tet2 in macrophages negatively regulated liver regeneration in the partial hepatectomy mice model. Tet2 interacted with Stat1, inhibiting the expression of proinflammatory factors and suppressing liver regeneration. The Tet2 inhibitor attenuated the interaction between Stat1 and Tet2, enhanced Stat1 phosphorylation, and promoted hepatocyte proliferation. The proliferative function of the Tet2 inhibitor relied on macrophages and did not affect hepatocytes directly.

CONCLUSION

Our findings underscore that Tet2 in macrophages negatively regulates liver regeneration by interacting with Stat1. Targeting Tet2 in macrophages promotes liver regeneration and function after a hepatectomy, presenting a novel target to promote liver regeneration and function.

Case report: A novel compound heterozygous variant in the TNXB gene causes single kidney agenesis and vesicoureteral reflux

Primary vesicoureteral reflux (VUR) is the prevailing congenital anomaly of the kidneys and urinary tract, posing a significant risk for pyelonephritis scarring and chronic renal insufficiency in pediatric patients. Nevertheless, the precise genetic etiology of VUR remains enigmatic. In this current investigation, we conducted whole-exome sequencing on a child exhibiting single kidney, devoid of any familial VUR background, along with both biological parents. Two missense variants (NM_019105.8: exon11: c.4111G>A and NM_019105.8: exon2: c.31A>T) in the TNXB gene were identified through whole-exome sequencing of the child. These variants were found to be inherited from the child's parents, with each parent carrying one of the variants. Molecular dynamics simulations were conducted to assess the impact of these variants on the tenascin XB proteins encoded by them, revealing varying degrees of impairment. Based on our findings, it is suggested that the TNXB compound heterozygous variant, consisting of c.4111G>A and c.31A>T, may be the underlying cause of right renal agenesis and left hydronephrosis in afflicted child. This discovery broadens the genetic range of the TNXB gene and establishes a genetic foundation for disease-specific preimplantation genetic diagnosis (PGD) in prospective pregnancies involving the parents of this afflicted child.

Understanding the Positive Role of Ionic Liquids in CO2 Capture by Poly(ethylenimine)

CO2 capture technology is one of the most important technical methods for significantly mitigating CO2 emissions in a low-carbon context. The recent invention of mixed absorbents comprising poly(ethylenimine) (PEI) and ionic liquids (ILs) provides a novel strategy for efficiently capturing CO2, and this has garnered widespread attention. However, the intermolecular interactions between the IL and other constituents during the CO2 absorption process remain unclear. In this present work, a series of density functional theory (DFT) calculations and molecular dynamics simulations were conducted to investigate the positive role of IL in CO2 capture by PEI. The results showed that the formation of hydrogen bonds between the IL anion and the amino groups of PEI primarily drives the addition of IL to PEI. During the CO2 absorption process, the IL anion not only can absorb CO2 but also exerts a dehydrogenation effect on the amino group of PEI, facilitating enhanced interaction between PEI and CO2. Additionally, the IL substantially reduces the viscosity of PEI, promoting the diffusion of CO2 within the system and enhancing the absorption rate. Based on the information on interaction energy and viscosity, we can easily make theoretical predictions for the optimal proportion of IL to be added. The above results provide fundamental insights to promote the industrial application of the PEI/IL system for CO2 capture.

Computing Accurate True Self-Diffusion Coefficients and Shear Viscosities Using the OrthoBoXY Approach

In a recent paper [Busch, J.; Paschek, D. J. Phys. Chem. B 2023, 127, 7983-7987], we have shown that for molecular dynamics (MD) simulations using orthorhombic periodic boundary conditions with "magic" box length ratios of Lz/Lx = Lz/Ly = 2.7933596497, the self-diffusion coefficients Dx and Dy in x- and y-directions are independent of the system size. They both represent the true self-diffusion coefficient D0 = (Dx + Dy)/2, while the shear viscosity can be calculated from diffusion coefficients in x-, y-, and z-directions, using η = kBT·8.1711245653/[3πLz(Dx + Dy - 2Dz)]. In this contribution, we test this "OrthoBoXY" approach by its application to a variety of different systems: liquid water, dimethyl ether, methanol, triglyme, water/methanol mixtures, water/triglyme mixtures, and imidazolium-based ionic liquids. The chosen systems range from small-sized molecular liquids to complex mixtures and ionic liquids, while spanning a viscosity range of almost 3 orders of magnitude. We assess the efficiency of the method for computing true self-diffusion and viscosity data and provide simple formulas for estimating the required MD simulation lengths and sizes for delivering reliable data with targeted uncertainty levels. Our analysis of the system size dependence of statistical uncertainties for both the viscosity and the self-diffusion coefficient leads us to the conclusion that it is preferable to extend the simulation length instead of increasing the system size. MD simulations consisting of 768 molecules or ion pairs seem to be perfectly adequate.

A rationally engineered small antimicrobial peptide with potent antibacterial activity

Antimicrobial resistance (AMR) is a silent pandemic declared by the WHO that requires urgent attention in the post-COVID world. AMR is a critical public health concern worldwide, potentially affecting people at different stages of life, including the veterinary and agriculture industries. Notably, very few new-age antimicrobial agents are in the current developmental pipeline. Thus, the design, discovery, and development of new antimicrobial agents are required to address the menace of AMR. Antimicrobial peptides (AMPs) are an important class of antimicrobial agents for combating AMR due to their broad-spectrum activity and ability to evade AMR through a multimodal mechanism of action. However, molecular size, aggregability, proteolytic degradation, cytotoxicity, and hemolysis activity significantly limit the clinical application of natural AMPs. The de novo design and engineering of a short synthetic amphipathic AMP (≤16 aa, Mol. Wt. ≤ 2 kDa) with an unusual architecture comprised of coded and noncoded amino acids (NCAAs) is presented here, which demonstrates potent antibacterial activity against a few selected bacterial strains mentioned in the WHO priority list. The designer AMP is conformationally ordered in solution and effectively permeabilizes the outer and inner membranes, leading to bacterial growth inhibition and death. Additionally, the peptide is resistant to proteolysis and has negligible cytotoxicity and hemolysis activity up to 150 μM toward cultured human cell lines and erythrocytes. The designer AMP is unique and appears to be a potent therapeutic candidate, which can be subsequently subjected to preclinical studies to explicitly understand and address the menace of AMR.

Ebola virus VP35 interacts non-covalently with ubiquitin chains to promote viral replication

Ebolavirus (EBOV) belongs to a family of highly pathogenic viruses that cause severe hemorrhagic fever in humans. EBOV replication requires the activity of the viral polymerase complex, which includes the cofactor and Interferon antagonist VP35. We previously showed that the covalent ubiquitination of VP35 promotes virus replication by regulating interactions with the polymerase complex. In addition, VP35 can also interact non-covalently with ubiquitin (Ub); however, the function of this interaction is unknown. Here, we report that VP35 interacts with free (unanchored) K63-linked polyUb chains. Ectopic expression of Isopeptidase T (USP5), which is known to degrade unanchored polyUb chains, reduced VP35 association with Ub and correlated with diminished polymerase activity in a minigenome assay. Using computational methods, we modeled the VP35-Ub non-covalent interacting complex, identified the VP35-Ub interacting surface, and tested mutations to validate the interface. Docking simulations identified chemical compounds that can block VP35-Ub interactions leading to reduced viral polymerase activity. Treatment with the compounds reduced replication of infectious EBOV in cells and in vivo in a mouse model. In conclusion, we identified a novel role of unanchored polyUb in regulating Ebola virus polymerase function and discovered compounds that have promising anti-Ebola virus activity.

Inhibitor design for TMPRSS2: insights from computational analysis of its backbone hydrogen bonds using a simple descriptor

Transmembrane protease serine 2 (TMPRSS2) is an important drug target due to its role in the infection mechanism of coronaviruses including SARS-CoV-2. Current understanding regarding the molecular mechanisms of known inhibitors and insights required for inhibitor design are limited. This study investigates the effect of inhibitor binding on the intramolecular backbone hydrogen bonds (BHBs) of TMPRSS2 using the concept of hydrogen bond wrapping, which is the phenomenon of stabilization of a hydrogen bond in a solvent environment as a result of being surrounded by non-polar groups. A molecular descriptor which quantifies the extent of wrapping around BHBs is introduced for this. First, virtual screening for TMPRSS2 inhibitors is performed by molecular docking using the program DOCK 6 with a Generalized Born surface area (GBSA) scoring function. The docking results are then analyzed using this descriptor and its relationship to the solvent-accessible surface area term ΔGsa of the GBSA score is demonstrated with machine learning regression and principal component analysis. The effect of binding of the inhibitors camostat, nafamostat, and 4-guanidinobenzoic acid (GBA) on the wrapping of important BHBs in TMPRSS2 is also studied using molecular dynamics. For BHBs with a large increase in wrapping groups due to these inhibitors, the radial distribution function of water revealed that certain residues involved in these BHBs, like Gln438, Asp440, and Ser441, undergo preferential desolvation. The findings offer valuable insights into the mechanisms of these inhibitors and may prove useful in the design of new inhibitors.

Transformation of Diffusion and Local Structure of CH4 , CO2 , SO2 and H2 O Mixtures in Graphene Under Wide Temperature and Pressure Range: A Molecular Dynamics Simulation Study

As a material with high specific surface area and excellent chemical stability, graphene exhibited remarkable adsorption and separation performance as well as a wide range of potential applications. The graphene layer played a significant role in influencing gas transmission. In this study, we employed molecular dynamics simulation to investigate the diffusion characteristics and local structures of a mixed system consisting of CH4 , CO2 , SO2 and H2 O. Additionally, we further examined the transformation of the behavior of these mixtures within graphene layers. The order of diffusion coefficients of the four molecules without graphene was H2 O>SO2 >CO2 ≫CH4 . However, in the double-layer graphene, the order changed to CH4 >CO2 ≫H2 O>SO2 . Higher temperatures and lower pressures were found to facilitate gas diffusion. Temperature and pressure had great effects on the local structures of CH4 , CO2 and SO2 , while their impact on H2 O was limited due to the extensive network of hydrogen bonds formed by H2 O molecules. The statistical results of average coordination number revealed that CH4 tended to aggregate with itself, whereas CO2 and SO2 exhibited a tendency to aggregate with H2 O. The graphene structure enhanced the separation and transportation of CH4 from mixed systems.

In-silico study MM/GBSA binding free energy and molecular dynamics simulation of some designed remdesivir derivatives as the inhibitory potential of SARS-CoV-2 main protease

Background and purpose

Coronavirus disease (COVID-19) is one of the greatest challenges of the twentieth century. Recently, in silico tools help to predict new inhibitors of SARS-CoV-2. In this study, the new compounds based on the remdesivir structure (12 compounds) were designed.

Experimental approach

The main interactions of remdesivir and designed compounds were investigated in the 3CLpro active site. The binding free energy of compounds by the MM-GBSA method was calculated and the best compound (compound 12 with the value of -88.173 kcal/mol) was introduced to the molecular dynamics simulation study.

Findings/Results

The simulation results were compared with the results of protein simulation without the presence of an inhibitor and in the presence of remdesivir. Additionally, the RMSD results for the protein backbone showed that compound 12 in the second 50 nanoseconds has less fluctuation than the protein alone and in the presence of remdesivir, which indicates the stability of the compound in the active site of the Mpro protein. Furthermore, protein compactness was investigated in the absence of compounds and the presence of compound 12 and remdesivir. The Rg diagram shows a fluctuation of approximately 0.05 A, which indicates the compressibility of the protein in the presence and absence of compounds. The results of the RMSF plot also show the stability of essential amino acids during protein binding.

Conclusion and implications

Supported by the theoretical results, compound 12 could have the potential to inhibit the 3CLpro enzyme, which requires further in vitro studies and enzyme inhibition must also be confirmed at protein levels.

Structure-based virtual screening and molecular dynamics approaches to identify new inhibitors of Staphylococcus aureus sortase A

Staphylococcus aureus is a prevalent Gram-positive bacteria leading cause of a wide range of human pathologies. Moreover, antibiotic résistance of pathogenesis bacteria is one of the worldwide health problems. In Gram-positive bacteria, the enzyme of SrtA, is responsible for the anchoring of surface-exposed proteins to the cell wall peptidoglycan. Because of its critical role in Gram-positive bacterial pathogenesis, SrtA is an attractive target for anti-virulence during drug development. To date, some SrtA inhibitors have been discovered most of them being derived from flavonoid compounds, like Myricetin. In order to provide potential hit molecules against SrtA for clinical use, we obtained a total of 293 compounds by performing in silico shape-based screening of compound libraries against Myristin as a reference structure. Employing molecular docking and scoring functions, the top 3 compounds Apigenin, Efloxate, and Compound 8261032 were screened by comparing their docking scores with Myricetin. Furthermore, MD simulations and MM-PBSA binding energy calculation studies revealed that only Compound 8261032 strongly binds to the catalytic core of the SrtA enzyme than Myricetin, and stable behavior was consistently observed in the docking complex. Compound 8261032 showed a good number of hydrogen bonds with SrtA and higher MM-PBSA binding energy when compared to all three molecules. Also, it makes strength interactions with Arg139 and His62, which are critical for SrtA biological activity. This study showed that the development of this inhibitor could be a fundamental strategy against resistant bacteria, but further studies in vitro are needed to confirm this claim.Communicated by Ramaswamy H. Sarma.

Structural and functional interactions between the EF hand domain and S2-S3 loop in the type-1 ryanodine receptor ion channel

Previous cryo-electron micrographs suggested that the skeletal muscle Ca2+ release channel, ryanodine receptor (RyR)1, is regulated by intricate interactions between the EF hand Ca2+ binding domain and the cytosolic loop (S2-S3 loop). However, the precise molecular details of these interactions and functional consequences of the interactions remain elusive. Here, we used molecular dynamics simulations to explore the specific amino acid pairs involved in hydrogen bond interactions within the EF hand-S2-S3 loop interface. Our simulations unveiled two key interactions: (1) K4101 (EF hand) with D4730 (S2-S3 loop) and (2) E4075, Q4078, and D4079 (EF hand) with R4736 (S2-S3 loop). To probe the functional significance of these interactions, we constructed mutant RyR1 complementary DNAs and expressed them in HEK293 cells for [3H]ryanodine binding assays. Our results demonstrated that mutations in the EF hand, specifically K4101E and K4101M, resulted in reduced affinities for Ca2+/Mg2+-dependent inhibitions. Interestingly, the K4101E mutation increased the affinity for Ca2+-dependent activation. Conversely, mutations in the S2-S3 loop, D4730K and D4730N, did not significantly change the affinities for Ca2+/Mg2+-dependent inhibitions. Our previous finding that skeletal disease-associated RyR1 mutations, R4736Q and R4736W, impaired Ca2+-dependent inhibition, is consistent with the current results. In silico mutagenesis analysis aligned with our functional data, indicating altered hydrogen bonding patterns upon mutations. Taken together, our findings emphasize the critical role of the EF hand-S2-S3 loop interaction in Ca2+/Mg2+-dependent inhibition of RyR1 and provide insights into potential therapeutic strategies targeting this domain interaction for the treatment of skeletal myopathies.

The Stokes-Einstein-Sutherland Equation at the Nanoscale Revisited

The Stokes-Einstein-Sutherland (SES) equation is at the foundation of statistical physics, relating a particle's diffusion coefficient and size with the fluid viscosity, temperature, and the boundary condition for the particle-solvent interface. It is assumed that it relies on the separation of scales between the particle and the solvent, hence it is expected to break down for diffusive transport on the molecular scale. This assumption is however challenged by a number of experimental studies showing a remarkably small, if any, violation, while simulations systematically report the opposite. To understand these discrepancies, analytical ultracentrifugation experiments are combined with molecular simulations, both performed at unprecedented accuracies, to study the transport of buckminsterfullerene C60 in toluene at infinite dilution. This system is demonstrated to clearly violate the conditions of slow momentum relaxation. Yet, through a linear response to a constant force, the SES equation can be recovered in the long time limit with no more than 4% uncertainty both in experiments and in simulations. This nonetheless requires partial slip on the particle interface, extracted consistently from all the data. These results, thus, resolve a long-standing discussion on the validity and limits of the SES equation at the molecular scale.

Genomic surveillance of omicron B.1.1.529 SARS-CoV-2 and its variants between December 2021 and March 2023 in Tamil Nadu, India-A state-wide prospective longitudinal study

A state-wide prospective longitudinal investigation of the genomic surveillance of the omicron B.1.1.529 SARS-CoV-2 variant and its sublineages in Tamil Nadu, India, was conducted between December 2021 and March 2023. The study aimed to elucidate their mutational patterns and their genetic interrelationship in the Indian population. The study identified several unique mutations at different time-points, which likely could attribute to the changing disease characteristics, transmission, and pathogenicity attributes of omicron variants. The study found that the omicron variant is highly competent in its mutating potentials, and that it continues to evolve in the general population, likely escaping from natural as well as vaccine-induced immune responses. Our findings suggest that continuous surveillance of viral variants at the global scenario is warranted to undertake intervention measures against potentially precarious SARS-CoV-2 variants and their evolution.

Investigation of Zika virus methyl transferase inhibitors using steered molecular dynamics

Zika virus (ZIKV) spread is considered a major public health threat by the World Health Organization (WHO). There are no vaccines or drugs available to control the infection of the Zika virus, therefore a highly effective medicinal molecule is urgently required. In this study, a computationally intensive investigation was performed to identify a potent natural compound that could inhibit the ZIKV NS5 methyltransferase. This research approach is based on target-based drug identification principles where the native inhibitor SAH (S-adenosylhomocysteine) of ZIKV NS5 methyltransferase was selected as a reference. High-throughput virtual screening and tanimoto similarity coefficient were applied to the natural compound library for ranking the potential candidates. The top five compounds were selected for interaction analysis, MD simulation, total binding free energy through MM/GBSA, and steered MD simulation. Among these compounds, Adenosine 5'-monophosphate monohydrate, Tubercidin, and 5-Iodotubercidin showed stable binding to the protein compared to the native compound, SAH. These three compounds also showed less fluctuations in RMSF in contrast to native compound. Additionally, the same interacting residues observed in SAH also made strong interactions with these three compounds. Adenosine 5'-monophosphate monohydrate and 5-Iodotubercidin had greater total binding free energies than the reference ligand. Moreover, the dissociation resistance of all three compounds was equivalent to that of the reference ligand. This study suggested binding properties of three-hit compounds that could be used to develop drugs against Zika virus infections.Communicated by Ramaswamy H. Sarma.

Separation of protein corona from nanoparticles under intracellular acidic conditions: effect of protonation on nanoparticle-protein and protein-protein interactions

Protein coronas separate from nanoparticles under intracellular acidic conditions however, competitive adsorption of multiple proteins and their protein network formation under different pH conditions have not yet been systematically studied at the atomic scale. Herein, we report all-atom molecular dynamics simulations of plasma proteins (human serum albumin and immunoglobulin gamma-1 chain C) adsorbed to 10 nm-sized carboxyl-terminated polystyrene (PS) nanoparticles at different protonation states that mimic extracellular and intracellular pH conditions of 7, 6-5, and 4.5. Binding free energies are calculated from umbrella sampling simulations, showing the significantly weakened binding between PS particles and proteins at the protonation state at pH 4.5, in agreement with experiments showing the separation of protein corona from nanoparticles at pH 4.5. Mixtures of multiple proteins and PS particles are also simulated, showing much less protein adsorption and protein cluster formation at the protonation state at pH 4.5 than that at higher pH values, which are further confirmed by calculating the diffusivities and hydrodynamic radii of individual proteins. In particular, electrostatic particle-protein and protein-protein interactions are significantly weakened by a combination of particle and protein protonation rather than by particle protonation alone, to an extent dependent on different proteins. These findings help explain the experimental observations regarding separation of protein corona from nanoparticles under intracellular acidic conditions at pH 4.5 but not at higher pH, supporting that acidification cannot be the only reason for this separation during the process of endosome maturation.

Stimuli-responsive and self-healing supramolecular Zn(II)-guanosine metal-organic gel for Schottky barrier diode application

Spontaneous formation of a supramolecular metal-organic hydrogel using unsubstituted guanosine as a ligand and Zn2+ ions is reported. Guanosine, in the presence of NaOH, self-assembled into a stable G-quadruplex structure, which underwent crosslinking through Zn2+ ions to afford a stable hydrogel. The gel has been characterized using several spectroscopic as well as microscopic studies. The hydrogel demonstrated excellent stimuli responsiveness towards various chemicals and pH. Furthermore, the gel exhibited intrinsic thixotropic behavior and showed self-healing and injectable properties. The optical properties of the Zn-guanosine metallo-hydrogel suggested a semiconducting nature of the gel, which has been exploited for fabricating a thin film device based on a Schottky diode interface between metal and a semiconductor. The fabricated device shows excellent charge transport characteristics and linear rectifying behavior. The findings are likely to pave the way for newer research in the area of soft electronic devices fabricated using materials synthesized by employing simple biomolecules.

Protein-activated atomic layer deposition for robust crude-oil-repellent hierarchical nano-armored membranes

Atomic layer deposition (ALD) offers unique capabilities to fabricate atomically engineered porous materials with precise pore tuning and multi-functionalization for diverse applications like advanced membrane separations towards sustainable energy-water systems. However, current ALD technique is inhibited on most non-polar polymeric membranes due to lack of accessible nucleation sites. Here, we report a facile method to efficiently promote ALD coating on hydrophobic surface of polymeric membranes via novel protein activation/sensitization. As a proof of concept, TiO2 ALD-coated membranes activated by bovine serum albumin exhibit remarkable superhydrophilicity, ultralow underwater crude oil adhesion, and robust tolerance to rigorous environments including acid, alkali, saline, and ethanol. Most importantly, excellent cyclable crude oil-in-water emulsion separation performance can be achieved. The mechanism for activation/sensitization is rooted in reactivity for a particular set of amino acids. Furthermore, the universality of protein-sensitized ALD is demonstrated using common egg white, promising numerous potential usages in biomedical engineering, environmental remediation, low-carbon manufacturing, catalysis, and beyond.

Immunoinformatics-driven approach for development of potential multi-epitope vaccine against the secreted protein FlaC of Campylobacter jejuni

Campylobacter jejuni causes a leading human gastrointestinal infection which is associated with foodborne diarrhea, stomach cramping, and fever. In the recent years, numerous multidrug-resistant strains of C. jejuni has evolved and is considered in the priority pathogens category. Therefore, an increasing demand exists to develop an effective vaccine against Campylobacteriosis. The T cell and B cell epitopes from the FlaC protein were predicted using comprehensive immunoinformatics tools. The predicted epitopes were chosen based on their antigenicity, allergenicity, and toxicity profiles. Using the bioinformatics approach various physicochemical properties of the constructed vaccine were determined. The molecular docking analysis of the vaccine with the TLRs demonstrated that TLR5 has a higher binding affinity of -1159.0 kcal/mol. Molecular dynamics simulation has confirmed the stable association of the vaccine with TLR5. The immune response of the constructed vaccine was validated using immunostimulation. Based on this study, we recommend the formulation of a multi-epitope vaccine as a promising agent to effectively combat the dreadful campylobacteriosis infection.Communicated by Ramaswamy H. Sarma.