MD investigation on the differences in the dynamic interactions between the specific ligand azamulin and two CYP3A isoforms, 3A4 and 3A5.

The unmarked potential drug molecule azamulin has been found to be a specific inhibitor of CYP3A4 and CYP3A5 in recent years, but this molecule also shows different binding ability and affinity to the two CYP3A isoforms. In order to explore the microscopic mechanism, conventional molecular dynamics (MD) simulation methods were performed to study the dynamic interactions between two isoforms and azamulin. The simulation results show that the binding of the ligand leads to different structural properties of two CYP3A proteins. First of all, compared with apo-CYP3A4, the binding of the ligand azamulin can lead to changes in the structural flexibility of CYP3A4, i.e., holo-CYP3A4 is more flexible than apo-CYP3A4. The structural changes of CYP3A5 are just the opposite. The ligand binding increases the rigidity of CYP3A5. Furthermore, the representative structures of the production phase in the MD simulation were in details analyzed to obtain the microscopic interactions between the ligand azamulin and two CYP3A isoforms at the atomic level. It is speculated that the difference of composition and interaction of the active sites is the fundamental cause of the change of structural properties of the two proteins.Communicated by Ramaswamy H. Sarma.

Molecular dynamics investigation on the interaction of human angiotensin-converting enzyme with tetrapeptide inhibitors.

Angiotensin-converting enzyme (ACE) is a well-known zinc metalloenzyme whose physiological functions are vital to blood pressure regulation and management of hypertension. The development of more efficient peptide inhibitors is of great significance for the prevention and treatment of hypertension. In this research, molecular dynamics (MD) simulations were implemented to study the specific binding mechanism and interaction between human ACE (hACE) and tetrapeptides, YIHP, YKHP, YLVR, and YRHP. The calculation of relative binding free energy on the one hand verified that YLVR, an experimentally identified inhibitor, has a stronger inhibitory effect and, on the other hand, indicated that YRHP is the "best" inhibitor with the strongest binding affinity. Inspection of atomic interactions discriminated the specific binding mode of each tetrapeptide inhibitor with hACE and explained the difference of their affinity. Moreover, in-depth analysis of the MD production trajectories, including clustering, principal component analysis, and dynamic network analysis, determined the dynamic correlation between tetrapeptides and hACE and obtained the communities' distribution of a protein-ligand complex. The present study provides essential insights into the binding mode and interaction mechanism of the hACE-peptide complex, which paves a path for designing effective anti-hypertensive peptides.

Interaction Mechanism of the Germination Stimulants Karrikins and Their Receptor ShKAI2iB.

The significance of karrikins (KARs) in plant physiology opens a door for their application in the agricultural production. As the first event of the whole signaling pathway, the binding of smoke-derived signal molecules KARs to the receptor protein KAI2 triggers the germination of the primary dormant seeds of all angiosperms, not just the "fire-prone" taxa. In the present study, all-atom molecular dynamics simulations, along with the accurate estimation of the ligand-receptor binding free energy, were used to investigate the atomic level interaction of all the members of the KARs family (from KAR1 to KAR6) with the receptor ShKAI2iB, an intermediate-evolving KAI2 from Striga hermonthica. The calculated binding energy value of KAR1 to ShKAI2iB, -5.64 kcal/mol, is in good agreement with the available experimental data, -5.67 kcal/mol. The further analysis of the detailed interaction between each KAR and the protein reveals the primary reasons for the difference of the affinity of the diverse ligands with the receptor and displays the regional characteristics of the protein's active site. Our research will not only provide clues for the study of equivalent endogenous phytohormone, but also contribute to the development of synthetic germinating chemicals.

Biotreatment of restaurant wastewater with an oily high concentration by newly isolated bacteria from oily sludge.

