Steered Molecular Dynamics Research Articles

Effects of Substitution Ratios of Zinc-Substituted Hydroxyapatite on Adsorption and Desorption Behaviors of Bone Morphogenetic Protein-2.

Understanding interactions between bone morphogenetic proteins (BMPs) and biomaterials is of great significance in preserving the structure and bioactivity of BMPs when utilized in clinical applications. Currently, bone morphogenetic protein-2 (BMP-2) is one of the most important growth factors in bone tissue engineering; however, atomistic interactions between BMP-2 and zinc-substituted hydroxyapatite (Zn-HAP, commonly used in artificial bone implants) have not been well clarified until now. Thus, in this work, the interaction energies, binding/debinding states, and molecular structures of BMP-2 upon a series of Zn-HAP surfaces (Zn-HAPs, 1 at%, 2.5 at%, 5 at%, and 10 at% substitution) were investigated by hybrid molecular dynamics (MD) and steered molecular dynamics (SMD) simulations. Meanwhile, cellular studies including alkaline phosphatase (ALP) activity and reverse transcription-polymerase chain reaction (RT-PCR) assay were performed to verify the theoretical modeling findings. It was found that, compared to pure HAP, Zn-HAPs exhibited a higher binding affinity of BMP-2 at the adsorption process; meanwhile, the detachment of BMP-2 upon Zn-HAPs was more difficult at the desorption process. In addition, molecular structures of BMP-2 could be well stabilized upon Zn-HAPs, especially for Zn10-HAP (with a 10 at% substitution), which showed both the higher stability of cystine-knots and less change in the secondary structures of BMP-2 than those upon HAP. Cellular studies confirmed that higher ALP activity and osteogenic marker gene expression were achieved upon BMP-2/Zn-HAPs than those upon BMP-2/HAP. These findings verified that Zn-HAPs favor the adsorption of BMP-2 and leverage the bioactivity of BMP-2. Together, this work clarified the interaction mechanisms between BMP-2 and Zn-HAPs at the atom level, which could provide new molecular-level insights into the design of BMP-2-loaded biomaterials for bone tissue engineering.

Calculation of binding affinity of JAK1 inhibitors via accurately computational estimation

Janus kinase 1 (JAK1) is a tyrosine kinase that is involved in the initiation of responses to a number of different cytokine receptor families. The JAK1-dependent pathway is a therapeutic target, and several JAK inhibitors have been developed thanks to intensive research. However, since the ATP binding sites of JAK family members are quite alike, JAK1 inhibitors can thus be less selective, resulting in unanticipated adverse effects. Despite this, minor variations in the ATP-binding site have been extensively used to find a variety of small compounds with different inhibitory properties. Stronger binding affinity of JAK1 inhibitors is believed to be able to reduce the negative effects, leading to better treatment results. Therefore, a thorough computational search that can effectively identify ligands with extremely high binding affinity for JAK1 to serve as promising inhibitors is required. Here, a method combining steered-molecular dynamic (SMD) simulations with a modified linear interaction energy (LIE) model has been developed to evaluate the binding affinities of known JAK1 inhibitors. The correlation coefficient between the estimated and experimental values was 0.72 and a root-mean-square error was 0.97 kcal•mol−1, revealing that the SMD/LIE method can precisely and quickly predict the binding free energies of JAK1 inhibitors. Furthermore, three marine fungus-derived compounds, namely hansforesters E, hansforesters G and tetroazolemycins B, were identified to be particularly promising JAK1 inhibitors, accordingly. These findings show that the SMD/LIE method has a lot of promise for in silico screening of possible JAK1 inhibitors from a vast number of compounds that are now accessible. Communicated by Ramaswamy H. Sarma

Comparison of the unbinding process of RBD-ACE2 complex between SARS-CoV-2 variants (Delta, delta plus, and Lambda): A steered molecular dynamics simulation

