Conformational selection in silico: loop latching motions and ligand binding in enzymes

Ligand binding frequently induces significant conformational changes in a protein receptor. Understanding and predicting such conformational changes represent an important challenge for computational biology, including applications to structure-based drug design. We describe an approach to this problem based on the assumption that the holo state is at least transiently populated in the absence of a ligand; this hypothesis has been referred to as "conformational selection." Here, we apply a method that tests this hypothesis on a challenging class of ligand-induced conformational changes, which we refer to as loop latching: the closing of a loop around an active site that sequesters the ligand from solvent. The method uses a combination of replica exchange molecular dynamics and a loop prediction algorithm to generate low-energy loop structures, and docking to select the conformation appropriate for binding a particular ligand. On a test set of six proteins, it yields loop structures including hololike conformations, generally below 2 A RMSD from the liganded structure, for loops that span up to 15 residues. Docking serves as a stringent test of the predictions. In five of the six cases, the predicted loop conformations improve the ranks of cognate ligands relative to using the apo structure, although the results remain, in most cases, significantly worse than using a holo structure. The poses of the cognate ligands are correct in four of the six test cases, while they are correct for five of the six using a holo structure.     

Essential motions in a fungal lipase with bound substrate,covalently attached inhibitor and product

As an aid to understanding the influence of dynamic fluctuations during esterolytic catalysis, we follow protein flexibility at three different steps along the catalytic pathway from substrate binding to product clearance via a covalently attached inhibitor, which represents a transition-state mimic. We have applied a classical approach, using molecular dynamics simulations to monitor protein dynamics in the nanosecond regime. We filter out small amplitude fluctuations and focus on the anharmonic contributions to the overall dynamics. This 'essential dynamics' analysis reveals different modes of response along the pathway suggesting that binding, catalysis and product clearance occur along different energy surfaces. Motions in the enzyme with a covalently attached ligand are more complex and occur along several eigenvectors. The magnitudes of the fluctuations in these individual subspaces are significantly smaller than those observed for the substrate and product molecules, indicating that the energy surface is shallow and that a relatively large number of conformational substates are accessible. On the other hand, substrate binding and product release occur at distinct modes of the protein flexibility suggesting that these processes occur along rough energy surfaces with only a few minima. Detailed energetic analyses along the trajectories indicated that in all cases binding is dominated by van der Waals interactions. The carboxylate form of the product is stabilized by a tight hydrogen bond network involving in particular Ser82, which may be a potential cause of product inhibition. Considerations such as these should aid the understanding of mechanisms of substrate, inhibitor or product recognition and could become of importance in the design of new substrates or inhibitors for enzymes.     

Insight derived from molecular dynamics simulation into dynamics and molecular motions of cuticle-degrading serine protease Ver112

Cuticle-degrading serine protease Ver112, which derived from a nematophagous fungus Lecanicillium psalliotae, has been exhibited to have high cuticle-degrading and nematicidal activities. We have performed molecular dynamics (MD) simulation based on the crystal structure of Ver112 to investigate its dynamic properties and large-scale concerted motions. The results indicate that the structural core of Ver112 shows a small fluctuation amplitude, whereas the substrate binding sites, and the regions close to and opposite the substrate binding sites experience significant conformational fluctuations. The large concerted motions obtained from essential dynamics (ED) analysis of MD trajectory can lead to open or close of the substrate binding sites, which are proposed to be linked to the functional properties of Ver112, such as substrate binding, orientation, catalytic, and release. The significant motion in the loop regions that is located opposite the binding sites are considered to play an important role in modulating the dynamics of the substrate binding sites. Furthermore, the bottom of free energy landscape (FEL) of Ver112 are rugged, which is mainly caused by the fluctuations of substrate binding regions and loops located opposite the binding site. In addition, the mechanism underlying the high flexibility and catalytic activity of Ver112 was also discussed. Our simulation study complements the biochemical and structural studies, and provides insight into the dynamics-function relationship of cuticle-degrading serine protease Ver112.     

Subtle functional collective motions in pancreatic-like ribonucleases: from ribonuclease A to angiogenin

