CHARMM: the biomolecular simulation program

Electrostatic interactions of enzymes in non-aqueous conditions: insights from molecular dynamics simulations

Electrostatic interactions of enzymes and their effects on enzyme activity and stability are poorly understood in non-aqueous conditions. Here, we investigate the contribution of the electrostatic interactions on the stability and activity of enzymes in the non-aqueous environment using molecular dynamics simulations. Lipase was selected as active and lysozyme as inactive model enzymes in non-aqueous media. Hexane was used as a common non-aqueous solvent model. In agreement with the previous experiments, simulations show that lysozyme has more structural instabilities than lipase in hexane. The number of hydrogen bonds and salt bridges of both enzymes is dramatically increased in hexane. In contrast to the other opinions, we show that the increase of the electrostatic interactions in non-aqueous media is not so favorable for enzymatic function and stability. In this condition, the newly formed hydrogen bonds and salt bridges can partially denature the local structure of the enzymes. For lysozyme, the changes in electrostatic interactions occur in all domains including the active site cleft, which leads to enzyme inactivation and destabilization. Interestingly, most of the changes in electrostatic interactions of lipase occur far from the active site regions. Therefore, the active site entrance regions remain functional in hexane. The results of this study reveal how the changes in electrostatic interactions can affect enzyme stability and activity in non-aqueous conditions. Moreover, we show for the first time how some enzymes, such as lipase, remain active in a non-aqueous environment.Communicated by Ramaswamy H. Sarma.

Interactions of Li+ ions with NCS1: A potential mechanism of Li+ neuroprotective action against psychotic disorders

Li+ based drugs have been used for the treatment of psychiatric disorders due to their mood stabilizing role for decades. Recently, several studies reported the protective effect of Li+ against severe neuropathologies such as Parkinson's, Alzheimer's, and Huntington's disease. Surprisingly, despite a broad range of Li+ effects on neurological conditions, little is known about its molecular mechanism. In this study, we propose that neuronal calcium sensor 1 (NCS1), can be an effective molecular target for Li+ action. Here we show that the EF-hands in ApoNCS1 have submillimolar affinity for Li+ with Kd = 223 ± 19 μM. Li+ binding to ApoNCS1 quenches Trp emission intensity, suggesting distinct Trp sidechains environment in Li+NCS1 compared to ApoNCS1 and Ca2+NCS1. Li+ association also stabilizes the protein α-helical structure, in a similar way to Ca2+. Li+ association does not promote NCS1 dimerization. Association of Li+ increases NCS1 affinity for the D2R receptor binding peptide, in a similar way to Ca2+, however, the affinity of NCS1 for chlorpromazine is reduced with respect to Ca2+NCS1, possibly due to a decrease in solvent exposed hydrophobic area on the NCS1 surface in the presence of Li+. MD simulation data suggests that Li+ ions are coordinated by four oxygens from Asp and Glu sidechains and one carbonyl oxygen, in a similar way as reported previously for Li+ binding to DREAM. Overall, the data shows that Li+ binds to EF-hands of NCS1 and Li+NCS1 interactions may be involved in the potential neuroprotective role of Li+ against psychotic disorders.

Crystal structures of Aspergillus oryzae exo-β-(1,3)-glucanase reveal insights into oligosaccharide binding, recognition, and hydrolysis

Exo-β-(1,3)-glucanases are promising enzymes for use in the biofuel industry as they hydrolyse sugars such as laminarin, a major constituent of the algal cell wall. This study reports structural and biochemical characterizations of Aspergillus oryzae exo-β-(1,3)-glucanase (AoBgl) belonging to the GH5 family. Purified AoBgl hydrolyses β-(1,3)-glycosidic linkages of the oligosaccharide laminaritriose and the polysaccharide laminarin effectively. We have determined three high-resolution structures of AoBgl: (a) the apo form at 1.75 Å, (b) the complexed form with bound cellobiose at 1.73 Å and (c) the glucose-bound form at 1.20 Å. The crystal structures, molecular dynamics simulation studies and site-directed mutagenesis reveal the mode of substrate binding and interactions at the active site. The results also indicate that AoBgl effectively hydrolyses trisaccharides and higher oligosaccharides. The findings from our structural and biochemical studies would aid in rational engineering efforts to generate superior AoBgl variants and similar GH5 enzymes for their industrial use.

Insight into Protein Engineering: From In silico Modelling to In vitro Synthesis

Protein engineering alters the polypeptide chain to obtain a novel protein with improved functional properties. This field constantly evolves with advanced in silico tools and techniques to design novel proteins and peptides. Rational incorporating mutations, unnatural amino acids, and post-translational modifications increases the applications of engineered proteins and peptides. It aids in developing drugs with maximum efficacy and minimum side effects. Currently, the engineering of peptides is gaining attention due to their high stability, binding specificity, less immunogenic, and reduced toxicity properties. Engineered peptides are potent candidates for drug development due to their high specificity and low cost of production compared with other biologics, including proteins and antibodies. Therefore, understanding the current perception of designing and engineering peptides with the help of currently available in silico tools is crucial. This review extensively studies various in silico tools available for protein engineering in the prospect of designing peptides as therapeutics, followed by in vitro aspects. Moreover, a discussion on the chemical synthesis and purification of peptides, a case study, and challenges are also incorporated.

A novel synthetic peptide analog enhanced antibacterial activity of the frog-derived skin peptide wuchuanin-A1

In recent years, there has been a growing focus on the development of novel antibacterial compounds for clinical applications, such as antimicrobial peptide (AMP). Among the developed AMP, wuchuanin-A1, a coil-shaped bioactive peptide derived from Odorrana wuchuanensis frog skin, has been reported to exhibit antibacterial, antifungal, and antioxidant activity, but there are limited studies on its potential as an antibacterial agent. Therefore, this study aims to molecularly modify the sequence of wuchuanin-A1 to enhance its antibacterial properties. The interaction of both the native and analog peptide with bacterial inner membranes was initially assessed using computational methods. Specific amino acid substitutions were then used to enhance the modified peptide's antibacterial efficacy, followed by several preliminary tests to evaluate its activity. This study bridges the gap in exploring the potential of wuchuanin-A1 for antibacterial purposes, providing insights into the design of effective antimicrobial agents.Communicated by Ramaswamy H. Sarma.

In silico docking and molecular dynamics for the discovery of inhibitors of enteric methane production in ruminants - A review

The increase in methane emissions, a major greenhouse gas, threatens human well-being and global ecosystems due to its contribution to global warming. Livestock, particularly ruminants, have been a major research topic in recent decades due to their methane production. Therefore, the objective of the current review was to comprehensively discuss the in silico techniques used to mitigate methane production from ruminants. The review covers the principles of in silico docking and molecular dynamics, which can be used to develop methanogenesis inhibitors. It also discusses specific methanogen enzymes as potential targets for inhibitor development. Furthermore, in silico-based methanogenesis inhibitor development studies have been reviewed with the authors' opinions. The further use of in silico-based research techniques, including artificial intelligence-based systems, is encouraged to help reduce methane production from livestock more efficiently and costeffectively.

Development and characterization of fluorescent cholesteryl probes with enhanced solvatochromic and pH-sensitive properties for live-cell imaging

We present novel fluorescent cholesteryl probes (CNDs) with a modular design based on the solvatochromic 1,8-phthalimide scaffold. We have explored how different modules-linkers and head groups-affect the ability of these probes to integrate into lipid membranes and how they distribute intracellularly in mouse astrocytes and fibroblasts targeting lysosomes and lipid droplets. Each compound was assessed for its solvatochromic behavior in organic solvents and model membranes. Molecular dynamics simulations and lipid partitioning using giant unilamellar vesicles showed how these analogs behave in model membranes compared to cholesterol. Live-cell imaging demonstrated distinct staining patterns and cellular uptake behaviors, further validating the utility of these probes in biological systems. We compared the empirical results with those of BODIPY-cholesterol, a well-regarded fluorescent cholesterol analog. The internalization efficiency of fluorescent CND probes varies in different cell types and is affected mainly by the head groups. Our results demonstrate that the modular design significantly simplifies the creation of fluorescent cholesteryl probes bearing distinct spectral, biophysical, and cellular targeting features. It is a valuable toolkit for imaging in live cells, measuring cellular membrane dynamics, and studying cholesterol-related processes.

Structure-guided design of a peripherally restricted chemogenetic system

Designer receptors exclusively activated by designer drugs (DREADDs) are chemogenetic tools for remotely controlling cellular signaling, neural activity, behavior, and physiology. Using a structure-guided approach, we provide a peripherally restricted Gi-DREADD, hydroxycarboxylic acid receptor DREADD (HCAD), whose native receptor is minimally expressed in the brain, and a chemical actuator that does not cross the blood-brain barrier (BBB). This was accomplished by combined mutagenesis, analoging via an ultra-large make-on-demand library, structural determination of the designed DREADD receptor via cryoelectron microscopy (cryo-EM), and validation of HCAD function. Expression and activation of HCAD in dorsal root ganglion (DRG) neurons inhibit action potential (AP) firing and reduce both acute and tissue-injury-induced inflammatory pain. The HCAD chemogenetic system expands the possibilities for studying numerous peripheral systems with little adverse effects on the central nervous system (CNS). The structure-guided approach used to generate HCAD also has the potential to accelerate the development of emerging chemogenetic tools for basic and translational sciences.

