Catalytic mechanism of butane anaerobic oxidation for alkyl-coenzyme M reductase.

Methane is among the most potent of the greenhouse gases, which plays a key role in global climate change. As an excellent carbon and energy source, methane can be utilized by anaerobic methane oxidizing archaea and aerobic methane oxidizing bacteria. The previous work shows that an anaerobic thermophilic enrichment culture composed of dense consortia of archaea and bacteria apparently uses partly similar pathways to oxidize the C4 hydrocarbon butane. However, the catalytic mechanism of butane anaerobic oxidation for alkyl-coenzyme M reductase is still unknown. Therefore, molecular dynamics (MD) simulation was used to investigate the dynamics differences of catalytic mechanism between methane coenzyme M reductase (MCR) and alkyl-coenzyme M reductase (ACR). At first, the binding pocket of ACR is larger than that of MCR. Then, the complex of butane and ACR is more stable than that of methane and ACR. Protein conformation cloud suggests that the position of methane is dynamics and methane escapes from the binding pocket of ACR during most of the simulation time, while butane tightly binds in the pocket of ACR. The hydrophobic interactions between butane and ACR are more and stronger than those between methane and ACR. At the same time, the binding free energy between butane and ACR is significantly lower than that between methane and ACR. The dynamics correlation network indicates that the transformation of information flow for ACR-butane is smoother than that for ACR-methane. The shortest pathway for ACR-butane is from Gln144, Ala141, Hie135, Ile133, Ala160, Arg206, Asp97, Met94, Tyr347 to Phe345 with synergistic effect for two butane molecules. This study can insight into the catalytic mechanism for butane/ACR complex.

Insights from molecular dynamics simulations and steered molecular dynamics simulations to exploit new trends of the interaction between HIF-1α and p300.

Hypoxia-inducible factor-1 (HIF-1) is a transcription factor that plays an important role in the expression of genes, whose function is exerted through protein-protein interactions (PPIs), such as the transcriptional co-activator (CREB)-binding protein (CBP) and p300. Under hypoxic conditions, HIF-1is stabilized and translocated to CBP or p300, leading to the hypoxic response cascade. Furthermore, the PPI between HIF and p300/CBP is a potential cancer target for their role in the hypoxic response. In this study, molecular dynamics (MD) simulation was used to explore the conformational change for the p300 binding to one subunit of HIF-1, namely HIF-1α. Results indicated that HIF-1α-p300 complex was stable during MD simulation. New H-bonds were made in the intra-chain of p300 with HIF-1α binding. Inhibiting the HIF-1α-p300 interaction modulated the HIF-1α identification of selective molecules, which may indicate the target metabolic and cellular processes that enable the survival and growth of tumors in cancer chemotherapy. CAVER 3.0 results suggested that three main tunnels were present, according to helices 1, 2 and 3 of p300. To explore the unbinding pathway for HIF-1α via p300, we selected helices 1, 2 and 3 on the HIF-1α as a new ligand to explore the unbinding pathway via its own tunnel. For helix 1, R368 in p300 formed a H-bond with E816 in HIF1-α. A345 and D346 in p300 formed H-bonds with N803 in HIF-1α. A H-bond existed between K351(p300) and E789 (Hif1-α). These molecules may be the key residues in the unbinding pathway via its tunnel.Communicated by Ramaswamy H. Sarma.

Probing inhibition mechanisms of adenosine deaminase by using molecular dynamics simulations.

