Search

Pim-1 kinase, a serine/threonine kinase, is often overexpressed in various cancers, contributing to disease progression and poor prognosis. In this study, we explored the potential of flavonoids as inhibitors of Pim-1 kinase using a combination of molecular docking and steered molecular dynamics (SMD) simulations. Our docking studies revealed two main binding orientations for the flavonoid molecules. The SMD simulations showed that the binding mode with higher pulling forces was linked to stronger inhibitory activity, with a strong positive correlation (R2 = 0.92) between pulling forces and IC50 values. Quercetin stood out as the most potent inhibitor, showing a pulling force of about 820 pN and an IC_(5) 0 of less than 6 µM. Further dynamic simulations indicated that quercetin's hydroxyl groups at the C3, C-5 and C-7 positions formed stable hydrogen bonds with key residues GLU-121, Leu-44 and Val-126, respectively enhancing its binding stability and effectiveness. Our results emphasized the critical role of the hydroxyl group at the C-3 position, which plays a pivotal function in effectively anchoring these molecules in the active site of Pim-1 kinase. Principal component analysis (PCA) of Pim-1 kinase's conformational changes revealed that potent inhibitors like quercetin, galangin, and kaempferol significantly restricted the enzyme's flexibility, suggesting potential inhibitory effect. These findings provide insights into the structural interactions between flavonoids and Pim-1 kinase, offering a foundation for future experimental investigations. However, further studies, including in vitro and in vivo validation, are necessary to assess the pharmacological relevance and specificity of flavonoids in cancer therapy. [ABSTRACT FROM AUTHOR]

π–π stacking are omnipresent interactions, crucial in many areas of chemistry, and often studied using quantum chemical methods. Here, we report a simple and computationally efficient method of estimating the binding energies of stacked polycyclic aromatic hydrocarbons based on steered molecular dynamics. This method leverages the force field parameters for accurate calculation. The presented results show good agreement with those obtained through DFT at the ωB97X-D3/cc-pVQZ level of theory. It is demonstrated that this force field-driven SMD method can be applied to other aromatic molecules, allowing insight into the complexity of the stacking interactions and, more importantly, reporting π–π stacking energy values with reasonable precision. [ABSTRACT FROM AUTHOR]

Pim-1 kinase, a serine/threonine kinase, is often overexpressed in various cancers, contributing to disease progression and poor prognosis. In this study, we explored the potential of flavonoids as inhibitors of Pim-1 kinase using a combination of molecular docking and steered molecular dynamics (SMD) simulations. Our docking studies revealed two main binding orientations for the flavonoid molecules. The SMD simulations showed that the binding mode with higher pulling forces was linked to stronger inhibitory activity, with a strong positive correlation (R2 ≈ 0.92) between pulling forces and IC50 values. Quercetin stood out as the most potent inhibitor, showing a pulling force of about 820 pN and an IC_(5) 0 of less than 6 µM. Further dynamic simulations indicated that quercetin’s hydroxyl groups at the C3, C-5 and C-7 positions formed stable hydrogen bonds with key residues GLU-121, Leu-44 and Val-126, respectively enhancing its binding stability and effectiveness. Our results emphasized the critical role of the hydroxyl group at the C-3 position, which plays a pivotal function in effectively anchoring these molecules in the active site of Pim-1 kinase. Principal component analysis (PCA) of Pim-1 kinase’s conformational changes revealed that potent inhibitors like quercetin, galangin, and kaempferol significantly restricted the enzyme’s flexibility, suggesting potential inhibitory effect. These findings provide insights into the structural interactions between flavonoids and Pim-1 kinase, offering a foundation for future experimental investigations. However, further studies, including in vitro and in vivo validation, are necessary to assess the pharmacological relevance and specificity of flavonoids in cancer therapy.

Molecular dynamics (MD) simulations can provide detailed insights into complex molecular systems, such as DNA, at high resolution in space and time. Using current computer architectures, time scales of tens of microseconds are feasible with contemporary all-atom force fields. However, these timescales are insufficient to accurately characterize large conformational transitions in DNA and compare calculations to experimental data. This review discusses the advantages and draw-backs of two simulation approaches to overcome the timescale challenge. The first approach is based on adding biasing potentials to the system to drive transitions. Umbrella sampling, steered MD, and metadynamics are examples of these methods. A key challenge of such methods is the necessity of selecting one or a few efficient coordinates, commonly referred to as collective variables (CVs), along which to apply the biasing potential. The path-metadynamics methodology addresses this issue by finding the optimal route(s) between states in a multi-dimensional CV space. The second strategy is path sampling, which focuses MD simulations on the transitions. The assumption is that even though transitions between states are rare, they are generally fast. Stopping the simulations as soon as they reach a stable state can significantly increase simulation efficiency. We introduce these methods on the two-dimensional Müller--Brown potential. DNA applications are featured for two different processes: the Watson--Crick--Franklin to Hoogsteen transition in adenine--thymine base pairs and the binding of a DNA-binding protein domain to DNA. This article is categorized under: Molecular and Statistical Mechanics Molecular Dynamics and Monte-Carlo Methods Molecular and Statistical Mechanics Free Energy Methods Software Simulation Methods [ABSTRACT FROM AUTHOR]

Protein folding is a prominent area of research in the life sciences. Exendin-4, a 39-amino acid peptide hormone, is of particular interest due to its potential therapeutic applications. In this study, we employed a steered molecular dynamics method to investigate the folding of Exendin-4. Our simulations reveal an intermediate state during the folding process, suggesting a more complex three-state mode of folding. Structural analysis indicates that Exendin-4 folds through two distinct pathways: pathway 1 involves gradual growth of a helical structure after the collapse of the hydrophobic core in the Trp-cage, while pathway 2 involves initial formation of local microdomains (helical structure and Trp-cage) which then combine with each other to form a stable native structure. We found that the folding along these pathways follows the hydrophobic collapse mechanism and the diffusion-collision mechanism, respectively. We believe that these two mechanisms together govern the folding of Exendin-4. These results significantly differ from those observed in most small protein folding. Our study on the folding of Exendin-4 will advance our understanding of the protein folding mechanism. [ABSTRACT FROM AUTHOR]

An atomic-scale mechanism of T Cell Receptor (TCR) mechanosensing of peptides in the binding groove of the peptide-major histocompatibility complex (pMHC) may inform the design of novel TCRs for immunotherapies. Using steered molecular dynamics simulations, our study demonstrates that mutations to peptides in the binding groove of the pMHC - which are known to discretely alter the T cell response to an antigen - alter the MHC conformation at equilibrium. This subsequently impacts the overall strength (duration and length) of the TCR-pMHC bond under constant load. Moreover, physiochemical features of the TCR-pMHC dynamic bond strength, such as hydrogen bonds and Lennard-Jones contacts, correlate with the immunogenic response elicited by the specific peptide in the MHC groove. Thus, formation of transient TCR-pMHC bonds is characteristic of immunogenic peptides, and steered molecular dynamics simulations can be used in the overall design strategy of TCRs for immunotherapies.

