Single-Molecule Investigation of the Binding Interface Stability of SARS-CoV-2 Variants with ACE2

The SARS-CoV-2 pandemic spurred numerous research endeavors to comprehend the virus and mitigate its global severity. Understanding the binding interface between the virus and human receptors is pivotal to these efforts and paramount to curbing infection and transmission. Here we employ atomic force microscopy and steered molecular dynamics simulation to explore SARS-CoV-2 receptor binding domain (RBD) variants and angiotensin-converting enzyme 2 (ACE2), examining the impact of mutations at key residues upon binding affinity. Our results show that the Omicron and Delta variants possess strengthened binding affinity in comparison to the Mu variant. Further, using sera from individuals either vaccinated or with acquired immunity following Delta strain infection, we assess the impact of immunity upon variant RBD/ACE2 complex formation. Single-molecule force spectroscopy analysis suggests that vaccination before infection may provide stronger protection across variants. These results underscore the need to monitor antigenic changes in order to continue developing innovative and effective SARS-CoV-2 abrogation strategies.

Exploration of novel hydroxamate zinc binding group inhibitors against HDAC-1-3 enzymes by AI-based virtual screening: atomistic insights from steered molecular dynamics

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 IC50 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 mechanism and energetics of the dynein priming stroke

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.

Investigation of Zika virus methyl transferase inhibitors using steered molecular dynamics

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

Application of temperature-dependent and steered molecular dynamics simulation to screen anti-dengue compounds against Marburg virus

Marburg virus infections are extremely fatal with a fatality range of 23% to 90%, therefore there is an urgent requirement to design and develop efficient therapeutic molecules. Here, a comprehensive temperature-dependent molecular dynamics (MD) simulation method was implemented to identify the potential molecule from the anti-dengue compound library that can inhibit the function of the VP24 protein of Marburg. Virtual high throughput screening identified five effective binders of VP24 after screening 484 anti-dengue compounds. These compounds were treated in MD simulation at four different temperatures: 300, 340, 380, and 420 K. Higher temperatures showed dissociation of hit compounds from the protein. Further, triplicates of 100 ns MD simulation were conducted which showed that compounds ID = 118717693, and ID = 5361 showed strong stability with the protein molecule. These compounds were further validated using Δ G binding free energies and they showed: -30.38 kcal/mol, and -67.83 kcal/mol binding free energies, respectively. Later, these two compounds were used in steered MD simulation to detect its dissociation. Compound ID = 5361 showed the maximum pulling force of 199.02 kcal/mol/nm to dissociate the protein-ligand complex while ID = 118717693 had a pulling force of 101.11 kcal/mol/nm, respectively. This ligand highest number of hydrogen bonds with varying occupancies at 89.93%, 69.80%, 57.93%, 52.33%, and 50.63%. This study showed that ID = 5361 can bind with the VP24 strongly and has the potential to inhibit its function which can be validated in the in-vitro experiment.Communicated by Ramaswamy H. Sarma.

Identification of synthetically tractable MERS-CoV main protease inhibitors using structure-based virtual screening and molecular dynamics potential of mean force (PMF) calculations

The Middle East Respiratory Syndrome Coronavirus (MERS-CoV) is a potentially lethal infection that presents a substantial threat to health, especially in Middle East nations. Given that no FDA-approved specific therapy for MERS infection exists, designing and discovering a potent antiviral therapy for MERS-CoV is crucial. One pivotal strategy for inhibiting MERS replication is to focus on the viral main protease (Mpro). In this study, we identify potential novel Mpro inhibitors employing structure-based virtual screening of our recently reported Ugi reaction-derived library (URDL) consisting of cherry-picked molecules from the literature. The key features of the URDL library include synthetic tractability (1-2 pot synthesis) of the molecules scaffold and unexplored chemical space. The hits were ranked based on the docking score, MM-GBSA free energy of binding, and the interaction pattern with the active site residues. A molecular dynamics (MD) simulation study was performed for the first two top-ranked compounds to analyze the stability and free binding energy based on the molecular mechanics Poisson-Boltzmann surface area. The potential mean force calculated from the steered molecular dynamics (SMD) simulations of the hits indicates improved H-bond potential, enhanced conformational stability, and binding affinity toward the target, compared to the cocrystallized ligand. The discovered hits represent novel synthetically tractable scaffolds as potential MERS-CoV Mpro inhibitors.Communicated by Ramaswamy H. Sarma.

A multi-tier computational screening framework to effectively search the mutational space of SARS-CoV-2 receptor binding motif to identify mutants with enhanced ACE2 binding abilities

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.

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

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

A Study of a Protein-Folding Machine: Transient Rotation of the Polypeptide Backbone Facilitates Rapid Folding of Protein Domains in All-Atom Molecular Dynamics Simulations

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.

Three structural solutions for bacterial adhesion pilus stability and superelasticity

Bacterial adhesion pili are key virulence factors that mediate host-pathogen interactions in diverse epithelial environments. Deploying a multimodal approach, we probed the structural basis underpinning the biophysical properties of pili originating from enterotoxigenic (ETEC) and uropathogenic bacteria. Using cryo-electron microscopy we solved the structures of three vaccine target pili from ETEC bacteria, CFA/I, CS17, and CS20. Pairing these and previous pilus structures with force spectroscopy and steered molecular dynamics simulations, we find a strong correlation between subunit-subunit interaction energies and the force required for pilus unwinding, irrespective of genetic similarity. Pili integrate three structural solutions for stabilizing their assemblies: layer-to-layer interactions, N-terminal interactions to distant subunits, and extended loop interactions from adjacent subunits. Tuning of these structural solutions alters the biophysical properties of pili and promotes the superelastic behavior that is essential for sustained bacterial attachment.

Steered Molecular Dynamics Simulations Study on FABP4 Inhibitors

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.

Static Binding and Dynamic Transporting-Based Design of Specific Ring-Chain-Ring Acetylcholinesterase Inhibitor: From Galantamine to Natural Product

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.

