Targeted molecular dynamics simulation studies of binding and conformational changes in E. coli MurD

Enzymes involved in the biosynthesis of bacterial peptidoglycan, an essential cell wall polymer unique to prokaryotic cells, represent a highly interesting target for antibacterial drug design. Structural studies of E. coli MurD, a three-domain ATP hydrolysis driven muramyl ligase revealed two inactive open conformations of the enzyme with a distinct C-terminal domain position. It was hypothesized that the rigid body rotation of this domain brings the enzyme to its closed active conformation, a structure, which was also determined experimentally. Targeted molecular dynamics 1 ns-length simulations were performed in order to examine the substrate binding process and gain insight into structural changes in the enzyme that occur during the conformational transitions into the active conformation. The key interactions essential for the conformational transitions and substrate binding were identified. The results of such studies provide an important step toward more powerful exploitation of experimental protein structures in structure-based inhibitor design.     

MD simulations reveal alternate conformations of the oxyanion hole in the Zika virus NS2B/NS3 protease

Recent crystallography studies have shown that the binding site oxyanion hole plays an important role in inhibitor binding, but can exist in two conformations (active/inactive). We have undertaken molecular dynamics (MD) calculations to better understand oxyanion hole dynamics and thermodynamics. We find that the Zika virus (ZIKV) NS2B/NS3 protease maintains a stable closed conformation over multiple 100-ns conventional MD simulations in both the presence and absence of inhibitors. The S1, S2, and S3 pockets are stable as well. However, in two of eight simulations, the A132-G133 peptide bond in the binding pocket of S1' spontaneously flips to form a 3₁₀-helix that corresponds to the inactive conformation of the oxyanion hole, and then maintains this conformation until the end of the 100-ns conventional MD simulations without inversion of the flip. This conformational change affects the S1' pocket in ZIKV NS2B/NS3 protease active site, which is important for small molecule binding. The simulation results provide evidence at the atomic level that the inactive conformation of the oxyanion hole is more favored energetically when no specific interactions are formed between substrate/inhibitor and oxyanion hole residues. Interestingly, however, transition between the active and inactive conformation of the oxyanion hole can be observed by boosting the valley potential in accelerated MD simulations. This supports a proposed induced-fit mechanism of ZIKV NS2B/NS3 protease from computational methods and provides useful direction to enhance inhibitor binding predictions in structure-based drug design.     

Dynamic flaps in HIV-1 protease adopt unique ordering at different stages in the catalytic cycle

The flexibility of HIV-1 protease flaps is known to be essential for the enzymatic activity. Here we attempt to capture a multitude of conformations of the free and substrate-bound HIV-1 protease that differ drastically in their flap arrangements. The substrate binding process suggests the opening of active site gate in conjunction with a reversal of flap tip ordering, from the native semiopen state. The reversed-flap, open-gated enzyme readily transforms to a closed conformation after proper placement of the substrate into the binding cleft. After substrate processing, the closed state protease which possessed opposite flap ordering relative to the semiopen state, encounters another flap reversal via a second open conformation that facilitates the evolution of native semiopen state of correct flap ordering. The complicated transitional pathway, comprising of many high and low energy states, is explored by combining standard and activated molecular dynamics (MD) simulation techniques. The study not only complements the existing findings from X-ray, NMR, EPR, and MD studies but also provides a wealth of detailed information that could help the structure-based drug design process.     

Naphthalene-based RNA editing inhibitor blocks RNA editing activities and editosome assembly in Trypanosoma brucei