The rational design of T Cell Receptors (TCRs) for immunotherapy has stagnated due to a limited understanding of the dynamic physiochemical features of the TCR that elicit an immunogenic response. The physiochemical features of the TCR-peptide major histocompatibility complex (pMHC) bond dictate bond lifetime which, in turn, correlates with immunogenicity. Here, we: i) characterize the force-dependent dissociation kinetics of the bond between a TCR and a set of pMHC ligands using Steered Molecular Dynamics (SMD); and ii) implement a machine learning algorithm to identify which physiochemical features of the TCR govern dissociation kinetics. Our results demonstrate that the total number of hydrogen bonds between the CDR2β-MHC⍺(β), CDR1α-Peptide, and CDR3β-Peptide are critical features that determine bond lifetime.

SARS‐CoV‐2 gained crucial mutations at the receptor binding domain (RBD) that often changed the course of the pandemic leading to new waves with increased case fatality. Variants are observed with enhanced transmission and immune invasion abilities. Thus, predicting future variants with enhanced transmission ability is a problem of utmost research interest. Here, we have developed a multi‐tier exhaustive SARS‐CoV‐2 mutation screening platform combining MM/GBSA, extensive molecular dynamics simulations, and steered molecular dynamics to identify RBD mutants with enhanced ACE2 binding capability. We have identified four RBM mutations (F490K, S494K, G504F, and the P499L) with significantly higher ACE2 binding abilities than wild‐type RBD. Compared to wild‐type RBD, they all form stable complexes with more hydrogen bonds and salt‐bridge interactions with ACE2. Our simulation data suggest that these mutations allosterically alter the packing of the RBM interface of the RBD‐ACE2 complex. As a result, the rupture force required to break the RBD‐ACE2 contacts is significantly higher for these mutants. [ABSTRACT FROM AUTHOR]

π–π stacking are omnipresent interactions, crucial in many areas of chemistry, and often studied using quantum chemical methods. Here, we report a simple and computationally efficient method of estimating the binding energies of stacked polycyclic aromatic hydrocarbons based on steered molecular dynamics. This method leverages the force field parameters for accurate calculation. The presented results show good agreement with those obtained through DFT at the ωB97X-D3/cc-pVQZ level of theory. It is demonstrated that this force field-driven SMD method can be applied to other aromatic molecules, allowing insight into the complexity of the stacking interactions and, more importantly, reporting π–π stacking energy values with reasonable precision.

The COVID-19 pandemic raised an unprecedented race in biotechnology in search for effective therapies and a preventive vaccine. Scientists worldwide have been attempting to stop the viral infection by interfering with the biological function of the SARS-CoV-2 main protease (Mpro), a critical protein required for viral transcription and replication during infection. In this study, we employed an effective approach integrating deep learning model calculations and steered molecular dynamic simulations to generate highly promising inhibitors of SARS-CoV-2 Mpro. First, using deep learning calculations, a natural molecule that was identified as a potential inhibitor of SARS-CoV-2 Mpro was chemically altered to boost its ligand-binding affinity to the Mpro protease. The proposed compounds were then verified using steered molecular dynamic simulations to estimate their binding free energies to SARS-CoV-2 Mpro. The procedure was repeated until the binding free energies of the proposed compounds did not improve further. Overall, one proposed compound was shown to have a high nanomolar affinity, and two others were estimated to possess nanomolar affinities for SARS-CoV-2 Mpro, indicating that they are highly promising inhibitors of the protease. Absorption, distribution, metabolism, and excretion and toxicity analysis show that all three chemicals are drug-like compounds following the MACCS-II Drug Data Report database, orally absorbed, tightly attached to the plasma membrane, and noncarcinogenic in rats. The results obtained potentially support COVID-19 treatment. Cannabisin A, a natural molecule that was found to potentially inhibit SARS-CoV-2 Mpro, was modified using deep-learning calculations to improve its binding to the protease. The resulting compounds were verified using steered-molecular dynamics simulations, and three were identified with high binding affinities. All three proposed compounds are considered drug-like, absorbed orally, attached to the plasma membrane, and non-carcinogenic in rats according to ADMET analysis. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations of protein folding typically consider the polypeptide chain at equilibrium and in isolation from the cellular components. We argue that in order to understand protein folding as it occurs in vivo, it should be modeled as an active, energy-dependent process, in which the cellular protein-folding machine directly manipulates the polypeptide. We conducted all-atom molecular dynamics simulations of four protein domains, whose folding from the extended state was augmented by the application of rotational force to the C-terminal amino acid, while the movement of the N-terminal amino acid was restrained. We have shown earlier that such a simple manipulation of peptide backbone facilitated the formation of native structures in diverse α-helical peptides. In this study, the simulation protocol was modified, to apply the backbone rotation and movement restriction only for a short time at the start of simulation. This transient application of a mechanical force to the peptide is sufficient to accelerate, by at least an order of magnitude, the folding of four protein domains from different structural classes to their native or native-like conformations. Our in silico experiments show that a compact stable fold may be attained more readily when the motions of the polypeptide are biased by external forces and constraints. [ABSTRACT FROM AUTHOR]

Mechanics of polysaccharides, the most diverse class of biomolecules, is explored in this work for a variety of strands bearing α- or β-D-glucopyranose residues and linked by different glycosidic linkages. The purpose is to understand how far the glycosidic linkage type and anomericity of the monomers influence the elasticity of the polysaccharides. We have used steered molecular dynamics (SMD) simulations to understand the mechanical strain under a pulling stress. The force-extension (F-E) curves thus generated are compared to the results of atomic force microscopy (AFM) for validation. In line with expectations, our SMD results indicate that in general, the α D-glucopyranosides show conformational switching from the most stable 4C1 chair to a half-chair or envelope and even to a skew-boat or boat form in some cases. The β-series, however, show minimal distortion, except β D-1-2 polysaccharide, which shows some transition. The exceptional cases of 1-6 linked polysaccharides and their extra elasticity owing to the presence of an additional –CH2 group in the glycosidic linkage are also evident. In general, amongst all the polysaccharides studied, α D-1-6 glucopyranoside shows a maximum elasticity of ~ 37% per ring when exposed to a pulling stress of ~ 2100 pN. Density functional theory (DFT) calculations help estimate the single point energies of the conformers as the polysaccharide stretches in the SMD simulations. A comparative understanding of polysaccharide mechanics, as obtained from this work, would set the platform for modeling and simulation of innovative polysaccharide-based membranes down the line. [ABSTRACT FROM AUTHOR]

