Understanding the effect of tert-butanol on Candida antarctica lipase B using molecular dynamics simulations

The use of organic solvents as reaction media for enzymatic reactions has many advantages. Several organic solvents have been proposed as reaction media, especially for transesterifications using Candida antarctica lipase B (CalB). Among organic solvents, tert-butanol is associated with an enhanced conversion rate in bio-diesel production. Thus, it is necessary to understand the effect of tert-butanol on CalB to explain the high-catalytic efficiency compared with the reaction in other hydrophilic organic solvents. In this study, the effects of tert-butanol on the structure of CalB were investigated by MD simulations. The overall flexibility was increased in the presence of tert-butanol. The substrate entrance and the binding pocket size of CalB in tert-butanol were maintained as in TIP3P water. The distance between the catalytic residues of CalB in tert-butanol indicated a higher likelihood of forming hydrogen bonds. These structural analyses could be useful for understanding the effect of tert-butanol on lipase transesterification.     

Sequence-dependent DNA deformability studied using molecular dynamics simulations

Proteins recognize specific DNA sequences not only through direct contact between amino acids and bases, but also indirectly based on the sequence-dependent conformation and deformability of the DNA (indirect readout). We used molecular dynamics simulations to analyze the sequence-dependent DNA conformations of all 136 possible tetrameric sequences sandwiched between CGCG sequences. The deformability of dimeric steps obtained by the simulations is consistent with that by the crystal structures. The simulation results further showed that the conformation and deformability of the tetramers can highly depend on the flanking base pairs. The conformations of xATx tetramers show the most rigidity and are not affected by the flanking base pairs and the xYRx show by contrast the greatest flexibility and change their conformations depending on the base pairs at both ends, suggesting tetramers with the same central dimer can show different deformabilities. These results suggest that analysis of dimeric steps alone may overlook some conformational features of DNA and provide insight into the mechanism of indirect readout during protein–DNA recognition. Moreover, the sequence dependence of DNA conformation and deformability may be used to estimate the contribution of indirect readout to the specificity of protein–DNA recognition as well as nucleosome positioning and large-scale behavior of nucleic acids.     

Analysis of ligand binding to proteins using molecular dynamics simulations

This work aims to explore theoretically the molecular mechanisms of ligand binding to proteins through the use of molecular dynamics simulations. The binding of sodium dodecyl sulfate (SDS) to cobra cardio toxin A3 (CTX A3) and thiourea (TOU) to lysozyme have been chosen as the two model systems. Data acquisitions were made by Gromacs software. To begin with, the collisions of ligand molecules with every residue of CTX A3 and lysozyme were evaluated. With this information in hand, the average numbers of collisions with each residue was defined and then assessed. Next, a measure of the affinity of a residue, P_i, referred to as conformational factor, toward a ligand molecule was established. Based on the results provided, all site-making residues for CTX A3 and lysozyme were identified. The results are in good agreement with the experimental data. Finally, based on this method, all site-making residues of bovine carbonic anhydrase (BCA) toward the SDS ligand were predicted.     

Characterization of diadzein-hemoglobin binding using optical spectroscopy and molecular dynamics simulations

The present study establishes the effectiveness of natural drug delivery mechanisms and investigates the interactions between drug and its natural carrier. The binding between the isoflavone diadzein (DZN) and the natural carrier hemoglobin (HbA) was studied using optical spectroscopy and molecular dynamics simulations. The inherent fluorescence emission characteristics of DZN along with that of tryptophan (Trp) residues of the protein HbA were exploited to elucidate the binding location and other relevant parameters of the drug inside its delivery vehicle HbA. Stern-Volmer studies at different temperatures indicate that static along with collisional quenching mechanisms are responsible for the quenching of protein fluorescence by the drug. Molecular dynamics and docking studies supported the hydrophobic interactions between ligand and protein, as was observed from spectroscopy. DZN binds between the subunits of HbA, ～15 ? away from the closest heme group of chain α1, emphasizing the fact that the drug does not interfere with oxygen binding site of HbA.     