Quercetin inhibits mitophagy-mediated apoptosis and inflammatory response by targeting the PPARγ/PGC-1α/NF-κB axis to improve acute liver failure

BACKGROUND

Reactive oxygen species (ROS) from mitochondrial dysfunction are critical in triggering apoptosis and inflammation in acute liver failure (ALF). Quercetin (QUE), an antioxidant, is renowned for its therapeutic effects onliverdiseases. There are no studies on whether QUE regulates mitophagy level in hepatocytes to inhibit ALF.

OBJECTIVE

This study investigates QUE's protective effects on ALF and elucidates the mechanisms involved.

METHODS

The ALF and hepatocyte inflammatory injury model was established using LPS and D-Galn. To predict potential targets and mechanisms of QUE in ALF treatment, transcriptomics, network pharmacology, molecular docking techniques, and ChIP were employed. The expression level related to mitophagy, apoptosis, and signaling pathways were detected by CCK8, IHC, IF staining, TUNEL, RT-qPCR, TEM, Western blotting, ELISA, and flow cytometry.

RESULTS

Network pharmacology and transcriptomics revealed common targets between QUE and ALF. Enrichment analysis showed that the anti-ALF targets of QUE were significantly associated with mitochondria and NF-κB-related pathways. Subsequent experiments showed that QUE pretreatment significantly alleviated the loss of hepatocyte viability, enhanced mitochondrial membrane potential, activated mitophagy, and promoted the clearance of damaged mitochondria, thereby reducing ROS accumulation, significantly reducing cell apoptosis and inflammatory responses, reducing ALT and AST levels, and improving liver tissue pathology. Mechanistically, molecular docking, DARTS, and CETSA analyses confirmed that QUE directly binds to the PPARγ molecule, which reduced binding to IκB and significantly inhibit the NF-κB pathway to exert its protective effects.

CONCLUSION

In short, our results provide the first evidence that QUE improves acute liver failure by promoting mitophagy through regulating the PPARγ/PGC-1α/NF-κB axis and inhibiting apoptosis and inflammatory responses mediated by mitochondrial dysfunction, which provides evidence for the potential of QUE in the treatment of ALF.

Extended Sampling of Macromolecular Conformations from Uniformly Distributed Points on Multidimensional Normal Mode Hyperspheres

Proteins are dynamic entities that adopt diverse conformations, which play a pivotal role in their function. Understanding these conformations is essential, and protein collective motions, particularly those captured by normal mode (NM) and their linear combinations, provide a robust means for conformational sampling. This work introduces a novel approach to obtaining a uniformly oriented set of a given number of lowest frequency NM combined vectors and generating harmonically equidistant restrained structures along them. They are all thus uniformly located on concentric hyperspheres, systematically covering the defined NM space fully. The generated structures are further relaxed with standard molecular dynamics (MD) simulations to explore the conformational space. The efficiency of the approach we termed "distributed points Molecular Dynamics using Normal Modes" (dpMDNM) was assessed by applying it to hen egg-white lysozyme and human cytochrome P450 3A4 (CYP3A4). To this purpose, we compared our new approach with other methods and analyzed the sampling of existing experimental structures. Our results demonstrate the efficacy of dpMDNM in extensive conformational sampling, particularly as more NMs are considered. Ensembles generated by dpMDNM exhibited a broad coverage of the experimental structures, providing valuable insights into the functional aspects of lysozyme and CYP3A4. Furthermore, dpMDNM also covered lysozyme structures with relatively elevated energies corresponding to transient states not easily obtained by standard MD simulations, in conformity with nuclear magnetic resonance structural indications. This method offers an efficient and rational framework for comprehensive protein conformational sampling, contributing significantly to our understanding of protein dynamics and function.

drMD: Molecular Dynamics for Experimentalists

Molecular dynamics (MD) simulations can be used by protein scientists to investigate a wide array of biologically relevant properties such as the effects of mutations on a protein's structure and activity, or probing intermolecular interactions with small molecule substrates or other macromolecules. Within the world of computational structural biology, several programs have become popular for running these simulations, but each of these programs requires a significant time investment from the researcher to run even simple simulations. Even after learning how to run and analyse simulations, many elements remain a "black box." This greatly limits the accessibility of molecular dynamics simulations for non-experts. Here we present drMD, an automated pipeline for running MD simulations using the OpenMM molecular mechanics toolkit. We have created drMD with non-experts in computational biology in mind. The drMD codebase has several functions that automatically handle routine procedures associated with running MD simulations. This greatly reduces the expertise required to run MD simulations. We have also introduced a series of quality-of-life features to make the process of running MD simulations both easier and more pleasant. Finally, drMD explains the steps it is taking interactively and, where useful, provides relevant references so the user can learn more. All these features make drMD an effective tool for learning MD while running publication-quality simulations. drMD is open source and can be found on GitHub: https://github.com/wells-wood-research/drMD.

Structure-Function Relationships of the CMP-Sialic Acid Transporter through Analysis of a Pathogenic Variant in an Alternatively Spliced Functional Isoform

The human CMP-sialic acid transporter (hCST) is a mammalian highly conserved type III antiporter that translocates CMP-sialic acid into the Golgi lumen, supporting sialylation. Although different works have focused on elucidating structure-function relationships in the hCST, this is the first study to address them in an alternatively spliced isoform. We have previously reported the expression of a functional human del177 isoform that has skipping of exon 6, resulting in a loss of 59 amino acids, without change in the open reading frame and conserving its C-terminal region. To elucidate structure-function relationships, we interrogated this isoform with a known pathogenic variant c.303C>T (p.Q101H) for the wt isoform, showing that its pathogenicity is significatively reduced in the mutated del177 isoform (del177Q101H). This is further explained by using a homology model based on previously reported mouse and maize crystal structures.

Monomer size effect in inelastic collisional dynamics of non-equilibrium soot nucleation

Molecular dynamics (MD) simulations of the collisional dynamics of the coronene-acepyrene and coronene radical-acepyrene pairs have been carried out to investigate the size effect of monomers of polycyclic aromatic hydrocarbons (PAH) on their non-equilibrium dimerization. The results compared to the previous MD simulations of the smaller pyrene-acepyrene and pyrenyl-acepyrene systems corroborate the non-equilibrium hypothesis of crosslinking PAH dimerization enhanced by physical interaction between the monomers. The phenomenon of inelastic collisional dynamics responsible for non-equilibrium van der Waals dimerization, which fosters a covalent bond formation between the monomers, amplifies with increasing PAH size. The increase in the size of the colliding monomers enhances the non-equilibrium effects as the growing pool of low-frequency modes provides a larger sink for the energy of the colliding PAH monomers. Based on the direct count of the crosslinking reaction events observed in the MD simulations, the forward rate constant for the coronene radical-acepyrene association is estimated at ∼10-11 cm3 molecule-1 s-1, showing a 15-fold increase with respect to the value from the statistical Rice-Ramsperger-Kassel-Marcus calculations. A comparison with the eightfold increase reported previously for the pyrenyl-acepyrene system shows that the statistical (equilibrium-based) calculations increasingly underestimate the reaction rate with the increasing size of the interacting PAHs from pyrene to coronene. The total increase of the MD-assessed rate constant for the coronene radical-acepyrene dimerization reaction as compared to pyrenyl-acepyrene is a factor of 2.4, with the overall collision efficiency to produce dimerized products growing by 30%.

Extracellular Peptide-Ligand Dimerization Actuator Receptor Design for Reversible and Spatially Dosed 3D Cell-Material Communication

Transmembrane receptors that endow mammalian cells with the ability to sense and respond to biomaterial-bound ligands will prove instrumental in bridging the fields of synthetic biology and biomaterials. Materials formed with thiol-norbornene chemistry are amenable to thiol-peptide patterning, and this study reports the rational design of synthetic receptors that reversibly activate cellular responses based on peptide-ligand recognition. This transmembrane receptor platform, termed Extracellular Peptide-ligand Dimerization Actuator (EPDA), consists of stimulatory or inhibitory receptor pairs that come together upon extracellular peptide dimer binding with corresponding monobody receptors. Intracellularly, Stimulatory EPDAs phosphorylate a substrate that merges two protein halves, whereas Inhibitory EPDAs revert split proteins back to their unmerged, inactive state via substrate dephosphorylation. To identify ligand-receptor pairs, over 2000 candidate monobodies were built in silico using PETEI, a novel computational algorithm we developed. The top 30 monobodies based on predicted peptide binding affinity were tested experimentally, and monobodies that induced the highest change in protein merging (green fluorescent protein, GFP) were incorporated in the final EPDA receptor design. In soluble form, stimulatory peptides induce intracellular GFP merging in a time- and concentration-dependent manner, and varying levels of green fluorescence were observed based on stimulatory and inhibitory peptide-ligand dosing. EPDA-programmed cells encapsulated in thiol-norbornene hydrogels patterned with stimulatory and inhibitory domains exhibited 3D activation or deactivation based on their location within peptide-patterned hydrogels. EPDA receptors can recognize a myriad of peptide-ligands bound to 3D materials, can reversibly induce cellular responses beyond fluorescence, and are widely applicable in biological research and regenerative medicine.