Acetylcholinesterase (AChE) is a key target for the current symptomatic treatment of Alzheimer's disease, and galantamine is a clinical anticholinesterase drug with transiently acting characteristic and good selectivity for AChE. The present theoretical‐experimental work improves the drug's residence time without reducing the inhibition effect, thus providing a crucial breakthrough for modifying the inhibitor of AChE with better kinetic behavior. The static binding and dynamic delivery properties acquired from atomic view reveal that the galantamine simply occupies a catalytic anionic site, and its release from AChE needs only ∼8.6 kcal/mol. Both of these may cause the short residence time of galantamine. The hotspots and most favorable transport mechanism are identified, and the hydrogen bond and aromatic stacking interactions are observed to play crucial roles for galantamine binding and release in AChE. The typical peripheral anionic site arisen at the delivery process would provide another key occupation to enhance the anti‐release ability for inhibitors. The compound with "specific‐ring‐chain‐ring" framework with detailed beneficial modification scheme is summarized, which may improve the residence time of the inhibitor in AChE. The thermodynamic and dynamic properties of galantamine derivatives are also studied. Based on dictamnine, a natural alkaloid, two novel eligible derivatives are designed, synthesized and evaluated, which verifies our prediction. Multiple computational approaches and experimental combinations probably provide a train of thought from both static and dynamic views to modify or design appropriate inhibitors on the basis of specific binding and transportation features. [ABSTRACT FROM AUTHOR]

Ordinary small molecule de novo drug design is time-consuming and expensive. Recently, computational tools were employed and proved their efficacy in accelerating the overall drug design process. Molecular dynamics (MD) simulations and a derivative of MD, steered molecular dynamics (SMD), turned out to be promising rational drug design tools. In this paper, we report the first application of SMD to evaluate the binding properties of small molecules toward FABP4, considering our recent interest in inhibiting fatty acid binding protein 4 (FABP4). FABP4 inhibitors (FABP4is) are small molecules of therapeutic interest, and ongoing clinical studies indicate that they are promising for treating cancer and other diseases such as metabolic syndrome and diabetes. [ABSTRACT FROM AUTHOR]

This study aims to understand the mechanical properties of cellulose nanofibers (CNFs), a nano-sized material element of woods or plants. We develop all-atom (AA) molecular dynamics models of cellulose microfibrils (CMFs), which are the smallest constituent of CNFs. The models were designed for the process of structural failure or the degradation of a hierarchical material of multiple CMF fibers, due to shear deformation. It was assumed that two CMFs were arranged in parallel and in close contact, either in a vacuum or in water. The CMF models in water were built by surrounding AA-modeled water molecules with a few nanometers. Shear deformation was applied in the axial direction of the CMF or in the direction parallel to molecular sheets. Shear moduli were measured, and they agree with previous experimental and computational values. The presence of water molecules reduced the elastic modulus, because of the behavior of water molecules at the interface between CMFs as a function of temperature. In the inelastic region, the CMF often broke down inside CMFs in a vacuum condition. However, in water environments, two CMFs tend to slip away from each other at the interface. Water molecules act like a lubricant between multiple CMFs and promote smooth sliding. [ABSTRACT FROM AUTHOR]

Directed self-assembly of nanoparticles (DSA-n) over a templated surface is of wide interest in many technological applications such as semiconducting, plasmonic, and photovoltaic devices and, in particular, molecular-specific imaging and sensing. Understanding the complex energetics of nanoparticle deposition will aid in improving the yield and scalability of DSA-n. In this study, we report steered molecular dynamics (SMD) simulations to predict the free-energy contribution to DSA-n over a templated surface. We utilize many-body dissipative particle dynamics along with Jarzynski's equality for both stagnant and receding liquid film, as a nonequilibrium system, for DSA-n. We apply the harmonic spring approximation in SMD simulations and calculate the average work, internal energy, entropy, and free energy for various nanoparticle densities in the liquid film. The results of SMD simulations in the stagnant and receding liquid film suggest that the free energy is a concave function of the nanoparticle density. The simulation results also indicate that the impact of the receding liquid film on the free energy associated with DSA of a nanoparticle is similar to that of the stagnant liquid film. The understanding gained from this study will shed some light on the mechanism of DSA-n, and the ways to increase the yield and decrease the defect density. [ABSTRACT FROM AUTHOR]

Non-equilibrium molecular dynamics (NEMD) simulation has been recognized as a powerful tool for examining biomolecules and provides fruitful insights into not only non-equilibrium but also equilibrium processes. We review recent advances in NEMD simulation and relevant, fundamental results of non-equilibrium statistical mechanics. We first introduce Crooks fluctuation theorem and Jarzynski equality that relate free energy difference to work done on a physical system during a non-equilibrium process. The theorems are beneficial for the analysis of NEMD trajectories. We then describe rate theory, a framework to calculate molecular kinetics from a non-equilibrium process; this theoretical framework enables us to calculate a reaction time—mean-first passage time—from NEMD trajectories. We, in turn, present recent NEMD techniques that apply an external force to a system to enhance molecular dissociation and introduce their application to biomolecules. Lastly, we show the current status of an appropriate selection of reaction coordinates for NEMD simulation. [ABSTRACT FROM AUTHOR]

In this work, the effect of common functional groups, namely hydroxyl, formyl, carboxyl, and amine groups on the interfacial behavior of surface-modified glass fiber-reinforced epoxy is investigated at molecular scale. The interfacial properties of the epoxy/silica coated with different functional group systems are quantified by performing pulling test using the steered molecular dynamics simulations. It is found that the system with hydroxyl groups has a relatively lower interfacial interaction, exhibiting an adhesive failure mode. When partial hydroxyl groups are replaced by carboxyl, amine, and formyl groups, respectively, the interfacial interactions are increased and these systems exhibit a cohesive failure mode where failure happens in the epoxy close to interface. A relatively higher force is required for the adhesive debonding, while more energy can be dissipated for the cohesive debonding. Because the increased interfacial interactions can prevent the mobility of polymer chains, and delay the propagation of micropores in the matrix, leading to the epoxy matrix with a high ability of energy absorption. Our work provides an insight into how functional groups affect the interface debonding behavior of glass fiber-reinforced epoxy, offering a guideline for control of the interfacial properties of such composites through surface modification techniques.