Adenosine deaminase (ADA) catalyzes the deamination of adenosine, which is important in purine metabolism. ADA is ubiquitous to almost all human tissues, and ADA abnormalities have been reported in various diseases, including rheumatoid arthritis. ADA can be divided into two conformations based on the inhibitor that it binds to: open and closed forms. Here, we chose three ligands, namely, FR117016 (FR0), FR221647 (FR2) (open form), and HDPR (PRH, closed form), to investigate the inhibition mechanism of ADA and its effect on ADA through molecular dynamics simulations. In open forms, Egap and electrostatic potential (ESP) indicated that electron transfer might occur more easily in FR0 than in FR2. Binding free energy and hydrogen bond occupation revealed that the ADA-FR0 complex had a more stable structure than ADA-FR2. The probability of residues Pro159 to Lys171 of ADA-FR0 and ADA-FR2 to form a helix moderately increased compared with that in nonligated ADA. In comparison with FR0 and FR2 PRH could maintain ADA in a closed form to inhibit the function of ADA. The α7 helix (residues Thr57 to Ala73) of ADA in the closed form was mostly unfastened because of the effect of PRH. The number of H bonds and the relative superiority of the binding free energy indicated that the binding strength of PRH to ADA was significantly lower than that of an open inhibitor, thereby supporting the comparison of the inhibitory activities of the three ligands. Alanine scanning results showed that His17, Gly184, Asp295, and Asp296 exerted the greatest effects on protein energy, suggesting that they played crucial roles in binding to inhibitors. This study served as a theoretical basis for the development of new ADA inhibitors.

Theoretical Study on Zearalenol Compounds Binding with Wild Type Zearalenone Hydrolase and V153H Mutant.

Zearalenone hydrolase (ZHD) is the only reported α/ß-hydrolase that can detoxify zearalenone (ZEN). ZHD has demonstrated its potential as a treatment for ZEN contamination that will not result in damage to cereal crops. Recent researches have shown that the V153H mutant ZHD increased the specific activity against α-ZOL, but decreased its specific activity to ß-ZOL. To understand whyV153H mutation showed catalytic specificity for α-ZOL, four molecular dynamics simulations combining with protein network analysis for wild type ZHD α-ZOL, ZHD ß-ZOL, V153H α-ZOL, and V153H ß-ZOL complexes were performed using Gromacs software. Our theoretical results indicated that the V153H mutant could cause a conformational switch at the cap domain (residues Gly161⁻Thr190) to affect the relative position catalytic residue (H242). Protein network analysis illustrated that the V153H mutation enhanced the communication with the whole protein and residues with high betweenness in the four complexes, which were primarily assembled in the cap domain and residues Met241 to Tyr245 regions. In addition, the existence of α-ZOL binding to V153H mutation enlarged the distance from the OAE atom in α-ZOL to the NE2 atom in His242, which prompted the side chain of H242 to the position with catalytic activity, thereby increasing the activity of V153H on the α-ZOL. Furthermore, α-ZOL could easily form a right attack angle and attack distance in the ZHD and α-ZOL complex to guarantee catalytic reaction. The alanine scanning results indicated that modifications of the residues in the cap domain produced significant changes in the binding affinity for α-ZOL and ß-ZOL. Our results may provide useful theoretical evidence for the mechanism underlying the catalytic specificity of ZHD.

Exploration of binding and inhibition mechanism of a small molecule inhibitor of influenza virus H1N1 hemagglutinin by molecular dynamics simulation.

Influenza viruses are a major public health threat worldwide. The influenza hemagglutinin (HA) plays an essential role in the virus life cycle. Due to the high conservation of the HA stem region, it has become an especially attractive target for inhibitors for therapeutics. In this study, molecular simulation was applied to study the mechanism of a small molecule inhibitor (MBX2329) of influenza HA. Behaviors of the small molecule under neutral and acidic conditions were investigated, and an interesting dynamic binding mechanism was found. The results suggested that the binding of the inhibitor with HA under neutral conditions facilitates only its intake, while it interacts with HA under acidic conditions using a different mechanism at a new binding site. After a series of experiments, we believe that binding of the inhibitor can prevent the release of HA1 from HA2, further maintaining the rigidity of the HA2 loop and stabilizing the distance between the long helix and short helices. The investigated residues in the new binding site show high conservation, implying that the new binding pocket has the potential to be an effective drug target. The results of this study will provide a theoretical basis for the mechanism of new influenza virus inhibitors.