Effect of mutation at the interface of Trp-repressor dimeric protein: a steered molecular dynamics simulation

Serine 477 plays a crucial role in the interaction of the SARS-CoV-2 spike protein with the human receptor ACE2

Since the worldwide outbreak of the infectious disease COVID-19, several studies have been published to understand the structural mechanism of the novel coronavirus SARS-CoV-2. During the infection process, the SARS-CoV-2 spike (S) protein plays a crucial role in the receptor recognition and cell membrane fusion process by interacting with the human angiotensin-converting enzyme 2 (hACE2) receptor. However, new variants of these spike proteins emerge as the virus passes through the disease reservoir. This poses a major challenge for designing a potent antigen for an effective immune response against the spike protein. Through a normal mode analysis (NMA) we identified the highly flexible region in the receptor binding domain (RBD) of SARS-CoV-2, starting from residue 475 up to residue 485. Structurally, the position S477 shows the highest flexibility among them. At the same time, S477 is hitherto the most frequently exchanged amino acid residue in the RBDs of SARS-CoV-2 mutants. Therefore, using MD simulations, we have investigated the role of S477 and its two frequent mutations (S477G and S477N) at the RBD during the binding to hACE2. We found that the amino acid exchanges S477G and S477N strengthen the binding of the SARS-COV-2 spike with the hACE2 receptor.

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.

Steered Molecular Dynamics Simulations Predict Conformational Stability of Glutamate Receptors

The stability of protein-protein interfaces can be essential for protein function. For ionotropic glutamate receptors, a family of ligand-gated ion channels vital for normal function of the central nervous system, such an interface exists between the extracellular ligand binding domains (LBDs). In the full-length protein, the LBDs are arranged as a dimer of dimers. Agonist binding to the LBDs opens the ion channel, and briefly after activation the receptor desensitizes. Several residues at the LBD dimer interface are known to modulate desensitization, and conformational changes around these residues are believed to be involved in the state transition. The general hypothesis is that the interface is disrupted upon desensitization, and structural evidence suggests that the disruption might be substantial. However, when cross-linking the central part of this interface, functional data suggest that the receptor can still undergo desensitization, contradicting the hypothesis of major interface disruption. Here, we illustrate how opening the dimer interface using steered molecular dynamics (SMD) simulations, and analyzing the work values required, provides a quantitative measure for interface stability. For one subtype of glutamate receptors, which is regulated by ion binding to the dimer interface, we show that opening the interface without ions bound requires less work than with ions present, suggesting that ion binding indeed stabilizes the interface. Likewise, for interface mutants with longer-lived active states, the interface is more stable, while the work required to open the interface is reduced for less active mutants. Moreover, a cross-linked mutant can still undergo initial interface opening motions similar to the native receptor and at similar energetic cost. Thus, our results support that interface opening is involved in desensitization. Furthermore, they provide reconciliation of apparently opposing data and demonstrate that SMD simulations can give relevant biological insight into longer time scale processes without the need for expensive calculations.

The molecular mechanism of ligand unbinding from the human telomeric G-quadruplex by steered molecular dynamics and umbrella sampling simulations

G-quadruplexes are attractive drug targets in cancer therapy. Understanding the mechanisms of the binding-unbinding processes involving biomolecules and molecular recognition is essential for designing new drugs of G-quadruplexes. We performed steered molecular dynamics and umbrella sampling simulations to investigate the molecular mechanism and kinetics of ligand unbinding processes of the basket, propeller and hybrid G-quadruplex structures. Our studies of the ligand charge effect showed that Coulomb interaction plays a significant role in stabilizing the G-quadruplex structure in the unbinding process. The free energy profiles were carried out and the free energy changes associated with the unbinding process were computed quantitatively, whereas these results could help to identify accessible binding sites and transient interactions. The dynamics of the hydration shell water molecules around the G-quadruplex exhibits an abnormal Brownian motion, and the thickness and free energy of the hydration shell were estimated. A two-step relaxation scheme was theoretically developed to describe the kinetic reaction of BMVC and G-quadruplex interactions. Our computed results fall in a reasonable range of experimental data. The present investigation could be helpful in the structure-based drug design.

Exploration of the chlorpyrifos escape pathway from acylpeptide hydrolases using steered molecular dynamics simulations

