Use of normal modes for structural modeling of proteins: the case study of rat heme oxygenase 1

We present an original approach based on full-atom normal mode analysis (NMA) aimed to expand the general framework of homology modeling. Using the rat heme-free oxygenase 1 as a case system, we show how NMA can be used to model different physiologically relevant conformations of the same protein. Starting from a unique heme-bound X-ray structure, and using two structural templates corresponding to a human and an incomplete rat heme-free structures, we generate models of the rat unbound species with open and closed conformations. Less than 100 lowest frequency modes of the target were sufficient to obtain the heme-free conformations, the closest to the templates. The rat HO-1 model built for the open form shows features similar to the open form of the human heme-free oxygenase, and the one built for the closed form was similar to the incompletely resolved X-ray structure of the same protein available in the Protein DataBank. In the latter case, the use of NMA was particularly useful since it allowed to build a complete structure and therefore to discuss on the reason of the structural differences between open and closed forms. This study shows that the amount of main chain flexibility provided by the normal modes can lead to major improvements in homology modeling approaches. Such applications will allow the characterization of alternative conformations of a target protein with respect to the templates and/or the construction of good quality 3D models based on existing templates with unresolved parts in their tertiary structure.     

Domain motions and the open-to-closed conformational transition of an enzyme: a normal mode analysis of S-adenosyl-L-homocysteine hydrolase

The structure and fluctuations of the enzyme S-adenosyl-L-homocysteine hydrolase (SAHH) are analyzed in an effort to explain its biological function. Besides the previously identified open structure, characteristic of the substrate-free enzyme, we find two distinct structures in enzyme-inhibitor complexes, the closed and closed-twisted conformers. Both closed conformers differ from the open form by a hinge bending motion of two large domains within each subunit, which isolate the inhibitor bound in the active site from the bulk solvent. The closed-twisted form further differs from the closed form by a rigid body twist of the two-subunit dimers. The local structural fluctuations of SAHH are analyzed by performing block normal mode analysis of the tetrameric enzyme in its three forms. For the open form, we find that the four lowest-frequency normal modes, corresponding to the collective motions of the protein with the largest amplitudes, are essentially combinations of the hinge bending deformations of the individual subunits. Thus, the mechanical properties of the open structure of SAHH lead to the presence of structural fluctuations in the direction of the open-to-closed conformational transition. A candidate for such a motion has been observed in previous fluorescence depolarization studies of the enzyme. Both structural and normal mode analyses indicate that residues 180-190 and 350-356 form hinge regions, connecting large domains which tend to move as rigid bodies in response to interactions with substrate, intermediates, and the product of the enzymatic reactions. We propose that these hinge regions play a crucial role in the enzymatic mechanism of SAHH. In contrast to the open form, normal mode calculations for the closed conformations show strong coupling of the hinge bending motions of the individual subunits to each other and to other low-frequency vibrations. Thus, information about structural changes related to reaction progress in one active site may be mechanically transmitted to other subunits of the protein, explaining the cooperativity found in the enzyme kinetics.     

Refinement of docked protein-ligand and protein-DNA structures using low frequency normal mode amplitude optimization

Prediction of structural changes resulting from complex formation, both in ligands and receptors, is an important and unsolved problem in structural biology. In this work, we use all-atom normal modes calculated with the Elastic Network Model as a basis set to model structural flexibility during formation of macromolecular complexes and refine the non-bonded intermolecular energy between the two partners (protein-ligand or protein-DNA) along 5-10 of the lowest frequency normal mode directions. The method handles motions unrelated to the docking transparently by first applying the modes that improve non-bonded energy most and optionally restraining amplitudes; in addition, the method can correct small errors in the ligand position when the first six rigid-body modes are switched on. For a test set of six protein receptors that show an open-to-close transition when binding small ligands, our refinement scheme reduces the protein coordinate cRMS by 0.3-3.2 A. For two test cases of DNA structures interacting with proteins, the program correctly refines the docked B-DNA starting form into the expected bent DNA, reducing the DNA cRMS from 8.4 to 4.8 A and from 8.7 to 5.4 A, respectively. A public web server implementation of the refinement method is available at http://lorentz.immstr.pasteur.fr.     

Conformational Equilibrium of CDK/Cyclin Complexes by Molecular Dynamics with Excited Normal Modes