High concentration restaurant oily wastewater from restaurants and food processing industries discharged into water environment usually results in environment pollution and inhibits the activity of microorganisms in biological wastewater treatment systems. In this study, 75 strains from oily sludge were isolated with oil degradation activity for edible oil-contained wastewater. Eight isolates were able to grow well in liquid cultures with edible oil as the sole carbon source and discovered with high efficient oil-degrading ability. Seven out of eight isolates were identified as Acinetobacter and one isolate as Kluyvera cryocrescens, based on their 16S rRNA gene sequences. Three highly efficient oil degrading bacteria (Acinetobacter dijkshoorniae LYC46-2, Kluyvera cryocrescens LYC50-1a and Acinetobacter pittii LYC73-4b) were selected and their degradation characteristic were examined, the results showed that the three isolates were effective under pH range from 7.0 to 10.0, and temperature from 25 to 35 °C. For degradation of 2-4% (v/v) of vegetable oil, > 85% degradation percentage were obtained within 30 h. Degradation of the higher concentration oil (6-8%, v/v) result in 50-70% degradation percentage within 72 h, and the degradation percentage for the isolated strains were decreased about 50% for the degradation of 10% oil (< 45%) compared to 2% oil. Different type of oils were also tested, > 90% of degradation percentage were obtained by the three isolates, implied that these strains are capable of removing various oils efficiently. These results suggested that Acinetobacter dijkshoorniae LYC46-2, Kluyvera cryocrescens LYC50-1a and Acinetobacter pittii LYC73-4b are potential species could be efficiently used for high concentration restaurant oily wastewater treatment and might be applicable to a wastewater treatment system for the removal of oil.

MD Simulation Investigation on the Binding Process of Smoke-Derived Germination Stimulants to Its Receptor.

Karrikins (KARs) are a class of smoke-derived seed germination stimulants with great significance in both agriculture and plant biology. By means of direct binding to the receptor protein KAI2, the compounds can initiate the KAR signal transduction pathway, hence triggering germination of the dormant seeds in the soil. In the research, several molecular dynamics (MD) simulation techniques were properly integrated to investigate the binding process of KAR1 to KAI2 and reveal the details of the whole binding event. The calculated binding free energy, -7.00 kcal/mol, is in good agreement with the experimental measurement, -6.83 kcal/mol. The obtained PMF profile indicates the existence of three intermediate states in the binding process. The analysis of the simulation trajectories demonstrates that, in the intermediate structures, KAR1 is stabilized by some hydrophobic residues (Phe26, Phe134, Leu142, Trp153, Phe157, Leu160, Phe194), along with several bridging water molecules, and meanwhile, the significant shifting occurs in the local conformation of the protein as the ligand's binding. A series of the residues (Gln141-Phe157) on the so-called "cap domain" are proposed to be responsible for capturing the ligand at the initial stage of the binding. Besides, the changes of the ligand's poses are also quantitatively characterized by the proper choice of the coordinate system. Our work will contribute to the more penetrating understanding of the ligand binding process and the receptor affinity difference between several members in the KAR family and help design new, more effective germination stimulants.

The molecular mechanisms of plant plasma membrane intrinsic proteins trafficking and stress response.

Plasma membrane intrinsic proteins (PIPs) are plant channel proteins located on the plasma membrane. PIPs transfer water, CO2 and small uncharged solutes through the plasma membrane. PIPs have high selectivity to substrates, suggestive of a central role in maintaining cellular water balance. The expression, activity and localization of PIPs are regulated at the transcriptional and post-translational levels, and also affected by environmental factors. Numerous studies indicate that the expression patterns and localizations of PIPs can change in response to abiotic stresses. In this review, we summarize the mechanisms of PIP trafficking, transcriptional and post-translational regulations, and abiotic stress responses. Moreover, we also discuss the current research trends and future directions on PIPs.

Effect of External Electric Field on Substrate Transport of a Secondary Active Transporter.

Substrate transport across a membrane accomplished by a secondary active transporter (SAT) is essential to the normal physiological function of living cells. In the present research, a series of all-atom molecular dynamics (MD) simulations under different electric field (EF) strengths was performed to investigate the effect of an external EF on the substrate transport of an SAT. The results show that EF both affects the interaction between substrate and related protein's residues by changing their conformations and tunes the timeline of the transport event, which collectively reduces the height of energy barrier for substrate transport and results in the appearance of two intermediate conformations under the existence of an external EF. Our work spotlights the crucial influence of external EFs on the substrate transport of SATs and could provide a more penetrating understanding of the substrate transport mechanism of SATs.

Exploring the interaction between human focal adhesion kinase and inhibitors: a molecular dynamic simulation and free energy calculations.