RNA editing, catalyzed by the multiprotein editosome complex, is an essential step for the expression of most mitochondrial genes in trypanosomatid pathogens. It has been shown previously that Trypanosoma brucei RNA editing ligase 1 (TbREL1), a core catalytic component of the editosome, is essential in the mammalian life stage of these parasitic pathogens. Because of the availability of its crystal structure and absence from human, the adenylylation domain of TbREL1 has recently become the focus of several studies for designing inhibitors that target its adenylylation pocket. Here, we have studied new and existing inhibitors of TbREL1 to better understand their mechanism of action. We found that these compounds are moderate to weak inhibitors of adenylylation of TbREL1 and in fact enhance adenylylation at higher concentrations of protein. Nevertheless, they can efficiently block deadenylylation of TbREL1 in the editosome and, consequently, result in inhibition of the ligation step of RNA editing. Further experiments directly showed that the studied compounds inhibit the interaction of the editosome with substrate RNA. This was supported by the observation that not only the ligation activity of TbREL1 but also the activities of other editosome proteins such as endoribonuclease, terminal RNA uridylyltransferase, and uridylate-specific exoribonuclease, all of which require the interaction of the editosome with the substrate RNA, are efficiently inhibited by these compounds. In addition, we found that these compounds can interfere with the integrity and/or assembly of the editosome complex, opening the exciting possibility of using them to study the mechanism of assembly of the editosome components.     

In silico analyses of substrate interactions with human serum paraoxonase 1

Human paraoxonase (HuPON1) is a serum enzyme that exhibits a broad spectrum of hydrolytic activities, including the hydrolysis of various organophosphates, esters, and recently identified lactone substrates. Despite intensive site-directed mutagenesis and other biological studies, the structural basis for the specificity of substrate interactions of HuPON1 remains elusive. In this study, we apply homology modeling, docking, and molecular dynamic (MD) simulations to probe the binding interactions of HuPON1 with representative substrates. The results suggest that the active site of HuPON1 is characterized by two distinct binding regions: the hydrophobic binding site for arylesters/lactones, and the paraoxon binding site for phosphotriesters. The unique binding modes proposed for each type of substrate reveal a number of key residues governing substrate specificity. The polymorphic residue R/Q192 interacts with the leaving group of paraoxon, suggesting it plays an important role in the proper positioning of this substrate in the active site. MD simulations of the optimal binding complexes show that residue Y71 undergoes an "open-closed" conformational change upon ligand binding, and forms strong interactions with substrates. Further binding free energy calculations and residual decomposition give a more refined molecular view of the energetics and origin of HuPON1/substrate interactions. These studies provide a theoretical model of substrate binding and specificity associated with wild type and mutant forms of HuPON1, which can be applied in the rational design of HuPON1 variants as bioscavengers with enhanced catalytic activity.     

Microsomal prostaglandin E synthase-1 exhibits one-third-of-the-sites reactivity

mPGES-1 (microsomal prostaglandin E synthase-1) is a newly recognized target for the treatment of inflammatory diseases. As the terminal enzyme of the prostaglandin production pathway, mPGES-1 inhibition may have a low risk of side effects. Inhibitors of mPGES-1 have attracted considerable attention as next-generation anti-inflammatory drugs. However, as mPGES-1 is a membrane protein, its enzymatic mechanism remains to be disclosed fully. We used MD (molecular dynamics) simulations, mutation analysis, hybrid experiments and co-IP (co-immunoprecipitation) to investigate the conformation transitions of mPGES-1 during catalysis. mPGES-1 forms a homotrimer with three substrate-binding sites (pockets). In the MD simulation, only one substrate molecule could bind to one of the pockets and form the active complex, suggesting that the mPGES-1 trimer has only one pocket active at any given time. This one-third-of-the-sites reactivity enzyme mechanism was verified further by hybridization experiments and MD simulations. The results of the present study revealed for the first time a novel one-third-of-the-sites reactivity enzyme mechanism for mPGES-1, and the unique substrate-binding pocket in our model constituted an active conformation that was suitable for further enzymatic mechanism study and structural-based drug design against mPGES-1.     

Cofactor assisted gating mechanism in the active site of NADH oxidase from Thermus thermophilus