Identification of novel inhibitors of cancer target telomerase using a dual structure-based pharmacophore approach to virtually screen libraries, molecular docking and validation by molecular dynamics simulations

In about 85% of cancer malignancies, replicative immortality caused by increased telomerase activity makes it an attractive target for developing anticancer therapeutics. However, the lack of approved small-molecule inhibitors rooted in the structural ambiguity of telomerase has impeded drug development for decades. In this study, we have exploited the FVYL pocket in the thumb domain, which plays a key role in the enzyme's processivity. Due to the unavailability of a co-crystalized structure of BIBR1532 with the catalytic hTERT thumb domain, we utilized the molecular dynamics method to identify the precise binding site of the inhibitor. Two pharmacophore models were generated and validated for the putative (Site-I) and newly identified (Site-II) binding pockets which were screened virtually through the ChemDiv anticancer library, Otava drug-like green collection to identify novel lead compounds, and Binding database to screen out thumb domain-specific telomerase inhibitors. The top hits obtained were filtered using drug-likeliness parameters followed by redocking using a three-level screening strategy into their binding site. The structural investigation, molecular docking studies, and confirmatory molecular dynamics revealed that the exact binding site of BIBR1532 is away from the reported FVYL pocket with characteristic interactions conserved. Subsequently, the lead compounds with the highest docking scores and significant interactions in the newly discovered extended FVYL pocket were validated using 100 ns MD simulations. Additionally, cross-validated binding free energy calculations were performed using MM-PB(GB)SA methods followed by PCA and FEL characterization. The identified top lead compounds can be validated in vitro and taken forward for anticancer drug development.

Organophosphate Hydrolysis by a Designed Metalloenzyme: Impact of Mutations Explained

The efficient design of novel enzymes has been attainable only by a combination of theoretical approaches and experimental refinement, suggesting inadequate performance of de novo design protocols. Based on the analysis of the evolutionary trajectory of a designed organophosphate hydrolase, this work aimed at developing and validating the improved theoretical models describing the catalytic activity of five enzyme variants (including wild-type as well as theoretically derived and experimentally refined enzymes) performing the hydrolysis of diethyl 7-hydroxycoumarinyl phosphate. The following aspects possibly important for enzyme design were addressed: the level of theory sufficient for a reliable description of enzyme-reactant interactions, the issue of ground state (GS) destabilization versus transition state (TS) stabilization, and the derivation of the proper side chain rotamers of amino acid residues. For enzyme variants analyzed herein, differential transition state stabilization (DTSS, i.e., preferential TS binding by an enzyme over the GS binding) calculated with a non-empirical model of the interaction energy (i.e., multipole electrostatic plus approximate dispersion terms, MED) displayed a superior performance in ranking the enzyme catalytic activity. The MED DTSS-based systematic rotamer refinement performed with an efficient scanning procedure and accounting for long-range interaction energy terms is an important step capable of unlocking the full potential impact of the given residue that could be otherwise overlooked with a conventional static approach featuring optimization to the nearby minimum. While TS stabilization is the main factor contributing to the increased catalytic activity of the de novo-designed variant studied in this work, directed evolution refinement appears to impact the catalytic activity of another enzyme variant analyzed herein via GS destabilization.

Functional Roles of the Charged Residues of the C- and M-Gates in the Yeast Mitochondrial NAD+ Transporter Ndt1p

Mitochondrial carriers transport organic acids, amino acids, nucleotides and cofactors across the mitochondrial inner membrane. These transporters consist of a three-fold symmetric bundle of six transmembrane α-helices that encircle a pore with a central substrate binding site, whose alternating access is controlled by a cytoplasmic and a matrix gate (C- and M-gates). The C- and M-gates close by forming two different salt-bridge networks involving the conserved motifs [YF][DE]XX[KR] on the even-numbered and PX[DE]XX[KR] on the odd-numbered transmembrane α-helices, respectively. We have investigated the effects on transport of mutating the C-gate charged residues of the yeast NAD+ transporter Ndt1p and performed molecular docking with NAD+ and other substrates into structural models of Ndt1p. Double-cysteine substitutions and swapping the positions of the C-gate charged-pair residues showed that all of them contribute to the high transport rate of wild-type Ndt1p, although no single salt bridge is essential for activity. The in silico docking results strongly suggest that both the C-gate motif mutations and our previously reported M-gate mutations affect gate closing, whereas those of the M-gate also affect substrate binding, which is further supported by molecular dynamics. In particular, NAD+ most likely interferes with the cation-π interaction between R303-W198, which has been proposed to exist in the Ndt1p M-gate in the place of one of the salt bridges. These findings contribute to understanding the roles of the charged C- and M-gate residues in the transport mechanism of Ndt1p.

Structural mechanism of CB1R binding to peripheral and biased inverse agonists

The cannabinoid receptor 1 (CB1R) regulates synaptic transmission in the central nervous system, but also has important roles in the peripheral organs controlling cellular metabolism. While earlier generations of brain penetrant CB1R antagonists advanced to the clinic for their effective treatment of obesity, such molecules were ultimately shown to exhibit negative effects on central reward pathways that thwarted their further therapeutic development. The peripherally restricted CB1R inverse agonists MRI-1867 and MRI-1891 represent a new generation of compounds that retain the metabolic benefits of CB1R inhibitors while sparing the negative psychiatric effects. To understand the mechanism of binding and inhibition of CB1R by peripherally restricted antagonists, we developed a nanobody/fusion protein strategy for high-resolution cryo-EM structure determination of the GPCR inactive state, and used this method to determine structures of CB1R bound to either MRI-1867 or MRI-1891. These structures reveal how these compounds retain high affinity and specificity for CB1R's hydrophobic orthosteric site despite incorporating polar functionalities that lead to peripheral restriction. Further, the structure of the MRI-1891 complex along with accompanying molecular dynamics simulations shows how differential engagement with transmembrane helices and the proximal N-terminus can propagate through the receptor to contribute to biased inhibition of β-arrestin signaling.

Interior pH-sensing residue of human voltage-gated proton channel Hv1 is histidine 168

The molecular mechanisms governing the human voltage-gated proton channel hHv1 remain elusive. Here, we used membrane-enabled hybrid-solvent continuous constant pH molecular dynamics (CpHMD) simulations with pH replica exchange to further evaluate the structural models of hHv1 in the closed (hyperpolarized) and open (depolarized) states recently obtained with MD simulations and explore potential pH-sensing residues. The CpHMD titration at a set of symmetric pH conditions revealed three residues that can gain or lose protons upon channel depolarization. Among them, residue H168 at the intracellular end of the S3 helix switches from the deprotonated to the protonated state and its protonation is correlated with the increased tilting of the S3 helix during the transition from the closed to the open state. Thus, the simulation data suggest H168 as an interior pH sensor, in support of a recent finding based on electrophysiological experiments of Hv1 mutants. We propose that protonation of H168 acts as a key that unlocks the closed channel configuration by increasing the flexibility of the S2-S3 linker, which increases the tilt angle of S3 and enhances the mobility of the S4 helix, thus promoting channel opening. Our work represents an important step toward deciphering the pH-dependent gating mechanism of hHv1.

The C-terminal activating domain promotes pannexin 1 channel opening

Pannexin 1 (Panx1) constitutes a large pore channel responsible for the release of adenosine triphosphate (ATP) from apoptotic cells. Strong evidence indicates that caspase-mediated cleavage of the C-terminus promotes the opening of the Panx1 channel by unplugging the pore. However, this simple pore-plugging mechanism alone cannot account for the observation that a Panx1 construct ending before the caspase cleavage site remains closed. Here, we show that a helical region located immediately before the caspase cleavage site, referred to as the "C-terminal activating domain (CAD)", plays a pivotal role in facilitating Panx1 activation. Electrophysiology and mutagenesis studies uncovered that two conserved leucine residues within the CAD play a pivotal role. Cryoelectron microscopy (Cryo-EM) analysis of the construct ending before reaching the CAD demonstrated that the N terminus extends into an intracellular pocket. In contrast, the construct including the CAD revealed that this domain occupies the intracellular pocket, causing the N terminus to flip upward within the pore. Analysis of electrostatic free energy landscape in the closed conformation indicated that the intracellular side of the ion permeation pore may be occupied by anions like ATP, creating an electrostatic barrier for anions attempting to permeate the pore. When the N terminus flips up, it diminishes the positively charged surface, thereby reducing the drive to accumulate anions inside the pore. This dynamic change in the electrostatic landscape likely contributes to the selection of permeant ions. Collectively, these experiments put forth a mechanism in which C-terminal cleavage liberates the CAD, causing the repositioning of the N terminus to promote Panx1 channel opening.