Focal adhesion kinase is an important target for the treatment of many kinds of cancers. Inhibitors of FAK are proposed to be the anticancer agents for multiple tumors. The interaction characteristic between FAK and its inhibitors is crucial to develop new inhibitors. In the present article, we used Molecular Dynamic (MD) simulation method to explore the characteristic of interaction between FAK and three inhibitors (PHM16, TAE226, and ligand3). The MD simulation results together with MM-GB/SA calculations show that the combinations are enthalpy-driven process. Cys502 and Asp564 are both essential residues due to the hydrogen bond interactions with inhibitors, which was in good agreement with experimental data. Glu500 can form a non-classical hydrogen bond with each inhibitor. Arg426 can form electrostatic interactions with PHM16 and ligand3, while weaker with TAE226. The electronic static potential was employed, and we found that the ortho-position methoxy of TAE226 has a weaker negative charge than the meta-position one in PHM16 or ligand3. Ile428, Val436, Ala452, Val484, Leu501, Glu505, Glu506, Leu553, Gly563 Leu567, Ser568 are all crucial residues in hydrophobic interactions. The key residues in this work will be available for further inhibitor design of FAK and also give assistance to further research of cancer.

Structural features and dynamic investigations of the membrane-bound cytochrome P450 17A1.

Cytochrome P450 (CYP) 17A1 is a dual-function monooxygenase with a critical role in the synthesis of many human steroid hormones. The enzyme is an important target for treatment of breast and prostate cancers that proliferate in response to estrogens and androgens. Despite the crystallographic structures available for CYP17A1, no membrane-bound structural features of this enzyme at atomic level are available. Accumulating evidence has indicated that the interactions between bounded CYPs and membrane could contribute to the recruitment of lipophilic substrates. To this end, we have investigated the effects on structural characteristics in the presence of the membrane for CYP17A1. The MD simulation results demonstrate a spontaneous insertion process of the enzyme to the lipid. Two predominant modes of CYP17A1 in the membrane are captured, characterized by the depths of insertion and orientations of the enzyme to the membrane surface. The measured heme tilt angles show good consistence with experimental data, thereby verifying the validity of the structural models. Moreover, conformational changes induced by the membrane might have impact on the accessibility of the active site to lipophilic substrates. The dynamics of internal aromatic gate formed by Trp220 and Phe224 are suggested to regulate tunnel opening motions. The knowledge of the membrane binding characteristics could guide future experimental and computational works on membrane-bound CYPs so that various investigations of CYPs in their natural, lipid environment rather than in artificially solubilized forms may be achieved.

Molecular basis of the recognition of arachidonic acid by cytochrome P450 2E1 along major access tunnel.

Cytochrome P450 2E1 is widely known for its ability to oxidize both low molecular weight xenobiotics and endogenous fatty acids (e.g., arachidonic acid (AA)). In this study, we investigated the structural features of the AA-bound CYP2E1 complex utilizing molecular dynamics (MD) and found that the distinct binding modes for both AA and fatty acid analog are conserved. Moreover, multiple random acceleration MD simulations and steered MD simulations uncovered the most possible tunnel for fatty acids. The main attractions are derived from three key residues, His107, Ala108, and His109, whose side chains reorient to keep ligands bound via hydrogen bonds during the initial unbinding process. More importantly, based on the calculated binding free energy results, we hypothesize that the hydrogen bonds between the receptor and the ligand are the most important contributors involved in the binding affinity. Thus, it is inferred that the hydrogen bonds between these three residues and the ligand may help offer insights into the structural basis of the different ligand egress mechanisms for fatty acids and small weight compounds. Our investigation provides detailed atomistic insights into the structural features of human CYP2E1-fatty acid complex structures. Furthermore, the ligand-binding characteristics obtained in the present study are helpful for both experimental and computational studies of CYPs and may allow future researchers to achieve desirable changes in enzymatic activities.

Heparin makes differences: a molecular dynamics simulation study on the human ßII-tryptase monomer.