Using steered molecular dynamic tension for assessing quality of computational protein structure models

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.

Coarse-Grained Molecular Dynamics Simulations of Paclitaxel-Loaded Polymeric Micelles

A continuous manufacturing technology based on coaxial turbulent jet in coflow was previously developed to produce paclitaxel-loaded polymeric micelles. Herein, coarse-grained molecular dynamics (CG-MD) simulations were implemented to better understand the effect of the material attributes (i.e., the drug-polymer ratio and the ethanol concentration) and process parameters (i.e., temperature) on the self-assembly process of polymeric micelles as well as to provide molecular details on micelle instability. An all-atom (AA) poly (ethylene glycol)-poly (lactic acid) (PEG-PLA) polymer model was developed as the reference for parameterizing a coarse-grained (CG) model, and the AA polymer model was further validated with experimental glass transition temperature (T_g). The model transferability was verified by comparing structural properties between the AA and CG models. The CG model was further validated with experimental data, including micelle particle size measurements and drug encapsulation efficiency. Furthermore, the encapsulation of paclitaxel into the polymeric micelles was included in the simulations, taking into consideration the interactions between the paclitaxel and the polymers. The results from various points of view demonstrated a strong dependence of the shape of the micelles on the drug encapsulation, with micelles transitioning from spherical to ellipsoidal structures with an increasing paclitaxel amount. Simulation data were also used to identify the critical aggregation number (i.e., the number of polymer and drug molecules required for transition from one shape to another). Improved micellar structural stability was found with a larger micellar size and less solvent accessibility. Lastly, an evaluation was performed on the micellar dissociation free energy using a steered molecular dynamics simulation over a range of temperatures and ethanol concentrations. These simulations revealed that at higher ethanol and temperature conditions, micelles become destabilized, resulting in greater paclitaxel release. The increased drug release was determined to originate from the solvation of the hydrophobic core, which promoted micellar swelling and an associated reduction in hydrophobic interactions, leading to a loosely packed micellar structure.

Steered Molecular Dynamics of Lipid Membrane Indentation by Carbon and Silicon-Carbide Nanotubes-The Impact of Indenting Angle Uncertainty

Due to the semi-liquid nature and uneven morphologies of biological membranes, indentation may occur in a range of non-ideal conditions. These conditions are relatively unstudied and may alter the physical characteristics of the process. One of the basic challenges in the construction of nanoindenters is to appropriately align the nanotube tip and approach the membrane at a perpendicular angle. To investigate the impact of deviations from this ideal, we performed non-equilibrium steered molecular dynamics simulations of the indentation of phospholipid membranes by homogeneous CNT and non-homogeneous SiCNT indenters. We used various angles, rates, and modes of indentation, and the withdrawal of the relative indenter out of the membrane in corresponding conditions was simulated.

Structural Basis of Salicylic Acid Decarboxylase Reveals a Unique Substrate Recognition Mode and Access Channel

Salicylic acid (SA) decarboxylase from Trichosporon moniliiforme (TmSdc), which reversibly catalyses the decarboxylation of SA to yield phenol, is of significant interest because of its potential for the production of benzoic acid derivatives under environmentally friendly reaction conditions. TmSdc showed a preference for C2 hydroxybenzoate derivatives, with k_cat/K_m of SA being 3.2 × 10³ M^-1 s^-1. Here, we presented the first crystal structures of TmSdc, including a complex with SA. The three conserved residues of Glu8, His169, and Asp298 are the catalytic residues within the TIM-barrel domain of TmSdc. Trp239 forms a unique hydrophobic recognition site by interacting with the phenyl ring of SA, while Arg235 is responsible for recognizing the hydroxyl group at the C2 of SA via hydrogen bond interactions. Using a semi-rational combinatorial active-site saturation test, we obtained the TmSdc mutant MT3 (Y64T/P191G/F195V/E302D), which exhibited a 26.4-fold increase in k_cat/K_m with SA, reaching 8.4 × 10⁴ M^-1 s^-1. Steered molecular dynamics simulations suggested that the structural changes in MT3 relieved the steric hindrance within the substrate access channel and enlarged the substrate-binding pocket, leading to the increased activity by improving substrate access. Our data elucidate the unique substrate recognition mode and the substrate entrance tunnel of SA decarboxylase.

Impact of an alpha helix and a cysteine-cysteine disulfide bond on the resistance of bacterial adhesion pili to stress

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.

Helix-coil transition and conformational deformity in Aβ42-monomer: a case study using the Zn2+ cation

The metal ions (like Fe2+, Zn²⁺, Cu2+) are known to influence the amyloid beta (Aβ) aggregation. In this study, we have examined the conformational and dynamical changes during the coordination of Aβ-monomer with the Zn2+ ion using all-atom molecular dynamics (MD) simulations using explicit solvent models. We have probed the unfolding of the full-length Aβ₄₂ monomer both inclusive and exclusive of the Zn2+ cation, with 1:1 ratio of the peptide and the Zn2+ cation. The inclusion of the Zn2+ cation shows differential intra-peptide interactions which has been probed using various analyses. The Helix - Coil transition of the wild type Aβ42 monomer is studied using the steered molecular dynamics simulations by taking the end-to-end C-α distance across the peptide. This gives an idea of the unequal intra - peptide and peptide - water interactions being found across the length of the Aβ monomer. The transition of an α-helix dominated wild-type (WT) Aβ structure to the unfolded coil structure gives significant evidence of the intra-peptide hydrogen bonding shifts in the presence of the Zn2+ cation. This accounts for the structural and the dynamical variations that take place in the Aβ monomer in the presence of the Zn2+ cation to mimic the conditions/environment at the onset of fibrillation.Communicated by Ramaswamy H. Sarma.

Qualitative Prediction of Ligand Dissociation Kinetics from Focal Adhesion Kinase Using Steered Molecular Dynamics

