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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Constant pH molecular dynamics (CpHMD) and mutation studies: insights into AaegOBP1 pH-induced ligand releasing mechanism.

AaegOBP1, isolated from the male and female antenna of yellow fever mosquitoes, may serve as crucial molecular targets for the development of mosquitoes' attractants and for the control of mosquito populations. Nowadays crystal structures of AaegOBP1 in the neutral environment have been obtained, whereas little research is focused on the conformational change of AaegOBP1 in the acid solution. In our study, the conformational change and the ligand bound poses in different solution pH were investigated using constant pH molecular dynamics (CpHMD) as well as mutation studies. Results demonstrate that the protein changes dramatically in low pH solution and that the pH-sensing triad (Arg23-Tyr54-Ile125) plays an indispensable role in the structural stability and ligand binding. In addition, we have proved that the residue Arg23 is the most important one of the pH-sensing triad. This work could provide more penetrating understanding of the pH-induced ligand-releasing mechanism.

Structural and dynamic basis of human cytochrome P450 7B1: a survey of substrate selectivity and major active site access channels.

Cytochrome P450 (CYP) 7B1 is a steroid cytochrome P450 7α-hydroxylase that has been linked directly with bile salt synthesis and hereditary spastic paraplegia type 5 (SPG5). The enzyme provides the primary metabolic route for neurosteroids dehydroepiandrosterone (DHEA), cholesterol derivatives 25-hydroxycholesterol (25-HOChol), and other steroids such as 5α-androstane-3ß,17ß-diol (anediol), and 5α-androstene-3ß,17ß-diol (enediol). A series of investigations including homology modeling, molecular dynamics (MD), and automatic docking, combined with the results of previous experimental site-directed mutagenesis studies and access channels analysis, have identified the structural features relevant to the substrate selectivity of CYP7B1. The results clearly identify the dominant access channels and critical residues responsible for ligand binding. Both binding free energy analysis and total interaction energy analysis are consistent with the experimental conclusion that 25-HOChol is the best substrate. According to 20 ns MD simulations, the Phe cluster residues that lie above the active site, particularly Phe489, are proposed to merge the active site with the adjacent channel to the surface and accommodate substrate binding in a reasonable orientation. The investigation of CYP7B1-substrate binding modes provides detailed insights into the poorly understood structural features of human CYP7B1 at the atomic level, and will be valuable information for drug development and protein engineering.

Exploring the molecular basis of dsRNA recognition by Mss116p using molecular dynamics simulations and free-energy calculations.

DEAD-box proteins are the largest family of helicase that are important in nearly all aspects of RNA metabolism. However, it is unclear how these proteins recognize and bind RNA. Here, we present a detailed analysis of the related DEAD-box protein Mss116p-RNA interaction, using molecular dynamics simulations with MM-GBSA calculations. The energetic analysis indicates that the two strands of double strands RNA (dsRNA) are recognized asymmetrically by Mss116p. The strand 1 of dsRNA provides the main binding affinity. Meanwhile, the nonpolar interaction provides the main driving force for the binding process. Although the contribution of polar interaction is small, it is vital in stabilizing the protein-RNA interaction. Compared with the wild type Mss116p, two studied mutants Q412A and D441A have obviously reduced binding free energies with dsRNA because of the decreasing of polar interaction. Three important residues Lys409, Arg415 and Arg438 lose their binding affinity significantly in mutants. In conclusion, these results complement previous experiments to advance comprehensive understanding of Mss116p-dsRNA interaction. The results also would provide support for the application of similar approaches to the understanding of other DEAD-box protein-RNA complexes.

Molecular dynamic investigations of the mutational effects on structural characteristics and tunnel geometry in CYP17A1.

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 the treatment of breast and prostate cancers that proliferate in response to estrogens and androgens. Despite the ample experimental mutagenesis data, the molecular origin and the structural motifs for the enzymatic activities deficiencies have not been rationalized at the atomic resolution. To this end, we have investigated the effects on structural characteristics and tunnel geometry upon single point mutations in CYP17A1. The MD simulation results combined with PMF calculations and MM-GBSA calculations render an "access mechanism" which encapsulates the effects of mutations on the changes in both structural flexibility and tunnel dynamics, bridging the gap between the theory and the experimentally observed results of enzymatic activity decrease. The underlying molecular mechanism of the heterogeneities in open/closed conformational changes, as well as the wider opening of their respective major tunnels between wt17A1 and two mutants, may be attributed to the closer distances of hydrophobic residues or the disruption of a hydrophobic core. The knowledge of ligand binding characteristics and key residues contributions could guide future experimental and computational work on CYPs so that desirable changes in their enzymatic activities may be achieved. The present study provides important insights into the structure-function relationships of CYP17A1 protein, which could contribute to further understanding about 17-hydroxylase deficiencies and may also improve the understanding of polycystic ovary disease.