ABSTRACT The Delta and Lambda SARS-CoV-2 variants are characterised with the double mutations T478K/L452R, and F490S/L452Q, respectively, in their receptor binding domains (RBDs) of the spike proteins. The Delta Plus variant is one of the Delta variant types with K417N substitution in RBD. Basically, SARS-CoV-2 variants need to bind to Angiotensin-Converting Enzyme 2(ACE2) to enter the host cell. Here, the effects of these mutations were studied on the binding affinity of RBD-ACE2 using steered molecular dynamics (SMD) simulation. According to the results, the rupture force of RBD-ACE2 for the wild type and Delta variant were close to each other, and the minimum and maximum rupture forces occurred in the Delta Plus and Lambda variants, respectively. While T478K and L452R occur adjacent to the RBD binding site and have no major effect on the rupture force of the RBD-ACE2 complex, K417N occurs exactly at the binding site of RBD-ACE2 and causes a decrease in the rupture force for the Delta Plus variant. The results of this study showed that T478K, L452R, and K417N could not justify the infection and transmissibility of the Delta and Delta Plus variants. However, L452Q and F490S were able to affect the binding affinity of the RBD-ACE2 complex.

Effect of mixed surfactants on foam stabilization: A molecular dynamics simulation

Mixed surfactants have antagonistic or synergistic effects on the air/liquid interface of foams, representing two opposite interactions about destroying or stabilizing foams. In this paper, mixed surfactant systems containing non-ionic surfactant polyoxyethylene alkyl ether (AEO) and anionic surfactant sodium dodecyl sulfate (SDS) were discussed by molecular dynamics simulation to study the mechanism of decreasing stability of SDS foam after AEO was added. The results showed that with the addition of AEO, the interaction between SDS molecules is weakened, and the thickness of hydrated shell around the polar head of SDS decreases. In order to simulate the process of foam drainage and rupture, steered molecular dynamics (SMD) containing two external forces on the foam film was performed. One is about the z-axis vertical external force, representing the stability of SDS hydrated layer at the interface; the other is about the external force in horizontal direction, representing the rupture process of the liquid film. The simulated results showed that the higher the concentration of AEO in the mixed surfactant interface layer, the smaller the external force to destroy the foam film interface. In addition, the free energy change of AEO molecules from the water phase through the SDS layer to the air phase further showed that the addition of AEO can reduce the stability of SDS foam films and play a defoaming role in the studied system. The simulated conclusions were consistent with the experimental results about the decreasing stability of SDS foam after AEO was added.

Mechanical Stability of Ribonuclease A Heavily Depends on the Redox Environment.

Disulfide bonds are covalent bonds that connect nonlocal fragments of proteins, and they are unique post-translational modifications of proteins. They require the oxidizing environment to be stable, which occurs for example during oxidative stress; however, in a cell the reductive environment is maintained, lowering their stability. Despite many years of research on disulfide bonds, their role in the protein life cycle is not fully understood and seems to strictly depend on a system or process in which they are involved. In this article, coarse-grained UNited RESidue (UNRES), and all-atom Assisted Model Building with Energy Refinement (AMBER) force fields were applied to run a series of steered molecular dynamics (SMD) simulations of one of the most studied, but still not fully understood, proteins—ribonuclease A (RNase A). SMD simulations were performed to study the mechanical stability of RNase A in different oxidative–reductive environments. As disulfide bonds (and any other covalent bonds) cannot break/form in any classical all-atom force field, we applied additional restraints between sulfur atoms of reduced cysteines which were able to mimic the breaking of the disulfide bonds. On the other hand, the coarse-grained UNRES force field enables us to study the breaking/formation of the disulfide bonds and control the reducing/oxidizing environment owing to the presence of the designed distance/orientation-dependent potential. This study reveals that disulfide bonds have a strong influence on the mechanical stability of RNase A only in a highly oxidative environment. However, the local stability of the secondary structure seems to play a major factor in the overall stability of the protein. Both our thermal unfolding and mechanical stretching studies show that the most stable disulfide bond is Cys65–Cys72. The breaking of disulfide bonds Cys26–Cys84 and Cys58–Cys110 is associated with large force peaks. They are structural bridges, which are mostly responsible for stabilizing the RNase A conformation, while the presence of the remaining two bonds (Cys65–Cys72 and Cys40–Cys95) is most likely connected with the enzymatic activity rather than the structural stability of RNase A in the cytoplasm. Our results prove that disulfide bonds are indeed stabilizing fragments of the proteins, but their role is strongly redox environment-dependent.