Most early-stage drug discovery projects focus on equilibrium binding affinity to the target alongside selectivity and other pharmaceutical properties. Since many approved drugs have nonequilibrium binding characteristics, there has been increasing interest in optimizing binding kinetics early in the drug discovery process. As focal adhesion kinase (FAK) is an important drug target, we examine whether steered molecular dynamics (SMD) can be useful for identifying drug candidates with the desired drug-binding kinetics. In simulating the dissociation of 14 ligands from FAK, we find an empirical power–law relationship between the simulated time needed for ligand unbinding and the experimental rate constant for dissociation, with a strong correlation depending on the SMD force used. To improve predictions, we further develop regression models connecting experimental dissociation rate with various structural and energetic quantities derived from the simulations. These models can be used to predict dissociation rates from FAK for related compounds.

Mechanochemical Cellular Membrane Internalization of Nanohydrogels: A Large-Scale Mesoscopic Simulation

By combining large-scale dissipative particle dynamics and steered molecular dynamics simulations, we investigate the mechanochemical cellular internalization pathways of homogeneous and heterogeneous nanohydrogels and demonstrate that membrane internalization is determined by the crosslink density and encapsulation ability of nanohydrogels. The homogeneous nanohydrogels with a high crosslink density and low encapsulation ability behave as soft nanoparticles partially wrapped by the membrane, while those with a low crosslink density and high encapsulation ability permeate into the membrane. Regardless of the crosslink density, the homogeneous nanohydrogels undergo typical dual morphological deformations. The local lipid nanodomains are identified at the contacting region between the membrane and nanohydrogels because of different diffusion behaviors between lipid and receptor molecules during the internalization process. The yolk@shell heterogeneous nanohydrogels present a different mechanochemical cellular internalization pathway. The yolk with strong affinity is directly in contact with the membrane, resulting in partial membrane wrapping, and the contacting area is much reduced when compared to homogenous nanohydrogels, leading to a smaller lipid nanodomain and thus avoiding related cellular toxicity. Our findings provide a critical mechanism understanding of the biological pathways of nanohydrogels and may guide the molecular design of the hydrogel-based materials for controlled release drug delivery, tissue engineering, and cell culture.

High-throughput virtual screening of drug databanks for potential inhibitors of SARS-CoV-2 spike glycoprotein

COVID-19, which is caused by a novel coronavirus known as SARS-CoV-2, has spread rapidly around the world, and it has infected more than 29 million individuals as recorded on 16 September 2020. Much effort has been made to stop the virus from spreading, and there are currently no approved pharmaceutical products to treat COVID-19. Here, we apply an in silico approach to investigate more than 3800 FDA approved drugs on the viral RBD S₁-ACE2 interface as a target. The compounds were investigated through flexible ligand docking, ADME property calculations and protein–ligand interaction maps. Molecular dynamics (MD) simulations were also performed on eleven compounds to study the stability and the interactions of the protein–ligand complexes. The MD simulations show that bagrosin, chidamide, ebastine, indacaterol, regorafenib, salazosulfadimidine, silodosin and tasosartan are relatively stable near the C terminal domain (CTD1) of the S₁ subunit of the viral S protein. The relative MMGBSA binding energies show that silodosin has the best binding to the target. The constant velocity steered molecular dynamics (SMD) simulations show that silodosin preferentially interacts with the RBD S₁ and has potential to act as an interfering compound between viral spike–host ACE2 interactions.

Communicated by Ramaswamy H. Sarma

Molecular dynamics simulations provide molecular insights into the role of HLA-B51 in Behçet's disease pathogenesis

Behçet's disease is an inflammatory disorder of unknown etiology. Genetic tendency has an important role in its pathogenesis, and HLA-B51, a class I MHC antigen, has been recognized as the strongest susceptibility factor for Behçet's disease. Despite the confirmation of the association of HLA-B51 with Behçet's disease in different populations, its pathogenic mechanisms remain elusive. HLA-B51 differs in only two amino acids from HLA-B52, other split antigen of HLA-B5, which is not associated with Behçet's disease. These two amino acids are located in the B pocket of the antigen-binding groove, which occupies the second amino acids of the bound peptides. To understand the nature of the HLA-peptide interactions, differences in structure and dynamics of two HLA alleles were investigated by molecular dynamics simulations using YAYDGKDYI, LPRSTVINI, and IPYQDLPHL peptides. For HLA-B51, all bound peptides fluctuated to larger extent than HLA-B52. Free energy profiles of unbinding process for YAYDGKDYI by steered molecular dynamics simulations showed that unbinding from HLA-B52 results in greater free energy differences than HLA-B51. These results suggest the possibility of an instability of HLA-B51 associated with the repertoire of peptides, and this finding may provide significant insight to its pathogenic role in Behçet's disease.

The Hydrophobic Ligands Entry and Exit from the GPCR Binding Site-SMD and SuMD Simulations

Most G protein-coupled receptors that bind the hydrophobic ligands (lipid receptors and steroid receptors) belong to the most populated class A (rhodopsin-like) of these receptors. Typical examples of lipid receptors are: rhodopsin, cannabinoid (CB), sphingosine-1-phosphate (S1P) and lysophosphatidic (LPA) receptors. The hydrophobic ligands access the receptor binding site from the lipid bilayer not only because of their low solubility in water but also because of a large N-terminal domain plug preventing access to the orthosteric binding site from the extracellular milieu. In order to identify the most probable ligand exit pathway from lipid receptors CB1, S1P1 and LPA1 orthosteric binding sites we performed at least three repeats of steered molecular dynamics simulations in which ligands were pulled in various directions. For specific ligands being agonists, the supervised molecular dynamics approach was used to simulate the ligand entry events to the inactive receptor structures. For all investigated receptors the ligand entry/exit pathway goes through the gate between transmembrane helices TM1 and TM7, however, in some cases it combined with a direction toward water milieu.

Entropic bonding of the type 1 pilus from experiment and simulation

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.

