Molecular dynamics study of the connection between flap closing and binding of fullerene-based inhibitors of the HIV-1 protease

The complementary spatial relationship between fullerene C(60) and the hydrophobic cavity region of the human immunodeficiency virus (HIV) protease, which houses the active site of the enzyme, has led to the suggestion that fullerene-based derivatives could have potential use as effective HIV protease inhibitors. The ability of such compounds to desolvate the cavity region leads to a strong hydrophobic interaction between the C(60) moiety and residues in the cavity region. In this study, the connection between the motion of the so-called flexible flaps of the cavity and favorable binding of a fullerene-based protease inhibitor is explored using multiple-time scale molecular dynamics simulations and free energy techniques. In addition, the effect of the interaction between the C(60) moiety and the residues in the cavity region on the water content of the cavity is also investigated. Conformational free energy profiles along a suitably chosen flap opening coordinate show a considerable barrier to flap opening in the presence of the inhibitor, while no such barrier exists for the protease alone. This result is interpreted in terms of a strong hydrophobic interaction between the C(60) moiety and the flexible flaps, which cause the latter to close tightly around the inhibitor, thereby expelling water from the cavity and leading to a favorable binding interaction. This interpretation is rationalized by direct analysis of the water content in the cavity in the presence and absence of the inhibitor.     

A dual target of Plasmepsin IX and X: Unveiling the atomistic superiority of a core chemical scaffold in malaria therapy

Plasmepsin IX and X, members of the prominent aspartic family of proteases whose function were hitherto unknown have only recently been established as key mediators of erythrocyte invasion and egress of the virulent malarial parasite. Inhibitor 49c, a potent antimalarial peptidomimetic inhibitor initially developed to target Plasmepsin II has lately been proven to exhibit potent inhibitory activity against Plasmepsin IX and X. However, the molecular and structural dynamics supporting its inhibitory activity remain inconclusive. Hindering the motion of the flap and hinge region of an aspartic protease remains essential for disabling the catalytic activity of the enzyme. Integrating molecular dynamic simulations coupled with other advanced biocomputational tools, we reveal the enhanced structural mechanistic competence of 49c in complex with Plasmepsin IX and X relative to Pepstatin. Pepstatin, a known aspartic protease inhibitor which actively hinders the opening and closing of the flap tip and flexible loop and consequently limits access to the catalytic aspartic residues, however, its administration has been related to elevated levels of toxicity. Thermodynamic calculations reveal a higher relative binding free energy associated with Plasmepsin IX and X in complex with 49c as opposed to Pepstatin. A relatively compact and structurally rigid 49c bound complexes sequel into the restriction of the flap and hinge residues by restraining cohesive movement, consequently hindering their “twisting motion” from transpiring. Findings unveil an atomistic perspective into the structural superiority of 49c in complex with Plasmepsin IX and X.     

Molecular dynamics of the long neurotoxin LSIII

Long neurotoxins bind tightly and specifically to the nicotinic acetylcholine receptor (AChR) in postsynaptic membranes and are useful for exploring the biology of synapses. In crystallographic studies of long neurotoxins the principal binding loop appears disordered, but the NMR solution structure of the long neurotoxin LSIII revealed significant local order, even though the loop is disordered with respect to the globular core. A possible mechanism for conferring global disorder while preserving local order is rigid-body motion of the loop about a hinge region. Here we report investigations of LSIII dynamics based on (13)C(alpha) magnetic relaxation rates and molecular dynamics simulation. The relaxation rates and MD simulation both confirm the hypothesis of rigid-body motion of the loop and place bounds on the extent and time scale of the motion. The bending motion of the loop is slow compared to the rapid fluctuations of individual dihedral angles, reflecting the collective nature and largely entropic free energy profile for hinge bending. The dynamics of the central binding loop in LSIII illustrates two distinct mechanisms by which molecular dynamics directly impacts biological activity. The relative rigidity of key residues involved in recognition at the tip of the central binding loop lowers the otherwise substantial entropic cost of binding. Large excursions of the loop hinge angle may endow the protein with structural plasticity, allowing it to adapt to conformational changes induced in the receptor.     

Molecular dynamics simulation of the interaction between the complex iron-sulfur flavoprotein glutamate synthase and its substrates