A conserved H-bond network in human aquaporin-1 is necessary for native folding and oligomerization

Aquaporins (AQPs) are α-helical transmembrane proteins that conduct water through membranes with high selectivity and permeability. For human AQP1, in addition to the functional Asn-Pro-Ala motifs and the aromatic/Arg selectivity filter within the pore, there are several highly conserved residues that form an expansive hydrogen-bonding network. Previous solid-state nuclear magnetic resonance studies and structural conservation analysis have detailed which residues may be involved in this network. We explored this network by mutating the side chains or backbones involved in hydrogen-bonding, generating the following mutants: N127A, V133P, E142A, T187A, R195A, and S196A. The fold and stability of these mutants were assessed with attenuated total reflection Fourier transform infrared spectroscopy coupled with hydrogen/deuterium exchange upon increasing temperature. We found that replacement of any of the chosen residues to alanine leads to either partial instability or outright misfolding at room temperature, with the latter being most pronounced for the N127A, V133P, T187A, and R195A mutants. Deconvolution analysis of the amide I band revealed considerable secondary structure deviations, with some mutants exhibiting new random coil and β sheet structures. We also found that some of these mutations potentially disrupt the oligomerization of human AQP1. BN-PAGE and DLS data provide evidence toward the loss of tetramers within most of the mutants, meanwhile only the S196A mutant retains tetrameric organization. The molecular dynamics simulation of the wild-type, and the N127A, E142A, and T187A mutants show that these mutations result in major rearrangements of intra- and intermonomer hydrogen-bond networks. Overall, we show that specific point mutations that perturb hydrogen-bonding clusters result in severe misfolding in hAQP1 and disruption of its oligomerization. These data provide valuable insight into the structural stability of human aquaporin-1 and have implications toward other members of the AQP family, as these networks are largely conserved among a variety of human and nonmammalian AQP homologs.

Structure-based virtual screening of Trachyspermum ammi metabolites targeting acetylcholinesterase for Alzheimer's disease treatment

In recent decades, Alzheimer's disease (AD) has garnered significant attention due to its rapid global prevalence. The cholinergic hypothesis posits that the degradation of acetylcholine by acetylcholinesterase (AChE) contributes to AD development. Despite existing anti-AChE drugs, their adverse side effects necessitate new agents. This study analyzed 150 bioactive phytochemicals from Trachyspermum ammi using structure-based drug design and various in-silico tools to identify potent anti-AChE compounds. Compounds were screened for drug-likeness (QEDw ≥50%) and bioavailability (≥55%) and underwent toxicity profiling via the ProTox-II server. Selected compounds were prepared for molecular docking with the human AChE protein as the receptor. Viridifloral, 2-Methyl-3-glucosyloxy-5-isopropyl phenol, Alpha-Curcumene, and Sterol emerged as top candidates with high AChE affinity. These results were validated by molecular dynamics simulations, confirming stable interactions. The hit compounds were further evaluated for drug-likeness using Lipinski's rule and ADMET properties, confirming favorable pharmacokinetic profiles. DFT optimization analyzed frontier molecular orbitals and electrostatic potential, demonstrating favorable chemical reactivity and stability. This study suggests that these identified compounds could be novel nature-derived AChE inhibitors, potentially contributing to AD treatment. However, further in-vitro and in-vivo studies are necessary to confirm their efficacy in biological systems. Future research will focus on developing these compounds into safe and effective drugs to combat Alzheimer's disease.

Fc-FcγRI Complexes: Molecular Dynamics Simulations Shed Light on Ectodomain D3's Potential Role in IgG Binding

FcγRI plays a crucial role in the effector function of IgG antibodies, interacting with the lower hinge region of IgG1 with nanomolar affinity. Binding occurs specifically in domain 2 (D2) of the FcγRI ectodomain, while domain 3 (D3) is a flexible linker. The D3 domain is positioned away from the IgG binding site on the FcγRI and does not directly contact the Fc region. This study investigates the structural and functional properties of FcγRI D3 using 200 ns classical MD simulations of two models: (1) a full FcγRI ectodomain complex with Fc and (2) a truncated model excluding D3. Our findings suggest that the D3 ectodomain provides additional structural flexibility to the FcγRI-Fc complex without altering the C backbone motion or flexibility of the KHR binding motif in the FG loop. Critical residues involved in binding and contributing to complex stability were evaluated regarding changes in intramolecular interactions and destabilization tendency upon D3 truncation. Truncation did not significantly alter interactions around glycan-interacting residues in Fc chains or FcγRI-Fc binding interfaces. These findings provide valuable insights into the role of FcγRI D3 in modulating the structural dynamics of the FcγRI-Fc complex. While D3 does not directly contact Fc, its mobility and positioning may modulate the receptor's affinity, accessibility, and ability to bind IgG immune complexes. We suggest that a truncated FcγRI construct lacking the D3 domain may be a promising candidate for biosensor or capturing agents' development and optimization, offering improved performance in IgG capture assays without compromising critical binding interactions.

Crosstalk of Nucleic Acid Mimics with Lipid Membranes: A Multifaceted Computational and Experimental Study