Modeling of DNA binding to the condensin hinge domain using molecular dynamics simulations guided by atomic force microscopy

The condensin protein complex compacts chromatin during mitosis using its DNA-loop extrusion activity. Previous studies proposed scrunching and loop-capture models as molecular mechanisms for the loop extrusion process, both of which assume the binding of double-strand (ds) DNA to the hinge domain formed at the interface of the condensin subunits Smc2 and Smc4. However, how the hinge domain contacts dsDNA has remained unknown. Here, we conducted atomic force microscopy imaging of the budding yeast condensin holo-complex and used this data as basis for coarse-grained molecular dynamics simulations to model the hinge structure in a transient open conformation. We then simulated the dsDNA binding to open and closed hinge conformations, predicting that dsDNA binds to the outside surface when closed and to the outside and inside surfaces when open. Our simulations also suggested that the hinge can close around dsDNA bound to the inside surface. Based on these simulation results, we speculate that the conformational change of the hinge domain might be essential for the dsDNA binding regulation and play roles in condensin-mediated DNA-loop extrusion.     

Accurate detection of protein:ligand binding sites using molecular dynamics simulations

Accurate prediction of location of cavities and surface grooves in proteins is important, as these are potential sites for ligand binding. Several currently available programs for cavity detection are unable to detect cavities near the surface or surface grooves. In the present study, an optimized molecular dynamics based procedure is described for detection and quantification of interior cavities as well as surface pockets. This is based on the observation that the mobility of water in such pockets is significantly lower than that of bulk water. The algorithm efficiently detects surface grooves that are sites of protein-ligand and protein-protein interaction. The algorithm was also used to substantially improve the performance of an automated docking procedure for docking monomers of nonobligate protein-protein complexes. In addition, it was applied to predict key residues involved in the binding of the E. coli toxin CcdB with its inhibitor. Predictions were subsequently validated by mutagenesis experiments.     

Understanding beta-hairpin formation by molecular dynamics simulations of unfolding

We have studied the mechanism of formation of a 16-residue beta-hairpin from the protein GB1 using molecular dynamics simulations in an aqueous environment. The analysis of unfolding trajectories at high temperatures suggests a refolding pathway consisting of several transient intermediates. The changes in the interaction energies of residues are related with the structural changes during the unfolding of the hairpin. The electrostatic energies of the residues in the turn region are found to be responsible for the transition between the folded state and the hydrophobic core state. The van der Waals interaction energies of the residues in the hydrophobic core reflect the behavior of the radius of gyration of the core region. We have examined the opposing influences of the protein-protein (PP) energy, which favors the native state, and the protein-solvent (PS) energy, which favors unfolding, in the formation of the beta-hairpin structure. It is found that the behavior of the electrostatic components of PP and PS energies reflects the structural changes associated with the loss of backbone hydrogen bonding. Relative changes in the PP and PS van der Waals interactions are related with the disruption of the hydrophobic core of a protein. The results of the simulations support the hydrophobic collapse mechanism of beta-hairpin folding.     

Absolute binding free energy calculations using molecular dynamics simulations with restraining potentials

The absolute (standard) binding free energy of eight FK506-related ligands to FKBP12 is calculated using free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent. A number of features are implemented to improve the accuracy and enhance the convergence of the calculations. First, the absolute binding free energy is decomposed into sequential steps during which the ligand-surrounding interactions as well as various biasing potentials restraining the translation, orientation, and conformation of the ligand are turned "on" and "off." Second, sampling of the ligand conformation is enforced by a restraining potential based on the root mean-square deviation relative to the bound state conformation. The effect of all the restraining potentials is rigorously unbiased, and it is shown explicitly that the final results are independent of all artificial restraints. Third, the repulsive and dispersive free energy contribution arising from the Lennard-Jones interactions of the ligand with its surrounding (protein and solvent) is calculated using the Weeks-Chandler-Andersen separation. This separation also improves convergence of the FEP/MD calculations. Fourth, to decrease the computational cost, only a small number of atoms in the vicinity of the binding site are simulated explicitly, while all the influence of the remaining atoms is incorporated implicitly using the generalized solvent boundary potential (GSBP) method. With GSBP, the size of the simulated FKBP12/ligand systems is significantly reduced, from approximately 25,000 to 2500. The computations are very efficient and the statistical error is small ( approximately 1 kcal/mol). The calculated binding free energies are generally in good agreement with available experimental data and previous calculations (within approximately 2 kcal/mol). The present results indicate that a strategy based on FEP/MD simulations of a reduced GSBP atomic model sampled with conformational, translational, and orientational restraining potentials can be computationally inexpensive and accurate.     