NADH oxidase (NOX) from Thermus thermophilus is a member of a structurally homologous flavoprotein family of nitroreductases and flavin reductases. The importance of local conformational dynamics in the active site of NOX has been recently demonstrated. The enzyme activity was increased by 250% in the presence of 1 M urea with no apparent perturbation of the native structure of the protein. The present in silico results correlate with the in vitro data and suggest the possible explanation about the effect of urea on NOX activity at the molecular level. Both, X-ray structure and molecular dynamics (MD) simulations, show open conformation of the active site represented by approximately 0.9 nm distance between the indole ring of Trp47 and the isoalloxazine ring of FMN412. In this conformation, the substrate molecule can bind in the active site without sterical restraints. MD simulations also indicate more stable conformation of the active site called "closed" conformation. In this conformation, Trp47 and the isoalloxazine ring of FMN412 are so close to each other (approximately 0.5 nm) that the substrate molecule is unable to bind between them without perturbing this conformation. The open/close transition of the active site between Trp47 and the flavin ring is accompanied by release of the "tightly" bound water molecule from the active site--cofactor assisted gating mechanism. The presence of urea in aqueous solutions of NOX prohibits closing of the active site and even unlocks the closed active site because of the concomitant binding of a urea molecule in the active site cavity. The binding of urea in the active site is stabilized by formation of one/two persistent hydrogen bonds involving the carbonyl group of the urea molecule. Our report represents the first MD study of an enzyme from the novel flavoprotein family of nitroreductases and flavin reductases. The common occurrence of aromatic residues covering the active sites in homologous enzymes suggests the possibility of a general gating mechanism and the importance of local dynamics within this flavoprotein family.     

Structural insight into Mycobacterium tuberculosis maltosyl transferase inhibitors: pharmacophore-based virtual screening,docking, and molecular dynamics simulations

Pharmacophore-based virtual screening, subsequent docking, and molecular dynamics (MD) simulations have been done to identify potential inhibitors of maltosyl transferase of Mycobacterium tuberculosis (mtb GlgE). Ligand and structure-based pharmacophore models representing its primary binding site (pbs) and unique secondary binding site 2 (sbs2), respectively, were constructed based on the three dimensional structure of mtb GlgE. These pharmacophore models were further used for screening of ZINC and antituberculosis compounds database (ATD). Virtually screened molecules satisfying Lipinski’s rule of five were then analyzed using docking studies and have identified 23 molecules with better binding affinity than its natural substrate, maltose. Four top scoring ligands from ZINC and ATD that either binds to pbs or sbs2 have been subjected to 10 ns each MD simulations and binding free energy calculations. Results of these studies have confirmed stable protein ligand binding. Results reported in the article are likely to be helpful in antitubercular therapeutic development research.     

Large-scale conformational changes of Trypanosoma cruzi proline racemase predicted by accelerated molecular dynamics simulation

Chagas' disease, caused by the protozoan parasite Trypanosoma cruzi (T. cruzi), is a life-threatening illness affecting 11-18 million people. Currently available treatments are limited, with unacceptable efficacy and safety profiles. Recent studies have revealed an essential T. cruzi proline racemase enzyme (TcPR) as an attractive candidate for improved chemotherapeutic intervention. Conformational changes associated with substrate binding to TcPR are believed to expose critical residues that elicit a host mitogenic B-cell response, a process contributing to parasite persistence and immune system evasion. Characterization of the conformational states of TcPR requires access to long-time-scale motions that are currently inaccessible by standard molecular dynamics simulations. Here we describe advanced accelerated molecular dynamics that extend the effective simulation time and capture large-scale motions of functional relevance. Conservation and fragment mapping analyses identified potential conformational epitopes located in the vicinity of newly identified transient binding pockets. The newly identified open TcPR conformations revealed by this study along with knowledge of the closed to open interconversion mechanism advances our understanding of TcPR function. The results and the strategy adopted in this work constitute an important step toward the rationalization of the molecular basis behind the mitogenic B-cell response of TcPR and provide new insights for future structure-based drug discovery.     

Unique MAP Kinase binding sites