• The bonding role of NPBA on different nitramine oxidant surfaces was investigated by MD simulation. • The peak pulling force of NPBA in the peeling process on HMX surface is higher than that of PEG. • A new peeling model is proposed to reveal the adhesion weakening effect in the nitrate ester environment. • The effect of NPBA on interfacial micro-peeling phenomenon and mechanical response is observed. • Efficient engineering application schemes of NPBA are recommended based on MD simulation results. Excellent mechanical properties are the lifeline for the application of solid propellants. Understanding the strengthening mechanism of bonding agent is of great significance for design and application of high-performance bonding agents. In this study, the strengthening mechanism of neutral polymer bonding agent (NPBA) for the NEPE propellant and dynamic mechanical response were explored through MD simulation. Significant difference in bonding role for RDX and HMX surface is first revealed with a more effective improvement in HMX surface. NPBA requires higher peak pulling force to achieve conformational transformation compared with binder throughout the peeling process. The chain orientation and movement undergo significant changes in nitroglycerin environment and an interface peeling model is proposed to reveal the potential mechanism. Finally, the effect of NPBA insertion on interfacial micro-peeling process indicates that no significant detachment of NPBA until the binder was severely deformed, resulting in a higher pulling force required. [ABSTRACT FROM AUTHOR]

This study investigates the role of the NH 2 -functionalized MIL-53(Fe) nanoparticle as a carrier for the anti-cancer drug 5-fluorouracil (5-Fu) through molecular dynamics (MD) simulation. The Atom in Molecules (AIM) method was used to determine the interaction type between 5-Fu and NH 2 -MIL-53(Fe). The Steered Molecular Dynamics (SMD) approach was employed to calculate the free energy of drug encapsulation. The simulation included the NH 2 -MIL-53(Fe) nanoparticle containing 5-FU molecules and water molecules at different temperatures and 1 atm pressure. The properties analyzed included density, radial distribution functions, displacement, diffusion, and binding energy of the drug with NH 2 -MIL-55(Fe). It is shown that the drug was encapsulated in the framework channels and the variation of density of the system with temperature in the presence of nanoparticles decreased. The drug's mean square displacement, total self-diffusion coefficient, and diffusion coefficient in channel direction increased with time as temperature rose. The drug was aligned to have its oxygen atoms towards the metal nodes of the framework. However, The nitrogen atom of the amine functional group in NH2-MIL-53(Fe) interacts more with the F atom of 5-Fu which shows the effect of this functional group on the adsorption of the 5-Fu. The drug's adsorption percentage increased from 298 K to 328 K. The calculated binding energy increased with temperatures and was desirable as it countered the drug's repulsive forces with atoms of NH 2 -MIL-53(Fe). The results of the molecular docking simulation confirmed MD simulation. [Display omitted] • Loading of 5-Fu into the NH 2 -MIL-53(Fe) was measured by MD method and compared with experiment. • Orientation of the 5-Fu into the NH 2 -MIL-53(Fe) was investigated by MD and Docking methods. • The NH 2 group of the NH 2 -MIL-53(Fe) has more interaction with the F atom of 5-Fu. • Binding energy, diffusion of 5-Fu and the RDF of 5-Fu atoms were investigated. • The free energy of drug encapsulation was calculated by Steered Molecular Dynamics approach. [ABSTRACT FROM AUTHOR]

The role of hydrophobic and polar interactions in providing thermodynamic stability to folded proteins has been intensively studied, but the relative contribution of these interactions to the mechanical stability is less explored. We used steered molecular dynamics simulations with constant-velocity pulling to generate force-extension curves of selected protein domains and monitor hydrophobic surface unravelling upon extension. Hydrophobic contribution was found to vary between one fifth and one third of the total force while the rest of the contribution is attributed primarily to hydrogen bonds. Moreover, hydrophobic force peaks were shifted towards larger protein extensions with respect to the force peaks attributed to hydrogen bonds. The higher importance of hydrogen bonds compared to hydrophobic interactions in providing mechanical resistance is in contrast with the relative importance of the hydrophobic interactions in providing thermodynamic stability of proteins. The different contributions of these interactions to the mechanical stability are explained by the steeper free energy dependence of hydrogen bonds compared to hydrophobic interactions on the relative positions of interacting atoms. Comparative analyses for several protein domains revealed that the variation of hydrophobic forces is modest, while the contribution of hydrogen bonds to the force peaks becomes increasingly important for mechanically resistant protein domains.

The rational design of T Cell Receptors (TCRs) for immunotherapy has stagnated due to a limited understanding of the dynamic physiochemical features of the TCR that elicit an immunogenic response. The physiochemical features of the TCR-peptide major histocompatibility complex (pMHC) bond dictate bond lifetime which, in turn, correlates with immunogenicity. Here, we: i) characterize the force-dependent dissociation kinetics of the bond between a TCR and a set of pMHC ligands using Steered Molecular Dynamics (SMD); and ii) implement a machine learning algorithm to identify which physiochemical features of the TCR govern dissociation kinetics. Our results demonstrate that the total number of hydrogen bonds between the CDR2β-MHC⍺(β), CDR1α-Peptide, and CDR3β-Peptide are critical features that determine bond lifetime.

The coronavirus E proteins are small membrane proteins found in the virus envelope of alpha and beta coronaviruses that have a high degree of overlap in their biochemical and functional properties despite minor sequence variations. The SARS-CoV-2 E is a 75-amino acid transmembrane protein capable of acting as an ion channel when assembled in a pentameric fashion. Various studies have found that hexamethylene amiloride (HMA) can inhibit the ion channel activity of the E protein in bilayers and also inhibit viral replication in cultured cells. Here, we use the available structural data in conjunction with homology modelling to build a comprehensive model of the E protein to assess potential binding sites and molecular interactions of HMA derivatives. Furthermore, we employed an iterative cycle of molecular modelling, extensive docking simulations, molecular dynamics and leveraging steered molecular dynamics to better understand the pore characteristics and quantify the affinity of the bound ligands. Results from this work highlight the potential of acylguanidines as blockers of the E protein and guide the development of subsequent small molecule inhibitors. [ABSTRACT FROM AUTHOR]

A variety of 10 polysaccharide motifs comprising naturally inspired sequences of monosaccharide building blocks are studied to understand the inter-relation between their structural and mechanical behaviors. Equilibrium and steered molecular dynamics (SMD) simulations are employed to investigate the stress-strain relationships and the associated conformational flips of the pyranose moieties along the polysaccharide chains. The presence of a variety of glycosidic linkages connecting the diverse monosaccharide units along with chain-branching in some cases induce wide diversity in the carbohydrate-Ramachandran plots of the glycosidic dihedrals. Similar variations are observed in the Cremer-Pople ring puckering patterns across the polysaccharide variants. The work further provides a comparison between the experimentally obtained atomic force microscopic data of mechanical stretching for some polysaccharides with the stress-strain curves generated from our SMD simulations. Out of all the systems studied, pectin having an axial-axial orientation of the glycosidic linkage showed maximum stretching potential, while acetan-M, with an equatorial-equatorial disposition of the glycosidic bond, stretched the least. The experimental Young's modulus of the corresponding natural polysaccharides could be reasonably compared to the values obtained from our simulation models. Force distribution analysis is done to understand the propagation of punctual stress in the polysaccharides under SMD conditions. Changes in local electrophilicity or nucleophilicity of atomic centers in puckered pyranose rings are estimated through the condensed Fukui functions. All of this information can help understand the physical behavior and chemical reactivity of complex polysaccharides in a complicated milieu of electronic and steric effects experienced by them.