Deformation Mechanisms of "Two-Part" Natural Adhesive in Bone Interfibrillar Nano-Interfaces

Noncollagenous proteins at nanoscale interfaces in bone are less than 2-3% of bone content by weight, while they contribute more than 30% to fracture toughness. Major gaps in quantitative understanding of noncollagenous proteins' role in the interfibrillar interfaces, largely because of the limitation of probing their nanoscale dimension, have resulted in ongoing controversies and several outstanding hypotheses on their role and function, arguably going back to centuries ago to the original work from Galileo. Our results from the first detailed computational model of the nano-interface in the bone reveal "synergistic" deformation mechanism of a "double-part" natural glue, that is, noncollagenous osteopontin and osteocalcin at the interfibrillar interface. Specifically, through strong anchoring and formation of dynamic binding sites on mineral nanoplatelets, the nano-interface can sustain a large nonlinear deformation with ductility approaching 5000%. This large deformation results in an outstanding specific energy to failure exceeding ∼350 J/g, which is larger than the most known tough materials (such as Kevlar, spider silk, and so forth.).

Molecular dynamics simulations to the bidirectional adhesion signaling pathway of integrin αV β3

The bidirectional force transmission process of integrin through the cell membrane is still not well understood. Several possible mechanisms have been discussed in literature on the basis of experimental data, and in this study, we investigate these mechanisms by free and steered molecular dynamics simulations. For the first time, constant velocity pulling on the complete integrin molecule inside a dipalmitoyl-phosphatidylcholine membrane is conducted. From the results, the most likely mechanism for inside-out and outside-in signaling is the switchblade model with further separation of the transmembrane helices.

[Steered Molecular Dynamics Simulation Study of Quantified Effects of Point Mutation Induced by Breast Cancer on Mechanical Behavior of E-Cadherin]

E-cadherin is a member of the cadherin family that plays a key role in the formation of cell-cell adhesion among epithelial tissues. Point mutations are one of the structural abnormalities of E-cadherin in human carcinomas. Such abnormalities can alter mechanical properties of proteins that play an important role in their biological activities. To determine the impact of point mutations on protein mechanical properties, the second fragment of extracellular domain of E-cadherin was modeled using steered molecular dynamics simulations. The molecular dynamics modeling included application of tensile forces in both constant velocity and constant force modes to examine the effects of Met282 to He and Asn315 to Ser mutations on mechanical behavior of protein structure. The stabilities of the wild type and mutant structures were also obtained by the protein design foldX algorithm. Results confirmed the lower stability of the mutant domains compared to the wild type. The mutated proteins displayed softer behavior than the reference protein and their stiffness decreased by up to 34%. Our findings suggest that local changes in molecular structure due to mutations may lead to noticeable alterations in mechanical properties within the entire domain. Since the function of protein is related to its structure, these changes may influence the function of the protein.

Connecting the dots in the mechanism of action of Cucurbitacin E (CurE) - path analysis and steered molecular dynamics reveal the precise site of entry and the passage of CurE in filamentous actin

Cucurbitacin E (CurE) modulates actin cytoskeleton by forming an irreversible covalent bond with Cys257 of actin. The reported binding conformation of CurE is deeply buried in the subdomain 4 of actin and is closely situated to the ATP-binding site. The entry and the path taken by CurE to reach this buried site remain a mystery. In this study, steered molecular dynamics (SMD) simulations were conducted to delineate the diffusion of CurE to its binding site. SMD simulations reveal that the distinctive entry site of CurE found in subdomain 4 is by itself a closed and compact region of two loops lying beside each other like a closed door and CurE induces it to open. From this point, CurE moves toward its binding site through a path facilitated by Thr188, Leu261, Ile267, Ile309, Tyr306. This study is also an insight into how CurE distinctly differentiates its tunnel to the binding site from the ATP-binding site. The conformational changes of CurE along the path to the binding site are surprisingly very minimal and closely resemble the attack conformation at the end of simulation. The study reveals that the little energy spent by the molecule is compensated by the enthalpic contribution to binding-free energy barrier making it undoubtedly the most preferred path of CurE. This study is the first of its kind in which the SMD was used to derive the complete and continuous translocation of one of the most potent phytochemicals, CurE through the binding site gorge of actin. [Formula: see text] Communicated by Ramaswamy H. Sarma.

Computational Study for the Unbinding Routes of β-N-Acetyl-d-Hexosaminidase Inhibitor: Insight from Steered Molecular Dynamics Simulations

β-N-Acetyl-d-hexosaminidase from Ostrinia furnacalis (OfHex1) is a new target for the design of insecticides. Although some of its inhibitors have been found, there is still no commercial drug available at present. The residence time of the ligand may be important for its pharmacodynamic effect. However, the unbinding routes of ligands from OfHex1 still remain largely unexplored. In the present study, we first simulated the six dissociation routes of N,N,N-trimethyl-d-glucosamine-chitotriomycin (TMG-chitotriomycin, a highly selective inhibitor of OfHex1) from the active pocket of OfHex1 by steered molecular dynamics simulations. By comparing the potential of mean forces (PMFs) of six routes, Route 1 was considered as the most possible route with the lowest energy barrier. Furthermore, the structures of six different states for Route 1 were snapshotted, and the key amino acid residues affecting the dissociated time were analyzed in the unbinding pathway. Moreover, we also analyzed the "open⁻close" mechanism of Glu368 and Trp448 and found that their conformational changes directly affected the dissociation of TMG-chitotriomycin. Our findings would be helpful to understanding and identifying novel inhibitors against OfHex1 from virtual screening or lead-optimization.

Lessons learned about steered molecular dynamics simulations and free energy calculations

