A histidine acid phosphatase (HAP) (PhySc) with 99.50% protein sequence similarity with PHO5 from Saccharomyces cerevisiae was expressed functionally with the molecular mass of ∼110 kDa through co-expression along with the set of molecular chaperones dnaK, dnaJ, GroESL. The purified HAP illustrated the optimum activity of 28.75 ± 0.39 U/mg at pH 5.5 and 40 ˚C. The Km and Kcat values towards calcium phytate were 0.608 ± 0.09 mM and 650.89 ± 3.6 s- 1. The half-lives (T1/2) at 55 and 60 ˚C were 2.75 min and 55 s, respectively. The circular dichroism (CD) demonstrated that PhySc includes 30.5, 28.1, 21.3, and 20.1% of random coils, α-Helix, β-Turns, and β-Sheet, respectively. The Tm recorded by CD for PhySc was 56.5 ± 0.34˚C. The molecular docking illustrated that His59 and Asp322 act as catalytic residues in the PhySc. MD simulation showed that PhySc at 40 ˚C has higher structural stability over those of the temperatures 60 and 80 ˚C that support the thermodynamic in vitro investigations. Secondary structure content results obtained from MD simulation indicated that PhySc consists of 34.03, 33.09, 17.5, 12.31, and 3.05% of coil, helix, turn, sheet, and helix310, respectively, which is almost consistent with the experimental results.
In this work, molecular docking, pharmacophore modeling and molecular dynamics (MD) simulation were rendered for the mouse P-glycoprotein (P-gp) (code: 4Q9H) and bioflavonoids; amorphigenin, chrysin, epigallocatechin, formononetin and rotenone including a positive control; verapamil to identify protein-ligand interaction features including binding affinities, interaction characteristics, hot-spot amino acid residues and complex stabilities. These flavonoids occupied the same binding site with high binding affinities and shared the same key residues for their binding interactions and the binding region of the flavonoids was revealed that overlapped the ATP binding region with hydrophobic and hydrophilic interactions suggesting a competitive inhibition mechanism of the compounds. Root mean square deviations (RMSDs) analysis of MD trajectories of the protein-ligand complexes and NBD2 residues, and ligands pointed out these residues were stable throughout the duration of MD simulations. Thus, the applied preliminary structure-based molecular modeling approach of interactions between NBD2 and flavonoids may be gainful to realize the intimate inhibition mechanism of P-gp at NBD2 level and on the basis of the obtained data, it can be concluded that these bioflavonoids have the potential to cause herb-drug interactions or be used as lead molecules for the inhibition of P-gp (as anti-multidrug resistance agents) via the NBD2 blocking mechanism in future.
Lipases are known for their versatility in addition to their ability to digest fat. They can be used for the formulation of detergents, as food ingredients and as biocatalysts in many industrial processes. Because conventional enzymes are frangible at high temperatures, the replacement of conventional chemical routes with biochemical processes that utilize thermostable lipases is vital in the industrial setting. Recent theoretical studies on enzymes have provided numerous fundamental insights into the structures, folding mechanisms and stabilities of these proteins. The studies corroborate the experimental results and provide additional information regarding the structures that were determined experimentally. In this paper, we review the computational studies that have described how temperature affects the structure and dynamics of thermoenzymes, including the thermoalkalophilic L1 lipase derived from Bacillus stearothermophilus. We will also discuss the potential of using pressure for the analysis of the stability of thermoenzymes because high pressure is also important for the processing and preservation of foods.
Mutations in human laforin lead to an autosomal neurodegenerative disorder Lafora disease. In N-terminal carbohydrate binding domain of laforin, two mutations W32G and K87A are reported as highly disease causing laforin mutants. Experimental studies reported that mutations are responsible for the abolishment of glycogen binding which is a critical function of laforin. Our current computational study focused on the role of conformational changes in human laforin structure due to existing single mutation W32G and prepared double mutation W32G/K87A related to loss of glycogen binding. We performed 10 ns molecular dynamics (MD) simulation studies in the Gromacs package for both mutations and analyzed the trajectories. From the results, the global properties like root mean square deviation, root mean square fluctuation, radius of gyration, solvent accessible surface area and hydrogen bonds showed structural changes in atomic level observed in W32G and W32G/K87A laforin mutants. The conformational change induced by mutants influenced the loss of the overall stability of the native laforin. Moreover, the change in overall motion of protein was analyzed by principal component analysis and results showed protein clusters expanded more than native and also change in direction in case of double mutant in conformational space. Overall, our report provides theoretical information on loss of structure-function relationship due to flexible nature of laforin mutants. In conclusion, comparative MD simulation studies support the experimental data on W32G and W32G/K87A related to the lafora disease mechanism on glycogen binding.
The limited sequence similarity of protein sequences with known structures has led to an indispensable need for computational technology to predict their structures. Structural bioinformatics (SB) has become integral in elucidating the sequence-structure-function relationship of a protein. This report focuses on the applications of SB within the context of protein engineering including its limitation and future challenges.
Molecular dynamics simulation was used to study the dynamic differences between native Aspergillus niger PhyA phytase and a mutant with 20 % greater thermostability. Atomic root mean square deviation, radius of gyration, and number of hydrogen bonds and salt bridges are examined to determine thermostability factors. The results suggest that, among secondary structure elements, loops have the most impact on the thermal stability of A. niger phytase. In addition, the location rather than the number of hydrogen bonds is found to have an important contribution to thermostability. The results also show that salt bridges may have stabilizing or destabilizing effect on the enzyme and influence its thermostability accordingly.
Glycolipids form materials of considerable potential for a wide range of surfactant and thin film applications. Understanding the effect of glycolipid covalent structure on the properties of their thermotropic and lyotropic assemblies is a key step toward rational design of new glycolipid-based materials. Here, we perform molecular dynamics simulations of anhydrous bilayers of dodecyl β-maltoside, dodecyl β-cellobioside, dodecyl β-isomaltoside, and a C(12)C(10) branched β-maltoside. Specifically, we examine the consequences of chain branching and headgroup identity on the structure and dynamics of the lamellar assemblies. Chain branching of the glycolipid leads to measurable differences in the dimensions and interactions of the lamellar assembly, as well as a more fluid-like hydrophobic chain region. Substitution of the maltosyl headgroup of βMal-C(12) by an isomaltosyl moiety leads to a significant decrease in bilayer spacing as well as a markedly altered pattern of inter-headgroup hydrogen bonding. The distinctive simulated structures of the two regioisomers provide insight into the difference of ~90 °C in their observed clearing temperatures. For all four simulated glycolipid systems, with the exception of the sn-2 chain of the branched maltoside, the alkyl chains are ordered and exhibit a distinct tilt, consistent with recent crystallographic analysis of a branched chain Guerbet glycoside. These insights into structure-property relationships from simulation provide an important molecular basis for future design of synthetic glycolipid materials.
In this study, the proton dynamics of hydrogen bonds for two forms of crystalline aspirin was investigated by the Born-Oppenheimer molecular dynamics (BOMD) method. Analysis of the geometrical parameters of hydrogen bonds using BOMD reveals significant differences in hydrogen bonding between the two crystalline forms of aspirin, Form I and Form II. Analysis of the trajectory for Form I shows spontaneous proton transfer in cyclic dimers, which is absent in Form II. Quantization of the O-H stretching modes allows a detailed discussion on the strength of hydrogen-bonding interactions. The focal point of our study is examination of the hydrogen bond characteristics in the crystal structure and clarification of the influence of hydrogen bonding on the presence of the two crystalline forms of aspirin. In the BOMD method, thermal motions were taken into account. Solving the Schrödinger equation for the snapshots of 2D proton potentials, extracted from MD, gives the best agreement with IR spectra. The character of medium-strong hydrogen bonds in Form I of aspirin was compared with that of weaker hydrogen bonds in aspirin Form II. Two proton minima are present in the potential function for the hydrogen bonds in Form I. The band contours, calculated by using one- and two-dimensional O-H quantization, reflect the differences in the hydrogen bond strengths between the two crystalline forms of aspirin, as well as the strong hydrogen bonding in the cyclic dimers of Form I and the medium-strong hydrogen bonding in Form II.
Herein, we detail an atomic-level investigation of the cutinase enzyme encapsulated within a model metal-organic framework (MOF) platform using quantum mechanics calculations and molecular dynamics simulations. Cutinase, when encapsulated in an isoreticularly expanded MOF-74 (cutinase@IRMOF-74-VI), was proven to maintain its structural stability at temperatures that would otherwise denature the enzyme in its unprotected native state. Hydrogen bonding and salt bridge interactions, most notably involving arginine residues at the surface of the enzyme, were critical for stabilizing cutinase within the pore channels of IRMOF-74-VI. The findings reported support the viability of enzyme encapsulation in a porous material by demonstrating that a model enzyme not only retains its structural integrity but also remains accessible and active under extreme and foreign conditions.
The stability of enzymes is critical for their application in industrial processes, which generally require different conditions from the natural enzyme environment. Both rational and random protein engineering approaches have been used to increase stability, with the latter requiring extensive experimental effort for the screening of variants. Moreover, some general rules addressing the molecular origin of protein thermostability have been established. Herein, we demonstrate the use of molecular dynamics simulations to gain molecular level understanding of protein thermostability and to engineer stabilizing mutations. Carbonic anhydrase (CA) is an enzyme with a high potential for biotechnological carbon capture applications, provided it can be engineered to withstand the high temperature process environments, inevitable in most gas treatment units. In this study, we used molecular dynamics simulations at 343, 353, and 363 K to study the relationship between structure flexibility and thermostability in bacterial α-CAs and applied this knowledge to the design of mutants with increased stability. The most thermostable α-CA known, TaCA from Thermovibrio ammonificans, had the most rigid structure during molecular dynamics simulations, but also showed regions with high flexibility. The most flexible amino acids in these regions were identified from root mean square fluctuation (RMSF) studies, and stabilizing point mutations were predicted based on their capacity to improve the calculated free energy of unfolding. Disulfide bonds were also designed at sites with suitable geometries and selected based on their location at flexible sites, assessed by B-factor calculation. Molecular dynamics simulations allowed the identification of five mutants with lower RMSF of the overall structure at 400 K, compared to wild-type TaCA. Comparison of free-energy landscapes between wild-type TaCA and the most promising mutants, Pro165Cys-Gln170Cys and Asn140Gly, showed an increased conformational stability of the mutants at 400 K.
Glycolipid, found commonly in membranes, is also a liquid crystal material which can self-assemble without the presence of a solvent. Here, the dielectric and conductivity properties of three synthetic glycolipid thin films in different thermotropic liquid crystal phases were investigated over a frequency and temperature range of (10(-2)-10(6) Hz) and (303-463 K), respectively. The observed relaxation processes distinguish between the different phases (smectic A, columnar/hexagonal, and bicontinuous cubic Q) and the glycolipid molecular structures. Large dielectric responses were observed in the columnar and bicontinuous cubic phases of the longer branched alkyl chain glycolipids. Glycolipids with the shortest branched alkyl chain experience the most restricted self-assembly dynamic process over the broad temperature range studied compared to the longer ones. A high frequency dielectric absorption (Process I) was observed in all samples. This is related to the dynamics of the hydrogen bond network from the sugar group. An additional low-frequency mechanism (Process II) with a large dielectric strength was observed due to the internal dynamics of the self-assembly organization. Phase sensitive domain heterogeneity in the bicontinuous cubic phase was related to the diffusion of charge carriers. The microscopic features of charge hopping were modelled using the random walk scheme, and two charge carrier hopping lengths were estimated for two glycolipid systems. For Process I, the hopping length is comparable to the hydrogen bond and is related to the dynamics of the hydrogen bond network. Additionally, that for Process II is comparable to the bilayer spacing, hence confirming that this low-frequency mechanism is associated with the internal dynamics within the phase.
Interaction of an anticancer drug, vandetanib (VDB) with a ligand transporter, lysozyme (LYZ) was explored using multispectroscopic techniques, such as fluorescence, absorption and circular dichroism along with computational analysis. Fluorescence data and absorption results confirmed VDB-LYZ complexation. VDB-induced quenching was characterized as static quenching based on inverse correlation of KSV with temperature as well as kq values. The complex was characterized by the weak binding constant (Ka=4.96-3.14×103M-1). Thermodynamic data (ΔS=+12.82Jmol-1K-1; ΔH=-16.73kJmol-1) of VDB-LYZ interaction revealed participation of hydrophobic and van der Waals forces along with hydrogen bonds in VDB-LYZ complexation. Microenvironmental perturbations around tryptophan and tyrosine residues as well as secondary and tertiary structural alterations in LYZ upon addition of VDB were evident from the 3-D fluorescence, far- and near-UV CD spectral analyses, respectively. Interestingly, addition of VDB to LYZ significantly increased protein's thermostability. Molecular docking results suggested the location of VDB binding site near the LYZ active site while molecular dynamics simulation results suggested stability of VDB-LYZ complex. Presence of Mg2+, Ba2+ and Zn2+ was found to interfere with VDB-LYZ interaction.
The ability of human serum albumin (HSA) to bind medium-sized hydrophobic molecules is important for the distribution, metabolism, and efficacy of many drugs. Herein, the interaction between pyrene, a hydrophobic fluorescent probe, and HSA was thoroughly investigated using steady-state and time-resolved fluorescence techniques, ligand docking, and molecular dynamics (MD) simulations. A slight quenching of the fluorescence signal from Trp214 (the sole tryptophan residue in the protein) in the presence of pyrene was used to determine the ligand binding site in the protein, using Förster's resonance energy transfer (FRET) theory. The estimated FRET apparent distance between pyrene and Trp214 was 27Å, which was closely reproduced by the docking analysis (29Å) and MD simulation (32Å). The highest affinity site for pyrene was found to be in subdomain IB from the docking results. The calculated equilibrium structure of the complex using MD simulation shows that the ligand is largely stabilized by hydrophobic interaction with Phe165, Phe127, and the nonpolar moieties of Tyr138 and Tyr161. The fluorescence vibronic peak ratio I1/I3 of bound pyrene inside HSA indicates the presence of polar effect in the local environment of pyrene which is less than that of free pyrene in buffer. This was clarified by the MD simulation results in which an average of 5.7 water molecules were found within 0.5nm of pyrene in the binding site. Comparing the fluorescence signals and lifetimes of pyrene inside HSA to that free in buffer, the high tendency of pyrene to form dimer was almost completely suppressed inside HSA, indicating a high selectivity of the binding pocket toward pyrene monomer. The current results emphasize the ability of HSA, as a major carrier of several drugs and ligands in blood, to bind hydrophobic molecules in cavities other than subdomain IIA which is known to bind most hydrophobic drugs. This ability stems from the nature of the amino acids forming the binding sites of the protein that can easily adapt their shape to accommodate a variety of molecular structures.
Topological defect nucleation and boundary branching in crystal growth on a curved surface are two typical elastic instabilities driven by curvature induced stress, and have usually been discussed separately in the past. In this work they are simultaneously considered during crystal growth on a sphere. Phase diagrams with respect to sphere radius, size, edge energy and stiffness of the crystal for the equilibrium crystal morphologies are achieved by theoretical analysis and validated by Brownian dynamics simulations. The simulation results further demonstrate the detail of morphological evolution governed by these two different stress relaxation modes. Topological defect nucleation and boundary branching not only compete with each other but also coexist in a range of combinations of factors. Clarification of the interaction mechanism provides a better understanding of various curved crystal morphologies for their potential applications.
The α-amylases from Anoxybacillus species (ASKA and ADTA), Bacillus aquimaris (BaqA) and Geobacillus thermoleovorans (GTA, Pizzo and GtamyII) were proposed as a novel group of the α-amylase family GH13. An ASKA yielding a high percentage of maltose upon its reaction on starch was chosen as a model to study the residues responsible for the biochemical properties. Four residues from conserved sequence regions (CSRs) were thus selected, and the mutants F113V (CSR-I), Y187F and L189I (CSR-II) and A161D (CSR-V) were characterised. Few changes in the optimum reaction temperature and pH were observed for all mutants. Whereas the Y187F (t1/2 43 h) and L189I (t1/2 36 h) mutants had a lower thermostability at 65°C than the native ASKA (t1/2 48 h), the mutants F113V and A161D exhibited an improved t1/2 of 51 h and 53 h, respectively. Among the mutants, only the A161D had a specific activity, k(cat) and k(cat)/K(m) higher (1.23-, 1.17- and 2.88-times, respectively) than the values determined for the ASKA. The replacement of the Ala-161 in the CSR-V with an aspartic acid also caused a significant reduction in the ratio of maltose formed. This finding suggests the Ala-161 may contribute to the high maltose production of the ASKA.
Immobilisation of aptameric ligands on solid stationary supports for effective binding of target molecules requires understanding of the relationship between aptamer-polymer interactions and the conditions governing the mass transfer of the binding process. Herein, key process parameters affecting the molecular anchoring of a thrombin-binding aptamer (TBA) onto polymethacrylate monolith pore surface, and the binding characteristics of the resulting macroporous aptasensor were investigated. Molecular dynamics (MD) simulations of the TBA-thrombin binding indicated enhanced Guanine 4 (G4) structural stability of TBA upon interaction with thrombin in an ionic environment. Fourier-transform infrared spectroscopy and thermogravimetric analyses were used to characterise the available functional groups and thermo-molecular stability of the immobilised polymer generated with Schiff-base activation and immobilisation scheme. The initial degradation temperature of the polymethacrylate stationary support increased with each step of the Schiff-base process: poly(Ethylene glycol Dimethacrylate-co-Glycidyl methacrylate) or poly(EDMA-co-GMA) [196.0 °C (±1.8)]; poly(EDMA-co-GMA)-Ethylenediamine [235.9 °C (±6.1)]; poly(EDMA-co-GMA)-Ethylenediamine-Glutaraldehyde [255.4 °C (±2.7)]; and aptamer-modified monolith [273.7 °C (±2.5)]. These initial temperature increments reflected in the associated endothermic energies were determined with differential scanning calorimetry. The aptameric ligand density obtained after immobilisation was 480 pmol/μL. Increase in pH and ionic concentration affected the surface charge distribution and the binding characteristics of the aptamer-modified disk-monoliths, resulting in the optimum binding pH and ionic concentration of 8.0 and 5 mM Mg2+, respectively. These results are critical in understanding and setting parametric constraints indispensable to develop and enhance the performance of aptasensors.
Lipid membranes serve as effective barriers allowing cells to maintain internal composition differing from that of extracellular medium. Membrane permeation, both natural and artificial, can take place via appearance of transversal pores. The rearrangements of lipids leading to pore formation in the intact membrane are not yet understood in details. We applied continuum elasticity theory to obtain continuous trajectory of pore formation and closure, and analyzed molecular dynamics trajectories of pre-formed pore reseal. We hypothesized that a transversal pore is preceded by a hydrophobic defect: intermediate structure spanning through the membrane, the side walls of which are partially aligned by lipid tails. This prediction was confirmed by our molecular dynamics simulations. Conversion of the hydrophobic defect into the hydrophilic pore required surmounting some energy barrier. A metastable state was found for the hydrophilic pore at the radius of a few nanometers. The dependence of the energy on radius was approximately quadratic for hydrophobic defect and small hydrophilic pore, while for large radii it depended on the radius linearly. The pore energy related to its perimeter, line tension, thus depends of the pore radius. Calculated values of the line tension for large pores were in quantitative agreement with available experimental data.
Lipid membranes are extremely stable envelopes allowing cells to survive in various environments and to maintain desired internal composition. Membrane permeation through formation of transversal pores requires substantial external stress. Practically, pores are usually formed by application of lateral tension or transmembrane voltage. Using the same approach as was used for obtaining continuous trajectory of pore formation in the stress-less membrane in the previous article, we now consider the process of pore formation under the external stress. The waiting time to pore formation proved a non-monotonous function of the lateral tension, dropping from infinity at zero tension to a minimum at the tension of several millinewtons per meter. Transmembrane voltage, on the contrary, caused the waiting time to decrease monotonously. Analysis of pore formation trajectories for several lipid species with different spontaneous curvatures and elastic moduli under various external conditions provided instrumental insights into the mechanisms underlying some experimentally observed phenomena.
Wax appearance temperature (WAT), defined as the temperature at which the first solid paraffin crystal appears in a crude oil, is one of the key flow assurance indicators in the oil industry. Although there are several commonly-used experimental techniques to determine WAT, none provides unambiguous molecular-level information to characterize the phase transition between the homogeneous fluid and the underlying solid phase. Molecular Dynamics (MD) simulations employing the statistical associating fluid theory (SAFT) force field are used to interrogate the incipient solidification states of models for long-chain alkanes cooled from a melt to an arrested state. We monitor the phase change of pure long chain n-alkanes: tetracosane (C24H50) and triacontane (C30H62), and an 8-component surrogate n-alkane mixture (C12-C33) built upon the compositional information of a waxy crude. Comparison to Diffusion Ordered Spectroscopy Nuclear Magnetic Resonance (DOSY NMR) results allows the assessment of the limitations of the coarse-grained models proposed. We show that upon approach to freezing, the heavier components restrict their motion first while the lighter ones retain their mobility and help fluidize the mixture. We further demonstrate that upon sub-cooling of long n-alkane fluids and mixtures, a discontinuity arises in the slope of the self-diffusion coefficient with decreasing temperature, which can be employed as a marker for the appearance of an arrested state commensurate with conventional WAT measurements.
Fibrinogen is one of the key proteins that participate in the protein corona composition of many types of nanoparticles (NPs), and its conformational changes are crucial for activation of immune systems. Recently, we demonstrated that the fibrinogen highly contributed in the protein corona composition at the surface of zeolite nanoparticles. Therefore, understanding the interaction of fibrinogen with zeolite nanoparticles in more details could shed light of their safe applications in medicine. Thus, we probed the molecular interactions between fibrinogen and zeolite nanoparticles using both experimental and simulation approaches. The results indicated that fibrinogen has a strong and thermodynamically favorable interaction with zeolite nanoparticles in a non-cooperative manner. Additionally, fibrinogen experienced a substantial conformational change in the presence of zeolite nanoparticles through a concentration-dependent manner. Simulation results showed that both E- and D-domain of fibrinogen are bound to the EMT zeolite NPs via strong electrostatic interactions, and undergo structural changes leading to exposing normally buried sequences. D-domain has more contribution in this interaction and the C-terminus of γ chain (γ377-394), located in D-domain, showed the highest level of exposure compared to other sequences/residues.