Validation of the mechanism of cholesterol binding by StAR using short molecular dynamics simulations

We previously proposed an original two-state cholesterol binding mechanism by StAR, in which the C-terminal α-helix of StAR gates the access of cholesterol to its binding site cavity. This cavity, which can accommodate one cholesterol molecule, was proposed to promote the reversible unfolding of the C-terminal α-helix and allow for the entry and dissociation of cholesterol. In our molecular model of the cholesterol–StAR complex, the hydrophobic moiety of cholesterol interacts with hydrophobic amino acid side-chains located in the C-terminal α-helix and at the bottom of the cavity. In this study, we present a structural in silico analysis of StAR. Molecular dynamics simulations showed that point mutations of Phe²⁶⁷, Leu²⁷¹ or Leu²⁷⁵ at the α-helix 4 increased the gyration radius (more flexibility) of the protein's structure, whereas the salt bridge double mutant E169M/R188M showed a decrease in flexibility (more compactness). Also, in the latter case, an interaction between Met¹⁶⁹ and Phe²⁶⁷ disrupted the hydrophobic cavity, rendering it impervious to ligand binding. These obtained results are in agreement with previous in vitro experiments, and provide further validation of the two-state binding mode of action.     

Interdomain dynamics and ligand binding: molecular dynamics simulations of glutamine binding protein

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.     

Characterizing the binding of annexin V to a lipid bilayer using molecular dynamics simulations

Annexins play critical roles in membrane organization, membrane trafficking and vesicle transport. The family members share the ability to bind to membranes with high affinities, but the interactions between annexins and membranes remain unclear. Here, using long‐time molecular dynamics simulations, we provide detailed information for the binding of an annexin V trimer to a POPC/POPS lipid bilayer. Calcium ions function as bridges between several negatively charged residues of annexin V and the oxygen atoms of lipids. The preferred calcium‐bridges are those formed via the carboxyl oxygen atoms of POPS lipids. H‐bonds and hydrophobic interactions formed by several critical residues have also been observed in the annexin‐membrane interface. The annexin‐membrane binding causes small changes of annexin trimer structures, while has significant effects on lipid bilayer structures. The lipid bilayer shows a bent shape and forms a concave region in the annexin‐membrane interaction interface, which provides an atomic‐level evidence to support the view that annexins could disturb the stability of lipids and bend membranes. This study provides insights into the commonly occurring PS‐dependent and calcium‐dependent binding of proteins to membranes. Proteins 2014; 82:312–322. © 2013 Wiley Periodicals, Inc.     

HLA-DP2 binding prediction by molecular dynamics simulations

Major histocompatibility complex (MHC) II proteins bind peptide fragments derived from pathogen antigens and present them at the cell surface for recognition by T cells. MHC proteins are divided into Class I and Class II. Human MHC Class II alleles are grouped into three loci: HLA-DP, HLA-DQ, and HLA-DR. They are involved in many autoimmune diseases. In contrast to HLA-DR and HLA-DQ proteins, the X-ray structure of the HLA-DP2 protein has been solved quite recently. In this study, we have used structure-based molecular dynamics simulation to derive a tool for rapid and accurate virtual screening for the prediction of HLA-DP2-peptide binding. A combinatorial library of 247 peptides was built using the "single amino acid substitution" approach and docked into the HLA-DP2 binding site. The complexes were simulated for 1 ns and the short range interaction energies (Lennard-Jones and Coulumb) were used as binding scores after normalization. The normalized values were collected into quantitative matrices (QMs) and their predictive abilities were validated on a large external test set. The validation shows that the best performing QM consisted of Lennard-Jones energies normalized over all positions for anchor residues only plus cross terms between anchor-residues.     