The calculation of free energy profiles is central in understanding differential enzymatic activity, for instance, involving chemical reactions that require QM-MM tools, ligand migration, and conformational rearrangements that can be modeled using classical potentials. The use of steered molecular dynamics (sMD) together with the Jarzynski equality is a popular approach in calculating free energy profiles. Here, we first briefly review the application of the Jarzynski equality to sMD simulations, then revisit the so-called stiff-spring approximation and the consequent expectation of Gaussian work distributions and, finally, reiterate the practical utility of the second-order cumulant expansion, as it coincides with the parametric maximum-likelihood estimator in this scenario. We illustrate this procedure using simulations of CO, both in aqueous solution and in a carbon nanotube as a model system for biologically relevant nanoheterogeneous environments. We conclude the use of the second-order cumulant expansion permits the use of faster pulling velocities in sMD simulations, without introducing bias due to large dispersion in the non-equilibrium work distribution.

BAR-based optimum adaptive steered MD for configurational sampling

The equilibrium and nonequilibrium adaptive alchemical free energy simulation methods optimum Bennett's acceptance ratio and optimum crooks' equation (OCE), based on the statistically optimal bidirectional reweighting estimator named Bennett's Acceptance Ratio or Crooks' equation, perform initial sampling in the staging alchemical transformation and then determine the importance rank of different states via the time-derivative of the variance. The method is proven to give speedups compared with the equal time rule. In the current work, we extend the time derivative of variance guided adaptive sampling method to the configurational space, falling in the term of steered MD (SMD). The SMD approach biasing physically meaningful collective variable (CV) such as one dihedral or one distance to pulling the system from one conformational state to another. By minimizing the variance of the free energy differences along the pathway in an optimized way, a new type of adaptive SMD (ASMD) is introduced. As exhibits in the alchemical case, this adaptive sampling method outperforms the traditional equal-time SMD in nonequilibrium stratification. Also, the method gives much more efficient calculation of potential of mean force than the selection criterion-based ASMD scheme, which is proven to be more efficient than traditional SMD. The OCE workflow is periodicity-of-CV dependent while ASMD is not. The performance is demonstrated in a dihedral flipping case and two distance pulling cases, accounting for periodic and nonperiodic CVs, respectively. © 2019 Wiley Periodicals, Inc.

Role of Resultant Dipole Moment in Mechanical Dissociation of Biological Complexes

Protein-peptide interactions play essential roles in many cellular processes and their structural characterization is the major focus of current experimental and theoretical research. Two decades ago, it was proposed to employ the steered molecular dynamics (SMD) to assess the strength of protein-peptide interactions. The idea behind using SMD simulations is that the mechanical stability can be used as a promising and an efficient alternative to computationally highly demanding estimation of binding affinity. However, mechanical stability defined as a peak in force-extension profile depends on the choice of the pulling direction. Here we propose an uncommon choice of the pulling direction along resultant dipole moment (RDM) vector, which has not been explored in SMD simulations so far. Using explicit solvent all-atom MD simulations, we apply SMD technique to probe mechanical resistance of ligand-receptor system pulled along two different vectors. A novel pulling direction—when ligand unbinds along the RDM vector—results in stronger forces compared to commonly used ligand unbinding along center of masses vector. Our observation that RDM is one of the factors influencing the mechanical stability of protein-peptide complex can be used to improve the ranking of binding affinities by using mechanical stability as an effective scoring function.

Docking, steered molecular dynamics, and QSAR studies as strategies for studying isoflavonoids as 5-, 12-, and 15-lipoxygenase inhibitors

Lipoxygenases (LOX) are enzymes that catalyze polyunsaturated fatty acid peroxidation and have a non-heme iron atom located in their active site. They are implicated in the arachidonic acid pathway and involved in inflammation, fever, pain production, and in the origins of several diseases such as cancer, asthma, and psoriasis. The search for inhibitors of these enzymes has emerged in the last years, and isoflavonoids have a broad spectrum of biological activity with low cytotoxicity. Our previous results have shown that isoflavonoids inhibited different LOX isoforms in vitro. For this reason, we studied the most important interactions that govern the potency and selectivity of some isoflavones and isoflavans toward different LOX isoforms using computational methods. The docking results have shown that all the molecules can be located in different zones in the LOX active site. Steered molecular dynamics indicated that selectivity was present at the cavity entry, but not at its exit. We also observed the correlation between the potential mean force and the best (HIR-303) and worst inhibitors (IR-213) in 5-LOX. Finally, structure-activity relationship (QSAR) studies showed a good correlation between theoretical IC₅₀ values and experimental data for 5-LOX and 12-LOX with 96 and 95%, respectively, and a lower correlation for 15-LOX (79%). Conclusively, pharmacophore analysis showed that our proposed molecules should possess a donor-acceptor and aromatic centers to encourage interactions in the active site.

Chloride Ion Transport by the E. coli CLC Cl-/H+ Antiporter: A Combined Quantum-Mechanical and Molecular-Mechanical Study

We performed steered molecular dynamics (SMD) and umbrella sampling simulations of Cl^- ion migration through the transmembrane domain of a prototypical E. coli CLC Cl^-/H⁺ antiporter by employing combined quantum-mechanical (QM) and molecular-mechanical (MM) calculations. The SMD simulations revealed interesting conformational changes of the protein. While no large-amplitude motions of the protein were observed during pore opening, the side chain rotation of the protonated external gating residue Glu148 was found to be critical for full access of the channel entrance by Cl^-. Moving the anion into the external binding site (S_ext) induced small-amplitude shifting of the protein backbone at the N-terminal end of helix F. As Cl^- traveled through the pore, rigid-body swinging motions of helix R separated it from helix D. Helix R returned to its original position once Cl^- exited the channel. Population analysis based on polarized wavefunction from QM/MM calculations discovered significant (up to 20%) charge loss for Cl^- along the ion translocation pathway inside the pore. The delocalized charge was redistributed onto the pore residues, especially the functional groups containing π bonds (e.g., the Tyr445 side chain), while the charges of the H atoms coordinating Cl^- changed almost negligibly. Potentials of mean force computed from umbrella sampling at the QM/MM and MM levels both displayed barriers at the same locations near the pore entrance and exit. However, the QM/MM PMF showed higher barriers (~10 kcal/mol) than the MM PMF (~2 kcal/mol). Binding energy calculations indicated that the interactions between Cl^- and certain pore residues were overestimated by the semi-empirical PM3 Hamiltonian and underestimated by the CHARMM36 force fields, both of which were employed in the umbrella sampling simulations. In particular, CHARMM36 underestimated binding interactions for the functional groups containing π bonds, missing the stabilizations of the Cl^- ion due to electron delocalization. The results suggested that it is important to explore these quantum effects for accurate descriptions of the Cl^- transport.