Acylpeptide hydrolases (APH) catalyze the removal of an N-acylated amino acid from blocked peptides. APH is significantly more sensitive than acetylcholinesterase, a target of Alzheimer's disease, to inhibition by organophosphorus (OP) compounds. Thus, OP compounds can be used as a tool to probe the physiological functions of APH. Here, we report the results of a computational study of molecular dynamics simulations of APH bound to the OP compounds and an exploration of the chlorpyrifos escape pathway using steered molecular dynamics (SMD) simulations. In addition, we apply SMD simulations to identify potential escape routes of chlorpyrifos from hydrolase hydrophobic cavities in the APH-inhibitor complex. Two previously proposed APH pathways were reliably identified by CAVER 3.0, with the estimated relative importance of P1 > P2 for its size. We identify the major pathway, P2, using SMD simulations, and Arg526, Glu88, Gly86, and Asn65 are identified as important residues for the ligand leaving via P2. These results may help in the design of APH-targeting drugs with improved efficacy, as well as in understanding APH selectivity of the inhibitor binding in the prolyl oligopeptidase family.

Insight into the interactive residues between two domains of human somatic Angiotensin-converting enzyme and Angiotensin II by MM-PBSA calculation and steered molecular dynamics simulation

Angiotensin-converting enzyme (ACE), a membrane-bound zinc metallopeptidase, catalyzes the formation of Angiotensin-II (AngII) and the deactivation of bradykinin in the renin-angiotensin-aldosterone and kallikrein-kinin systems. As a hydrolysis product of ACE, AngII is regarded as an inhibitor and displays stronger competitive inhibition in the C-domain than the N-domain of ACE. However, the AngII binding differences between the two domains and the mechanisms behind AngII dissociation from the C-domain are rarely explored. In this work, molecular docking, Molecular Mechanics/Poisson-Boltzmann Surface Area calculation, and steered molecular dynamics (SMD) are applied to explore the structures and interactions in the binding or unbinding of AngII with the two domains of human somatic ACE. Calculated free energy values suggest that the C-domain-AngII complex is more stable than the N-domain-AngII complex, consistent with available experimental data. SMD simulation results imply that electrostatic interaction is dominant in the dissociation of AngII from the C-domain. Moreover, Gln106, Asp121, Glu123, and Tyr213 may be the key residues in the unbinding pathway of AngII. The simulation results in our work provide insights into the interactions between the two domains of ACE and its natural peptide inhibitor AngII at a molecular level. Moreover, the results provide theoretical clues for the design of new inhibitors.

Unbinding pathways of VEGFR2 inhibitors revealed by steered molecular dynamics

A detailed atomistic description of the unbinding process of sorafenib and sunitinib, two known VEGFR2 inhibitors clinically used to treat renal cell carcinoma, was unraveled by using steered molecular dynamics (SMD) simulations. While sunitinib is a fast-dissociating binder, sorafenib exhibits quite a long residence time at this enzyme, which might impact its duration of action in vivo. In order to gain insights into the kinetically different behaviors of the two inhibitors, an SMD study was carried out, which involved a careful optimization of the force and velocity parameters. We were able to identify two different binding pathways for the two inhibitors, as sunitinib exited the ATP binding site from the cavity entrance without a rupture point while sorafenib moved opposite to the ATP binding site entrance. Furthermore, the calculated ΔGoff values clearly reflect on a qualitative level the distinct off-rates of the two inhibitors, thus suggesting that this protocol could be tried on other VEGFR2 ligands to assess its robustness and then used to rank structural analogues of these derivatives.

Analyzing machupo virus-receptor binding by molecular dynamics simulations

In many biological applications, we would like to be able to computationally predict mutational effects on affinity in protein-protein interactions. However, many commonly used methods to predict these effects perform poorly in important test cases. In particular, the effects of multiple mutations, non alanine substitutions, and flexible loops are difficult to predict with available tools and protocols. We present here an existing method applied in a novel way to a new test case; we interrogate affinity differences resulting from mutations in a host-virus protein-protein interface. We use steered molecular dynamics (SMD) to computationally pull the machupo virus (MACV) spike glycoprotein (GP1) away from the human transferrin receptor (hTfR1). We then approximate affinity using the maximum applied force of separation and the area under the force-versus-distance curve. We find, even without the rigor and planning required for free energy calculations, that these quantities can provide novel biophysical insight into the GP1/hTfR1 interaction. First, with no prior knowledge of the system we can differentiate among wild type and mutant complexes. Moreover, we show that this simple SMD scheme correlates well with relative free energy differences computed via free energy perturbation. Second, although the static co-crystal structure shows two large hydrogen-bonding networks in the GP1/hTfR1 interface, our simulations indicate that one of them may not be important for tight binding. Third, one viral site known to be critical for infection may mark an important evolutionary suppressor site for infection-resistant hTfR1 mutants. Finally, our approach provides a framework to compare the effects of multiple mutations, individually and jointly, on protein-protein interactions.