The analysis of the dynamic behavior of enzymes is fundamental to structural biology. A direct relationship between protein flexibility and biological function has been shown for bovine pancreatic ribonuclease (RNase A) (Rasmussen et al., Nature 1992;357:423-424). More recently, crystallographic studies have shown that functional motions in RNase A involve the enzyme beta-sheet regions that move concertedly on substrate binding and release (Vitagliano et al., Proteins 2002;46:97-104). These motions have been shown to correspond to intrinsic dynamic properties of the native enzyme by molecular dynamics (MD) simulations. To unveil the occurrence of these collective motions in other members of pancreatic-like superfamily, we carried out MD simulations on human angiogenin (Ang). Essential dynamics (ED) analyses performed on the trajectories reveal that Ang exhibits collective motions similar to RNase A, despite the limited sequence identity (33%) of the two proteins. Furthermore, we show that these collective motions are also present in ensembles of experimentally determined structures of both Ang and RNase A. Finally, these subtle concerted beta-sheet motions were also observed for other two members of the pancreatic-like superfamily by comparing the ligand-bound and ligand-free structures of these enzymes. Taken together, these findings suggest that pancreatic-like ribonucleases share an evolutionary conserved dynamic behavior consisting of subtle beta-sheet motions, which are essential for substrate binding and release.     

Computational insights into the interaction of the anthrax lethal factor with the N-terminal region of its substrates

The anthrax lethal factor is a zinc metalloprotease toxin secreted by Bacillus anthracis which cleaves at the N-terminal region of six mitogen activated protein kinase kinases (MEKs) in the cell. Additionally, it is known to cleave a nine residue peptide "LF10," 50-fold more efficiently than nine residues of MEK1. There is very little sequence similarity between the MEK N-termini, thus, it is unclear how the lethal factor can accommodate and cleave the diverse N-termini of the MEKs and whether there is a hierarchy in this interaction, as there is between LF10 and MEK1. To investigate this problem, we carried out multiple molecular dynamics simulations of the lethal factor with nine residues of each of the substrates. Our simulations reveal that like LF10, certain MEK substrates have residue compositions that favor beta-sheet formation with the lethal factor over others. The formation of this secondary structure maintains a catalytic conformation. Binding energetics using the MM-PBSA method was used to rank-order the substrates for their affinity to LF (K(M)). On the basis of the results, we conclude that the LF does not equally accommodate the MEK substrates and further predict that there will be differences between rates of cleavage among the nine residue MEK N-termini.     

In silico screening,molecular docking,and molecular dynamics studies of SNP-derived human P5CR mutants

Pyrroline-5-carboxylate reductase (P5CR) encoded by PYCR1 gene is a housekeeping enzyme that catalyzes the reduction of P5C to proline using NAD(P)H as the cofactor. In this study, we used in silico approaches to examine the role of nonsynonymous single-nucleotide polymorphisms in the PYCR1 gene and their putative functions in the pathogenesis of Cutis Laxa. Among the 348 identified SNPs, 15 were predicted to be potentially damaging by both SIFT and PolyPhen tools; of them two SNP‐derived mutations, R119G and G206W, have been previously reported to correlate with Cutis Laxa. These two mutations were therefore selected to be mapped to the wild‐type (WT) P5CR structure for further structural and functional analyses. The results of comparative computational analyses using I-Mutant and Autodock reveal reductions in both stability and cofactor binding affinity of these two mutants. Comparative molecular dynamics (MD) simulations were performed to evaluate the changes in dynamic properties of P5CR upon mutations. The results reveal that the two mutations enhance the rigidity of P5CR structure, especially that of cofactor binding site, which could result in decreased kinetics of cofactor entrance and egress. Comparison between the structural properties of the WT and mutants during MD simulations shows that the enhanced rigidity of mutants results most likely from the increased number of inter‐atomic interactions and the decreased number of dynamic hydrogen bonds. Our study provides novel insight into the deleterious effects of the R119G and G206W mutations on P5CR, and sheds light on the mechanisms by which these mutations mediate Cutis Laxa.     

Effect of pH on structure,function, and stability of mitochondrial carbonic anhydrase VA

Mitochondrial carbonic anhydrase VA (CAVA) catalyzes the hydration of carbon dioxide to produce proton and bicarbonate which is primarily expressed in the mitochondrial matrix of liver, and involved in numerous physiological processes including lipogenesis, insulin secretion from pancreatic cells, ureagenesis, gluconeogenesis, and neuronal transmission. To understand the effect of pH on the structure, function, and stability of CAVA, we employed spectroscopic techniques such as circular dichroism, fluorescence, and absorbance measurements in wide range of pH (from pH 2.0 to pH 11.5). CAVA showed an aggregation at acidic pH range from pH 2.0 to pH 5.0. However, it remains stable and maintains its secondary structure in the pH range, pH 7.0–pH 11.5. Furthermore, this enzyme has an appreciable activity at more than pH 7.0 (7.0 < pH ≤ 11.5) with maximum activity at pH 9.0. The maximal values of k_cat and k_cat/K_m at pH 9.0 are 3.7?×?10⁶ s^?1 and 5.5?×?10⁷ M^?1 s^?1, respectively. However, this enzyme loses its activity in the acidic pH range. We further performed 20-ns molecular dynamics simulation of CAVA to see the dynamics at different pH values. An excellent agreement was observed between in silico and in vitro studies. This study provides an insight into the activity of CAVA in the pH range of subcellular environment.     