Since the dawn of the computer age, scientists have designed devices to represent molecular structures and developed tools to simulate their dynamic behavior in silico. To this day, these tools remain central to our understanding of biomolecular phenomena. In contrast to other fields such as fluid mechanics or meteorology, the observation of molecular motions at the atomic level remains a major experimental challenge. Continuous advances in computer graphics and numerical computation, combined with the emergence of human–computer interaction approaches, led to the methodology of so‐called "interactive molecular simulations", characterized by two main features. First, the possibility to visualize a running simulation in interactive time, that is, compatible with human perception. Second, the possibility to manipulate the simulation interactively by imposing a force, changing a biophysical property, or editing runtime parameters on the fly. Such simulations are still little used in computational biology, where it is more common to run a series of offline simulations and then visualize and analyze the results. However, interactive molecular simulation tools promise to handle time‐consuming tasks such as the modeling of particularly complex biomolecular structures more efficiently or to support approaches such as Rational Drug Design with regard to pharmaceutical applications. This article is categorized under:Computer and Information Science > Computer Algorithms and ProgrammingComputer and Information Science > VisualizationSoftware > Simulation Methods [ABSTRACT FROM AUTHOR]

The native structures of proteins, except for notable exceptions of intrinsically disordered proteins, in general take their most stable conformation in the physiological condition to maintain their structural framework so that their biological function can be properly carried out. Experimentally, the stability of a protein can be measured by several means, among which the pulling experiment using the atomic force microscope (AFM) stands as a unique method. AFM directly measures the resistance from unfolding, which can be quantified from the observed force‐extension profile. It has been shown that key features observed in an AFM pulling experiment can be well reproduced by computational molecular dynamics simulations. Here, we applied computational pulling for estimating the accuracy of computational protein structure models under the hypothesis that the structural stability would positively correlated with the accuracy, i.e. the closeness to the native, of a model. We used in total 4929 structure models for 24 target proteins from the Critical Assessment of Techniques of Structure Prediction (CASP) and investigated if the magnitude of the break force, that is, the force required to rearrange the model's structure, from the force profile was sufficient information for selecting near‐native models. We found that near‐native models can be successfully selected by examining their break forces suggesting that high break force indeed indicates high stability of models. On the other hand, there were also near‐native models that had relatively low peak forces. The mechanisms of the stability exhibited by the break forces were explored and discussed. [ABSTRACT FROM AUTHOR]

N‐acetylglucosamine containing compounds acting as pathogenic or symbiotic signals are perceived by plant‐specific Lysin Motif Receptor‐Like Kinases (LysM‐RLKs). The molecular mechanisms of this perception are not fully understood, notably those of lipo‐chitooligosaccharides (LCOs) produced during root endosymbioses with nitrogen‐fixing bacteria or arbuscular mycorrhizal fungi. In Medicago truncatula, we previously identified the LysM‐RLK LYR3 (MtLYR3) as a specific LCO‐binding protein. We also showed that the absence of LCO binding to LYR3 of the non‐mycorrhizal Lupinus angustifolius, (LanLYR3), was related to LysM3, which differs from that of MtLYR3 by several amino acids and, particularly, by a critical tyrosine residue absent in LanLYR3. Here, we aimed to define the LCO binding site of MtLYR3 by using molecular modelling and simulation approaches, combined with site‐directed mutagenesis and LCO binding experiments. 3D models of MtLYR3 and LanLYR3 ectodomains were built, and homology modelling and molecular dynamics (MD) simulations were performed. Molecular docking and MD simulation on the LysM3 identified potential key residues for LCO binding. We highlighted by steered MD simulations that in addition to the critical tyrosine, two other residues were important for LCO binding in MtLYR3. Substitution of these residues in LanLYR3‐LysM3 by those of MtLYR3‐LysM3 allowed the recovery of high‐affinity LCO binding in experimental radioligand‐binding assays. An analysis of selective constraints revealed that the critical tyrosine has experienced positive selection pressure and is absent in some LYR3 proteins. These findings now pave the way to uncover the functional significance of this specific evolutionary pattern. [ABSTRACT FROM AUTHOR]

YidC is a membrane protein that facilitates the insertion of newly synthesized proteins into lipid membranes. Through YidC, proteins are inserted into the lipid bilayer via the SecYEG-dependent complex. Additionally, YidC functions as a chaperone in protein folding processes. Several studies have provided evidence of its independent insertion mechanism. However, the mechanistic details of the YidC SecY-independent protein insertion mechanism remain elusive at the molecular level. This study elucidates the insertion mechanism of YidC at an atomic level through a combination of equilibrium and non-equilibrium molecular dynamics (MD) simulations. Different docking models of YidC-Pf3 in the lipid bilayer were built in this study to better understand the insertion mechanism. To conduct a complete investigation of the conformational difference between the two docking models developed, we used classical molecular dynamics simulations supplemented with a non-equilibrium technique. Our findings indicate that the YidC transmembrane (TM) groove is essential for this high-affinity interaction and that the hydrophilic nature of the YidC groove plays an important role in protein transport across the cytoplasmic membrane bilayer to the periplasmic side. At different stages of the insertion process, conformational changes in YidC’s TM domain and membrane core have a mechanistic effect on the Pf3 coat protein. Furthermore, during the insertion phase, the hydration and dehydration of the YidC’s hydrophilic groove are critical. These results demonstrate that Pf3 coat protein interactions with the membrane and YidC vary in different conformational states during the insertion process. Finally, this extensive study directly confirms that YidC functions as an independent insertase.

Ordinary small molecule de novo drug design is time-consuming and expensive. Recently, computational tools were employed and proved their efficacy in accelerating the overall drug design process. Molecular dynamics (MD) simulations and a derivative of MD, steered molecular dynamics (SMD), turned out to be promising rational drug design tools. In this paper, we report the first application of SMD to evaluate the binding properties of small molecules toward FABP4, considering our recent interest in inhibiting fatty acid binding protein 4 (FABP4). FABP4 inhibitors (FABP4is) are small molecules of therapeutic interest, and ongoing clinical studies indicate that they are promising for treating cancer and other diseases such as metabolic syndrome and diabetes.