Prospecting and Structural Insight into the Binding of Novel Plant-Derived Molecules of Leea indica as Inhibitors of BACE1

BACKGROUND

Alzheimer disease (AD) can be considered as the most common age related neurodegenerative disorder and also an important cause of death in elderly patients. A number of studies showed the correlation of this disease pathology with BACE1 inhibitor and it is also evident that BACE1 inhibitor can function as a very potent strategy in treating AD.

METHODS

In this present study, we aimed to prospect for novel plant-derived BACE1 inhibitors from Leea indica and to realise structural basis of their interactions and mechanisms using combined molecular docking and molecular dynamics based approaches. An extensive library of Leea indica plant derived molecule was compiled and computationally screened for inhibitory action against BACE1 by using virtual screening approaches. Furthermore, induced fit docking and classical molecular dynamics along with steered molecular dynamics simulations were employed to get insight of the binding mechanisms.

RESULTS

Two triterpenoids, ursolic acid and lupeol were identified through virtual screening; wherein, lupeol showed better binding free energy in MM/GBSA, MM/PBSA and MM/GBVI approaches. Furthermore classical and steered dynamics revealed the favourable hydrophobic interactions between the lupeol and the residues of flap or catalytic dyadof BACE1; however, ursolic acid showed disfavorable interactions with the BACE1.

CONCLUSION

This study therefore unveiled the potent BACE1 inhibitor from a manually curated dataset of Leea indica molecules, which may provide a novel dimension of designing novel BACE1 inhibitors for AD therapy.

Combining in Vitro and in Silico Single-Molecule Force Spectroscopy to Characterize and Tune Cellulosomal Scaffoldin Mechanics

Cellulosomes are polyprotein machineries that efficiently degrade cellulosic material. Crucial to their function are scaffolds consisting of highly homologous cohesin domains, which serve a dual role by coordinating a multiplicity of enzymes as well as anchoring the microbe to its substrate. Here we combined two approaches to elucidate the mechanical properties of the main scaffold ScaA of Acetivibrio cellulolyticus. A newly developed parallelized one-pot in vitro transcription-translation and protein pull-down protocol enabled high-throughput atomic force microscopy (AFM)-based single-molecule force spectroscopy (SMFS) measurements of all cohesins from ScaA with a single cantilever, thus promising improved relative force comparability. Albeit very similar in sequence, the hanging cohesins showed considerably lower unfolding forces than the bridging cohesins, which are subjected to force when the microbe is anchored to its substrate. Additionally, all-atom steered molecular dynamics (SMD) simulations on homology models offered insight into the process of cohesin unfolding under force. Based on the differences among the individual force propagation pathways and their associated correlation communities, we designed mutants to tune the mechanical stability of the weakest hanging cohesin. The proposed mutants were tested in a second high-throughput AFM SMFS experiment revealing that in one case a single alanine to glycine point mutation suffices to more than double the mechanical stability. In summary, we have successfully characterized the force induced unfolding behavior of all cohesins from the scaffoldin ScaA, as well as revealed how small changes in sequence can have large effects on force resilience in cohesin domains. Our strategy provides an efficient way to test and improve the mechanical integrity of protein domains in general.

Modeling the human Nav1.5 sodium channel: structural and mechanistic insights of ion permeation and drug blockade

Abnormalities in the human Nav1.5 (hNav1.5) voltage-gated sodium ion channel (VGSC) are associated with a wide range of cardiac problems and diseases in humans. Current structural models of hNav1.5 are still far from complete and, consequently, their ability to study atomistic interactions of this channel is very limited. Here, we report a comprehensive atomistic model of the hNav1.5 ion channel, constructed using homology modeling technique and refined through long molecular dynamics simulations (680 ns) in the lipid membrane bilayer. Our model was comprehensively validated by using reported mutagenesis data, comparisons with previous models, and binding to a panel of known hNav1.5 blockers. The relatively long classical MD simulation was sufficient to observe a natural sodium permeation event across the channel's selectivity filters to reach the channel's central cavity, together with the identification of a unique role of the lysine residue. Electrostatic potential calculations revealed the existence of two potential binding sites for the sodium ion at the outer selectivity filters. To obtain further mechanistic insight into the permeation event from the central cavity to the intracellular region of the channel, we further employed "state-of-the-art" steered molecular dynamics (SMD) simulations. Our SMD simulations revealed two different pathways through which a sodium ion can be expelled from the channel. Further, the SMD simulations identified the key residues that are likely to control these processes. Finally, we discuss the potential binding modes of a panel of known hNav1.5 blockers to our structural model of hNav1.5. We believe that the data presented here will enhance our understanding of the structure-property relationships of the hNav1.5 ion channel and the underlying molecular mechanisms in sodium ion permeation and drug interactions. The results presented here could be useful for designing safer drugs that do not block the hNav1.5 channel.

Unraveling Hydrophobic Interactions at the Molecular Scale Using Force Spectroscopy and Molecular Dynamics Simulations