Protein aggregation rate depends on mechanical stability of fibrillar structure.

The formation of the fibrillar structure of amyloid proteins/peptides is believed to be associated with neurodegenerative diseases, such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis. Since the rate of aggregation can influence neurotoxicity, finding the key factors that control this rate is of paramount importance. It was recently found that the rate of protein aggregation is related to the mechanical stability of the fibrillar structure such that the higher the mechanical stability, the faster the fibril is formed. However, this conclusion was supported by a limited dataset. In this work, we expand the previous study to a larger dataset, including the wild type of Aβ42 peptide and its 20 mutants, the aggregation rate of which was measured experimentally. By using all-atom steered molecular dynamics (SMD) simulations, we can assess the mechanical stability of the fibril structure, which is characterized by the rupture force, pulling work, and unbinding free energy barrier. Our result confirms that mechanical stability is indeed related to the aggregation rate. Since the estimation of the aggregation rate using all-atom simulations is almost forbidden by the current computational capabilities, our result is useful for predicting it based on information obtained from fast SMD simulations for fibrils.

Implementation of Telescoping Boxes in Adaptive Steered Molecular Dynamics.

Long-time dynamical processes, such as those involving protein unfolding and ligand interactions, can be accelerated and realized through steered molecular dynamics (SMD). The challenge has been the extraction of information from such simulations that generalize for complex nonequilibrium processes. The use of Jarzynski's equality opened the possibility of determining the free energy along the steered coordinate, but sampling over the nonequilibrium trajectories is slow to converge. Adaptive steered molecular dynamics (ASMD) and other related techniques have been introduced to overcome this challenge through the use of stages. Here, we take advantage of these stages to address the numerical cost that arises from the required use of very large solvent boxes. We introduce telescoping box schemes within adaptive steered molecular dynamics (ASMD) in which we adjust the solvent box between stages and thereby vary (and optimize) the required number of solvent molecules. We have benchmarked the method on a relatively long α-helical peptide, Ala30, with respect to the potential of mean force and hydrogen bonds. We show that the use of telescoping boxes introduces little numerical error while significantly reducing the computational cost.

Critical Extracellular Ca2+ Dependence of the Binding between PTH1R and a G-Protein Peptide Revealed by MD Simulations.

The parathyroid hormone type 1 receptor (PTH1R), a canonical class B GPCR, is regulated by a positive allosteric modulator, extracellular Ca2+. Calcium ions prolong the residence time of PTH on the PTH1R, leading to increased receptor activation and duration of cAMP signaling. But the essential mechanism of the allosteric behavior of PTH1R is not fully understood. Here, extensive molecular dynamics (MD) simulations are performed for the PTH1R-G-protein combinations with and without Ca2+ to describe how calcium ions allosterically engage receptor-G-protein coupling. We find that the binding of Ca2+ stabilizes the conformation of the PTH1R-PTH-spep (the α5 helix of Gs protein) complex, especially the extracellular loop 1 (ECL1). Moreover, the MM-GBSA result indicates that Ca2+ allosterically promotes the interaction between PTH1R and spep, consistent with the observation of steered molecular dynamics (SMD) simulations. We further illuminate the possible allosteric signaling pathway from the stable Ca2+-coupling site to the intracellular G-protein binding site. These results unveil structural determinants for Ca2+ allosterism in the PTH1R-PTH-spep complex and give insights into pluridimensional GPCR signaling regulated by calcium ions.

Prediction of COMT Inhibitors Using Machine Learning and Molecular Dynamics Methods.