Human ß-tryptase, an enzyme with trypsin-like activity in mast cells, is an important target for the treatment of inflammatory and allergy related diseases. Heparin has been inferred to play a vital role in the stabilization of the tryptase structure and the maintenance of its active form. Up to now, the structure-function relationship between heparin and the ßII-tryptase monomer has not been studied with atomic resolution due to the lack of a complex structure of tryptase and heparin. To this end, the exact effect of heparin bonding to the ßII-tryptase monomer structure has been investigated using molecular docking and molecular dynamics (MD) simulation. The MD simulation results combined with MM-GB/SA calculations showed that heparin stabilized the ß-tryptase structure mainly through salt bridge interaction. The averaged noncovalent interaction (aNCI) method was employed for the visualization of nonbonding interactions. A crucial loop, which is located in the core region of ßII-tryptase monomer structure, has been found. Arg188 and Asp189 from this loop act as a salt bridge intermediary between 4-mer heparin and 0GX. The observation of a salt bridge between Asp189 and P1 groups of 0GX confirms the supposed interaction between these two groups. These two residues have been proved to be responsible for the direction of the P1 group of 0GX. Our study revealed that how heparin affected the activity of the human ßII-tryptase monomer (hBTM) through salt bridge interactions. The knowledge of heparin binding characteristics and the key residue contributions in this study may enlighten further the inhibitor design of this enzyme and may also improve our understanding of inflammatory and allergy related diseases.

Molecular simulation investigation on the interaction between barrier-to-autointegration factor dimer or its Gly25Glu mutant and LEM domain of emerin.

The interaction between barrier-to-autointegration factor dimer (BAF2) and LEM domain of emerin (EmLEM) was studied by molecular simulation methods. Nonspecific fragment of double-strand DNA molecule was docked with each chain of BAF2 by ZDOCK program. The model of DNA2:BAF2:EmLEM was thus constructed. The mutant Gly25Glu of BAF2 was manually constructed to explore the detailed effect of the mutation on the binding of BAF2 and EmLEM. It has been experimentally suggested that point mutation Gly25Glu can disturb the binding between BAF2 and EmLEM. Then, molecular dynamics (MD) simulations were performed on DNA2:BAF2(WT):EmLEM and DNA2:BAF2(MT):EmLEM complexes. 30ns trajectories revealed that the trajectory fluctuations of MT complex are more violent than that of the WT complex. Further, the binding free energy analysis showed that the electronegative residues Asp57, Glu61 and Asp65 from chain A, glu36 from chain B of BAF2 mainly contribute to interact with EmLEM. Besides, a stable π-π stack between trp62 and phe39 from BAF2(WT) chain B is destroyed by Glu25 in BAF2(MT). As a result, trp62 forms an interaction with glu25, and phe39 converts to strengthen affinity to EmLEM. On the other hand, Trp62 from chain A also forms a strong interaction with MT Glu25. Thus, with the docking of DNA, BAF2(MT) has higher affinity with EmLEM than BAF2(WT).

Mechanism of a pH-induced peptide inserting into a POPC bilayer: a molecular dynamic study.

Membrane insertion peptides have been developed in recent years and serve as cargos to deliver small molecules into cells. A class of membrane insertion peptides is the so called pH-induced peptides (pHLIPs), which are able to insert into membrane when the environment pH is acidic. Despite a number of experimental studies, the insertion process as well as the penetration mechanism is still worth study with computational methods. Thus, we performed molecular dynamics simulations in this study to elucidate the detailed penetration process and mechanism. Both protonated and unprotonated peptides are employed to interact with a POPCs bilayer. By analyzing the trajectory of the simulation, the peptide travelling across membrane is expected to take milliseconds or seconds. While the peptide penetrating through the POPC bilayer boundary is much faster (several nanoseconds). More importantly, the elaborate energies between a peptide and water molecules, the energies between a peptide and POPCs have been analyzed throughout the simulation time correspondingly. A constant decrease of interaction energies have been observed for peptide-water interaction in the protonated condition. At last, we employ the statistics of hydrogen bonds to explain the penetration mechanism tentatively. For the protonated system, the decrease of hydrogen bonds of peptide-water and the increase of hydrogen bonds of peptide- POPCs have been considered as the main driven force for the peptide insertion.

Molecular dynamics simulation of the processive endocellulase Cel48F from Clostridium cellulolyticum: a novel "water-control mechanism" in enzymatic hydrolysis of cellulose.