Interactions between hydrophobic moieties steer ubiquitous processes in aqueous media, including the self-organization of biologic matter. Recent decades have seen tremendous progress in understanding these for macroscopic hydrophobic interfaces. Yet, it is still a challenge to experimentally measure hydrophobic interactions (HIs) at the single-molecule scale and thus to compare with theory. Here, we present a combined experimental-simulation approach to directly measure and quantify the sequence dependence and additivity of HIs in peptide systems at the single-molecule scale. We combine dynamic single-molecule force spectroscopy on model peptides with fully atomistic, both equilibrium and nonequilibrium, molecular dynamics (MD) simulations of the same systems. Specifically, we mutate a flexible (GS)₅ peptide scaffold with increasing numbers of hydrophobic leucine monomers and measure the peptides' desorption from hydrophobic self-assembled monolayer surfaces. Based on the analysis of nonequilibrium work-trajectories, we measure an interaction free energy that scales linearly with 3.0-3.4 k_BT per leucine. In good agreement, simulations indicate a similar trend with 2.1 k_BT per leucine, while also providing a detailed molecular view into HIs. This approach potentially provides a roadmap for directly extracting qualitative and quantitative single-molecule interactions at solid/liquid interfaces in a wide range of fields, including interactions at biointerfaces and adhesive interactions in industrial applications.

Characterizing Structural Stability of Amyloid Motif Fibrils Mediated by Water Molecules

In biological systems, structural confinements of amyloid fibrils can be mediated by the role of water molecules. However, the underlying effect of the dynamic behavior of water molecules on structural stabilities of amyloid fibrils is still unclear. By performing molecular dynamics simulations, we investigate the dynamic features and the effect of interior water molecules on conformations and mechanical characteristics of various amyloid fibrils. We find that a specific mechanism induced by the dynamic properties of interior water molecules can affect diffusion of water molecules inside amyloid fibrils, inducing their different structural stabilities. The conformation of amyloid fibrils induced by interior water molecules show the fibrils' different mechanical features. We elucidate the role of confined and movable interior water molecules in structural stabilities of various amyloid fibrils. Our results offer insights not only in further understanding of mechanical features of amyloids as mediated by water molecules, but also in the fine-tuning of the functional abilities of amyloid fibrils for applications.

Acidity Constant (pKa ) Calculation of Large Solvated Dye Molecules: Evaluation of Two Advanced Molecular Dynamics Methods

pH‐Sensitive dyes are increasingly applied on polymer substrates for the creation of novel sensor materials. Recently, these dye molecules were modified to form a covalent bond with the polymer host. This had a large influence on the pH‐sensitive properties, in particular on the acidity constant (pK_a). Obtaining molecular control over the factors that influence the pK_a value is mandatory for the future intelligent design of sensor materials. Herein, we show that advanced molecular dynamics (MD) methods have reached the level at which the pK_a values of large solvated dye molecules can be predicted with high accuracy. Two MD methods were used in this work: steered or restrained MD and the insertion/deletion scheme. Both were first calibrated on a set of phenol derivatives and afterwards applied to the dye molecule bromothymol blue. Excellent agreement with experimental values was obtained, which opens perspectives for using these methods for designing dye molecules.

All Subdomains of the Talin Rod Are Mechanically Vulnerable and May Contribute To Cellular Mechanosensing

Although the relevance of mechanotransduction in cell signaling is currently appreciated, the mechanisms that drive this process remain largely unknown. Mechanical unfolding of proteins may trigger distinct downstream signals in cells, providing a mechanism for cellular mechanotransduction. Force-induced unfolding of talin, a prominent focal adhesion protein, has been demonstrated previously for a small portion of its rod domain. Here, using single-molecule atomic force microscopy (smAFM), we show that the entire talin rod can be unfolded by mechanical extension, over a physiological range of forces between 10 and 40 pN. We also demonstrate, through a combination of smAFM and steered molecular dynamics, that the different bundles within the talin rod exhibit a distinct hierarchy of mechanical stability. These results provide a mechanism by which different force conditions within the cell control a graduated unfolding of the talin rod. Mechanical unfolding of the rod subdomains, and the subsequent effect on talin's binding interactions, would allow for a finely tuned cellular response to internally or externally applied forces.

Ranking protein-protein docking results using steered molecular dynamics and potential of mean force calculations

Crystallization of protein-protein complexes can often be problematic and therefore computational structural models are often relied on. Such models are often generated using protein-protein docking algorithms, where one of the main challenges is selecting which of several thousand potential predictions represents the most near-native complex. We have developed a novel technique that involves the use of steered molecular dynamics (sMD) and umbrella sampling to identify near-native complexes among protein-protein docking predictions. Using this technique, we have found a strong correlation between our predictions and the interface RMSD (iRMSD) in ten diverse test systems. On two of the systems, we investigated if the prediction results could be further improved using potential of mean force calculations. We demonstrated that a near-native (<2.0 Å iRMSD) structure could be identified in the top-1 ranked position for both systems. © 2016 Wiley Periodicals, Inc.

Atomic-Scale Sliding Friction on Graphene in Water

The sliding of a sharp nanotip on graphene completely immersed in water is investigated by molecular dynamics (MD) and atomic force microscopy. MD simulations predict that the atomic-scale stick-slip is almost identical to that found in ultrahigh vacuum. Furthermore, they show that water plays a purely stochastic role in sliding (solid-to-solid) friction. These observations are substantiated by friction measurements on graphene grown on Cu and Ni, where, oppositely of the operation in air, lattice resolution is readily achieved. Our results promote friction force microscopy in water as a robust alternative to ultra-high-vacuum measurements.

Steered molecular dynamics study reveals insights into the function of the repetitive B region of collagen- and fibrinogen-binding MSCRAMMs