Catechol O-methyltransferase (COMT) plays a vital role in deactivating neurotransmitters like dopamine, norepinephrine, etc., by methylating those compounds. However, the deactivation of an excess amount of neurotransmitters leads to serious mental ailments such as Parkinson's disease. Molecules that bind inside the enzyme's active site inhibit this methylation mechanism by methylating themselves, termed COMT inhibitors. Our study is focused on designing these inhibitors by various machine learning methods. First, we have developed a classification model with experimentally available COMT inhibitors, which helped us generate a new data set of small inhibitor-like molecules. Then, to predict the activity of the new molecules, we have applied regression techniques such as Random Forest, AdaBoost, gradient boosting, and support vector machines. Each of the regression models yielded an R2 value > 70% for both training and test data sets. Finally, to validate our models, 200 ns long molecular dynamics (MD) simulations of the two known inhibitors with known IC50 values and the resultant inhibitors were performed inside the binding pockets to check their stability within. The free energy barrier of the methyl transfer from S-adenosyl-l-methionine (SAM) to each inhibitor was determined by combining steered molecular dynamics (SMD) and umbrella sampling using the quantum mechanics/molecular mechanics (QM/MM) method.

Determination of Multidirectional Pathways for Ligand Release from the Receptor: A New Approach Based on Differential Evolution.

Steered molecular dynamics (SMD) simulation is a powerful method in computer-aided drug design as it can be used to access the relative binding affinity with high precision but with low computational cost. The success of SMD depends on the choice of the direction along which the ligand is pulled from the receptor-binding site. In most simulations, the unidirectional pathway was used, but in some cases, this choice resulted in the ligand colliding with the complex surface of the exit tunnel. To overcome this difficulty, several variants of SMD with multidirectional pulling have been proposed, but they are not completely devoid of disadvantages. Here, we have proposed to determine the direction of pulling with a simple scoring function that minimizes the receptor–ligand interaction, and an optimization algorithm called differential evolution is used for energy minimization. The effectiveness of our protocol was demonstrated by finding expulsion pathways of Huperzine A and camphor from the binding site of Torpedo California acetylcholinesterase and P450cam proteins, respectively, and comparing them with the previous results obtained using memetic sampling and random acceleration molecular dynamics. In addition, by applying this protocol to a set of ligands bound with LSD1 (lysine specific demethylase 1), we obtained a much higher correlation between the work of pulling force and experimental data on the inhibition constant IC50 compared to that obtained using the unidirectional approach based on minimal steric hindrance.

Analysis of the Passage Times for Unfolding/Folding of the Adenine Riboswitch Aptamer.

The conformational transitions of the adenosine deaminase A-riboswitch aptamer both with and without ligand binding are investigated within the tenets of the generalized Langevin equation in a complex viscoelastic cellular environment. Steered molecular dynamics (SMD) simulations are performed to evaluate and compare the results of the first passage times (FPTs) with those obtained from the theory for the unfold and fold transitions of the aptamer. The results of the distribution of Kramers's FPT reveal that the unfold-fold transitions are faster and hence more probable as compared to the fold-unfold transitions of the riboswitch aptamer for both ligand-bound and -unbound states. The transition path time is lower than Kramers's FPT for the riboswitch aptamer as the transition path times for the unfold-fold transition of both without and with ligand binding are insensitive to the details of the exact mechanism of the transition events. However, Kramers's FPTs show varied distributions which correspond to different transition pathways, unlike the transition path times. The mean FPT increases with an increase in the complexity of the cellular environment. The results of Kramers's FPT, transition path time distribution, and mean FPT obtained from our calculations qualitatively match with those obtained from the SMD simulations. Analytically derived values of the mean transition path time show good quantitative agreement with those estimated from the single-molecule force spectroscopy experiments for higher barrier heights.

Molecular properties of Ca2+ transport through TRPV2 channel: a molecular dynamics simulations study

TRPV channels are a category of nonselective cation channels that are activated by heat and ligands and permeate monovalent and divalent ions. The mechanism of Ca2+ transfer through TRPV2 channel is not well known. Here, we investigated the reaction coordination and energy fluctuation of Ca2+ transition in TRPV2 channel by steered molecular dynamics (SMD) simulations and potential of mean force (PMF) calculation. Results showed that electrostatic interactions between Ca2+ and residues of the first and second gates had main roles in ions transfer through the channel. Also, we recognized important amino acids in this path. Moreover, results indicated that enter and exit of calcium ions need to overcome barrier energies in the first and second gates. Communicated by Ramaswamy H. Sarma