Map kinases are drug targets for autoimmune disease, cancer, and apoptosis-related diseases. Drug discovery efforts have developed MAP kinase inhibitors directed toward the ATP binding site and neighboring "DFG-out" site, both of which are targets for inhibitors of other protein kinases. On the other hand, MAP kinases have unique substrate and small molecule binding sites that could serve as inhibition sites. The substrate and processing enzyme D-motif binding site is present in all MAP kinases, and has many features of a good small molecule binding site. Further, the MAP kinase p38alpha has a binding site near its C-terminus discovered in crystallographic studies. Finally, the MAP kinases ERK2 and p38alpha have a second substrate binding site, the FXFP binding site that is exposed in active ERK2 and the D-motif peptide induced conformation of MAP kinases. Crystallographic evidence of these latter two binding sites is presented.     

Regulatory phosphorylation of cyclin-dependent kinase 2: insights from molecular dynamics simulations

The structures of fully active cyclin-dependent kinase-2 (CDK2) complexed with ATP and peptide substrate, CDK2 after the catalytic reaction, and CDK2 inhibited by phosphorylation at Thr14/Tyr15 were studied using molecular dynamics (MD) simulations. The structural details of the CDK2 catalytic site and CDK2 substrate binding box were described. Comparison of MD simulations of inhibited complexes of CDK2 was used to help understand the role of inhibitory phosphorylation at Thr14/Tyr15. Phosphorylation at Thr14/Tyr15 causes ATP misalignment for the phosphate-group transfer, changes in the Mg²⁺ coordination sphere, and changes in the H-bond network formed by CDK2 catalytic residues (Asp127, Lys129, Asn132). The inhibitory phosphorylation causes the G-loop to shift from the ATP binding site, which leads to opening of the CDK2 substrate binding box, thus probably weakening substrate binding. All these effects explain the decrease in kinase activity observed after inhibitory phosphorylation at Thr14/Tyr15 in the G-loop. Interaction of the peptide substrate, and the phosphorylated peptide product, with CDK2 was also studied and compared. These results broaden hypotheses drawn from our previous MD studies as to why a basic residue (Arg/Lys) is preferred at the P₊₂ substrate position. Figure View of the substrate binding site of the fully active cyclin-dependent kinase-2 (CDK2) (pT160-CDK2/cyclin A/ATP). The pThr160 activation site is located in the T-loop (yellow secondary structure). The G-loop, which partly forms the ATP binding site, is shown in blue. The Thr14 and Tyr15 inhibitory phosphorylation sites located in the G-loop are shown in licorice representation     

Genetic and structural characterization of PvSERA4: potential implication as therapeutic target for Plasmodium vivax malaria

Plasmodium vivax malaria is geographically the most widely distributed and prevalent form of human malaria. The development of drug resistance by the parasite to existing drugs necessitates higher focus to explore and identify new drug targets. Plasmodial proteases have key roles in parasite biology and are involved in nutritional uptake, egress from infected reticulocytes, and invasion of the new target erythrocytes. Serine repeat antigens (SERA) of Plasmodium are parasite proteases that remain attractive drug targets and are important vaccine candidates due to their high expression profiles in the blood stages. SERA proteins have a unique putative papain-like cysteine protease motif that has either serine or cysteine in its active site. In P. vivax, PvSERA4 is the highest transcribed member of this multigene family. In this study, we have investigated the genetic polymorphism of PvSERA4 central protease domain and deduced its 3D model by homology modeling and also performed MD simulations to acquire refined protein structure. Sequence analysis of protease domain of PvSERA4 from Indian field isolates reveals that the central domain is highly conserved. The high sequence conservation of the PvSERA4 enzyme domain coupled with its high expression raises the possibility of it having a critical role in parasite biology and hence, being a reliable target for new selective inhibitor-based antimalarial chemotherapeutics. The 3D model showed the presence of an unusual antiparallel Beta hairpin motif between catalytic residues similar to hemoglobin binding motif of Plasmodial hemoglobinases. Our PvSERA4 model will aid in designing structure-based inhibitors against this enzyme.     