Cyclin-dependent kinases (CDKs) and their associated regulatory cyclins are central for timely regulation of cell-cycle progression. They constitute attractive pharmacological targets for development of anticancer therapeutics, since they are frequently deregulated in human cancers and contribute to sustained, uncontrolled tumor proliferation. Characterization of their structural/dynamic features is essential to gain in-depth insight into structure-activity relationships. In addition, the identification of druggable pockets or key intermediate conformations yields potential targets for the development of novel classes of inhibitors. Structural studies of CDK2/cyclin A have provided a wealth of information concerning monomeric/heterodimeric forms of this kinase. There is, however, much less structural information for other CDK/cyclin complexes, including CDK4/cyclin D1, which displays an alternative (open) position of the cyclin partner relative to CDK, contrasting with the closed CDK2/cyclin A conformation. In this study, we carried out normal-mode analysis and enhanced sampling simulations with our recently developed method, molecular dynamics with excited normal modes, to understand the conformational equilibrium on these complexes. Interestingly, the lowest-frequency normal mode computed for each complex described the transition between the open and closed conformations. Exploration of these motions with an explicit-solvent representation using molecular dynamics with excited normal modes confirmed that the closed conformation is the most stable for the CDK2/cyclin A complex, in agreement with their experimentally available structures. On the other hand, we clearly show that an open↔closed equilibrium may exist in CDK4/cyclin D1, with closed conformations resembling that captured for CDK2/cyclin A. Such conformational preferences may result from the distinct distributions of frustrated contacts in each complex. Using the same approach, the putative roles of the Thr¹⁶⁰ phosphoryl group and the T-loop conformation were investigated. These results provide a dynamic view of CDKs revealing intermediate conformations not yet characterized for CDK members other than CDK2, which will be useful for the design of inhibitors targeting critical conformational transitions.     

The effect of physicochemical factors on the activity of DNA-dependent DNA polymerases in Acholeplasma laidlawii PG-8]

The biological and physico-chemical properties of DNA-dependent DNA-polymerases of Acholeplasma laidlawii PG-8 have been studied. The optimal parameters of maximal enzymatic activity are determined. It is stated that N-ethylmaleimide in concentration of 1 mM activated DNA-polymerase I by 52%, whereas DNA-polymerase II with reagent concentration of 0.5 mM demonstrated the peak of activity exceeding the control only by 10%. Spermidine in concentration of 1.5 mM for the first form of DNA-polymerase and 0.15 mM-for the second one increased the ability of both forms of polymerases to synthesize DNA by 10%. Aphidicolin added to the reaction medium up to concentration of 10 mg/ml decreased activity of forms I and II of enzymes by 83 and 68%, respectively. The presence of 0.6 mM of EDTA in the medium also negatively affected the activity of polymerases inhibiting it by 83% in form I and by 77%-in form II.     

Conformational transition in the substrate binding domain of β-secretase exploited by NMA and its implication in inhibitor recognition: BACE1-myricetin a case study

BACE1 is a key protease involved in the proteolysis of amyloid precursor protein (APP) that generates a toxic peptide amyloid beta (Aβ), a pathological feature of Alzheimer's disease (AD). The enzyme is believed to possess an open and a closed conformation that corresponds to its free and inhibitor-bound form respectively. Here, we study the dynamic transition of BACE1 employing normal mode analysis (NMA) using a simplified elastic network model (ENM). Estimation of the catalytic cavity volume on the structures of BACE1 encoded by the lowest frequency normal mode reveals the dynamical transition of the enzyme from the open to the closed conformer. Detailed analysis reveals that concerted movement of different loop segments in the active site of the protein, namely flap regions, 10s loop, A loop and F loop, squeeze the catalytic cavity between the N-terminal and C-terminal lobe of the substrate binding domain of BACE1. We also propose that the NMA encoded multiple receptor conformations (MRC) of BACE1 elucidate the pharmacophoric feature necessary to inhibit the enzyme by a polyphenol, myricetin. van der Waals interaction is found to be the main driving force that guides the ligand induced conformational switching to the closed conformer. We suggest that NMA derived MRC of BACE1 is an efficient way to treat the receptor flexibility in docking and thus can be further applied in virtual screening and structure based drug design.     

Mismatched dNTP incorporation by DNA polymerase beta does not proceed via globally different conformational pathways

Understanding how DNA polymerases control fidelity requires elucidation of the mechanisms of matched and mismatched dNTP incorporations. Little is known about the latter because mismatched complexes do not crystallize readily. In this report, we employed small-angle X-ray scattering (SAXS) and structural modeling to probe the conformations of different intermediate states of mammalian DNA polymerase β (Pol β) in its wild-type and an error-prone variant, I260Q. Our structural results indicate that the mismatched ternary complex lies in-between the open and the closed forms, but more closely resembles the open form for WT and the closed form for I260Q. On the basis of molecular modeling, this over-stabilization of mismatched ternary complex of I260Q is likely caused by formation of a hydrogen bonding network between the side chains of Gln²⁶⁰, Tyr²⁹⁶, Glu²⁹⁵ and Arg²⁵⁸, freeing up Asp¹⁹² to coordinate MgdNTP. These results argue against recent reports suggesting that mismatched dNTP incorporations follow a conformational path distinctly different from that of matched dNTP incorporation, or that its conformational closing is a major contributor to fidelity.     