In COVID-19 infection, the SARS-CoV-2 spike protein S1 interacts to the ACE2 receptor of human host, instigating the viral infection. To examine the competitive inhibitor efficacy of broad spectrum alpha helical AMPs extracted from frog skin, a comparative study of intermolecular interactions between viral S1 and AMPs was performed relative to S1-ACE2p interactions. The ACE2 binding region with S1 was extracted as ACE2p from the complex for ease of computation. Surprisingly, the Spike-Dermaseptin-S9 complex had more intermolecular interactions than the other peptide complexes and importantly, the S1-ACE2p complex. We observed how atomic displacements in docked complexes impacted structural integrity of a receptor-binding domain in S1 through conformational sampling analysis. Notably, this geometry-based sampling approach confers the robust interactions that endure in S1-Dermaseptin-S9 complex, demonstrating its conformational transition. Additionally, QM calculations revealed that the global hardness to resist chemical perturbations was found more in Dermaseptin-S9 compared to ACE2p. Moreover, the conventional MD through PCA and the torsional angle analyses indicated that Dermaseptin-S9 altered the conformations of S1 considerably. Our analysis further revealed the high structural stability of S1-Dermaseptin-S9 complex and particularly, the trajectory analysis of the secondary structural elements established the alpha helical conformations to be retained in S1-Dermaseptin-S9 complex, as substantiated by SMD results. In conclusion, the functional dynamics proved to be significant for viral Spike S1 and Dermaseptin-S9 peptide when compared to ACE2p complex. Hence, Dermaseptin-S9 peptide inhibitor could be a strong candidate for therapeutic scaffold to prevent infection of SARS-CoV-2. [ABSTRACT FROM AUTHOR]

Ivermectin (D1), an FDA-approved drug, is mainly used as an inhibitor for the treatment of parasitic infections. This drug was drawn attention due to its effectiveness in COVID-19 treatment. In this work, we investigated how ivermectin and its analogues are interacted with Impα/β1 heterodimer by molecular docking, steered molecular dynamics (SMD), classical molecular dynamics (MD), and MM/PBSA binding free energy analysis. Docking results showed that ivermectin dimer showed the highest binding affinity of −12.2 kcal/mol compared to its monomer. In SMD, the highest acceleration of 600 pm/ps2 is noticed for D1-Impα/β1 complex in the distance from 18.27 to 18.36 Å. The pulling force of 1745.86 pN is also detected for D1-Impα/β1. To validate the docking and SMD results, 100 ns molecular dynamics (MD) simulation is performed on the D1-Impα/β1 complex. The average RMSD value indicates a good structural stability of the complex despite some significant changes at the beginning. It is noted that most of the residues are stable over the simulation time with an average RMSF value of less than 2.65 Å. The MM/PBSA free energy of ivermectin also shows strong and spontaneous binding with the Impα/β1. [ABSTRACT FROM AUTHOR]

Recent studies have showed that five flavonoids (Epigallocatechin gallate [EGCG], Rutin, Silibinin, Genistein, and Quercetin) are potential therapeutic agents for type II diabetes mellitus (T2DM). Our previous study revealed the mechanism of these five flavonoids with human islet amyloid peptide, but the mechanism of action with other T2DM‐related target proteins is unknown. In this study, we investigate the competitive inhibition mechanism of N‐terminal catalytic domain maltase glucose‐amylase (NtMGAM), which is another important target of T2DM, with the five types of flavonoids through molecular docking, molecular dynamics, and steered molecular dynamics. Our results show that the NtMGAM active pocket tends to bind to the aromatic nucleus. Moreover, EGCG and Rutin can potentially inhibit NtMGAM, and both molecules exhibit a larger binding capability with the active pocket than Acarbose. Compared with the hexasaccharide group, the pentasaccharide group has the advantage to bind with the NtMGAM active pocket, suggesting that the structure embedded into the active pocket should be inadvisable to oversize. Meanwhile, solvent molecules are a crucial factor affecting the binding stability of inhibitor with NtMGAM. In conclusion, this work is helpful to provide a theoretical basis for screening, modifying, and designing new therapeutic or preventive drugs targeting T2DM. [ABSTRACT FROM AUTHOR]

The weak interfacial adhesion has significantly affected the durability, long-term reliability, and performance of glass fiber–reinforced epoxy composites. The coating of graphene and carbon nanotubes on the glass fiber can have a positive effect on the strength, toughness, and thermal insulation performance of glass fiber-reinforced composites. However, the strengthening mechanism of carbon nanomaterial coating on the interfacial adhesion between glass fiber and epoxy has not been fully explored. In this work, the effect of graphene and single-walled carbon nanotubes (SWCNTs) on the interfacial properties of the glass fiber–reinforced epoxy has been investigated at atomistic scale. The graphene and SWCNTs are sandwiched between epoxy and silica to study the debonding behavior of the sandwiched structures. It is found that the interfacial energy is significantly improved with the incorporation of graphene and SWCNTs between epoxy and silica, causing an obvious improvement in adhesion stress for graphene coating and in debonding displacement for SWCNT coating. Compared with the epoxy/silica without coatings where the silica and epoxy detach from the contact surface, the sandwiched structures display different failure modes. The sandwiched structure with graphene coating fails at the epoxy matrix close to the interface, exhibiting a cohesive failure mode because of the relatively stronger interfacial interactions. The structures with SWCNTs fail at the interface between silica and SWCNTs, representing an adhesive failure mode due to the interlocking between SWCNTs and polymer chains. This work provides a theoretical guideline to optimize the interface adhesion of coated glass fiber–reinforced epoxy via structure design and surface modification of coating materials.

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

The present article aims to investigate the mechanical properties of a new DNA origami nano-joint using the steered molecular dynamics (SMD) simulation. Since the analysis of mechanical properties is of great importance in bending conditions for a nano-joint, the forces are selected to achieve angular changes in the joint by the resultant torque. In this study, the nano-joint is considered as a beam in order to use mechanical equations to extract the mechanical properties of the designed nano-joint. In addition, the bending stiffness of the beam is investigated in different modes of deflection using the Euler–Bernoulli beam theory. The results revealed that the value of bending stiffness increases with increasing deflection, and the changes in the bending stiffness relative to the deflection is linear. The proposed DNA origami nano-joint can be used as a joint in nanorobots and can be effectively applied in nanorobotic systems to move different components. [ABSTRACT FROM AUTHOR]

Continuing our search to find more potent and selective 5-LOX inhibitors, we present now the enzymatic evaluation of seventeen isoflavones (IR) and nine isoflavans (HIR), and their in vitro and in cellulo potency against human leukocyte 5-LOX. Of the 26 compounds tested, 10 isoflavones and 9 isoflavans possessed micromolar potency, but only three were selective against 5-LOX (IR-2, HIR-303, and HIR-309), with IC50 values at least 10 times lower than those of 12-LOX, 15-LOX-1, and 15-LOX-2. Of these three, IR-2 (6,7-dihydroxy-4-methoxy-isoflavone, known as texasin) was the most selective 5-LOX inhibitor, with over 80-fold potency difference compared with other isozymes; Steered Molecular Dynamics (SMD) studies supported these findings. The presence of the catechol group on ring A (6,7-dihydroxy versus 7,8-dihydroxy) correlated with their biological activity, but the reduction of ring C, converting the isoflavones to isoflavans, and the substituent positions on ring B did not affect their potency against 5-LOX. Two of the most potent/selective inhibitors (HIR-303 and HIR-309) were reductive inhibitors and were potent against 5-LOX in human whole blood, indicating that isoflavans can be potent and selective inhibitors against human leukocyte 5-LOX in vitro and in cellulo.