Glutamate synthase (GltS) is a complex iron-sulfur flavoprotein that catalyzes the reductive transfer of L-glutamine amide group to the C2 carbon of 2-oxoglutarate yielding two molecules of L-glutamate. Molecular dynamics calculations in explicit solvent were carried out to gain insight into the conformational flexibility of GltS and into the role played by the enzyme substrates in regulating the catalytic cycle. We have modelled the free (unliganded) form of Azospirillum brasilense GltS alpha subunit and the structure of the reduced enzyme in complex with the L-glutamine and 2-oxoglutarate substrates starting from the crystallographically determined coordinates of the GltS alpha subunit in complex with L-methionine sulphone and 2-oxoglutarate. The present 4-ns molecular dynamics calculations reveal that the GltS glutaminase site may exist in a catalytically inactive conformation unable to bind glutamine, and in a catalytically competent conformation, which is stabilized by the glutamine substrate. Substrates binding also induce (1) closure of the loop formed by residues 263-271 with partial shielding of the glutaminase site from solvent, and (2) widening of the ammonia tunnel entrance at the glutaminase end to allow for ammonia diffusion toward the synthase site. The Q-loop of glutamate synthase, which acts as an active site lid in other amidotransferases, seems to maintain an open conformation. Finally, binding of L-methionine sulfone, a glutamine analog that mimics the tetrahedral transient species occurring during its hydrolysis, causes a coordinated rigid-body motion of segments of the glutaminase domain that results in the inactive conformation observed in the crystal structure of GltS alpha subunit.     

Effects of ligand binding on the dynamics of rice nonspecific lipid transfer protein 1: a model from molecular simulations

Plant nonspecific lipid transfer proteins (nsLTPs) are small, basic proteins constituted mainly of alpha-helices and stabilized by four conserved disulfide bridges. They are characterized by the presence of a tunnel-like hydrophobic cavity, capable of transferring various lipid molecules between lipid bilayers in vitro. In this study, molecular dynamics (MD) simulations were performed at room temperature to investigate the effects of lipid binding on the dynamic properties of rice nsLTP1. Rice nsLTP1, either in the free form or complexed with one or two lipids was subjected to MD simulations. The C-terminal loop was very flexible both before and after lipid binding, as revealed by calculating the root-mean-square fluctuation. After lipid binding, the flexibility of some residues that were not in direct contact with lipid molecules increased significantly, indicating an increase of entropy in the region distal from the binding site. Essential dynamics analysis revealed clear differences in motion between unliganded and liganded rice nsLTP1s. In the free form of rice nsLTP1, loop1 exhibited the largest directional motion. This specific essential motion mode diminished after binding one or two lipid molecules. To verify the origin of the essential motion observed in the free form of rice nsLTP1, we performed multiple sequence alignments to probe the intrinsic motion encoded in the primary sequence. We found that the amino acid sequence of loop1 is highly conserved among plant nsLTP1s, thus revealing its functional importance during evolution. Furthermore, the sequence of loop1 is composed mainly of amino acids with short side chains. In this study, we show that MD simulations, together with essential dynamics analysis, can be used to determine structural and dynamic differences of rice nsLTP1 upon lipid binding.     

Investigation of the mechanism of domain closure in citrate synthase by molecular dynamics simulation

Six, 2 ns molecular dynamics simulations have been performed on the homodimeric enzyme citrate synthase. In three, both monomers were started from the open, unliganded X-ray conformation. In the remaining three, both monomers started from a closed, liganded X-ray conformation, with the ligands removed. Projecting the motion from the simulations onto the experimental domain motion revealed that the free-energy profile is rather flat around the open conformation, with steep sides. The most closed conformations correspond to hinge-bending angles of 12-14 compared to the 20 degrees that occurs upon the binding of oxaloacetate. It is also found that the open, unliganded X-ray conformation is situated at the edge of the steep rise in free energy, although conformations that are about 5 degrees more open were sampled. A rigid-body essential dynamics analysis of the combined open trajectories has shown that domain motions in the direction of the closed X-ray conformation are compatible with the natural domain motion of the unliganded protein, which has just two main degrees of freedom. The simulations starting from the closed conformation suggest a free-energy profile with a small barrier in going from the closed to open conformation. A combined essential dynamics and hinge-bending analysis of a trajectory that spontaneously converts from the closed to open state shows an almost exact correspondence to the experimental transition that occurs upon ligand binding. The simulations support the conclusion from an earlier analysis of the experimental transition that the beta-hairpin acts as a mechanical hinge by attaching the small domain to the large domain through a conserved main-chain hydrogen bond and salt-bridges, and allowing rotation to occur via its two flexible termini. The results point to a mechanism of domain closure in citrate synthase that has analogy to the process of closing a door.     