Discrimination between closed and open forms of lipases using electrophoretic techniques

The enhanced catalytic activity of lipases is often associated with structural changes. The three-dimensional (3D) structures showed that the covalently inhibited lipases exist under their open conformations, in contrast to their native closed forms. We studied the inhibition of various lipases--human and dog gastric lipases, human pancreatic lipase, and Humicola lanuginosa lipase--by the octyl-undecyl phosphonate inhibitor, and we measured the subsequent modifications of their respective electrophoretic mobility. Furthermore, the experimental values of the isoelectric points found for the native (closed) and inhibited (open) lipases are in agreement with theoretical calculations based on the electrostatic potential. We concluded that there is a significant difference in the isoelectric points between the closed (native) and open (inhibited) conformations of the four lipases investigated. Thus, analysis of the electrophoretic pattern is proposed as an easy experimental tool to differentiate between a closed and an open form of a given lipase.     

Conformational transition pathways explored by Monte Carlo simulation integrated with collective modes

Conformational transitions between open/closed or free/bound states in proteins possess functional importance. We propose a technique in which the collective modes obtained from an anisotropic network model (ANM) are used in conjunction with a Monte Carlo (MC) simulation approach, to investigate conformational transition pathways and pathway intermediates. The ANM-MC technique is applied to adenylate kinase (AK) and hemoglobin. The iterative method, in which normal modes are continuously updated during the simulation, proves successful in accomplishing the transition between open-closed conformations of AK and tense-relaxed forms of hemoglobin (C^α− root mean square deviations between two end structures of 7.13 Å and 3.55 Å, respectively). Target conformations are reached by root mean-square deviations of 2.27 Å and 1.90 Å for AK and hemoglobin, respectively. The intermediate conformations overlap with crystal structures from the AK family within a 3.0-Å root mean-square deviation. In the case of hemoglobin, the transition of tense-to-relaxed passes through the relaxed state. In both cases, the lowest-frequency modes are effective during transitions. The targeted Monte Carlo approach is used without the application of collective modes. Both the ANM-MC and targeted Monte Carlo techniques can explore sequences of events in transition pathways with an efficient yet realistic conformational search.     

Ligand-induced conformational change of a protein reproduced by a linear combination of displacement vectors obtained from normal mode analysis

The conformational change of a protein upon ligand binding was examined by normal mode analysis (NMA) based on an elastic-network model (ENM) for a full-atom system using dihedral angles as independent variables. Specifically, we investigated the extent to which conformational change vectors of atoms from an apo form to a holo form of a protein can be represented by a linear combination of the displacement vectors of atoms in the apo form calculated for the lowest-frequency m normal modes (m=1, 2,…, 20). In this analysis, the latter vectors were best fitted to the former ones by the least-squares method. Twenty-two paired proteins in the holo and apo forms, including three dimer pairs, were examined. The results showed that, in most cases, the conformational change vectors were reproduced well by a linear combination of the displacement vectors of a small number of low-frequency normal modes. The conformational change around an active site was reproduced as well as the entire conformational change, except for some proteins that only undergo significant conformational changes around active sites. The weighting factors for 20 normal modes optimized by the least-squares fitting characterize the conformational changes upon ligand binding for these proteins. The conformational changes sampled around the apo form of a protein by the linear combination of the displacement vectors obtained by ENM-based NMA may help solve the flexible-docking problem of a protein with another molecule because the results presented herein suggest that they have a relatively high probability of being involved in an actual conformational change.     

The Artemis:DNA-PKcs endonuclease cleaves DNA loops, flaps, and gaps

In eukaryotic cells, nonhomologous DNA end joining (NHEJ) is a major pathway for repair of double-strand DNA breaks (DSBs). Artemis and the 469kDa DNA-dependent protein kinase (DNA-PKcs) together form a key nuclease for NHEJ in vertebrate organisms. The structure-specific endonucleolytic activity of Artemis is activated by binding to and phosphorylation by DNA-PKcs. We tested various DNA structures in order to understand the range of structural features that are recognized by the Artemis:DNA-PKcs complex. We find that all tested substrates that contain single-to-double-strand transitions can be cleaved by the Artemis:DNA-PKcs complex near the transition region. The cleaved substrates include heterologous loops, stem-loops, flaps, and gapped substrates. Such versatile activity on single-/double-strand transition regions is important in understanding how reconstituted NHEJ systems that lack DNA polymerases can join incompatible DNA ends and yet preserve 3' overhangs. Additionally, the flexibility of the Artemis:DNA-PKcs nuclease may be important in removing secondary structures that hinder processing of DNA ends during NHEJ.     

Hinge-bending motion in citrate synthase arising from normal mode calculations