Mapping the phase diagram of the writhe of DNA nanocircles using atomistic molecular dynamics simulations

We have investigated the effects of duplex length, sequence, salt concentration and superhelical density on the conformation of DNA nanocircles containing up to 178 base pairs using atomistic molecular dynamics simulation. These calculations reveal that the partitioning of twist and writhe is governed by a delicate balance of competing energetic terms. We have identified conditions which favour circular, positively or negatively writhed and denatured DNA conformations. Our simulations show that AT-rich DNA is more prone to denaturation when subjected to torsional stress than the corresponding GC containing circles. In contrast to the behaviour expected for a simple elastic rod, there is a distinct asymmetry in the behaviour of over and under-wound DNA nanocircles. The most biologically relevant negatively writhed state is more elusive than the corresponding positively writhed conformation, and is only observed for larger circles under conditions of high electrostatic screening. The simulation results have been summarised by plotting a phase diagram describing the various conformational states of nanocircles over the range of circle sizes and experimental conditions explored during the study. The changes in DNA structure that accompany supercoiling suggest a number of mechanisms whereby changes in DNA topology in vivo might be used to influence gene expression.     

Understanding the interactions of different substrates with wild-type and mutant acylaminoacyl peptidase using molecular dynamics simulations

Acylaminoacylpeptidase (AAP) belongs to peptidase protein family, which can degrade amyloid β-peptide forms in the brains of patients, and hence leads to Alzheimer’s disease. And so, AAP is considered to be a novel target in the design of drugs against Alzheimer’s disease. In this investigation, six molecular dynamics simulations were used to find that the interaction between the wild-type and R526V AAP with two different substrates (p-nitrophenylcaprylate and Ac-Leu-p-nitroanilide). Our results were as follows: firstly, Ac-Leu-p-nitroanilide bound to R526V AAP to form a more disordered loop (residues 552–562) in the α/β-hydrolase fold like of AAP, which caused an open and inactive AAP domain form, secondly, binding p-nitrophenylcaprylate and Ac-Leu-p-nitroanilide to AAP can decrease the flexibility of residues 225–250, 260–270, and 425–450, in which the ordered secondary structures may contain the suitable geometrical structure and so it is useful to serine attack. Our theoretical results showed that the binding of the two substrates can induce specific conformational changes responsible for the diverse AAP catalytic specificity. These theoretical substrate-induced structural diversities can help explain the abilities of AAPs to recognize and hydrolyze extremely different substrates.     

Non-B DNA conformations analysis through molecular dynamics simulations

BackgroundNon-B DNA conformations are molecular structures that do not follow the canonical DNA double helix. Mutagenetic instability in nuclear and mitochondrial DNA (mtDNA) genomes has been associated with simple non-B DNA conformations, as hairpins or more complex structures, as G-quadruplexes. One of these structures is Structure A, a cloverleaf-like non-B conformation predicted for a 93-nt (nucleotide) stretch of the mtDNA control region 5′-peripheral domain. Structure A is embedded in a hot spot for the 3′ end of human mtDNA deletions revealing its importance in influencing the mutational instability of the mtDNA genome.MethodsTo better characterize Structure A, we predicted its 3D conformation using state-of-art methods and algorithms. The methodologic workflow consisted in the prediction of non-B conformations using molecular dynamics simulations. The conservation scores of alignments of the Structure A region in humans, primates, and mammals, was also calculated.ResultsOur results show that these computational methods are able to measure the stability of non-B conformations by using the level of base pairing during molecular dynamics. Structure A showed high stability and low flexibility correlated with high conservation scores in mammalian, more specifically in primate lineages.ConclusionsWe showed that 3D non-B conformations can be predicted and characterized by our methodology. This allowed the in-depth analysis of the structure A, and the main results showed the structure remains stable during the simulations.General significanceThe fine-scale atomic molecular determination of this type of non-B conformation opens the way to perform computational molecular studies that can show their involvement in mtDNA cellular mechanisms.     