MSCRAMMs (microbial surface components recognizing adhesive matrix molecules) are modular proteins covalently anchored in the bacterial cell wall of many Gram-positive bacteria. The N-terminal region of most MSCRAMMs carries the ligand-binding domains (A region) which specifically target the host extracellular matrix (ECM) proteins such as collagen, fibrinogen and fibronectin. In Staphylococcus aureus Cna, the prototype collagen-binding MSCRAMM, the A region is followed by a repetitive B region which is found to be conserved among many Gram-positive bacteria. This conservation signifies an important functional role for the B region which is made of repetitive domains. It was suggested that this region could act as a 'stalk' as well as a 'spring' to present the ligand-binding A region, away from the bacterial surface. But there is no clear functional implication of this region available till date. Each repetitive domain in the B region possesses a variant of the Ig fold called the CnaB fold. Additionally, the B repeats are also paired and the pairs are clustered together. To investigate if the B domains have a function similar to the Ig domains in the I-band region of the giant muscle protein, titin, steered molecular dynamics simulations of one, two and four B repeats of Cna were carried out. The results of the simulations suggest that the B region could provide mechanical stability, extensibility and elasticity to Cna due to the CnaB fold as well as the clustered arrangement of their domains. This study thus provided further insights into the biological underpinnings of adhesin-host interaction.

Atomic resolution probe for allostery in the regulatory thin filament

Calcium binding and dissociation within the cardiac thin filament (CTF) is a fundamental regulator of normal contraction and relaxation. Although the disruption of this complex, allosterically mediated process has long been implicated in human disease, the precise atomic-level mechanisms remain opaque, greatly hampering the development of novel targeted therapies. To address this question, we used a fully atomistic CTF model to test both Ca(2+) binding strength and the energy required to remove Ca(2+) from the N-lobe binding site in WT and mutant troponin complexes that have been linked to genetic cardiomyopathies. This computational approach is combined with measurements of in vitro Ca(2+) dissociation rates in fully reconstituted WT and cardiac troponin T R92L and R92W thin filaments. These human disease mutations represent known substitutions at the same residue, reside at a significant distance from the calcium binding site in cardiac troponin C, and do not affect either the binding pocket affinity or EF-hand structure of the binding domain. Both have been shown to have significantly different effects on cardiac function in vivo. We now show that these mutations independently alter the interaction between the Ca(2+) ion and cardiac troponin I subunit. This interaction is a previously unidentified mechanism, in which mutations in one protein of a complex indirectly affect a third via structural and dynamic changes in a second to yield a pathogenic change in thin filament function that results in mutation-specific disease states. We can now provide atom-level insight that is potentially highly actionable in drug design.

Mapping Mechanical Force Propagation through Biomolecular Complexes

Here we employ single-molecule force spectroscopy with an atomic force microscope (AFM) and steered molecular dynamics (SMD) simulations to reveal force propagation pathways through a mechanically ultrastable multidomain cellulosome protein complex. We demonstrate a new combination of network-based correlation analysis supported by AFM directional pulling experiments, which allowed us to visualize stiff paths through the protein complex along which force is transmitted. The results implicate specific force-propagation routes nonparallel to the pulling axis that are advantageous for achieving high dissociation forces.

A Self-Adaptive Steered Molecular Dynamics Method Based on Minimization of Stretching Force Reveals the Binding Affinity of Protein-Ligand Complexes

Binding affinity prediction of protein–ligand complexes has attracted widespread interest. In this study, a self-adaptive steered molecular dynamics (SMD) method is proposed to reveal the binding affinity of protein–ligand complexes. The SMD method is executed through adjusting pulling direction to find an optimum trajectory of ligand dissociation, which is realized by minimizing the stretching force automatically. The SMD method is then used to simulate the dissociations of 19 common protein–ligand complexes which are derived from two homology families, and the binding free energy values are gained through experimental techniques. Results show that the proposed SMD method follows a different dissociation pathway with lower a rupture force and energy barrier when compared with the conventional SMD method, and further analysis indicates the rupture forces of the complexes in the same protein family correlate well with their binding free energy, which reveals the possibility of using the proposed SMD method to identify the active ligand.

Influence of Solvent in Controlling Peptide-Surface Interactions

Protein binding to surfaces is an important phenomenon in biology and in modern technological applications. Extensive experimental and theoretical research has been focused in recent years on revealing the factors that govern binding affinity to surfaces. Theoretical studies mainly focus on examining the contribution of the individual amino acids or, alternatively, the binding potential energies of the full peptide, which are unable to capture entropic contributions and neglect the dynamic nature of the system. We present here a methodology that involves the combination of nonequilibrium dynamics simulations with strategic mutation of polar residues to reveal the different factors governing the binding free energy of a peptide to a surface. Using a gold-binding peptide as an example, we show that relative binding free energies are a consequence of the balance between strong interactions of the peptide with the surface and the ability for the bulk solvent to stabilize the peptide.

Energetic evaluation of binding modes in the C3d and Factor H (CCP 19-20) complex

As a part of innate immunity, the complement system relies on activation of the alternative pathway (AP). While feed-forward amplification generates an immune response towards foreign surfaces, the process requires regulation to prevent an immune response on the surface of host cells. Factor H (FH) is a complement protein secreted by native cells to negatively regulate the AP. In terms of structure, FH is composed of 20 complement-control protein (CCP) modules that are structurally homologous but vary in composition and function. Mutations in these CCPs have been linked to states of autoimmunity. In particular, several mutations in CCP 19-20 are correlated to atypical hemolytic uremic syndrome (aHUS). From crystallographic structures there are three putative binding sites of CCP 19-20 on C3d. Since there has been some controversy over the primary mode of binding from experimental studies, we approach characterization of binding using computational methods. Specifically, we compare each binding mode in terms of electrostatic character, structural stability, dissociative and associative properties, and predicted free energy of binding. After a detailed investigation, we found two of the three binding sites to be similarly stable while varying in the number of contacts to C3d and in the energetic barrier to complex dissociation. These sites are likely physiologically relevant and may facilitate multivalent binding of FH CCP 19-20 to C3b and either C3d or host glycosaminoglycans. We propose thermodynamically stable binding with modules 19 and 20, the latter driven by electrostatics, acting synergistically to increase the apparent affinity of FH for host surfaces.

Enzymatic Studies of Isoflavonoids as Selective and Potent Inhibitors of Human Leukocyte 5-Lipo-Oxygenase

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.