Glycoside hydrolase of Cel48F from Clostridium cellulolyticum is an important processive cellulose, which can hydrolyze cellulose into cellobiose. Molecular dynamics simulations were used to investigate the hydrolysis mechanism of cellulose. The two conformations of the Cel48F-cellotetrose complex in which the cellotetroses are bound at different sites (known as the sliding conformation and the hydrolyzing conformation) were simulated. By comparing these two conformations, a water-control mechanism is proposed, in which the hydrolysis proceeds by providing a water molecule for every other glucosidic linkage. The roles of certain key residues are determined: Glu55 and Asp230 are the most probable candidates for acid and base, respectively, in the mechanism of inverting anomeric carbon. Met414 and Trp417 constitute the water-control system. Glu44 might keep the substrate at a certain location within the active site or help the substrate chain to move from the sliding conformation to the hydrolyzing conformation. The other hydrophobic residues around the substrate can decrease the sliding energy barrier or provide a hydrophobic environment to resist entry of the surrounding water molecules into the active site, except for those coming from a specific water channel.

Molecular simulation investigation on the interaction between barrier-to-autointegration factor or its Gly25Glu mutant and DNA.

In order to understand the binding mechanism between Barrier-to-autointegration factor (BAF) and DNA, two DNA:BAF complexes with wild type (WT) BAF and its Gly25Glu point mutate type (MT) were generated by molecular docking on the basis of the crystal structures of BAF (PDB code: 2ODG, chain A) and DNA (PDB code: 2BZF, chain B and C). Then, molecular dynamics (MD) simulations were performed on the two docked structures, as well as BAF (WT) and BAF (MT). The results show that monomer BAF is more flexible than BAF in DNA:BAF complex, suggesting that DNA is effective to stabilize conformation of BAF, which is in good agreement with the experimental results. Besides, the mutated Glu25 in DNA:BAF (MT) can change the BAF conformation to some extent. With deeper investigation on the DNA:BAF structures, the hydrogen bonds are found to make great contribution to the interaction between DNA and BAF. The hydrogen bonds in DNA:BAF (MT) are fewer than those in DNA:BAF (WT), indicating that the Gly25Glu mutation in BAF has an important effect on the hydrogen bonds in the DNA:BAF complex. Besides, the binding free energy in DNA:BAF (MT) is also higher than that in DNA:BAF (WT). It results from the influence of Glu25 side chain on the orientation of Lys6 and Lys33 in the interface between DNA and BAF. The binding free energy of Lys72, another key residue, decreases a lot in DNA:BAF (MT) anomalously. The decreasing energy causes the destruction of hydrophobic pocket in the binding site between DNA and BAF (MT). Our results are helpful for further experimental investigations.

Fosfomycin induced structural change in fosfomycin resistance kinases FomA: molecular dynamics and molecular docking studies.

Fosfomycin resistance kinases FomA, one of the key enzymes responsible for bacterial resistances to fosfomycin, has gained much attention recently due to the raising public concern for multi-drug resistant bacteria. Using molecular docking followed by molecular dynamics simulations, our group illustrated the process of fosfomycin induced conformational change of FomA. The detailed roles of the catalytic residues (Lys18, His58 and Thr210) during the formation of the enzyme-substrate complex were shown in our research. The organization functions of Gly53, Gly54, Ile61 and Leu75 were also highlighted. Furthermore, the cation-π interaction between Arg62 and Trp207 was observed and speculated to play an auxiliary role in the conformation change process of the enzyme. This detailed molecular level illustration of the formation of FomA·ATP·Mg·Fosfomycin complex could provide insight for both anti-biotic discovery and improvement of fosfomycin in the future.

Mutation and low pH effect on the stability as well as unfolding kinetics of transthyretin dimer.