Impact of tetramerization on the ligand recognition of N1 influenza neuraminidase via MMGBSA approach

Influenza virus neuraminidase (NA) is a homotetrameric surface protein that, in contrast to other non-influenza NAs, requires a quaternary assembly to exhibit enzymatic activity, suggesting that the oligomeric state significantly impacts the active site of influenza NA. Nevertheless, most structure-based drug design studies have been reported by employing the monomeric state in the closed or open-loop due to the computational cost of employing the tetrameric NA. In this work, we present MD simulations coupled to the MMGBSA approach of avian N1 type NA in its monomeric and tetrameric closed and open-loop state both with and without the inhibitor oseltamivir and its natural substrate, sialic acid. Structural and energetic analyses revealed that the tetrameric state impacts flexibility as well as the map of interactions participating in stabilizing the protein–ligand complexes with respect to the monomeric state. It was observed that the tetrameric state exerts dissimilar effects in binding affinity, characteristic of positive and negative cooperativity for oseltamivir and sialic acid, respectively. Based on our results, to perform a confident structure-based drug design, as well as to evaluate the impact of key mutations through MD simulations, it is important to consider the tetrameric state closed-loop state.     

How an Inhibitor Bound to Subunit Interface Alters Triosephosphate Isomerase Dynamics

The tunnel region at triosephosphate isomerase (TIM)’s dimer interface, distant from its catalytic site, is a target site for certain benzothiazole derivatives that inhibit TIM’s catalytic activity in Trypanosoma cruzi, the parasite that causes Chagas disease. We performed multiple 100-ns molecular-dynamics (MD) simulations and elastic network modeling (ENM) on both apo and complex structures to shed light on the still unclear inhibitory mechanism of one such inhibitor, named bt10. Within the time frame of our MD simulations, we observed stabilization of aromatic clusters at the dimer interface and enhancement of intersubunit hydrogen bonds in the presence of bt10, which point to an allosteric effect rather than destabilization of the dimeric structure. The collective dynamics dictated by the topology of TIM is known to facilitate the closure of its catalytic loop over the active site that is critical for substrate entrance and product release. We incorporated the ligand’s effect on vibrational dynamics by applying mixed coarse-grained ENM to each one of 54,000 MD snapshots. Using this computationally efficient technique, we observed altered collective modes and positive shifts in eigenvalues due to the constraining effect of bt10 binding. Accordingly, we observed allosteric changes in the catalytic loop’s dynamics, flexibility, and correlations, as well as the solvent exposure of catalytic residues. A newly (to our knowledge) introduced technique that performs residue-based ENM scanning of TIM revealed the tunnel region as a key binding site that can alter global dynamics of the enzyme.     

Role of a mutated residue at the entrance of the substrate access channel in cytochrome p450 engineered for vitamin D(3) hydroxylation activity

The cytochrome P450 enzyme engineered for enhancement of vitamin D(3) (VD(3)) hydroxylation activity, Vdh-K1, includes four mutations (T70R, V156L, E216M, and E384R) compared to the wild-type enzyme. Plausible roles for V156L, E216M, and E384R have been suggested by crystal structure analysis (Protein Data Bank 3A50 ), but the role of T70R, which is located at the entrance of the substrate access channel, remained unclear. In this study, the role of the T70R mutation was investigated by using computational approaches. Molecular dynamics (MD) simulations and steered molecular dynamics (SMD) simulations were performed, and differences between R70 and T70 were compared in terms of structural change, binding free energy change (PMF), and interaction force between the enzyme and substrate. MD simulations revealed that R70 forms a salt bridge with D42 and the salt bridge affects the locations and the conformations of VD(3) in the bound state. SMD simulations revealed that the salt bridge tends to be formed strongly when VD(3) passes through the binding pocket. PMFs showed that the T70R mutation leads to energetic stabilization of enzyme-VD(3) binding in the region near the heme active site. Interestingly, these results concluded that the D42-R70 salt bridge at the entrance of the substrate access channel affects the region near the heme active site where the hydroxylation of VD(3) occurs; i.e., it is thought that the T70R mutation plays an important role in enhancing VD(3) hydroxylation activity. A significant future challenge is to compare the hydroxylation activities of R70 and T70 directly by a quantum chemical calculation, and three-dimensional coordinates of the enzyme and VD(3) obtained from MD and SMD simulations will be available for the future challenge.     