The coronavirus E proteins are small membrane proteins found in the virus envelope of alpha and beta coronaviruses that have a high degree of overlap in their biochemical and functional properties despite minor sequence variations. The SARS-CoV-2 E is a 75-amino acid transmembrane protein capable of acting as an ion channel when assembled in a pentameric fashion. Various studies have found that hexamethylene amiloride (HMA) can inhibit the ion channel activity of the E protein in bilayers and also inhibit viral replication in cultured cells. Here, we use the available structural data in conjunction with homology modelling to build a comprehensive model of the E protein to assess potential binding sites and molecular interactions of HMA derivatives. Furthermore, we employed an iterative cycle of molecular modelling, extensive docking simulations, molecular dynamics and leveraging steered molecular dynamics to better understand the pore characteristics and quantify the affinity of the bound ligands. Results from this work highlight the potential of acylguanidines as blockers of the E protein and guide the development of subsequent small molecule inhibitors.

We applied a new large two membrane model of human urea transporter B (UT-B) with urea concentration gradient between intracellular and extracellular solution to explore the structure and dynamics of urea permeation. The first solved human UT-B X-ray crystal structure (PDB 6QD5) was studied using atomistic force field based large scale molecular dynamics (MD) simulations. The urea binding sites in UT-B were identified by the combination of first principles molecular dynamics docking and hybrid steered molecular dynamics (hSMD). In the simulations, urea molecules bound to the two binding sites although urea molecules did not go through the channels. The potential of mean force (PMF) of urea along the channel was calculated by hSMD to investigate the selectivity of the urea permeation across the channel pore. The simulations provide an opportunity to explore the binding sites spontaneously to help understand the urea permeation gating and transportation mechanism. A clearer analysis of the selectivity filter is shown in more detail. We also furthered the research about small organic molecules in the specific environment, such as channels in this work. [ABSTRACT FROM AUTHOR]

Dyneins are an AAA+ motor responsible for motility and force generation toward the minus end of microtubules. Dynein motility is powered by nucleotide-dependent transitions of its linker domain, which transitions between straight (post-powerstroke) and bent (pre-powerstroke) conformations. To understand the dynamics and energetics of the linker, we performed all-atom molecular dynamics simulations of human dynein-2 primed for its power stroke. Simulations revealed that the linker can adopt either a bent conformation or a semi-bent conformation, separated by a 5.7 kT energy barrier. The linker cannot switch back to its straight conformation in the pre-powerstroke state due to a steric clash with the AAA+ ring. Simulations also showed that an isolated linker has a free energy minimum near the semi-bent conformation in the absence of the AAA+ ring, indicating that the linker stores energy as it bends and releases this energy during the powerstroke. [Display omitted] • The priming stroke of dynein is modeled using all-atom MD simulations • The linker stores energy in its bent conformation • The linker can adopt a semi-bent conformation in the pre-powerstroke state • Straightening of the linker is gated due to a steric clash with the AAA+ ring Dynein motility is powered by the swinging motion of its linker domain at the surface of the AAA+ ring. Modeling studies by Golcuk et al. found that the linker can attain two distinct bent conformations and cannot straighten due to the steric clash with the pre-powerstroke conformation of the ring. [ABSTRACT FROM AUTHOR]

In this study, a reasonable 3D structure of enzyme IIC (EIIC) was established using ab initio methods. Curcumin, baicalein and phloretin were screened out from 127,695 natural products based on EIIC protein in Staphylococcus aureus (S. aureus) by using the virtual screening technology. The antibacterial assays results confirmed that these three compounds could inhibit S. aureus cell growth. The results of molecular docking, molecular dynamics and molecular mechanics/generalized born surface area calculation indicated that curcumin, baicalein and phloretin mainly bound to Glu275, His231, and Phe359 of EIIC through hydrogen bonds and van der Waals forces to obstruct the transport of glucose-based carbohydrates, and thereby inhibited the carbohydrate metabolism of S. aureus. The steered molecular dynamics analysis revealed the existence of transfer channel in EIIC and illustrated that curcumin, baicalein and phloretin could dissociate from EIIC when pulling force on them. Furthermore, alanine scanning and hydrogen bond occupancy analyses found that Glu275 played an important role in stabilizing EIIC-baicalein/curcumin/phloretin complexes, and 275Glu@OE1 was a stable hydrogen bond acceptor. Therefore, Glu275 may be a hot spot residue that can maintain the stable binding of agents to EIIC, which may be a direction for the development of new drugs against S. aureus. • Curcumin, baicalein and phloretin were screened out by virtual screening technology. • The compounds can bind to EIIC to hinder sugar intake inhibiting the bacteria growth. • H-bonds and van der Waals forces had contributions to the EIIC-ligand complexes. • When pulling force on these three compounds, they could dissociate from EIIC. [ABSTRACT FROM AUTHOR]

Escherichia coli express adhesion pili that mediate attachment to host cell surfaces and are exposed to body fluids in the urinary and gastrointestinal tracts. Pilin subunits are organized into helical polymers, with a tip adhesin for specific host binding. Pili can elastically unwind when exposed to fluid flow forces, reducing the adhesin load, thereby facilitating sustained attachment. Here we investigate biophysical and structural differences of pili commonly expressed on bacteria that inhabit the urinary and intestinal tracts. Optical tweezers measurements reveal that class 1a pili of uropathogenic E. coli (UPEC), as well as class 1b of enterotoxigenic E. coli (ETEC), undergo an additional conformational change beyond pilus unwinding, providing significantly more elasticity to their structure than ETEC class 5 pili. Examining structural and steered molecular dynamics simulation data, we find that this difference in class 1 pili subunit behavior originates from an α-helical motif that can unfold when exposed to force. A disulfide bond cross-linking β-strands in class 1 pili stabilizes subunits, allowing them to tolerate higher forces than class 5 pili that lack this covalent bond. We suggest that these extra contributions to pilus resiliency are relevant for the UPEC niche, since resident bacteria are exposed to stronger, more transient drag forces compared to those experienced by ETEC bacteria in the mucosa of the intestinal tract. Interestingly, class 1b ETEC pili include the same structural features seen in UPEC pili, while requiring lower unwinding forces that are more similar to those of class 5 ETEC pili. [ABSTRACT FROM AUTHOR]