Force-enhanced biophysical connectivity of platelet \u03b23 integrin signaling through Talin is predicted by steered molecular dynamics simulations

Platelet β3-integrin signaling through Talin is crucial in platelet transmembrane signaling, activation, adhesion, spreading and aggregation, and remains unclear in mechano-microenvironments. In order to examine Talin-β3 integrin biophysical connectivity, a series of “ramp-clamp” steered molecular dynamics (SMD) simulations were performed on complex of F3 domain of Talin and cytoplasmic tail of β3 integrin to imitate different force-loads in platelet. Pull-induced allostery of the hydrophobic pocket in F3 domain might markedly enhance complex rupture-force (> 150pN) and slow down breakage of the complex; the complex should mechano-stable for its conformational conservation under loads (≤ 80pN); increasing force below 60pN would decrease the complex dissociation probability, and force-induced extension of β5 strand on Talin and binding site residues, ASP740 and ALA742 as well as Asn744, on β3-integrin were responsible for the force-enhanced linkage of the Talin-β3 integrin. Force might enhance biophysical connectivity of β3-integrin signaling through Talin by a catch bond mechanism, which be mediated by the force-induced allostery of complex at clamped stage. This work provides a novel insight into the force-regulated transmembrane β3-integrin signaling and its molecular basis for platelet activation, and exhibited a potential power of the present computer strategy in predicting mechanical regulation on ligand-receptor interaction under loads.

Nanoscale Topographical Effects on the Adsorption Behavior of Bone Morphogenetic Protein-2 on Graphite.

The interaction between bone morphogenetic protein-2 (BMP-2) and the surface of biomaterials is essential for the restoration of bone and cartilage tissue, inducing cellular differentiation and proliferation. The properties of the surface, including topology features, regulate the conformation and bioactivity of the protein. In this research, we investigated the influence of nanopatterned surfaces on the interaction of a homodimer BMP-2 with graphite material by combining molecular dynamics (MD) and steered molecular dynamics (SMD) simulations. The graphite substrates were patterned as flat, linear grating, square, and circular profiles in combination with BMP-2 conformation in the side-on configuration. Ramachandran plots for the wrist and knuckle epitopes indicated no steric hindrances and provided binding sites to type I and type II receptors. Results showed two optimal patterns that increased protein adsorption of the lower monomer while preserving the secondary structure and leaving the upper monomer free to interact with the cells. Charged residues arginine and lysine and polar residues histidine and tyrosine were the main residues responsible for the strong interaction with the graphite surface. This research provides new molecular-level insights to further understand the mechanisms underlying protein adsorption on nanoscale patterned substrates.

Probing the Molecular Basis of Cofactor Affinity and Conformational Dynamics of Mycobacterium tuberculosis Elongation Factor Tu: An Integrated Approach Employing Steered Molecular Dynamics and Umbrella Sampling Simulations.