Nucleic acid mimics (NAMs) have demonstrated high potential as antibacterial drugs. However, very few studies have assessed their possible diffusion across the bacterial envelope. In this work, we studied NAMs' diffusion in lipid bilayer systems that mimic the bacterial outer membrane using molecular dynamics (MD) simulations. Additionally, we examined the interactions of a NAM sequence with lipid membranes and ascertained the partition constants (Kp) through MD and spectroscopic investigations. The NAM sequences were composed of locked nucleic acid (LNA) and 2'-O-methyl (2'-OMe) residues, whereas the membrane models were composed of 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) or 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (POPG) phospholipids. The parametrization protocol followed was validated against literature data and demonstrated the reliability of our approach for simulating NAM sequences. Investigation into the interaction of the sequences with zwitterionic and anionic membranes revealed a preference for hydrogen bond formation with the anionic model over the zwitterionic one. Additionally, potential of mean force (PMF) calculations unveiled a lower free energy barrier for translocation across the zwitterionic bilayer model. Contrarily, the partition constants derived suggested a slightly higher partitioning within the anionic membrane, emphasizing a nuanced interplay of factors. Finally, spectroscopic partition measurements with liposomes presented challenges in quantifying the partition of NAMs due to minimal signal variations. However, a tendency for quenching in anionic vesicles suggested a potential, albeit small, partitioning effect that warrants further investigation. In summary, our study revealed that NAMs will not, in principle, be able to cross an intact bacterial outer membrane by passive diffusion.

Dissolution of denim waste in N-methyl morpholine-N-oxide monohydrate for fabrication of regenerated cellulosic film: Experimental and simulation study

Despite the significant amount of denim waste and its potential as a cellulose source, its use has been neglected. This study uses N-methyl morpholine-N-oxide, an eco-friendly solvent, to dissolve denim (including 100 % cotton) and create a denim film. Achieving a 10 % denim record solubility, a cellulosic film was also fabricated for comparison. Characterisation techniques were applied, and molecular dynamics simulations explored intramolecular interactions and the influence of indigo dye on dissolution process. FTIR spectra indicated no chemical reactions during dissolution and regeneration, though a shift in OH stretching suggested a change in crystallinity, confirmed by XRD results showing decreased crystallinity and a structural shift from cellulose I to cellulose II. 13C NMR analysis revealed disruptions in interchain hydrogen bonds after regeneration. TGA results showed lower decomposition temperatures for both films compared to the powders. Testing mechanical properties showed the denim film had higher elongation at break but lower tensile strength than the cellulose film. MD simulations indicated indigo dye did not significantly affect fundamental interactions but decreased denim solubility by reducing the diffusion coefficient. Rheological tests supported the simulation results, showing higher viscosity and molecular weight for the denim solution compared to cellulose.

In Silico Analysis of Triamterene as a Potential Dual Inhibitor of VEGFR-2 and c-Met Receptors

The vascular endothelial growth factor receptor 2 (VEGFR2) and the hepatocyte growth factor receptor (C-Met) are critical receptors for signaling pathways controlling crucial cellular processes such as cell growth, angiogenesis and tissue regeneration. However, dysregulation of these proteins has been reported in different diseases, particularly cancer, where these proteins promote tumour growth, invasiveness, metastasis and resistance to conventional therapies. The identification of dual inhibitors targeting both VEGFR-2 and c-Met has emerged as a strategic therapeutic approach to overcome the limitations and resistance mechanisms associated with single-target therapies in clinical settings. Through molecular dynamics simulations and comparative docking analysis, we tested the inhibitory potential of 2,016 Food and Drug Administration (FDA)-approved drugs targeting VEGFR-2 and/or c-Met receptors. The results revealed that entacapone and telmisartan are potent and selective inhibitors for c-Met and VEGFR-2, respectively. Interestingly, triamterene was identified as a promising dual inhibitor, demonstrating specific and significant binding affinity to both proteins. Molecular dynamics simulations revealed key interactions between the identified compounds and critical residues in the catalytic domains of both VEGFR-2 (e.g., Lys868, Asp1028, Asp1046) and c-Met (e.g., Asp1204, His1202, Asp1222), providing insights into their mechanism of action. These findings underscore the therapeutic potential of triamterene in targeting multiple signaling pathways involved in cancer progression, metastasis and poor prognosis in patients. Our study provides a foundational framework for the development of novel anticancer compounds able to target multiple pathways in cancer. Further preclinical and clinical investigations are needed to validate the efficacy of these compounds in clinical settings and to test their ability to overcome resistance and improve patient outcome.

Grand canonical Monte Carlo and deep learning assisted enhanced sampling to characterize the distribution of Mg2+ and influence of the Drude polarizable force field on the stability of folded states of the twister ribozyme

Molecular dynamics simulations are crucial for understanding the structural and dynamical behavior of biomolecular systems, including the impact of their environment. However, there is a gap between the time scale of these simulations and that of real-world experiments. To address this problem, various enhanced simulation methods have been developed. In addition, there has been a significant advancement of the force fields used for simulations associated with the explicit treatment of electronic polarizability. In this study, we apply oscillating chemical potential grand canonical Monte Carlo and machine learning methods to determine reaction coordinates combined with metadynamics simulations to explore the role of Mg2+ distribution and electronic polarizability in the context of the classical Drude oscillator polarizable force field on the stability of the twister ribozyme. The introduction of electronic polarizability along with the details of the distribution of Mg2+ significantly stabilizes the simulations with respect to sampling the crystallographic conformation. The introduction of electronic polarizability leads to increased stability over that obtained with the additive CHARMM36 FF reported in a previous study, allowing for a distribution of a wider range of ions to stabilize twister. Specific interactions contributing to stabilization are identified, including both those observed in the crystal structures and additional experimentally unobserved interactions. Interactions of Mg2+ with the bases are indicated to make important contributions to stabilization. Notably, the presence of specific interactions between the Mg2+ ions and bases or the non-bridging phosphate oxygens (NBPOs) leads to enhanced dipole moments of all three moieties. Mg2+-NBPO interactions led to enhanced dipoles of the phosphates but, interestingly, not in all the participating ions. The present results further indicate the importance of electronic polarizability in stabilizing RNA in molecular simulations and the complicated nature of the relationship of Mg2+-RNA interactions with the polarization response of the bases and phosphates.

Balancing Group 1 Monoatomic Ion-Polar Compound Interactions in the Polarizable Drude Force Field: Application in Protein and Nucleic Acid Systems

An accurate force field (FF) is the foundation of reliable results from molecular dynamics (MD) simulations. In our recently published work, we developed a protocol to generate atom pair-specific Lennard-Jones (known as NBFIX in CHARMM) and through-space Thole dipole screening (NBTHOLE) parameters in the context of the Drude polarizable FF based on readily accessible quantum mechanical (QM) data to fit condensed phase experimental thermodynamic benchmarks, including the osmotic pressure, diffusion coefficient, ionic conductivity, and solvation free energy, when available. In the present work, the developed protocol is applied to generate NBFIX and NBTHOLE parameters for interactions between monatomic ions (specifically Li+, Na+, K+, Rb+, Cs+, and Cl-) and common functional groups found in proteins and nucleic acids. The parameters generated for each ion-functional group pair were then applied to the corresponding functional groups within proteins or nucleic acids followed by MD simulations to analyze the distribution of ions around these biomolecules. The modified FF successfully addresses the issue of overbinding observed in a previous iteration of the Drude FF. Quantitatively, the model accurately reproduces the effective charge of proteins and demonstrates a level of charge neutralization for a double-helix B-DNA in good agreement with the counterion condensation theory. Additionally, simulations involving ion competition correlate well with experimental results, following the trend Li+ > Na+ ≈ K+ > Rb+. These results validate the refined model for group 1 ion-biomolecule interactions that will facilitate the application of the polarizable Drude FF in systems in which group 1 ions play an important role.

Accurate Predictions of Molecular Properties of Proteins via Graph Neural Networks and Transfer Learning

Machine learning has emerged as a promising approach for predicting molecular properties of proteins, as it addresses limitations of experimental and traditional computational methods. Here, we introduce GSnet, a graph neural network (GNN) trained to predict physicochemical and geometric properties including solvation free energies, diffusion constants, and hydrodynamic radii, based on three-dimensional protein structures. By leveraging transfer learning, pre-trained GSnet embeddings were adapted to predict solvent-accessible surface area (SASA) and residue-specific pKa values, achieving high accuracy and generalizability. Notably, GSnet outperformed existing protein embeddings for SASA prediction, and a locally charge-aware variant, aLCnet, approached the accuracy of simulation-based and empirical methods for pKa prediction. Our GNN framework demonstrated robustness across diverse datasets, including intrinsically disordered peptides, and scalability for high-throughput applications. These results highlight the potential of GNN-based embeddings and transfer learning to advance protein structure analysis, providing a foundation for integrating predictive models into proteome-wide studies and structural biology pipelines.

Elucidating the mechanisms of a herbal compound fumarine and its modulation on the estrogen receptor 1

Stroke-related numbness and weakness (SRNW) are resultant disabilities following a stroke episode and may present with muscle weakness, numbness, tightness, spasticity, and pain in up to 85% of patients. Huangqi Guizhi Wuwu Decoction (HGWD) has been widely investigated to manage the sensorimotor deficiencies at the herb and formula level. However, detailed molecular mechanisms of its constituents are presently lacking. This project employed computational molecular modelling and docking methods to identify candidate compounds of HGWD which may serve as effective modulators of target proteins involved in SRNW. Estrogen Receptor 1 was identified as a promising target for HGWD compounds, while the herbal compound fumarine, a constituent of Jujubae Fructus, was predicted to exhibit high binding affinity and favourable ligand-receptor interactions with ESR1. There is currently a lack of scientific evidence for specific atomic-level interactions between ESR1 and this compound. Therefore, molecular docking and molecular dynamics simulations were used to elucidate the interaction mechanisms of fumarine with ESR1; and the molecular-level structural and functional consequences of ligand binding. Ligand-receptor contact analysis and free energy decomposition calculations identified Glu419 and Leu38 as stable hydrogen bond partners, while favourable contributions to the binding free energy include in Met421 (-10.74 kJ/mol) and Leu525 (-10.02 kJ/mol). This work provides the basis for further studies on discovering lead compounds which modulate the activity of ESR1.

Accurate in silico predictions of modified RNA interactions to a prototypical RNA-binding protein with λ-dynamics

RNA-binding proteins shape biology through their widespread functions in RNA biochemistry. Their function requires the recognition of specific RNA motifs for targeted binding. These RNA binding elements can be composed of both unmodified and chemically modified RNAs, of which over 170 chemical modifications have been identified in biology. Unmodified RNA sequence preferences for RNA-binding proteins have been widely studied, with numerous methods available to identify their preferred sequence motifs. However, only a few techniques can detect preferred RNA modifications, and no current method can comprehensively screen the vast array of hundreds of natural RNA modifications. Prior work demonstrated that λ-dynamics is an accurate in silico method to predict RNA base binding preferences of an RNA-binding antibody. This work extends that effort by using λ-dynamics to predict unmodified and modified RNA binding preferences of human Pumilio, a prototypical RNA binding protein. A library of RNA modifications was screened at eight nucleotide positions along the RNA to identify modifications predicted to affect Pumilio binding. Computed binding affinities were compared with experimental data to reveal high predictive accuracy. In silico force field accuracies were also evaluated between CHARMM and Amber RNA force fields to determine the best parameter set to use in binding calculations. This work demonstrates that λ-dynamics can predict RNA interactions to a bona fide RNA-binding protein without the requirements of chemical reagents or new methods to experimentally test binding at the bench. Advancing in silico methods like λ-dynamics will unlock new frontiers in understanding how RNA modifications shape RNA biochemistry.

Dynamic interplay between a TonB-dependent heme transporter and a TonB protein in a membrane environment

The envelope of Gram-negative bacteria is composed of two membranes separated by the periplasmic space. This organization imposes geometrical and distance constraints that are key for the mechanism of action of multicomponent systems spanning the envelope. However, consideration of all three compartments by experimental approaches is still elusive. Here, we have used the state-of-the-art molecular dynamics simulation in an Escherichia coli envelope model to obtain a dynamic view of molecular interactions between the outer membrane heme transporter HasR and the inner membrane TonB-like protein HasB. Their interaction allows the transfer of the inner membrane proton-motive force derived energy to the transporter for heme internalization. The simulations that incorporate both membranes show the key role of periplasmic domains of both proteins and their dynamics in complex formation and stability. They revealed a previously unidentified network of HasR-HasB protein-protein interactions in the periplasm. Experimental validation (mutations, in vivo phenotypic and biophysical assays) provides support for the simulation-predicted interactions. Based on structural and sequence conservation, the network of interaction revealed in this study is expected to occur in other nutrient import systems.

IMPORTANCE

Gram-negative bacteria import scarce nutrients such as metals and vitamins by an energized mechanism involving a multicomponent protein system that spans the cell envelope. It consists of an outer membrane TonB-dependent transporter (TBDT) and a TonB complex in the inner membrane that provides the proton motive force energy for the nutrient entry. Despite the intense research efforts focused on this system (a) from structural and fundamental microbiology perspectives and (b) for the interest in the development of new antibacterial strategies, the molecular mechanism of the system is not at all well understood. The lack of understanding comes from incomplete structural data and the experimental difficulties of studying an inherently flexible multicomponent complex that resides within the heterogeneous environment of the double membrane bacterial cell envelope. To address these challenges and obtain a comprehensive view of the molecular interactions at atomic level, here, we have used the combined power of advanced molecular simulations and complementary microbiology and biochemical experiments. Our results represent a significant step forward in understanding the structural and molecular bases of this vital mechanism.

Respiratory complex II acting as a homeostatic regulatory sensor

The succinate-ubiquinone oxidoreductase (SQR) complex connects two of the cell's most vital energy-producing metabolic processes: the tricarboxylic acid cycle and the electron transport chain. Hence, the SQR complex is essential in cell metabolism, and its malfunction leads to the progression of multiple metabolic disorders and other diseases, such as cancer. In the current study, we calculated the electron tunneling (ET) pathways between the different redox systems in the SQR complex, including the SQR ligands and the distant heme b redox center, using the broken-symmetry semi-empirical ZINDO method. Interestingly, we discovered a water channel running from the mitochondrial matrix, filling the space between Fe3S4 and heme b redox centers. To investigate the physiological function of the water channel, we performed extensive molecular dynamics (MD) simulations of the membrane-embedded SQR complex in small and large water boxes, representing regular (MDA) and extended (MDB) volume states, respectively. We found that under regular volume conditions (MDA), the ET reaction is conducted through both the iron-sulfur cluster chain (i.e., pathway A) and through heme b (i.e., pathway B). Hence, the SQR complex encompasses an internal interferometer similar to the Mach-Zender interferometer, such that the tunneling electron experiences a self-interference effect through pathways A and B, enhancing the SQR complex's overall ET thermodynamics and favoring the forward ET direction of oxidizing succinate to fumarate and reducing ubiquinone to ubiquinol. On the other hand, we found that under extended volume conditions (MDB), the internal water channel of the SQR complex "senses" the expansion in the mitochondrial volume, pushing the heme b and Fe4S3 redox centers apart and hence lowering the SQR equilibrium constant to almost unity. Therefore, the SQR complex could be driven to work in the reverse direction, catalyzing the production of ubiquinone molecules essential for the physiological function of respiratory complexes I and III and restoring the inner-mitochondrial membrane potential, which leads to restoring the function of the H-K anti-porter, pumping K+ outward from the matrix and restoring the regular mitochondrial volume.

Mechanism of Dimer Selectivity and Binding Cooperativity of BRAF Inhibitors

Aberrant signaling of BRAFV600E is a major cancer driver. Current FDA-approved RAF inhibitors selectively inhibit the monomeric BRAFV600E and suffer from tumor resistance. Recently, dimer-selective and equipotent RAF inhibitors have been developed; however, the mechanism of dimer selectivity is poorly understood. Here, we report extensive molecular dynamics (MD) simulations of the monomeric and dimeric BRAFV600E in the apo form or in complex with one or two dimer-selective (PHI1) or equipotent (LY3009120) inhibitor(s). The simulations uncovered the unprecedented details of the remarkable allostery in BRAFV600E dimerization and inhibitor binding. Specifically, dimerization retrains and shifts the αC helix inward and increases the flexibility of the DFG motif; dimer compatibility is due to the promotion of the αC-in conformation, which is stabilized by a hydrogen bond formation between the inhibitor and the αC Glu501. A more stable hydrogen bond further restrains and shifts the αC helix inward, which incurs a larger entropic penalty that disfavors monomer binding. This mechanism led us to propose an empirical way based on the co-crystal structure to assess the dimer selectivity of a BRAFV600E inhibitor. Simulations also revealed that the positive cooperativity of PHI1 is due to its ability to preorganize the αC and DFG conformation in the opposite protomer, priming it for binding the second inhibitor. The atomically detailed view of the interplay between BRAF dimerization and inhibitor allostery as well as cooperativity has implications for understanding kinase signaling and contributes to the design of protomer selective RAF inhibitors.

Hemoglobin targeting potential of aminocarb pesticide: Investigation into dynamics, conformational stability, and energetics in solvent environment

Aminocarb (AMC), a carbamate pesticide, due to its prevalent usage exhibits increased accumulation in the environment affecting both insects and humans. It enters the human body via food grains and be transported through bloodstream. AMC's chemical structure, containing specific molecular frameworks and functional groups, enables it to bind with proteins like albumin and hemoglobin. Given that molecules with similar architecture are known to bind with hemoglobin, we aimed to explore Aminocarb's binding capability and the potential mechanism or mode of its interaction with hemoglobin. Hb being a tetramer with a profound interface between amino acid chains offers multiple binding sites. It is therefore important to investigate the structural aspects of binding of AMC by employing various spectroscopic and in-silico methods. The surface of the α1 chain near the α1β2 interface emerges as the preferred binding site for AMC, primarily due to its conformational restrictions. In its bound state, AMC tends to maintain a relaxed conformation, closely resembling its globally optimized geometry, and resides in close proximity to the α1 chain via multiple hydrophobic contacts and water bridge as observed in molecular dynamics (MD) simulations. Fluorescence quenching experiments showed moderate binding strength (7.7 × 10⁴ L M⁻1 at 288 K, 7.8 × 10⁴ L M⁻1 at 298 K, 7.9 × 10⁴ L M⁻1 at 308 K) and spontaneous binding, driven by hydrophobic and van der Waals interactions, as indicated by enthalpy (0.80-0.91 kJ mol⁻1), entropy (0.0970-0.0974 kJ mol⁻1), and Gibbs free energy (-27.13 to - 29.08 kJ mol⁻1). Circular dichroism experiments reveal no major structural changes in Hb. Quantum chemical calculations and MD simulations reveal conformation-dependent energy differences, enhancing our understanding of AMC's binding mechanism to Hb.

Integrated study of Quercetin as a potent SARS-CoV-2 RdRp inhibitor: Binding interactions, MD simulations, and In vitro assays

To find an effective inhibitor for SARS-CoV-2, Quercetin's chemical structure was compared to nine ligands associated with nine key SARS-CoV-2 proteins. It was found that Quercetin closely resembles Remdesivir, the co-crystallized ligand of RNA-dependent RNA polymerase (RdRp). This similarity was confirmed through flexible alignment experiments and molecular docking studies, which showed that both Quercetin and Remdesivir bind similarly to the active site of RdRp. Molecular dynamics (MD) simulations over a 200 ns trajectory, analyzing various factors like RMSD, RG, RMSF, SASA, and hydrogen bonding were conducted. These simulations gave detailed insights into the binding interactions of Quercetin with RdRp compared to Remdesivir. Further analyses, including MM-GBSA, Protein-Ligand Interaction Fingerprints (ProLIF) and Profile PLIP studies, confirmed the stability of Quercetin's binding. Principal component analysis of trajectories (PCAT) provided insights into the coordinated movements within the systems studied. In vitro assays showed that Quercetin is highly effective in inhibiting RdRp, with an IC50 of 122.1 ±5.46 nM, which is better than Remdesivir's IC50 of 21.62 ±2.81 μM. Moreover, Quercetin showed greater efficacy against SARS-CoV-2 In vitro, with an IC50 of 1.149 μg/ml compared to Remdesivir's 9.54 μg/ml. The selectivity index (SI) values highlighted Quercetin's safety margin (SI: 791) over Remdesivir (SI: 6). In conclusion, our comprehensive study suggests that Quercetin is a promising candidate for further research as an inhibitor of SARS-CoV-2 RdRp, providing valuable insights for developing an effective anti-COVID-19 treatment.

Phonons reveal coupled cholesterol-lipid dynamics in ternary membranes

Experimental studies of collective dynamics in lipid bilayers have been challenging due to the energy resolution required to observe these low-energy phonon-like modes. However, inelastic x-ray scattering (IXS) measurements-a technique for probing vibrations in soft and biological materials-are now possible with sub-meV resolution, permitting direct observation of low-energy, phonon-like modes in lipid membranes. Here, IXS measurements with sub-meV energy resolution reveal a low-energy optic-like phonon mode at roughly 3 meV in the liquid-ordered (Lo) and liquid-disordered phases of a ternary lipid mixture. This mode is only observed experimentally at momentum transfers greater than 5 nm-1 in the Lo system. A similar gapped mode is also observed in all-atom molecular dynamics (MD) simulations of the same mixture, indicating that the simulations accurately represent the fast, collective dynamics in the Lo phase. Its optical nature and the Q range of the gap together suggest that the observed mode is due to the coupled motion of cholesterol-lipid pairs, separated by several hydrocarbon chains within the membrane plane. Analysis of the simulations provides molecular insight into the origin of the mode in transient, nanoscale substructures of hexagonally packed hydrocarbon chains. This nanoscale hexagonal packing was previously reported based on MD simulations and, later, by NMR measurements. Here, however, the integration of IXS and MD simulations identifies a new signature of the Lo substructure in the collective lipid dynamics, thanks to the recent confluence of IXS sensitivity and MD simulation capabilities.

Spc2 modulates substrate- and cleavage site-selection in the yeast signal peptidase complex

Secretory proteins are critically dependent on the correct processing of their signal sequence by the signal peptidase complex (SPC). This step, which is essential for the proper folding and localization of proteins in eukaryotic cells, is still not fully understood. In eukaryotes, the SPC comprises four evolutionarily conserved membrane subunits (Spc1-3 and Sec11). Here, we investigated the role of Spc2, examining SPC cleavage efficiency on various models and natural signal sequences in yeast cells depleted of or with mutations in Spc2. Our data show that discrimination between substrates and identification of the cleavage site by SPC is compromised when Spc2 is absent or mutated. Molecular dynamics simulation of the yeast SPC AlphaFold2-Multimer model indicates that membrane thinning at the center of SPC is reduced without Spc2, suggesting a molecular explanation for the altered substrate recognition properties of SPC lacking Spc2. These results provide new insights into the molecular mechanisms by which SPC governs protein biogenesis.

Rational design of monomeric IL37 variants guided by stability and dynamical analyses of IL37 dimers

IL37 plays important roles in the regulation of innate immunity and its oligomeric status is critical to these roles. In its monomeric state, IL37 can effectively inhibit the inflammatory response of IL18 by binding to IL18Rα, a capacity lost in its dimeric form, underlining the pivotal role of the oligomeric status of IL37 in its anti-inflammatory action. Until now, two IL37 dimer structures have been deposited in PDB, reflecting a substantial difference in their dimer interfaces. Given this discrepancy, we analyzed the PDB structures of the IL37 dimer (PDB IDs: 6ncu, 5hn1) along with a AF2-multimer prediction by molecular dynamics (MD) simulations. Results showed that the 5hn1 and AF2-predicted dimers have the same interface and stably maintained their conformations throughout simulations, while the recent IL37 dimer (PDB ID: 6ncu) with a different interface did not, proposing a possible issue with the recent IL37 dimer structure (6ncu). Next, focusing on the stable dimer structures, we have identified five critical positions of V71/Y85/I86/E89/S114, three new positions compared to the literature, that would reduce dimer stability without affecting the monomer structure. Two quintuple mutants were tested by MD simulations and showed partial or complete dissociation of the dimer. Overall, the insights gained from this study reinforce the validity of the 5hn1 and AF2 multimer structures, while also advancing our understanding of the IL37 dimer interface through the generation of monomer-locked IL37 variants.

Discovering new mode-of-action pesticide leads inhibiting protein-protein interactions: example targeting plant O-acetylserine sulfhydrylase

BACKGROUND

The widespread evolution of pesticide resistance poses a significant challenge to current agriculture, necessitating the discovery of molecules with new modes of action. Despite extensive efforts, no major molecules with new modes of action have been commercialized for decades. Most pesticides function by binding to specific pockets on target enzymes, enabling a single target site mutation to confer resistance. An alternative approach is the disruption of protein-protein interactions (PPI), which require complementary mutations on both interacting partners for resistance to occur. Thus, our aim is the discovery and design of small-molecule inhibitors that target the interface of the PPI complex of O-acetylserine sulfhydrylase (OASS) and serine acetyltransferase (SAT), key obligatory interacting plant enzymes involved in the biosynthesis of the amino acid cysteine.

RESULTS

By employing in silico filtering techniques on a virtual library of 30 million small molecules, we identified initial hits capable of binding OASS and interfering with its interaction with a peptide derived from SAT with a half-maximal inhibitory concentration (IC50) of 34 μm. Subsequently, we conducted molecular chemical optimizations, generating an early lead molecule (PJ4) with an IC50 value of 4 μm. PJ4 successfully inhibited the germination of Arabidopsis thaliana seedlings and inhibited clover growth in a pre-emergence application at an effective concentration of 4.6 kg ha-1.

CONCLUSION

These new compounds described herein can serve as promising leads for further optimization as herbicides with a new mode-of-action. This technology can be used for discovering new modes of action chemicals inhibiting all pest groups. © 2024 The Author(s). Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Design, synthesis and evaluation of 3-(2-(substituted benzyloxy)benzylidene) pyrrolidine-2,5-dione derivatives for novel ATX inhibitor

Autotaxin (ATX) has emerged as a promising therapeutic target for liver diseases. In this study, we identified potential drug candidates through in silico high-throughput screening. Subsequently, we synthesized a series of small molecules, specifically KR-40795 (2c), a pyrrolidine-2,5-dione-based analogue that binds to the allosteric tunnel and hydrophobic pocket of ATX. This compound was designed to inhibit the enzymatic activity of ATX for the treatment of liver diseases. The inhibitory potency of KR-40795 was evaluated using a biochemical assay that measured the hydrolysis of a specific substrate (FS-3). Notably, KR-40795 demonstrated significant inhibition of both collagen formation and lipid accumulation in liver cells, suggesting its potential as a therapeutic agent for liver diseases, particularly fibrosis and steatosis.

Catalytic mechanism, computational design, and crystal structure of a highly specific and efficient benzoylecgonine hydrolase

Enzyme therapy for cocaine detoxification should break down both cocaine and its primary toxic metabolite, benzoylecgonine (BZE), which is also the main form of cocaine contaminant in the environment. An ideal BZE-metabolizing enzyme (BZEase) is expected to be highly efficient and selective in BZE hydrolysis. Here, BZEase4 was engineered from bacterial cocaine esterase (CocE) by our reactant state-based enzyme design theories (RED), which has a 34,977-fold improved substrate discrimination between BZE and the neurotransmitter acetylcholine (ACh), compared with wild-type CocE. Under the physiological concentrations of BZE and ACh, the reaction velocity of BZEase4 against BZE is 2.25 × 106-fold higher than it against ACh, suggesting BZEase4 has extremely high substrate selectivity for BZE over ACh to minimize the potential cholinergic side-effects. This study provides additional evidence supporting the further development of BZEase4 toward a promising therapeutic for cocaine overdose, a potentially effective and eco-friendly enzymatic method for BZE degradation in the environment.

Challenging the Biginelli scaffold to surpass the first line antitubercular drugs: Mycobacterium tuberculosis thymidine monophosphate kinase (TMPKmt) inhibition activity and molecular modelling studies

New Biginelli adducts were rationalised, via the introduction of selected anti-tubercular (TB) pharmacophores into the dihydropyrimidine (DHPM) ring of deoxythymidine monophosphate (dTMP), the natural substrate of Mycobacterium tuberculosis thymidine monophosphate kinase (TMPKmt). Repurposing was one of the design rationale strategies for some selected mimics of the designed compounds. The anti-TB activity was screened against the Mtb H37Rv strain where 11a was superior to ethambutol (EMB), and was 9-fold more potent than pyrazinamide (PZA). Additionally, compounds 11b, 4a, 4b, 13a, 13b and 14a elicited higher anti-TB activity than PZA, showing better safety profiles than EMB against RAW 264.7 cells' growth. The in vitro TMPKmt inhibition assay released compounds 11a, 11b and 13b as the most potent inhibitors. Docking studies presumed the binding modes and molecular dynamics (MD) simulation revealed the dynamic stability of 11a-TMPKmt complex over 100 ns. In silico prediction of the chemo-informatics properties of the most active compounds was conducted.

Membrane tension-dependent conformational change of Isoleucine 106 of loop B diminishes water permeability in FaPIP2;1

Aquaporins (AQPs) are membrane proteins specialized in facilitating water transport across membranes. Mechanical stress is one of the various stimuli that regulate AQPs. Briefly, there are several studies that report a decrease in permeability upon an increase in membrane tension. However, the molecular details of this mechanosensitive (MS) response are still a matter of debate. Our work attempts to close that gap in knowledge by providing evidence of a conformational change that occurs inside the pore of the strawberry aquaporin FaPIP2;1. Via osmotic shock experiments and molecular dynamics (MD) simulations, we found that a residue of loop B, I106, is key to the blocking of the permeation pathway and such a change is almost exclusively found under membrane tensile stress. In detail, osmotic shock experiments exhibited a nonlinear increment in water fluxes for increasing osmolarities, evidencing a decrease in the FaPIP2;1 permeability. MD simulations under membrane tension showed the same trend, with a significant increase in states with a low water permeability. The latter was correlated with a conformational change in I106 that generates a permeation barrier of around 18 kJ mol-1, effectively closing the pore. This work constitutes the first report of a PIP type aquaporin reacting to tensile stress in the membrane. Our findings could pave the way to test whether this conformational change is also responsible for mechanical gating in the other MS aquaporins, both those already reported and those still waiting to be found.

Integrated all-atom and coarse-grained simulations uncover structural, dynamics and energetic shifts in SARS-CoV-2 JN.1 and BA.2.86 variants

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for the COVID-19 pandemic, is an enveloped, positive-stranded RNA virus that enters human cells by using its spike protein to bind to the human angiotensin-converting enzyme 2 (ACE2) receptor. Since its emergence, the virus has mutated, producing variants with increased transmissibility, immune evasion, and infectivity. The JN.1 variant, detected in January 2024, features a single substitution mutation (Leu455Ser) in the receptor-binding domain (RBD) of its spike protein, setting it apart from its parent lineage, BA.2.86. This variant has rapidly become globally predominant due to its enhanced transmission and significant epidemiological impact. To understand the causes behind the dominance of the JN.1 variant, we conducted a comprehensive study using all-atom molecular dynamics (MD) and coarse-grained MD simulations. This allowed us to examine the structural, dynamic, energetics and binding properties of the wild-type (Wuhan strain), BA.2.86, and JN.1 variants. Principal component and free energy landscape analyses revealed enhanced structural stability in the JN.1 variant. Molecular Mechanics Poisson-Boltzmann Surface Area (MM/PBSA) assessments indicated lower binding affinity for JN.1 as compared to BA.2.86. Intermolecular interaction analyses further confirmed BA.2.86's superior binding affinity over JN.1 and wild-type. Additionally, we compared and validated our findings against experimentally determined cryo-electron microscopy (cryo-EM) structures of JN.1 and BA.2.86 variants, confirming the reliability of our simulation results. Overall, this study provides crucial insights into the structural-dynamics-energetics features and physicochemical properties that have contributed to the global prevalence of the JN.1 variant and sheds light on its potential to generate future subvariants.

Osmolyte-induced protein stability changes explained by graph theory

Enhanced stabilization of protein structures via the presence of inert osmolytes is a key mechanism adopted both by physiological systems and in biotechnological applications. While the intrinsic stability of proteins is ultimately fixed by their amino acid composition and organization, the interactions between osmolytes and proteins together with their concentrations introduce an additional layer of complexity and in turn, a method of modulating protein stability. Here, we combined experimental measurements with molecular dynamics simulations and graph-theory-based analyses to predict the stabilizing/destabilizing effects of different kinds of osmolytes on proteins during heat-mediated denaturation. We found that (i) proteins in solution with stability-enhancing osmolytes tend to have more compact interaction networks than those assumed in the presence of destabilizing osmolytes; (ii) a strong negative correlation (R = -0.85) characterizes the relationship between the melting temperatureT m and the preferential interaction coefficient defined by the radial distribution functions of osmolytes and water around the protein and (iii) a positive correlation exists between osmolyte-osmolyte clustering and the extent of preferential exclusion from the local domain of the protein, suggesting that exclusion may be driven by enhanced steric hindrance of aggregated osmolytes.

Why is binding of a divalent metal cation to a structural motif containing four carboxylate residues not accompanied by a conformational change?

We earlier showed that Torpedo californica acetylcholinesterase (AChE) contains a cluster of four conserved aspartates that can strongly bind divalent cations, which we named the 4D motif. Binding of the divalent metal cations greatly increases its thermal stability. Here we systematically examined all available crystallographic structures of T. californica AChE. Two additional metal-binding sites were identified, both composed of acidic and histidine residues. Relative binding to the 4D and additional sites was studied using metadynamics simulations. It was observed that in crystal structures devoid of metal ions in the 4D site, the conformation of T. californica AChE is almost identical to that in structures in which it is occupied by a divalent metal ion. Closer examination of the 4D motif reveals that three of the four acidic residues form ion pairs with conserved basic residues surrounding them. We named this new motif the 4A/3B motif. Molecular dynamics with quantum potential simulations was used to quantify the 4D motif's binding strength compared with that of the metal-binding site in the protein fXIIIa, which consists of four aspartates, but is devoid of adjacent cationic residues. Whereas fXIIIa's 4D site, in the absence of a metal cation, expanded significantly in the simulation, that of Torpedo AChE displayed only minor periodic changes in size. Furthermore, the energy of metal ion unbinding from the two sites differs by ca. 10 kcal/mol. We identified several other proteins in the PDB that contain the 4A/3B motif, whose conformations are identical in the presence or absence of a metal ion. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at https://proteopedia.org/w/Journal:Protein_Science:4.

Integrated computational characterization of valosin-containing protein double-psi β-barrel domain: Insights into structural stability, binding mechanisms, and evolutionary significance

Valosin-containing protein (VCP) plays a crucial role in various cellular processes, yet the molecular mechanisms and structural dynamics of its double-psi β-barrel (DPBB) domain, particularly in human, remain insufficiently explored. While previous studies have characterized the VCP_DPBB domain in other organisms, such as thermoplasma acidophilum and methanopyrus kandleri, its evolutionary conservation, binding potential, and stability in human require further investigation. To address this gap, we first employed all-atom molecular dynamics (AAMD) simulations to examine the structural dynamics of the human VCP_DPBB domain. We also assessed its amino acid interaction energies, stability, folding enthalpy, evolutionary conservation, solubility, and crystallizability using various computational frameworks. Additionally, to uncover the plausible biological function, protein-peptide docking was performed to evaluate the interactions between the DPBB domain and the C-terminal gp78 peptide of the E3 ubiquitin ligase. Further, AAMD and coarse-grained molecular dynamics (CGMD) simulations explored the binding preferences, fluctuations, and stability of human VCP_DPBB-gp78 complexes. Our findings indicate that, while thermoplasma acidophilum VCP_DPBB-gp78 showed stronger initial binding, the human VCP_DPBB-gp78 complex exhibited superior stability, binding affinity, and more stabilizing interactions. This integrated analysis provides valuable insights into the evolutionary significance and functionality of the DPBB domain, with potential therapeutic implications for VCP-related diseases.

Leukotriene B4 receptor 1 (BLT1) activation by leukotriene B4 (LTB4) and E resolvins (RvE1 and RvE2)

Leukotriene B4 (LTB4) is a lipid inflammatory mediator derived from arachidonic acid (AA). Leukotriene B4 receptor 1 (BLT1), a G protein-coupled receptor (GPCR), is a receptor of LTB4. Nonetheless, the resolution of inflammation is driven by specialized pro-resolving lipid mediators (SPMs) such as resolvins E1 (RvE1) and E2 (RvE2). Both resolvins are derived from omega-3 fatty acid eicosapentaenoic acid (EPA). Here, long-term molecular dynamics simulations (MD) were performed to investigate the activation of the BLT1 receptor using two pro-resolution agonists (RvE1 and RvE2) and an inflammatory agonist (LTB4). We have analyzed the receptor's activation state, electrostatic interactions, and the binding affinity the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach. The results showed that LTB4 and RvE1 have kept the receptor in an active state by higher simulation time. MD showed that the ligand-receptor interactions occurred mainly through residues H94, R156, and R267. The MMPBSA calculations showed residues R156 and R267 were the two mainly hotspots. Our MMPBSA results were compatible with experimental results from other studies. Overall, the results from this study provide new insights into the activation mechanisms of the BLT1 receptor, reinforcing the role of critical residues and interactions in the binding of pro-resolution and pro-inflammatory agonists.

Controlling the trans effect induced by nitric oxide and carbon monoxide: H93C myoglobin versus H-NOX sensors and soluble guanylate cyclase

Myoglobin (Mb) has been engineered to replace the proximal histidine (His93) with a cysteine in order to investigate the trans effect induced by diatomic ligands using time-resolved electronic absorption spectroscopy. This single mutation induces a change of heme coordination state and bonding character which change carbon monoxide (CO) and nitric oxide (NO) dynamics. In H93C Mb the increased Fe2+-S distance weakens this bond which is replaced with a distal Fe2+-His64 ligation. We measured dynamics very different from wild type Mb but similar with those measured in soluble guanylate cyclase (sGC). Whereas NO induces a direct negative trans effect, the strain on His64 ligation is sufficient to counteract the positive trans effect due to CO. After photodissociation, geminate recombination of NO to the transient 4-coordinate heme of H93C occurred with a fast time constant (6.9 ps) identical to that in sGC. Remarkably, we also observed picosecond geminate rebinding of CO to H93C Mb, similarly with sGC in the simultaneous presence of CO and an allosteric stimulator. This CO rebinding dynamics to the 4c-heme in H93C Mb was never measured in other Mb mutants and demonstrates the existence of 5-coordinate heme with CO, explaining the synergistic activation of sGC in presence of CO and a stimulator.