Computational,pulse-radiolytic,and structural investigations of lysine-136 and its role in the electrostatic triad of human C u,Z n superoxide dismutase

Key charged residues in Cu,Zn superoxide dismutase (Cu,Zn SOD) promote electrostatic steering of the superoxide substrate to the active site Cu ion, resulting in dismutation of superoxide to oxygen and hydrogen peroxide. Lys-136, along with the adjacent residues Glu-132 and Glu-133, forms a proposed electrostatic triad contributing to substrate recognition. Human Cu,Zn SODs with single-site replacements of Lys-136 by Arg, Ala, Gln, or Glu or with a triple-site substitution (Glu-132 and Glu-133 to Gln and Lys-136 to Ala) were made to test hypotheses regarding contributions of these residues to Cu,Zn SOD activity. The structural effects of these mutations were modeled computationally and validated by the X-ray crystallographic structure determination of Cu,Zn SOD having the Lys-136-to-Glu replacement. Brownian dynamics simulations and multiple-site titration calculations predicted mutant reaction rates as well as ionic strength and pH effects measured by pulse-radiolytic experiments. Lys-136-to-Glu charge reversal decreased dismutation activity 50% from 2.2 × 10⁹ to 1.2 × 10⁹ M⁻¹ s⁻¹ due to repulsion of negatively charged superoxide, whereas charge-neutralizing substitutions (Lys-136 to Gln or Ala) had a less dramatic influence. In contrast, the triple-mutant Cu,Zn SOD (all three charges in the electrostatic triad neutralized) surprisingly doubled the reaction rate compared with wild-type enzyme but introduced phosphate inhibition. Computational and experimental reaction rates decreased with increasing ionic strength in all of the Lys-136 mutants, with charge reversal having a more pronounced effect than charge neutralization, implying that local electrostatic effects still govern the dismutation rates. Multiple-site titration analysis showed that deprotonation events throughout the enzyme are likely responsible for the gradual decrease in SOD activity above pH 9.5 and predicted a pK_a value of 11.7 for Lys-136. Overall, Lys-136 and Glu-132 make comparable contributions to substrate recognition but are less critical to enzyme function than Arg-143, which is both mechanistically and electrostatically essential. Thus, the sequence-conserved residues of this electrostatic triad are evidently important solely for their electrostatic properties, which maintain the high catalytic rate and turnover of Cu,Zn SOD while simultaneously providing specificity by selecting against binding by other anions. Proteins 29:103–112, 1997. © 1997 Wiley-Liss, Inc.     

Clamping,bending, and twisting inter-domain motions in the misfold-recognizing portion of UDP-glucose: Glycoprotein glucosyltransferase

1. Download : Download high-res image (119KB)
2. Download : Download full-size image

     

Investigating the inter‐subunit/subdomain interactions and motions relevant to disease mutations in the N‐terminal domain of ryanodine receptors by molecular dynamics simulation

The ryanodine receptors (RyR) are essential to calcium signaling in striated muscles, and numerous disease mutations have been identified in two RyR isoforms, RyR1 in skeletal muscle and RyR2 in cardiac muscle. A deep understanding of the activation/regulation mechanisms of RyRs has been hampered by the shortage of high‐resolution structures and dynamic information for this giant tetrameric complex in different functional states. Toward elucidating the molecular mechanisms of disease mutations in RyRs, we performed molecular dynamics simulation of the N‐terminal domain (NTD) which is not only the best‐resolved structural component of RyRs, but also a hotspot of disease mutations. First, we simulated the tetrameric NTD of wild‐type RyR1 and three disease mutants (K155E, R157Q, and R164Q) that perturb the inter‐subunit interfaces. Our simulations identified a dynamic network of salt bridges involving charged residues at the inter‐subunit/subdomain interfaces and disease‐mutation sites. By perturbing this key network, the above three mutations result in greater flexibility with the highest inter‐subunit opening probability for R157Q. Next, we simulated the monomeric NTD of RyR2 in the presence or absence of a central Cl^– anion which is known to stabilize the interfaces between the three NTD subdomains (A, B, and C). We found that the loss of Cl^– restructures the salt‐bridge network near the Cl^–‐binding site, leading to rotations of subdomain A/B relative to subdomain C and enhanced mobility between the subdomains. This finding supports a mechanism for disease mutations in the NTD of RyR2 via perturbation of the Cl^– binding. The rich structural and dynamic information gained from this study will guide future mutational and functional studies of the NTD of RyRs. Proteins 2017; 85:1633–1644. © 2017 Wiley Periodicals, Inc.     