A normal mode analysis of the closed form of dimeric citrate synthase has been performed. The largest-amplitude collective motion predicted by this method compares well with the crystallographically observed hinge-bending motion. Such a result supports those obtained previously in the case of hinge-bending motions of smaller systems, such as lysozyme or hexokinase. Taken together, all these results suggest that low-frequency normal modes may become useful for determining a first approximation of the conformational path between the closed and open forms of these proteins. © 1995 Wiley-Liss, Inc.     

Toward decrypting the allosteric mechanism of the ryanodine receptor based on coarse‐grained structural and dynamic modeling

The ryanodine receptors (RyRs) are a family of calcium (Ca) channels that regulate Ca release by undergoing a closed‐to‐open gating transition in response to action potential or Ca binding. The allosteric mechanism of RyRs gating, which is activated/regulated by ligand/protein binding >200 Å away from the channel gate, remains elusive for the lack of high‐resolution structures. Recent solution of the closed‐form structures of the RyR1 isoform by cryo‐electron microscopy has paved the way for detailed structure‐driven studies of RyRs functions. Toward elucidating the allosteric mechanism of RyRs gating, we performed coarse‐grained modeling based on the newly solved closed‐form structures of RyR1. Our normal mode analysis captured a key mode of collective motions dominating the observed structural variations in RyR1, which features large outward and downward movements of the peripheral domains with the channel remaining closed, and involves hotspot residues that overlap well with key functional sites and disease mutations. In particular, we found a key interaction between a peripheral domain and the Ca‐binding EF hand domain, which may allow for direct coupling of Ca binding to the collective motions as captured by the above mode. This key mode was robustly reproduced by the normal mode analysis of the other two closed‐form structures of RyR1 solved independently. To elucidate the closed‐to‐open conformational changes in RyR1 with amino‐acid level of details, we flexibly fitted the closed‐form structures of RyR1 into a 10‐Å cryo‐electron microscopy map of the open state. We observed extensive structural changes involving the peripheral domains and the central domains, resulting in the channel pore opening. In sum, our findings have offered unprecedented structural and dynamic insights to the allosteric mechanism of RyR1 via modulation of the key collective motions involved in RyR1 gating. The predicted hotspot residues and open‐form conformation of RyR1 will guide future mutational and functional studies. Proteins 2015; 83:2307–2318. © 2015 Wiley Periodicals, Inc.     

DNA-dependent RNA polymerases from Drosophila melanogaster adults: Isolation and partial characterization

In preparation for the isolation and biochemical characterization of putative RNA polymerase mutants, DNA-dependent RNA polymerases of Drosophila melanogaster adults were isolated and partially characterized. Approximately 70% of the female adult RNA polymerase is located in ovaries. Multiple forms of ovarian RNA polymerases I and II are separable by DEAE-Sephadex chromatography. The two forms of RNA polymerase II differ in ammonium sulfate optima. RNA polymerase IIA is more active with double-stranded DNA as template, whereas RNA polymerase IIB transcribes single-stranded DNA most efficiently. Rechromatography of RNA polymerase IIA on DEAE-Sephadex results in the loss of ability of this form to transcribed double-stranded DNA most efficiently. Ovariectomized carcasses have two forms of RNA polymerase I and one form of RNA polymerase II and each transcribes single-stranded DNA most efficiently. As judged by gel filtration chromatography, female adult extracts have forms of RNA polymerase II that differ in molecular weight and template preference.Supported by Grants GM23456 from the NIH and 11259 from the City University Research Foundation.     

Novel conformational states of peptide deformylase from pathogenic bacterium Leptospira interrogans: implications for population shift

Peptide deformylase is an attractive target for developing novel antibiotics. Previous studies at pH 3.0 showed peptide deformylase from Leptospira interrogans (LiPDF) exists as a dimer in which one monomer is in a closed form and the other is in an open form, with different conformations of the CD-loop controlling the entrance to the active pocket. Here we present structures of LiPDF at its active pH range. LiPDF forms a similar dimer at pH values 6.5-8.0 as it does at pH 3.0. Interestingly, both of the monomers are almost in the same closed form as that observed at pH 3.0. However, when the enzyme is complexed with the natural inhibitor actinotin, the conformation of the CD-loop is half-open. Two pairs of Arg109-mediated cation-pi interactions, as well as hydrogen bonds, have been identified to stabilize the different CD-loop conformations. These results indicate that LiPDF may be found in different structural states, a feature that has never before been observed in the peptide deformylase family. Based on our results, a novel substrate binding model, featured by an equilibrium between the closed and the open forms, is proposed. Our results present crystallographic evidence supporting population shift theory, which is distinguished from the conventional lock-and-key or induced-fit models. These results not only facilitate the development of peptide deformylase-targeted drugs but also provide structural insights into the mechanism of an unusual type of protein binding event.