Building a foundation for structure‐based cellulosome design for cellulosic ethanol: Insight into cohesin‐dockerin complexation from computer simulation

The organization and assembly of the cellulosome, an extracellular multienzyme complex produced by anaerobic bacteria, is mediated by the high‐affinity interaction of cohesin domains from scaffolding proteins with dockerins of cellulosomal enzymes. We have performed molecular dynamics simulations and free energy calculations on both the wild type (WT) and D39N mutant of the C. thermocellum Type I cohesin‐dockerin complex in aqueous solution. The D39N mutation has been experimentally demonstrated to disrupt cohesin‐dockerin binding. The present MD simulations indicate that the substitution triggers significant protein flexibility and causes a major change of the hydrogen‐bonding network in the recognition strips—the conserved loop regions previously proposed to be involved in binding—through electrostatic and salt‐bridge interactions between β‐strands 3 and 5 of the cohesin and α‐helix 3 of the dockerin. The mutation‐induced subtle disturbance in the local hydrogen‐bond network is accompanied by conformational rearrangements of the protein side chains and bound water molecules. Additional free energy perturbation calculations of the D39N mutation provide differences in the cohesin‐dockerin binding energy, thus offering a direct, quantitative comparison with experiments. The underlying molecular mechanism of cohesin‐dockerin complexation is further investigated through the free energy profile, that is, potential of mean force (PMF) calculations of WT cohesin‐dockerin complex. The PMF shows a high‐free energy barrier against the dissociation and reveals a stepwise pattern involving both the central β‐sheet interface and its adjacent solvent‐exposed loop/turn regions clustered at both ends of the β‐barrel structure.     

Dynamics of the active site loops in catalyzing aminoacylation reaction in seryl and histidyl tRNA synthetases

Aminoacylation reaction is the first step of protein biosynthesis. The catalytic reorganization at the active site of aminoacyl tRNA synthetases (aaRSs) is driven by the loop motions. There remain lacunae of understanding concerning the catalytic loop dynamics in aaRSs. We analyzed the functional loop dynamics in seryl tRNA synthetase from Methanopyrus kandleri (^mkSerRS) and histidyl tRNA synthetases from Thermus thermophilus (^ttHisRS), respectively, using molecular dynamics. Results confirm that the motif 2 loop and other active site loops are flexible spots within the catalytic domain. Catalytic residues of the loops form a network of interaction with the substrates to form a reactive state. The loops undergo transitions between closed state and open state and the relaxation of the constituent residues occurs in femtosecond to nanosecond time scale. Order parameters are higher for constituent catalytic residues which form a specific network of interaction with the substrates to form a reactive state compared to the Gly residues within the loop. The development of interaction is supported from mutation studies where the catalytic domain with mutated loop exhibits unfavorable binding energy with the substrates. During the open-close motion of the loops, the catalytic residues make relaxation by ultrafast librational motion as well as fast diffusive motion and subsequently relax rather slowly via slower diffusive motion. The Gly residues act as a hinge to facilitate the loop closing and opening by their faster relaxation behavior. The role of bound water is analyzed by comparing implicit solvent-based and explicit solvent-based simulations. Loops fail to form catalytically competent geometry in absence of water. The present result, for the first time reveals the nature of the active site loop dynamics in aaRS and their influence on catalysis.     

Multiple Molecular Dynamics Simulations of TEM β-Lactamase: Dynamics and Water Binding of the Ω-Loop