The emergence of multidrug-resistant and extensively drug-resistant tuberculosis strains is the reason that the infectious tuberculosis pathogen is still the most common cause of death. The quest for new antitubercular drugs that can fit into multidrug regimens, function swiftly, and overcome the ever-increasing prevalence of drug resistance continues. The crucial role of MtbEF-Tu in translation and trans-translation processes makes it an excellent target for antitubercular drug design. In this study, the primary sequence of MtbEF-Tu was used to model the three-dimensional structures of MtbEF-Tu in the presence of GDP ("off" state) and GTP ("on" state). The binding free energy computed using both the molecular mechanics/Poisson-Boltzmann surface area and umbrella sampling approaches shows that GDP binds to MtbEF-Tu with an ∼2-fold affinity compared to GTP. The steered molecular dynamics (SMD) and umbrella sampling simulation also shows that the dissociation of GDP from MtbEF-Tu in the presence of Mg2+ is a thermodynamically intensive process, while in the absence of Mg2+, the destabilized GDP dissociates very easily from the MtbEF-Tu. Naturally, the dissociation of Mg2+ from the MtbEF-Tu is facilitated by the nucleotide exchange factor EF-Ts, and this prior release of magnesium makes the dissociation process of destabilized GDP easy, similar to that observed in the umbrella sampling and SMD study. The MD simulations of MtbEF-Tu's "on" state conformation in the presence of GTP reveal that the secondary structure of switch-I and Mg2+ coordination network remains similar to its template despite the absence of identity in the conserved region of switch-I. On the other hand, the secondary structure in the conserved region of the switch-I of MtbEF-Tu unwinds from a helix to a loop in the presence of GDP. The major conformational changes observed in switch-I and the movement of Thr64 away from Mg2+ mainly reflect essential conformational changes to make the shift of MtbEF-Tu's "on" state to the "off" state in the presence of GDP. These obtained structural and functional insights into MtbEF-Tu are pivotal for a better understanding of structural-functional linkages of MtbEF-Tu, and these findings may serve as a basis for the design and development of MtbEF-Tu-specific inhibitors.

In silico screening of potential β-secretase (BACE1) inhibitors from VIETHERB database.

For the purpose of discovering potential inhibitors of β-amyloid (BACE1), which is a crucial element in Alzheimer's disease (AD) pathogenesis, an in silico study of naturally occurring compounds was performed using precise computational approaches. Autodock4 package was preliminary used to predict the binding affinities to BACE1 of more than four thousand compounds presented in the Vietnamese plants (VIETHERB) database. Based on docking results, twenty top-lead compounds having the largest docking energy to BACE1 were rigorously examined using steered molecular dynamics (SMD) simulations. Interestingly, SMD results found that the binding affinity values of three compounds, including myricetin 3-O-(3''-galloylrhamnopyranoside), quercetin 3-O-neohesperidoside, and hydroxysafflor yellow A, are remarkably higher than that of the well-known BACE1 inhibitor, 23I, and these compounds can thus be considered the promising candidates for inhibitors of BACE1.

Identification of novel bioactive compounds from Olea europaea by evaluation of chemical compounds in the OliveNet™ library: in silico bioactivity and molecular modelling, and in vitro validation of hERG activity

BackgroundAs highlighted in the OliveNet™ library, Olea europaea consists of a diverse collection of chemical compounds. We have classified over 600 compounds into 13 main classes and 47 subclasses. Various compounds, including oleuropein and hydroxytyrosol, have been investigated for their potential beneficial effects in multiple human pathologies. However, the vast majority of compounds remain largely unexplored and approximately 50% are currently non-commercially available. MethodHere, we utilise conventional software to characterise the absorption, distribution, metabolism, excretion, and toxicity profile of OliveNet™ compounds. Molecular docking was performed for assessment of P-glycoprotein (P-gp) inhibition and interactions with the human ether-à-go-go-related gene (hERG) channel. Potential hERG ion channel inhibition was calibrated using in vitro patch clamp assays and steered molecular dynamics (SMD) simulations were used to examine membrane permeability of a subset of compounds. ResultsOur findings indicated that 313 out of 675 olive compounds were predicted to be absorbed by the gastrointestinal tract. Hydroxytyrosol required the least amount of force to pass through the lipid bilayer compared to elenolic acid diglucoside in SMD simulations. Based on the ADMET and molecular docking data, the hERG inhibitory activities of verbascoside, oleuropein, and hydroxytyrosol were investigated using patch clamp assays and they were identified as non-inhibitors. ConclusionsWhile the favourable properties of well-known compounds were confirmed, we identified oleuropein aglycone decarboxymethyl dialdehyde acetal form, decarboxymethyl elenolic acid dialdehyde, acetal of decarboxymethyl elenolic acid dialdehyde, methyl malate-β-hydroxytyrosol ester, hydroxytyrosil elenolate, D-(+)-erythro-1-(4-hydroxy-3-methoxy)-214-phenyl-1,2,3-propantriol, (+)-1-acetoxypinoresinol-4″-O-methyl ether, and 3-[1-(hydroxymethyl)-(E)-1-propenyl] glutaric acid as potential candidates for synthesis and further evaluation.

