Effects of Organic Solvent and Crystal Water on γ-Chymotrypsin in Acetonitrile Media: Observations from Molecular Dynamics Simulation and DFT Calculation

Solvent induced conformational changes for the altered activity of laccase: A molecular dynamics study

The growing demands of solvent-based industries like paint, pharmaceutical, petrochemical, paper and pulp, etc., have directly increased the release of effluents that are rich in hazardous aromatic compounds in the environment. A sustainable biotechnological approach utilizing laccases as biocatalyst enable in biodegradation of these aromatic toxin-rich effluents. However, this enzymatic process is ineffective as laccases lose their stability and catalytic activity at high organic solvent concentrations. In this study, molecular dynamic simulations of a novel solvent tolerant laccase, DLac from Cerrena sp. RSD1 was performed to explore the molecular-level understanding of DLac in 30%(v/v) acetone and acetonitrile. Solvent-induced conformational changes were analyzed via protein structure network, which was illustrated with respect to cliques and communities. In the presence of acetonitrile, the cliques around the active site and substrate-binding site were disjoined, thus the communities lost their network integrity. Whereas with acetone, the community near the substrate-binding site gained new residues and formed a rigidified network that corresponded to enhanced DLac's activity. Moreover, prominent solvent binding sites were speculated, which can be probable mutation targets to further improve solvent tolerance and catalytic activity. The molecular basis behind solvent induced catalytic activity will further aid in engineering laccase for its industrial application.

Molecular dynamics simulations of amino acid adsorption and transport at the acetonitrile–water–silica interface: the role of side chains

The solvation and transport of amino acid residues at liquid–solid interfaces have great importance for understanding the mechanism of separation of biomolecules in liquid chromatography. This study uses umbrella sampling molecular dynamics simulations to study the adsorption and transport of three amino acid molecules with different side chains (phenylalanine (Phe), leucine (Leu) and glutamine (Gln)) at the silica–water–acetonitrile interface in liquid chromatography. Free energy analysis shows that the Gln molecule has stronger binding affinity than the other two molecules, indicating the side chain polarity may play a primary role in adsorption at the liquid–solid interface. The Phe molecule with a phenyl side chain exhibits stronger adsorption free energy than Leu with a non-polar side chain, which can be ascribed to the better solvated configuration of Phe. Further analysis of molecular orientations found that the amino acid molecules with apolar side chains (Phe and Leu) have ‘standing up’ configurations at their stable adsorption state, where the polar functional groups are close to the interface and the side chain is far from the interface, whereas the amino acid molecule with a polar side chain (Gln) chooses the ‘lying’ configuration, and undergoes a sharp orientation transition when the molecule moves away from the silica surface. Extending our simulation studies to systems with different solute concentrations reveals that there is a decrease in the adsorption free energy as well as surface diffusion as the solute concentration increases, which is related to the crowding in the interfacial layers. This simulation study gives a detailed microscopic description of amino acid molecule solvation and transport at the acetonitrile–water–silica interface in liquid chromatography and will be helpful for understanding the retention mechanism for amino acid separation.

The solvation and transport of amino acid residues at liquid–solid interfaces have great importance for understanding the mechanism of separation of biomolecules in liquid chromatography.

In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function

Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.

Lid opening and conformational stability of T1 Lipase is mediated by increasing chain length polar solvents

The dynamics and conformational landscape of proteins in organic solvents are events of potential interest in nonaqueous process catalysis. Conformational changes, folding transitions, and stability often correspond to structural rearrangements that alter contacts between solvent molecules and amino acid residues. However, in nonaqueous enzymology, organic solvents limit stability and further application of proteins. In the present study, molecular dynamics (MD) of a thermostable Geobacillus zalihae T1 lipase was performed in different chain length polar organic solvents (methanol, ethanol, propanol, butanol, and pentanol) and water mixture systems to a concentration of 50%. On the basis of the MD results, the structural deviations of the backbone atoms elucidated the dynamic effects of water/organic solvent mixtures on the equilibrium state of the protein simulations in decreasing solvent polarity. The results show that the solvent mixture gives rise to deviations in enzyme structure from the native one simulated in water. The drop in the flexibility in H₂O, MtOH, EtOH and PrOH simulation mixtures shows that greater motions of residues were influenced in BtOH and PtOH simulation mixtures. Comparing the root mean square fluctuations value with the accessible solvent area (SASA) for every residue showed an almost correspondingly high SASA value of residues to high flexibility and low SASA value to low flexibility. The study further revealed that the organic solvents influenced the formation of more hydrogen bonds in MtOH, EtOH and PrOH and thus, it is assumed that increased intraprotein hydrogen bonding is ultimately correlated to the stability of the protein. However, the solvent accessibility analysis showed that in all solvent systems, hydrophobic residues were exposed and polar residues tended to be buried away from the solvent. Distance variation of the tetrahedral intermediate packing of the active pocket was not conserved in organic solvent systems, which could lead to weaknesses in the catalytic H-bond network and most likely a drop in catalytic activity. The conformational variation of the lid domain caused by the solvent molecules influenced its gradual opening. Formation of additional hydrogen bonds and hydrophobic interactions indicates that the contribution of the cooperative network of interactions could retain the stability of the protein in some solvent systems. Time-correlated atomic motions were used to characterize the correlations between the motions of the atoms from atomic coordinates. The resulting cross-correlation map revealed that the organic solvent mixtures performed functional, concerted, correlated motions in regions of residues of the lid domain to other residues. These observations suggest that varying lengths of polar organic solvents play a significant role in introducing dynamic conformational diversity in proteins in a decreasing order of polarity.