Activation and inhibition of cyclin-dependent kinase-2 by phosphorylation; a molecular dynamics study reveals the functional importance of the glycine-rich loop

Nanoseconds long molecular dynamics (MD) trajectories of differently active complexes of human cyclin-dependent kinase 2 (inactive CDK2/ATP, semiactive CDK2/Cyclin A/ATP, fully active pT160-CDK2/Cyclin A/ATP, inhibited pT14-; pY15-; and pT14,pY15,pT160-CDK2/Cyclin A/ATP) were compared. The MD simulations results of CDK2 inhibition by phosphorylation at T14 and/or Y15 sites provide insight into the structural aspects of CDK2 deactivation. The inhibitory sites are localized in the glycine-rich loop (G-loop) positioned opposite the activation T-loop. Phosphorylation of T14 and both inhibitory sites T14 and Y15 together causes ATP misalignment for phosphorylation and G-loop conformational change. This conformational change leads to the opening of the CDK2 substrate binding box. The phosphorylated Y15 residue negatively affects substrate binding or its correct alignment for ATP terminal phospho-group transfer to the CDK2 substrate. The MD simulations of the CDK2 activation process provide results in agreement with previous X-ray data.     

Long‐range molecular dynamics show that inactive forms of Protein Kinase A are more dynamic than active forms

Many protein kinases are characterized by at least two structural forms corresponding to the highest level of activity (active) and low or no activity, (inactive). Further, protein dynamics is an important consideration in understanding the molecular and mechanistic basis of enzyme function. In this work, we use protein kinase A (PKA) as the model system and perform microsecond range molecular dynamics (MD) simulations on six variants which differ from one another in terms of active and inactive form, with or without bound ligands, C‐terminal tail and phosphorylation at the activation loop. We find that the root mean square fluctuations in the MD simulations are generally higher for the inactive forms than the active forms. This difference is statistically significant. The higher dynamics of inactive states has significant contributions from ATP binding loop, catalytic loop, and αG helix. Simulations with and without C‐terminal tail show this differential dynamics as well, with lower dynamics both in the active and inactive forms if C‐terminal tail is present. Similarly, the dynamics associated with the inactive form is higher irrespective of the phosphorylation status of Thr 197. A relatively stable stature of active kinases may be better suited for binding of substrates and detachment of the product. Also, phosphoryl group transfer from ATP to the phosphosite on the substrate requires precise transient coordination of chemical entities from three different molecules, which may be facilitated by the higher stability of the active state.     

Structure and mechanism of yeast RNA triphosphatase: an essential component of the mRNA capping apparatus

RNA triphosphatase is an essential mRNA processing enzyme that catalyzes the first step in cap formation. The 2.05 A crystal structure of yeast RNA triphosphatase Cet1p reveals a novel active site fold whereby an eight-stranded beta barrel forms a topologically closed triphosphate tunnel. Interactions of a sulfate in the center of the tunnel with a divalent cation and basic amino acids projecting into the tunnel suggest a catalytic mechanism that is supported by mutational data. Discrete surface domains mediate Cet1p homodimerization and Cet1p binding to the guanylyltransferase component of the capping apparatus. The structure and mechanism of fungal RNA triphosphatases are completely different from those of mammalian mRNA capping enzymes. Hence, RNA triphosphatase presents an ideal target for structure-based antifungal drug discovery.