Transthyretin (TTR) dissociation and aggregation appear to cause several amyloid diseases. TTR dimer is an important intermediate that is hard to be observed from the biological experiments. To date, the molecular origin and the structural motifs for TTR dimer dissociation, as well as the unfolding process have not been rationalized at atomic resolution. To this end, we have investigated the effect of low pH and mutation L55P on stability as well as the unfolding pathway of TTR dimer using constant pH molecular dynamics simulations. The result shows that acidic environment results in loose TTR dimer structure. Mutation L55P causes the disruption of strand D and makes the CE-loop very flexible. In acidic conditions, dimeric L55P mutant exhibits notable conformation changes and an evident trend to separate. Our work shows that the movements of strand C and the loops nearby are the beginning of the unfolding process. In addition, hydrogen bond network at the interface of the two monomers plays a part in stabilizing TTR dimer. The dynamic investigation on TTR dimer provides important insights into the structure-function relationships of TTR, and rationalizes the structural origin for the tendency of unfolding and changes of structure that occur upon introduction of mutation and pH along the TTR dimer dissociation and unfolding process.

Ultrafast self-healing of polymer toward strength restoration.

Self-healing materials should take effect immediately following crack generation in principle, but the speed of autonomic recovery of mechanical properties through either extrinsic or intrinsic healing strategy reported so far is not that fast. Mostly, a couple of hours are taken for reaching steady state or maximum healing. To accelerate the healing process, the authors of this work make use of antimony pentafluoride as instant hardener of epoxy and successfully encapsulate the highly active antimony pentafluoride-ethanol complex in terms of hollow silica spheres. Accordingly, self-healing agent based on microencapsulated antimony pentafluoride-ethanol complex and epoxy monomer is developed. Epoxy material with the embedded healant capsules can thus be healed within a few seconds, as demonstrated by impact and fatigue tests. It is believed that the outcome presented here might help to move the self-healing technique closer to practical application, especially when the engineering significance of epoxy material is concerned.

Exploring the mechanism how Marburg virus VP35 recognizes and binds dsRNA by molecular dynamics simulations and free energy calculations.

Filoviruses often cause terrible infectious disease which has not been successfully dealt with pharmacologically. All filoviruses encode a unique protein termed VP35 which can mask doubled-stranded RNA to deactivate interferon. The interface of VP35-dsRNA would be a feasible target for structure-based antiviral agent design. To explore the essence of VP35-dsRNA interaction, molecular dynamics simulation combined with MM-GBSA calculations were performed on Marburg virus VP35-dsRNA complex and several mutational complexes. The energetic analysis indicates that nonpolar interactions provide the main driving force for the binding process. Although the intermolecular electrostatic interactions play important roles in VP35-dsRNA interaction, the whole polar interactions are unfavorable for binding which result in a low binding affinity. Compared with wild type VP35, the studied mutants F228A, R271A, and K298A have obviously reduced binding free energies with dsRNA reflecting in the reduction of polar or nonpolar interactions. The results also indicate that the loss of binding affinity for one dsRNA strand would abolish the total binding affinity. Three important residues Arg271, Arg294, and Lys298 which makes the largest contribution for binding in VP35 lose their binding affinity significantly in mutants. The uncovering of VP35-dsRNA recognition mechanism will provide some insights for development of antiviral drug.

Molecular dynamics simulations studies and free energy analysis on inhibitors of MDM2-p53 interaction.

The oncoprotein MDM2 (murine double minute 2) negatively regulates the activity and stability of tumor suppressor p53. Inactivation of the MDM2-p53 interaction by potent inhibitors offers new possibilities for anticancer therapy. Molecular dynamics (MD) simulations were performed on three inhibitors-MDM2 complexes to investigate the stability and structural transitions. Simulations show that the backbone of MDM2 maintains stable during the whole time. However, slightly structural changes of inhibitors and MDM2 are observed. Furthermore, the molecular mechanics generalized Born surface area (MM-GBSA) approach was introduced to analyze the interactions between inhibitors and MDM2. The results show that binding of inhibitor pDIQ to MDM2 is significantly stronger than that of pMI and pDI to MDM2. The side chains of residues have more contribution than backbone of residues in energy decomposition. The structure-affinity analyses show that L54, I61, M62, Y67, Q72, H73 and V93 produce important interaction energy with inhibitors. The residue W/Y22' is also very important to the interaction between the inhibitors and MDM2. The three-dimensional structures at different times indicate that the mobility of Y100 influences on the binding of inhibitors to MDM2, and its change has important role in conformations of inhibitors and MDM2.