The Ω-loop of TEM β-lactamase is involved in substrate recognition and catalysis. Its dynamical properties and interaction with water molecules were investigated by performing multiple molecular dynamics simulations of up to 50 ns. Protein flexibility was assessed by calculating the root mean-square fluctuations and the generalized order parameter, S². The residues in secondary structure elements are highly ordered, whereas loop regions are more flexible, which is in agreement with previous experimental observations. Interestingly, the Ω-loop (residues 161-179) is rigid with order parameters similar to secondary structure elements, with the exception of the tip of the loop (residues 173-177) that has a considerably higher flexibility and performs an opening and closing motion on the 50-ns timescale. The rigidity of the main part of the Ω-loop is mediated by stabilizing and highly conserved water bridges inside a cavity lined by the Ω-loop and residues 65-69 of the protein core. In contrast, the flexible tip of the Ω-loop lacks these interactions. Hydration of the cavity and exchange of the water molecules with the bulk solvent occurs via two pathways: the flexible tip that serves as a door to the cavity, and a temporary water channel involving the side chain of Arg¹⁶⁴.     

Molecular dynamics simulations of human carbonic anhydrase II: Insight into experimental results and the role of solvation

In this paper, the carbonic anhydrase II (CA II) enzyme active site is modeled using ab initio calculations and molecular dynamics simulations to examine a number of important issues for the enzyme function. It is found that the Zn²⁺ ion is dominantly tetrahedrally coordinated, which agrees with X-ray crystallographic studies. However, a transient five-fold coordination with an extra water molecule is also found. Studies of His64 conformations upon a change in the protonation states of the Zn-bound water and the His64 residue also confirm the results of an X-ray study which suggest that the His64 conformation is quite flexible. However, the degree of water solvation is found to affect this behavior. Water bridge formation between the Zn-bound water and the His64 residue was found to involve a free energy barrier of 2–3 kcal/mol and an average lifetime of several picoseconds, which supports the concept of a proton transfer mechanism through such a bridge. Mutations of various residues around the active site provide further insight into the corresponding experimental results and, in fact, suggest an important role for the solvent water molecules in the CA II catalytic mechanism. Proteins 33:119–134, 1998. © 1998 Wiley-Liss, Inc.     

Role of conformational fluctuations in the enzymatic reaction of HIV-1 protease

The emergence of compensatory drug-resistant mutations in HIV-1 protease challenges the common view of the reaction mechanism of this enzyme. Here, we address this issue by performing classical and ab initio molecular dynamics simulations (MD) on a complex between the enzyme and a peptide substrate. The classical MD calculation reveals large-scale protein motions involving the flaps and the cantilever. These motions modulate the conformational properties of the substrate at the cleavage site. The ab initio calculations show in turn that substrate motion modulates the activation free energy barrier of the enzymatic reaction dramatically. Thus, the catalytic power of the enzyme does not arise from the presence of a pre-organized active site but from the protein mechanical fluctuations. The implications of this finding for the emergence of drug-resistance are discussed.     

Docking study and free energy simulation of the complex between p53 DNA-binding domain and azurin

Molecular interaction between p53 tumor suppressor and the copper protein azurin (AZ) has been demonstrated to enhance p53 stability and hence antitumoral function, opening new perspectives in cancer treatment. While some experimental work has provided evidence for AZ binding to p53, no crystal structure for the p53-AZ complex was solved thus far. In this work the association between AZ and the p53 DNA-binding domain (DBD) was investigated by computational methods. Using a combination of rigid-body protein docking, experimental mutagenesis information, and cluster analysis 10 main p53 DBD-AZ binding modes were generated. The resulting structures were further characterized by molecular dynamics (MD) simulations and free energy calculations. We found that the highest scored docking conformation for the p53 DBD-AZ complex also yielded the most favorable free energy value. This best three-dimensional model for the complex was validated by using a computational mutagenesis strategy. In this structure AZ binds to the flexible L(1) and s(7)-s(8) loops of the p53 DBD and stabilizes them through protein-protein tight packing interactions, resulting in high degree of both surface matching and electrostatic complementarity.     

Calculations of free energy profiles for the staphylococcal nuclease catalyzed reaction