Polyvinylidene fluoride (PVDF) has several crystal forms of which the α‐form is nonpolar, while the β‐form is polar and has the highest piezoelectric constant. α PVDF, when stretched, transforms into the β form, which has wide applications in sensors and actuators. Steered molecular dynamics simulations are used to study the transformation of a single chain of PVDF from a trans–gauche conformation to an all trans one. The Helmholtz free energy change (∆F) is estimated using Jarzynski's equality. The transformation starts at the chain ends followed by the transformation of the remaining chain. The free energy change for the transformation is found to be always positive, indicating that the TGTG' form has higher thermodynamic stability than the all trans form throughout the studied temperature range. With increasing temperature, free energy change for the transformation increases monotonically. [ABSTRACT FROM AUTHOR]

Introduction: E-selectin is a member of the selectin family of cell adhesion molecules expressed on the plasma membrane of inflamed endothelium and facilitates initial leukocyte tethering and subsequent cell rolling during the early stages of the inflammatory response via binding to glycoproteins expressing sialyl LewisX and sialyl LewisA (sLeX/A). Existing crystal structures of the extracellular lectin/EGF-like domain of E-selectin complexed with sLeX have revealed that E-selectin can exist in two conformation states, a low affinity (bent) conformation, and a high affinity (extended) conformation. The differentiating characteristic of the two conformations is the interdomain angle between the lectin and the EGF-like domain. Methods: Using molecular dynamics (MD) simulations we observed that in the absence of tensile force E-selectin undergoes spontaneous switching between the two conformational states at equilibrium. A single amino acid substitution at residue 2 (serine to tyrosine) on the lectin domain favors the extended conformation. Results: Steered molecular dynamics (SMD) simulations of E-selectin and PSGL-1 in conjunction with experimental cell adhesion assays show a longer binding lifetime of E-selectin (S2Y) to PSGL-1 compared to wildtype protein. Conclusions: The findings in this study advance our understanding into how the structural makeup of E-selectin allosterically influences its adhesive dynamics. [ABSTRACT FROM AUTHOR]

Protein-protein complexes can vary in mechanical stability depending on the direction from which force is applied. Here we investigated the anisotropic mechanical stability of a molecular complex between a therapeutic non-immunoglobulin scaffold called Affibody and the extracellular domain of the immune checkpoint protein PD-L1. We used a combination of single-molecule AFM force spectroscopy (AFM-SMFS) with bioorthogonal clickable peptide handles, shear stress bead adhesion assays, molecular modeling, and steered molecular dynamics (SMD) simulations to understand the pulling point dependency of mechanostability of the Affibody:(PD-L1) complex. We observed diverse mechanical responses depending on the anchor point. For example, pulling from residue #22 on Affibody generated an intermediate unfolding event attributed to partial unfolding of PD-L1, while pulling from Affibody's N-terminus generated force-activated catch bond behavior. We found that pulling from residue #22 or #47 on Affibody generated the highest rupture forces, with the complex breaking at up to ~ 190 pN under loading rates of ~10 4 -10 5 pN/sec, representing a ~4-fold increase in mechanostability as compared with low force N-terminal pulling. SMD simulations provided consistent tendencies in rupture forces, and through visualization of force propagation networks provided mechanistic insights. These results demonstrate how mechanostability of therapeutic protein-protein interfaces can be controlled by informed selection of anchor points within molecules, with implications for optimal bioconjugation strategies in drug delivery vehicles.

Overexpression of histone deacetylase (HDAC) enzymes is linked to a wide variety of illnesses, including malignancies and neurological disorders, which makes HDAC inhibitors potentially therapeutic. However, most HDAC inhibitors lack subclass or isoform selectivity, which can be dangerous. Featuring both enhanced selectivity and toxicity profiles, slow-binding HDAC inhibitors offer promising treatment options for a variety of disorders. Diseases like cardiac, neurodegenerative disorders and diabetes are mainly associated with the HDAC1, HDAC2 and HDAC3 enzymes. The AI-based virtual screening tool PyRMD is implemented to identify the potential inhibitors from ∼2 million compounds. Based on the IC 50 values, the top 10 compounds were selected for molecular docking. From the docking and ADMET study, the top-ranked three compounds were selected for molecular dynamics (MD) simulations. Further, to get more insights into the binding/unbinding mechanism of the ligand, we have employed the steered molecular dynamics (SMD) simulations. This study assists in developing Amber force field parameters for the HDAC1, HDAC2 and HDAC3 proteins and sheds light on the discovery of a potent drug. Our study suggests that hydroxamic acid derivative (i.e. referred to as Comp-1, CHEMBL600072) is the potential inhibitor for the series of HDAC-related diseases.Communicated by Ramaswamy H. Sarma.

The type 1 pilus is a bacterial filament consisting of a long coiled proteic chain of subunits joined together by non-covalent bonding between complementing β-strands. Its strength and structural stability are critical for its anchoring function in uropathogenic Escherichia coli bacteria. The pulling and unravelling of the FimG subunit of the pilus was recently studied by atomic force microscopy experiments and steered molecular dynamics simulations (Alonso-Caballero et al. 2018 Nat. Commun. 9, 2758. (doi:10.1038/s41467-018-05107-6)). In this work, we perform a quantitative comparison between experiment and simulation, showing a good agreement in the underlying work values for the unfolding. The simulation results are then used to estimate the free energy difference for the detachment of FimG from the complementing strand of the neighbouring subunit in the chain, FimF. Finally, we show that the large free energy difference for the unravelling and detachment of the subunits which leads to the high stability of the chain is entirely entropic in nature.

Unfolding pathways of T3 liganded thyroid hormone receptor (THRT3) can be studied by using the protocols of steered molecular dynamics (SMD). Theory of constant velocity pulling has been implemented to the structure of THRT3 in a neutral water-ion solution equilibrated up to 20 ns. The globular form of THRT3 is completely unfolded extending N-C termini from 38 Å to 876 Å at a constant speed of 0.1 Å/ps by means of 8.5 ns long SMD simulations. The peak force measured in the intermediate conformations is related to a burst of backbone H-bonds among a-helices and b-hairpins. With decrease in H-bonds, electrostatic energy increases by losing gradually the secondary structure and separating a and b-strands in solution. The force at the end (t > 8.5 ns) increases steeply with the large increase in bond-angle and bond-length potentials when the system becomes completely unfolded. The hydrophobic ligand binding domain (LBD) of THR-b with load bearing H-bonds protects T3 from water attack. Even after complete unfolding of THR-b LBD, the position of T3 is not deviated more than 2.5 Å and a large number of water molecules remain in the surrounding of this domain area. This is a strong evidence for the mechanochemical stability of a receptor protein’s LBD towards hormone activated gene expressions followed by ligand binding and dissociation. BIBECHANA 17 (2020) 50-57