The interface of protein structure, protein biophysics, and molecular evolution

Abstract The interface of protein structural biology, protein biophysics, molecular evolution, and molecular population genetics forms the foundations for a mechanistic understanding of many aspects of protein biochemistry. Current efforts in interdisciplinary protein modeling are in their infancy and the state-of-the art of such models is described. Beyond the relationship between amino acid substitution and static protein structure, protein function, and corresponding organismal fitness, other considerations are also discussed. More complex mutational processes such as insertion and deletion and domain rearrangements and even circular permutations should be evaluated. The role of intrinsically disordered proteins is still controversial, but may be increasingly important to consider. Protein geometry and protein dynamics as a deviation from static considerations of protein structure are also important. Protein expression level is known to be a major determinant of evolutionary rate and several considerations including selection at the mRNA level and the role of interaction specificity are discussed. Lastly, the relationship between modeling and needed high-throughput experimental data as well as experimental examination of protein evolution using ancestral sequence resurrection and in vitro biochemistry are presented, towards an aim of ultimately generating better models for biological inference and prediction.     

Solution Conformation of [D-Pen2,D-Pen5]enkephalin in Water: A NMR and Molecular Dynamics Study

The solution conformation of [D -Pen²,D -Pen⁵] enkephalin (DPDPE), a highly potent δ-selective opioid agonist, was examined by means of NMR, molecular mechanics and molecular dynamics methods. The structural information in the solvent water was obtained employing one- and two-dimensional methods of ¹H and ¹³C-NMR spectroscopy. Based on the distance geometry technique using the ROE data as input, 400 conformers were obtained and considered in the structure analysis. Alternatively, about 2000 conformers were stochastically generated and related to the NMR data after energy minimization. The structure analysis provides one conformer in agreement with all NMR data, which belongs to the lowest energy conformation group. This structure may serve as a reference conformer for DPDPE analogues synthesized with the aim of activity increase.     

Intrinsic flexibility of NLRP pyrin domains is a key factor in their conformational dynamics,fold stability,and dimerization

Nucleotide‐binding domain leucine‐rich repeat‐containing receptors (NLRs) are key proteins in the innate immune system. The 14 members of the NLRP family of NLRs contain an N‐terminal pyrin domain which is central for complex formation and signal transduction. Recently, X‐ray structures of NLRP14 revealed an unexpected rearrangement of the α5/6 stem‐helix of the pyrin domain allowing a novel symmetric dimerization mode. We characterize the conformational transitions underlying NLRP oligomerization using molecular dynamics simulations. We describe conformational stability of native NLRP14 and mutants in their monomeric and dimeric states and compare them to NLRP4, a representative of a native pyrin domain fold. Thereby, we characterize the interplay of conformational dynamics, fold stability, and dimerization in NLRP pyrin domains. We show that intrinsic flexibility of NLRP pyrin domains is a key factor influencing their behavior in physiological conditions. Additionally, we provide further evidence for the crucial importance of a charge relay system within NLRPs that critically influences their conformational ensemble in solution.     

Evaluation of the binding effect and cytotoxicity assay of 2-Ethyl-5-(4-methylphenyl) pyramido pyrazole ophthalazine trione on calf thymus DNA: spectroscopic,calorimetric, and molecular dynamics approaches