Calculations of the free energy profile for the first two (rate-limiting) steps of the staphylococcal nuclease catalyzed reaction are reported. The calculations are based on the empirical valence bond method in combination with free energy perturbation molecular dynamics simulations. The calculated activation free energy is in good agreement with experimental kinetic data, and the catalytic effect of the enzyme is reproduced without any arbitrary adjustment of parameters. The enormous reduction of the activation barrier (relative to the reference reaction in water) appears to be largely associated with the strong electrostatic effect of the Ca2+ ion and the two arginine residues in the active site. This favorable electrostatic environment reduces the cost of the general-base catalysis step by almost 15 kcal/mol (by stabilizing the OH- nucleophile) and then stabilizes the developing negative charge on the 5'-phosphate group in the second step of the reaction by about 19 kcal/mol. The basic features of the originally postulated enzyme mechanism (Cotton et al., 1979) are found to be compatible with the observed activation free energy. However, the proposed modification of the mechanism (Sepersu et al., 1987), in which Arg 87 interacts only with the pentacoordinated transition state, is supported by the simulations. Further calculations on the D21E mutant also give results in good agreement with kinetic data.     

Protein dynamics. A time-resolved fluorescence, energetic and molecular dynamics study of ribonuclease T1

Studies using time-resolved fluorescence depolarization were performed on the internal motion of Trp 59 of ribonuclease T1 (EC 3.1.27.3) in the free enzyme, 2'-GMP-enzyme complex and 3'-GMP-enzyme complex. The Trp 59 motion was also studied in the free enzyme using molecular dynamics simulations. Energetic analysis of activation barriers to the Trp 59 motion was performed using both the transition state theory and Kramers' theory. The activation parameters showed a dependence on solvent viscosity indicating the transition state approach in aqueous solution to be inadequate. When taking solvent viscosity contributions into account agreement between the transition state and Kramers' theories was obtained. The results indicate the three enzyme forms to have different conformations with the free enzyme and 3'-GMP-enzyme complex being similar. Comparison of the experimental and theoretical results showed a good agreement on the Trp 59 motion in the free enzyme. Trp 59 appears to vibrate rapidly, with a relaxation time of the order of 1 ps, within free space in the protein matrix and to have a slower motion, with a relaxation time of the order of 100 ps, which is related to breathing of the surrounding protein matrix. Molecular dynamics results indicate high mobility in regions of the enzyme involved in the interaction with the guanine base of the inhibitor or substrate while much lower mobility occurred in residues involved in the catalytic mechanism of ribonuclease T1.     

Molecular insight into pseudolysin inhibition using the MM-PBSA and LIE methods

Pseudolysin, the extracellullar elastase of Pseudomonas aeruginosa (EC: 3.4.24.26) plays an important role in the pathogenesis of P. aeruginosa infections. In the present study, molecular dynamics simulations and theoretical affinity predictions were used to gain molecular insight into pseudolysin inhibition. Four low molecular weight inhibitors were docked at their putative binding sites and molecular dynamics (MD) simulations were performed for 5.0 ns, and the free energy of binding was calculated by the linear interaction energy method. The number and the contact surface area of stabilizing hydrophobic, aromatic, and hydrogen bonding interactions appears to reflect the affinity differences between the inhibitors. The proteinaceous inhibitor, Streptomyces metalloproteinase inhibitor (SMPI) was docked in three different binding positions and MD simulations were performed for 3.0 ns. The MD trajectories were used for molecular mechanics-Poisson-Boltzmann surface area analysis of the three binding positions. Computational alanine scanning of the average pseudolysin-SMPI complexes after MD revealed residues at the pseudolysin-SMPI interface giving the main contribution to the free energy of binding. The calculations indicated that SMPI interacts with pseudolysin via the rigid active site loop, but that also contact sites outside this loop contribute significantly to the free energy of association.     

On the unfolding of alpha-lytic protease and the role of the pro region

Molecular dynamics simulations of alpha-lytic protease (alphaLP) alone and complexed with its pro region (PRO) are performed to understand the origin of its high unfolding (and folding) barrier when it is alone and how the pro region lowers this barrier. At room temperature, alphaLP exhibits lower dynamic fluctuations than alpha-chymotrypsin. Simulation of PRO alone led to reorientation of its N terminal helix and collapse to a more compact state. A model for the uncleaved proenzyme was built and found to be stable in the time scale of the simulations. Energetic analysis suggests that the origin of strain in the uncleaved proenzyme compared with the cleaved complex is in the intramolecular backbone electrostatic interactions of the cleaved strand. In high temperature simulations, the interaction of the long beta hairpin of the enzyme with the C terminal beta sheet of PRO is among the most stable in the complex and a likely "nucleation site" for folding. In the course of unfolding, the C terminal tail of PRO is sometimes observed to intervene between the long hairpin and the aspartate loop of the enzyme, perhaps thereby lowering the energy barrier for separation of the two hairpins. Tighter interactions at the interface between the enzyme and its pro region are also occasionally observed, providing an additional mechanism for unfolding catalysis. Simulations of a mutant enzyme where the buried ion pair residues R102 and D142 were replaced by W and L, respectively, did not display any distinguishable behavior compared with the wild type.     