The effects of mutation on the drug binding affinity of Neuraminidase: case study of Capsaicin using steered molecular dynamics simulation.

The influenza virus is an important respiratory pathogen that causes many incidences of diseases and even death each year. One of the primary factors of this virus is the Neuraminidase surface protein, which causes the virus to leave the host cell and spread to new target cells. The main antiviral medication for influenza is designed as a protein inhibitor ligand that prevents further spread of the disease, and eventually relieves the emerged symptoms. The effectiveness of such inhibitory drugs is highly associated with their binding affinity. In this paper, the binding affinity of an herbal ligand of Capsaicin bound to Neuraminidase of the influenza virus is investigated using steered molecular dynamics (SMD) simulation. Since mutations of the virus directly impact the binding affinity of the inhibitory drugs, different mutations were generated by using Mutagenesis module. The rapid spread of infection during the avian influenza A/H5N1 epidemic has raised concerns about far more dangerous consequences if the virus becomes resistant to current drugs. Currently, oseltamivir (Tamiflu), zanamivir (Relenza), pramivir (Rapivab), and laninamivir (Inavir) are increasingly used to treat the flu. However, with the rapid evolution of the virus, some drug-resistant strains are emerging. Therefore, it is very important to seek alternative therapies and identify the roots of drug resistance. Obtained results demonstrated a reduced binding affinity for the applied mutations. This reduction in binding affinity will cause the virus mutation to become resistant to the drug, which will spread the disease and make it more difficult to treat. From a molecular prospect, this decrease in binding affinity is due to the loss of a number of effective bonds between the ligand and the receptor, which occurs with mutations of the wild-type (WT) species. The results of the present study can be used in the rational design of novel drugs that are compatible with specific mutations.

A computational algorithm to assess the physiochemical determinants of T cell receptor dissociation kinetics

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.

A computational study on strontium ion modified hydroxyapatite-fibronectin interactions.

Protein adsorption is the first key step in cell-material interactions. The initial phase of such an adsorption process can only be probed using modelling approaches like molecular dynamics (MD) simulations. Despite a large number of studies on the adsorption behaviour of proteins on different biomaterials including calcium phosphates (CaP), little attention has been paid towards the quantitative assessment of the effects of various physicochemical influencers like surface modification, pH, and ionic strength. In the case of doped CaPs, surface modification through isomorphic substitution of foreign ions inside the apatite structure is of particular interest in the context of protein-HA interactions, as it is widely used to tailor the biological response of HA. Given this background, we present here the molecular-level understanding of the fibronectin (FN) adsorption mechanism and kinetics on a Sr2+-doped hydroxyapatite, HA, (001) surface at 300 K by means of all-atom molecular dynamics simulations. Electrostatic interactions involved in the adsorption of FN on HA were found to be significantly modified due to Sr2+ doping into the apatite lattice. In harmony with the published experimental observations, the Sr-doped surfaces were found to better support FN adhesion compared to pure HA, with 10 mol% Sr-doped HA exhibiting the best FN adsorption. The observed altered adsorption behaviour of FN on Sr-doped HA was correlated with the Hofmeister effect. Moreover, the non-monotonous trend of the FN-material interaction energy can be attributed to the spatial rearrangement of the functional groups (PO43-, OH-) in the apatite crystal. Sr2+ ions also influence the stability of the secondary structure of FN, as observed from the root mean square deviation (RMSD) and root mean square fluctuation (RMSF) analysis. The presence of Sr2+ enhances the flexibility of specific residues (residue nos. 20-44, 74-88) of the FN module. Rupture forces to disentangle FN from the biomaterial surface, obtained from steered molecular dynamics (SMD) simulations, were found to corroborate well with the results of equilibrium MD simulations. One particular observation is that the availability of an RGD motif (Arginine-Glycine-aspartate sequence, which interacts with cell surface receptor integrin to form a focal adhesion complex) for the interaction with cell surface receptor integrin is not significantly influenced by Sr2+ substitution.