Understanding peptide competitive inhibition of botulinum neurotoxin a binding to SV2 protein via molecular dynamics simulations

Botulinum neurotoxins (BoNTs) are known as the most toxic natural substances. Synaptic vesicle protein 2 (SV2) has been proposed to be a protein receptor for BoNT/A. Recently, two short peptides (BoNT/A‐A2 and SV2C‐A3) were designed to inhibit complex formation between the BoNT/A receptor‐binding domain (BoNT/A‐RBD) and the synaptic vesicle protein 2C luminal domain (SV2C‐LD). In this article, the two peptide complex systems are studied by molecular dynamics (MD) simulations. The structural stability analysis indicates that BoNT/A‐A2 system is more stable than SV2C‐A3 system. The conformational analysis implies that the β‐sheet in BoNT/A‐A2 system maintains its secondary structure but the two β‐strands in SV2C‐A3 system have remarkable conformational changes. Based on the calculation of hydrogen bonds, hydrophobic interactions and cation‐π interactions, it is found that the internal hydrogen bonds play crucial roles in the structural stability of the peptides. Because of the stable secondary structure, the β‐sheet in BoNT/A‐A2 system establishes effective interactions at the interface and inhibits BoNT/A‐RBD binding to SV2C‐LD. In contrast, without other β‐strands forming internal hydrogen bonds, the two isolated β‐strands in SV2C‐A3 system become the random coil. This conformational change breaks important hydrogen bonds and weakens cation‐π interaction in the interface, so the complex formation is only partially inhibited by the two β‐strands. These results are consistent with experimental studies and may be helpful in understanding the inhibition mechanisms of peptide inhibitors. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 597–608, 2015.     

DNA basepair step deformability inferred from molecular dynamics simulations

The sequence-dependent DNA deformability at the basepair step level was investigated using large-scale atomic resolution molecular dynamics simulation of two 18-bp DNA oligomers: d(GCCTATAAACGCCTATAA) and d(CTAGGTGGATGACTCATT). From an analysis of the structural fluctuations, the harmonic potential energy functions for all 10 unique steps with respect to the six step parameters have been evaluated. In the case of roll, three distinct groups of steps have been identified: the flexible pyrimidine-purine (YR) steps, intermediate purine-purine (RR), and stiff purine-pyrimidine (RY). The YR steps appear to be the most flexible in tilt and partially in twist. Increasing stiffness from YR through RR to RY was observed for rise, whereas shift and slide lack simple trends. A proposed measure of the relative importance of couplings identifies the slide-rise, twist-roll, and twist-slide couplings to play a major role. The force constants obtained are of similar magnitudes to those based on a crystallographic ensemble. However, the current data have a less complicated and less pronounced sequence dependence. A correlation analysis reveals concerted motions of neighboring steps and thus exposes limitations in the dinucleotide model. The comparison of DNA deformability from this and other studies with recent quantum-chemical stacking energy calculations suggests poor correlation between the stacking and flexibility.     

Protein flexibility using constraints from molecular dynamics simulations

Proteins are held together in the native state by hydrophobic interactions, hydrogen bonds and interactions with the surrounding water, whose strength as well as spatial and temporal distribution affects protein flexibility and hence function. We study these effects using 10 ns molecular dynamics simulations of pure water and of two proteins, the glutamate receptor ligand binding domain and barnase. We find that most of the noncovalent interactions flicker on and off over typically nanoseconds, and so we can obtain good statistics from the molecular dynamics simulations. Based on this information, a topological network of rigid bonds corresponding to a protein structure with covalent and noncovalent bonds is constructed, with account being taken of the influence of the flickering hydrogen bonds. We define the duty cycle for the noncovalent interactions as the percentage of time a given interaction is present, which we use as an input to investigate flexibility/rigidity patterns, in the algorithm FIRST which constructs and analyses topological networks.