[Temperature compensation strategy and implementation for photoelectric modulation interferometer with large optical path difference].

For temperature drift in hypervelocity photoelectric modulation interferometer, a control model of temperature compensation is presented including voltage and phase compensation. First, according to the similar and modeling theory, an equivalent circuit model of mechanical properties of hypervelocity photoelectric modulation interferometer was established, the impact of temperature drift on its resonance frequency was analyzed, a mathematical model was set up, which contains drive voltage, frequency and resonance frequency, and the control method was determined for high optical path difference to get steady. Then, a digital method including voltage and phase compensation is given for optical path difference deviation control, which merges the DPLL and program of voltage and phase compensation. Finally, the control method was tested through experiment system. A test between drive control system including voltage and phase compensation and traditional drive control system was executed, using a laser doppler vibrometer to record the amount of change in optical path difference within 3 hours. Results show that the optical path difference deviation caused by temperature drift in long term is reduced by about 50%.

[Phase correction technology research and improvement based on single-sided interferograms in FTIR].

In the Fourier transform infrared spectrometers the interferograms are usually sampled by the means of the single-side beyond zero path difference and the interferograms are asymmetric. These phase-corrected technologies are researched and improved in the present paper. In the Mertz method which requires apodization to be higher, the interferograms are apodized by these asymmetric windows so that the width of the main lobe increases and spectral resolution reduced. In order to solve these problems, a method which combines the Mertz and Forman methods is presented in the paper. In the method a double-side interferogram is constructed by the mirrored processing of the single-side interferogram beyond zero path difference, followed by the symmetric apodization, fast Fourier transform (FFT) to rebuild spectrum. And in the method high resolution phase spectrum is computed by the symmetrization of double-side interfergram to improve the phase-correction accuracy. In the experiment, these asymmetric window functions have the lower resolution compared with these symmetric window functions. And the spectrum which is acquired by the improved phase corrected method is compared with the spectra which are acquired by the Mertz and Forman methods, and the improved method has less error than Mertz, less computation than Forman, and the spectral resolution achieved 2 cm(-1).

[Design of non-dispersed infrared (NDIR) methane gas sensor].

A non-dispersed infrared (NDIR) methane gas sensors system based on infrared absorption spectrum theory was designed according to single light beam and double wavelengths technology. In the system, an infrared LED IRL715 serving as the light power, a absorptive gas cell with the function of dust-proof and damp-proof and a pyroelectric detector LIM-262 are composed of optical probe. Signal condition uses active filter circuit and differential amplifier, and binomial expression fits the relation curve between methane concentration and voltage, which realizes accurate detection of gas concentration. Experiment approved that the sensor system with good consistency and applicability can detect the range of 5% methane reliably and have 0.5% of the sensitivity, possessing the conditions for industrial applications initially.

[Fiber monitoring network of methane concentration based on space division multiplexing in coal mine].

Methane explosion accidents occur frequently, and accurate and real-time detection and early warning of methane concentration are the effective means of preventing these accidents. The research was based on the spectrum absorption properties of methane, and a near-infrared tunable DFB laser diode of 1.65 microm wave band was used. With the mode-hopping features of laser diode, a differential absorption of double-wavelength and single-fiber optical sensor network was designed. Sixteen methane sensors were multiplexed in this system with space division multiple access technology and optical switch, and the key technologies of anti-dust in gas absorption cell were researched. All signals were gathered by the PCI data acquisition card, and information of each way was analyzed and displayed with virtual instrument. The results of experiment show that the method can reach the sensitivity of 0.05% even without using a phase-locked amplifier and the absorption light path is only 10 centimeters. Long-time accuracy and stability of all sensors could meet the practical demands, and the response time of each sensor was less than 1 seconds. With the replacement of lasers, the network can be used for the real-time detection of other gases.

[Research on concentration of multi-component pollution gas based on SVM with kernel optimized by rough set].

This paper introduced the application of support vector machines (SVM) regression method based on kernel function optimized by the rough set in the infrared spectrum quantitative calculation. According to kernel function with the rough set classification's method, the spectrum data (characteristic wavelength section) is optimized. The kernel function leads support vector machines, and the SVM project the two-dimensional room to the multi-dimensional room, and calculate the concentration of every kind of gas in multi-component pollution gas. By using two kinds of typical spectrum data processing algorithm to make the contrast, the comparison of five kinds of gaseous mixture various proximate analysis is carried out, and when the spectrum separable rate is high, the predicted values of the three methods approach the normal value, and the average error is smaller than 0.13; but when the spectrum separable rate is low, the RS-SVM predicted value is more precise than the first two kinds. Experimental data show that the consequence is better when there are more testing types, and the precision and operation of this method is of more remarkable superiority.

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.