With advances in new drug therapies, it is essential to understand the interactions between drugs and target molecules. In this study, we applied multiple spectroscopic techniques including absorbance, fluorescence, circular dichroism spectroscopy, viscosity, thermal melting, calorimetric, and molecular dynamics (MD) simulation to study the interaction between 2-Ethyl-5-(4-methylphenyl) pyramido pyrazole ophthalazine trione (PPF) and calf thymus DNA (ct DNA) in the absence or presence of histone H1. PPF exhibits a high binding affinity towards ct DNA in binary and ternary systems. In addition, the result for the binding constant was observed within the range 10⁴ M⁻¹ achieved through fluorescence quenching data, while the values for enthalpy and entropy changes for ct DNA–PPF and (ct DNA–H1) PPF complexes were measured to be −72.54 kJ.mol⁻¹, −161.14 J.mol⁻¹ K⁻¹, −85.34 kJ.mol⁻¹, and −19.023 J.mol⁻¹ K⁻¹, respectively. Furthermore, in accordance with circular dichroism spectra, the inducement of ct DNA structural changes was observed during binding of PPF and H1 in binary and ternary system forms. The essential roles of hydrogen bonding and van der Waals forces throughout the interaction were suggested using thermodynamic parameters. According to the obtained data, the interaction mode of ct DNA–PPF and (ct DNA–H1) PPF complexes was intercalation binding. Suggested by the MD simulation study, the ct DNA–H1 complex caused a reduction in the stability of the DNA structure in the presence or absence of ligand, which demonstrated that PPF as an intercalating agent can further distort the structure. The information achieved from this study will be very helpful in understanding the effects of PPF on the conformational state of ct DNA in the absence or presence of the H1 molecule, which seems to be quite significant for clarifying the mechanisms of action and its pharmacokinetics.     

CD4-gp120 interaction interface - a gateway for HIV-1 infection in human: molecular network,modeling and docking studies

The major causative agent for Acquired Immune Deficiency Syndrome (AIDS) is Human Immunodeficiency Virus-1 (HIV-1). HIV-1 is a predominant subtype of HIV which counts on human cellular mechanism virtually in every aspect of its life cycle. Binding of viral envelope glycoprotein-gp120 with human cell surface CD4 receptor triggers the early infection stage of HIV-1. This study focuses on the interaction interface between these two proteins that play a crucial role for viral infectivity. The CD4–gp120 interaction interface has been studied through a comprehensive protein–protein interaction network (PPIN) analysis and highlighted as a useful step towards identifying potential therapeutic drug targets against HIV-1 infection. We prioritized gp41, Nef and Tat proteins of HIV-1 as valuable drug targets at early stage of viral infection. Lack of crystal structure has made it difficult to understand the biological implication of these proteins during disease progression. Here, computational protein modeling techniques and molecular dynamics simulations were performed to generate three-dimensional models of these targets. Besides, molecular docking was initiated to determine the desirability of these target proteins for already available HIV-1 specific drugs which indicates the usefulness of these protein structures to identify an effective drug combination therapy against AIDS.     

Comprehensive modeling and functional analysis of Toll‐like receptor ligand‐recognition domains

Toll‐like receptors (TLRs) are innate immune pattern‐recognition receptors endowed with the capacity to detect microbial pathogens based on pathogen‐associated molecular patterns. The understanding of the molecular principles of ligand recognition by TLRs has been greatly accelerated by recent structural information, in particular the crystal structures of leucine‐rich repeat‐containing ectodomains of TLR2, 3, and 4 in complex with their cognate ligands. Unfortunately, for other family members such as TLR7, 8, and 9, no experimental structural information is currently available. Methods such as X‐ray crystallography or nuclear magnetic resonance are not applicable to all proteins. Homology modeling in combination with molecular dynamics may provide a straightforward yet powerful alternative to obtain structural information in the absence of experimental (structural) data, provided that the generated three‐dimensional models adequately approximate what is found in nature. Here, we report the development of modeling procedures tailored to the structural analysis of the extracellular domains of TLRs. We comprehensively compared secondary structure, torsion angles, accessibility for glycosylation, surface charge, and solvent accessibility between published crystal structures and independently built TLR2, 3, and 4 homology models. Finding that models and crystal structures were in good agreement, we extended our modeling approach to the remaining members of the TLR family from human and mouse, including TLR7, 8, and 9.     

Spectroscopy and computational studies on the interaction of octyl,dodecyl, and hexadecyl derivatives of anionic and cationic surfactants with adenosine deaminase

Effects of sodium (octyl, dodecyl, hexadecyl) sulfate and their cationic analogous on the structure of adenosine deaminase (ADA) were investigated by fluorescence and circular dichroism spectroscopy as well as molecular dynamics simulation and docking calculation. Root-mean-square derivations, radius of gyration, solvent accessible surface area, and radial distribution function were obtained. The results showed that anionic and cationic surfactants reduce protein stability. Cationic surfactants have more effect on the ADA structure in comparison with anionic surfactants. More concentration and longer surfactants are parallel to higher denaturation. Furthermore, aggregation in the presence of anionic surfactants is more than cationic surfactants. Docking data showed that longer surfactants have more interaction energy and smaller ones bound to the active site.