Crystal structure analysis,covalent docking,and molecular dynamics calculations reveal a conformational switch in PhaZ7 PHB depolymerase

An open and a closed conformation of a surface loop in PhaZ7 extracellular poly(3‐hydroxybutyrate) depolymerase were identified in two high‐resolution crystal structures of a PhaZ7 Y105E mutant. Molecular dynamics (MD) simulations revealed high root mean square fluctuations (RMSF) of the 281–295 loop, in particular at residue Asp289 (RMSF 7.62 Å). Covalent docking between a 3‐hydroxybutyric acid trimer and the catalytic residue Ser136 showed that the binding energy of the substrate is significantly more favorable in the open loop conformation compared to that in the closed loop conformation. MD simulations with the substrate covalently bound depicted 1 Å RMSF higher values for the residues 281–295 in comparison to the apo (substrate‐free) form. In addition, the presence of the substrate in the active site enhanced the ability of the loop to adopt a closed form. Taken together, the analysis suggests that the flexible loop 281–295 of PhaZ7 depolymerase can act as a lid domain to control substrate access to the active site of the enzyme. Proteins 2017; 85:1351–1361. © 2017 Wiley Periodicals, Inc.     

Flexibility and charge asymmetry in the activation loop of Src tyrosine kinases

Regulated activity of Src kinases is critical for cell growth. Src kinases can be activated by trans-phosphorylation of a tyrosine located in the central activation loop of the catalytic domain. However, because the required exposure of this tyrosine is not observed in the down-regulated X-ray structures of Src kinases, transient partial opening of the activation loop appears to be necessary for such processes. Umbrella sampling molecular dynamics simulations are used to characterize the free energy landscape of opening of the hydrophilic part of the activation loop in the Src kinase Hck. The loop prefers a partially open conformation where Tyr416 has increased accessibility, but remains partly shielded. An asymmetric distribution of the charged residues in the sequence near Tyr416, which contributes to shielding, is found to be conserved in Src family members. A conformational equilibrium involving exchange of electrostatic interactions between the conserved residues Glu310 and Arg385 or Arg409 affects activation loop opening. A mechanism for access of unphosphorylated Tyr416 into an external catalytic site is suggested based on these observations.     

Understanding the molecular basis of MK2-p38α signaling complex assembly: insights into protein-protein interaction by molecular dynamics and free energy studies

The formation of a p38 MAPK and MAPK-activated protein kinase 2 (MK2) signaling complex is physiologically relevant to cellular responses such as the proinflammatory cytokine production. The interaction between p38α isoform and MK2 is of great importance for this signaling. In this study, molecular dynamics simulation and binding free energy calculation were performed on the MK2-p38α signaling complex to investigate the protein-protein interaction between the two proteins. Dynamic domain motion analyses were performed to analyze the conformational changes between the unbound and bound states of proteins during the interaction. The activation loop, αF-I helices, and loops among α helices in the C-lobe of MK2 are found to be highly flexible and exhibit significant changes upon p38α binding. The results also show that after the binding of p38α, the N- and C-terminal domains of MK2 display an opening and twisting motion centered on the activation loop. The molecular mechanics Poisson-Boltzmann and generalized-Born surface area (MM-PB/GBSA) methods were used to calculate binding free energies between MK2 and p38α. The analysis of the components of binding free energy calculation indicates that the van der Waals interaction and the nonpolar solvation energy provide the driving force for the binding process, while the electrostatic interaction contributes critically to the specificity, rather than to MK2-p38α binding affinity. The contribution of each residue at the interaction interface to the binding affinity of MK2 with p38α was also analyzed by free energy decomposition. Several important residues responsible for the protein-protein interaction were also identified.