Electrostatics-mediated α-chymotrypsin inhibition by functionalized single-walled carbon nanotubes

Electrostatics-mediated α-chymotrypsin inhibition by functionalized single-walled carbon nanotubes.

To Keep or Not to Keep? The Question of Crystallographic Waters for Enzyme Simulations in Organic Solvent

The use of enzymes in non-aqueous solvents expands the use of biocatalysts to hydrophobic substrates, with the ability to tune selectivity of reactions through solvent selection. Non-aqueous enzymology also allows for fundamental studies on the role of water and other solvents in enzyme structure, dynamics, and function. Molecular dynamics simulations serve as a powerful tool in this area, providing detailed atomic information about the effect of solvents on enzyme properties. However, a common protocol for non-aqueous enzyme simulations does not exist. If you want to simulate enzymes in non-aqueous solutions, how many and which crystallographic waters do you keep? In the present work, this question is addressed by determining which crystallographic water molecules lead most quickly to an equilibrated protein structure. Five different methods of selecting and keeping crystallographic waters are used in order to discover which crystallographic waters lead the protein structure to reach an equilibrated structure more rapidly in organic solutions. It is found that buried waters contribute most to rapid equilibration in organic solvent, with slow-diffusing waters giving similar results.

Probing immobilization mechanism of alpha-chymotrypsin onto carbon nanotube in organic media by molecular dynamics simulation

The enzyme immobilization has been adopted to enhance the activity and stability of enzymes in non-aqueous enzymatic catalysis. However, the activation and stabilization mechanism has been poorly understood on experiments. Thus, we used molecular dynamics simulation to study the adsorption of α-chymotrypsin (α-ChT) on carbon nanotube (CNT) in aqueous solution and heptane media. The results indicate that α-ChT has stronger affinity with CNT in aqueous solution than in heptane media, as confirmed by more adsorption atoms, larger contact area and higher binding free energies. Although the immobilization causes significant structure deviations from the crystal one, no significant changes in secondary structure of the enzyme upon adsorption are observed in the two media. Different from aqueous solution, the stabilization effects on some local regions far from the surface of CNT were observed in heptane media, in particular for S1 pocket, which should contribute to the preservation of specificity reported by experiments. Also, CNT displays to some extent stabilization role in retaining the catalytic H-bond network of the active site in heptane media, which should be associated with the enhanced activity of enzymes. The observations from the work can provide valuable information for improving the catalytic properties of enzymes in non-aqueous media.

Understanding the effects on constitutive activation and drug binding of a D130N mutation in the β2 adrenergic receptor via molecular dynamics simulation

G-protein-coupled receptors (GPCRs) are currently one of the largest families of drug targets. The constitutive activation induced by mutation of key GPCR residues is associated closely with various diseases. However, the structural basis underlying such activation and its role in drug binding has remained unclear. Herein, we used all-atom molecular dynamics simulations and free energy calculations to study the effects of a D130N mutation on the structure of β2 adrenergic receptor (β2AR) and its binding of the agonist salbutamol. The results indicate that the mutation caused significant changes in some key helices. In particular, the mutation leads to the departure of transmembrane 3 (TM3) from transmembrane 6 (TM6) and marked changes in the NPxxY region as well as the complete disruption of a key ionic lock, all of which contribute to the observed constitutive activation. In addition, the D130N mutation weakens some important H-bonds, leading to structural changes in these regions. Binding free energy calculations indicate that van der Waals and electrostatic interactions are the main driving forces in binding salbutamol; however, binding strength in the mutant β2AR is significantly enhanced mainly through modifying electrostatic interactions. Further analysis revealed that the increase in binding energy upon mutation stems mainly from the H-bonds formed between the hydroxyl group of salbutamol and the serine residues of TM5. This observation suggests that modifications of the H-bond groups of this drug could significantly influence drug efficacy in the treatment of diseases associated with this mutation.

Synthesis, crystal structure and spectral properties of 2-[(1-Methyl-2-benzimidazolyl)azo]-p-cresol: an experimental and theoretical study

2-[(1-Methyl-2-benzimidazolyl)azo]-p-cresol (HL), containing phenolic-OH function and benzimidazole moiety has been synthesized and characterized. The chemical, electronic structure and photophysical properties have been studied by spectroscopic analysis abetted with DFT and TDDFT calculations. The change in electronic spectra of HL by titration with aq. NaOH is studied and well supported by TDDFT calculations. The structure is confirmed by single crystal X-ray study. In the unit cell, two HL molecules are H-bonded with H2O molecule and forms dimmeric structure. The molecule forms 2D-supramolecular structure by inter-molecular H-bonding and π-π interactions.

Effects of water content on the tetrahedral intermediate of chymotrypsin - trifluoromethylketone in polar and non-polar media: observations from molecular dynamics simulation

The work uses MD simulation to study effects of five water contents (3 %, 10 %, 20 %, 50 %, 100 % v/v) on the tetrahedral intermediate of chymotrypsin - trifluoromethyl ketone in polar acetonitrile and non-polar hexane media. The water content induced changes in the structure of the intermediate, solvent distribution and H-bonding are analyzed in the two organic media. Our results show that the changes in overall structure of the protein almost display a clear correlation with the water content in hexane media while to some extent U-shaped/bell-shaped dependence on the water content is observed in acetonitrile media with a minimum/maximum at 10-20 % water content. In contrast, the water content change in the two organic solvents does not play an observable role in the stability of catalytic hydrogen-bond network, which still exhibits high stability in all hydration levels, different from observations on the free enzyme system [Zhu L, Yang W, Meng YY, Xiao X, Guo Y, Pu X, Li M (2012) J Phys Chem B 116(10):3292-3304]. In low hydration levels, most water molecules mainly distribute near the protein surface and an increase in the water content could not fully exclude the organic solvent from the protein surface. However, the acetonitrile solvent displays a stronger ability to strip off water molecules from the protein than the hexane. In a summary, the difference in the calculated properties between the two organic solvents is almost significant in low water content (<10 %) and become to be small with increasing water content. In addition, some structural properties at 10 ~ 20 % v/v hydration zone, to large extent, approach to those in aqueous solution.

Calcium-ion-induced stabilization of the protease from Bacillus cereus WQ9-2 in aqueous hydrophilic solvents: effect of calcium ion binding on the hydration shell and intramolecular interactions

The neutral protease WQ from Bacillus cereus is stable in various aqueous organic mixtures, with the exception of those containing acetonitrile (ACN) and dimethylformamide (DMF). The stability of the enzyme in aqueous hydrophilic solvents was dramatically enhanced with the addition of calcium ions, with the degree of improvement in the half-life relative to different solutions ranging from fourfold to more than 70-fold. Studies of the kinetic constants showed that calcium ions induced slight conformational changes in the active site of the enzyme in aqueous ACN. We investigated the molecular mechanisms underlying this stabilizing effect by employing a combination of biophysical techniques and molecular dynamics simulation. In aqueous ACN, the intrinsic fluorescence and circular dichroism analysis demonstrated that the addition of calcium ions induced a relatively compact conformation and maintained both the native-like microenvironment near the tryptophan residues and the secondary structure. Alternatively, homology modeling confirmed the location of four calcium-ion-binding sites in the enzyme, and molecular dynamics simulation revealed that three other calcium ions were bound to the surface of the enzyme. Calcium ions, known as a type of kosmotrope, can strongly bond with water molecules, thus aiding in the formation of the regional hydration shell required for the maintenance of enzyme activity. In addition, the introduction of calcium ions resulted in the formation of additional ionic interactions, providing propitious means for protein stabilization. Thus, the stronger intramolecular interactions were also expected to contribute partially to the enhanced stability of the enzyme in an aqueous organic solvent.