Structural details of ribonuclease H from Escherichia coli as refined to an atomic resolution

A comparative analysis of mycobacterial ribonucleases: Towards a therapeutic novel drug target

Bacterial responses to continuously changing environments are addressed through modulation of gene expression at the level of transcription initiation, RNA processing and/or decay. Ribonucleases (RNases) are hydrolytic or phosphorolytic enzymes involved in a majority of RNA metabolism reactions. RNases play a crucial role in RNA degradation, either independently or in collaboration with various trans-acting regulatory factors. The genus Mycobacterium consists of five subgenera: Mycobacteroides, Mycolicibacterium, Mycobacterium, Mycolicibacter and Mycolicibacillus, which include 63 fully sequenced species (pathogenic/non-pathogenic) to date. These include 13 different RNases, among which 5 are exonucleases (RNase PH, PNPase, RNase D, nano-RNases and RNase AS) and 8 are endonucleases (RNase J, RNase H, RNase P, RNase III, RNase BN, RNase Z, RNase G and RNase E), although RNase J and RNase BN were later identified to have exoribonuclease functions also. Here, we provide a detailed comparative insight into the Escherichia coli and mycobacterial RNases with respect to their types, phylogeny, structure, function, regulation and mechanism of action, with the main emphasis on RNase E. Among these 13 different mycobacterial RNases, 10 are essential for cell survival and have diverse structures hence, they are promising drug targets. RNase E is also an essential endonuclease that is abundant in many bacteria, forms an RNA degradosome complex that controls central RNA processing/degradation and has a conserved 5' sensor domain/DNase-I like region in its RNase domain. The essential mycobacterial RNases especially RNase E provide a potential repertoire of drug targets that can be exploited for inhibitor/modulator screening against many deadly mycobacterial diseases.

Parsing Dynamics of Protein Backbone NH and Side-Chain Methyl Groups using Molecular Dynamics Simulations

Experimental NMR spectroscopy and theoretical molecular dynamics (MD) simulations provide complementary insights into protein conformational dynamics and hence into biological function. The present work describes an extensive set of backbone NH and side-chain methyl group generalized order parameters for the Escherichia coli ribonuclease HI (RNH) enzyme derived from 2-μs microsecond MD simulations using the OPLS4 and AMBER-FF19SB force fields. The simulated generalized order parameters are compared with values derived from NMR 15N and 13CH2D spin relaxation measurements. The squares of the generalized order parameters, S2 for the N-H bond vector and Saxis2 for the methyl group symmetry axis, characterize the equilibrium distribution of vector orientations in a molecular frame of reference. Optimal agreement between simulated and experimental results was obtained by averaging S2 or Saxis2 calculated by dividing the simulated trajectories into 50 ns blocks (∼five times the rotational diffusion correlation time for RNH). With this procedure, the median absolute deviations (MAD) between experimental and simulated values of S2 and Saxis2 are 0.030 (NH) and 0.061 (CH3) for OPLS4 and 0.041 (NH) and 0.078 (CH3) for AMBER-FF19SB. The MAD between OPLS4 and AMBER-FF19SB are 0.021 (NH) and 0.072 (CH3). The generalized order parameters for the methyl group symmetry axis can be decomposed into contributions from backbone fluctuations, between-rotamer dihedral angle transitions, and within-rotamer dihedral angle fluctuations. Analysis of the simulation trajectories shows that (i) backbone and side chain conformational fluctuations exhibit little correlation and that (ii) fluctuations within rotamers are limited and highly uniform with values that depend on the number of dihedral angles considered. Low values of Saxis2, indicative of enhanced side-chain flexibility, result from between-rotamer transitions that can be enhanced by increased local backbone flexibility.

The importance of input sequence set to consensus-derived proteins and their relationship to reconstructed ancestral proteins

A protein sequence encodes its energy landscape-all the accessible conformations, energetics, and dynamics. The evolutionary relationship between sequence and landscape can be probed phylogenetically by compiling a multiple sequence alignment of homologous sequences and generating common ancestors via Ancestral Sequence Reconstruction or a consensus protein containing the most common amino acid at each position. Both ancestral and consensus proteins are often more stable than their extant homologs-questioning the differences between them and suggesting that both approaches serve as general methods to engineer thermostability. We used the Ribonuclease H family to compare these approaches and evaluate how the evolutionary relationship of the input sequences affects the properties of the resulting consensus protein. While the consensus protein derived from our full Ribonuclease H sequence alignment is structured and active, it neither shows properties of a well-folded protein nor has enhanced stability. In contrast, the consensus protein derived from a phylogenetically-restricted set of sequences is significantly more stable and cooperatively folded, suggesting that cooperativity may be encoded by different mechanisms in separate clades and lost when too many diverse clades are combined to generate a consensus protein. To explore this, we compared pairwise covariance scores using a Potts formalism as well as higher-order sequence correlations using singular value decomposition (SVD). We find the SVD coordinates of a stable consensus sequence are close to coordinates of the analogous ancestor sequence and its descendants, whereas the unstable consensus sequences are outliers in SVD space.

Development of Inhibitory Compounds for Metallo-beta-lactamase through Computational Design and Crystallographic Analysis

Metallo-β-lactamases (MBL) deactivate β-lactam antibiotics through a catalytic reaction caused by two zinc ions at the active center. Since MBLs deteriorate a wide range of antibiotics, they are dangerous factors for bacterial multidrug resistance. In this work, organic synthesis, computational design, and crystal structure analysis were performed to obtain potent MBL inhibitors based on a previously identified hit compound. The hit compound comprised 3,4-dihydro-2(1H)-quinolinone linked with a phenyl-ether-methyl group via a thiazole ring. In the first step, the thiazole ring was replaced with a tertiary amine to avoid the planar structure. In the second step, we virtually modified the compound by keeping the quinolinone backbone. Every modified compound was bound to a kind of MBL, imipenemase-1 (IMP-1), and the binding pose was optimized by a molecular mechanics calculation. The binding scores were evaluated for the respective optimized binding poses. Given the predicted binding poses and calculated binding scores, candidate compounds were determined for organic syntheses. The inhibitory activities of the synthesized compounds were measured by an in vitro assay for two kinds of MBLs, IMP-1 and New Delhi metallo-β-lactamase (NDM-1). A quinolinone connected with an amine bound with methyl-phenyl-ether-propyl and cyclohexyl-ethyl showed a 50% inhibitory concentration of 4.8 μM. An X-ray crystal analysis clarified the binding structure of a synthesized compound to IMP-1. The δ-lactam ring of quinolinone was hydrolyzed, and the generated carboxyl group was coordinated with zinc ions. The findings on the chemical structure and binding pose are expected to be a base for developing MBL inhibitors.

Systematic enhancement of protein crystallization efficiency by bulk lysine-to-arginine (KR) substitution

Structural genomics consortia established that protein crystallization is the primary obstacle to structure determination using x-ray crystallography. We previously demonstrated that crystallization propensity is systematically related to primary sequence, and we subsequently performed computational analyses showing that arginine is the most overrepresented amino acid in crystal-packing interfaces in the Protein Data Bank. Given the similar physicochemical characteristics of arginine and lysine, we hypothesized that multiple lysine-to-arginine (KR) substitutions should improve crystallization. To test this hypothesis, we developed software that ranks lysine sites in a target protein based on the redundancy-corrected KR substitution frequency in homologs. This software can be run interactively on the worldwide web at https://www.pxengineering.org/. We demonstrate that three unrelated single-domain proteins can tolerate 5-11 KR substitutions with at most minor destabilization, and, for two of these three proteins, the construct with the largest number of KR substitutions exhibits significantly enhanced crystallization propensity. This approach rapidly produced a 1.9 Å crystal structure of a human protein domain refractory to crystallization with its native sequence. Structures from Bulk KR-substituted domains show the engineered arginine residues frequently make hydrogen-bonds across crystal-packing interfaces. We thus demonstrate that Bulk KR substitution represents a rational and efficient method for probabilistic engineering of protein surface properties to improve crystallization.

Selection pressures on evolution of ribonuclease H explored with rigorous free-energy-based design

Understanding natural protein evolution and designing novel proteins are motivating interest in development of high-throughput methods to explore large sequence spaces. In this work, we demonstrate the application of multisite λ dynamics (MSλD), a rigorous free energy simulation method, and chemical denaturation experiments to quantify evolutionary selection pressure from sequence-stability relationships and to address questions of design. This study examines a mesophilic phylogenetic clade of ribonuclease H (RNase H), furthering its extensive characterization in earlier studies, focusing on E. coli RNase H (ecRNH) and a more stable consensus sequence (AncCcons) differing at 15 positions. The stabilities of 32,768 chimeras between these two sequences were computed using the MSλD framework. The most stable and least stable chimeras were predicted and tested along with several other sequences, revealing a designed chimera with approximately the same stability increase as AncCcons, but requiring only half the mutations. Comparing the computed stabilities with experiment for 12 sequences reveals a Pearson correlation of 0.86 and root mean squared error of 1.18 kcal/mol, an unprecedented level of accuracy well beyond less rigorous computational design methods. We then quantified selection pressure using a simple evolutionary model in which sequences are selected according to the Boltzmann factor of their stability. Selection temperatures from 110 to 168 K are estimated in three ways by comparing experimental and computational results to evolutionary models. These estimates indicate selection pressure is high, which has implications for evolutionary dynamics and for the accuracy required for design, and suggests accurate high-throughput computational methods like MSλD may enable more effective protein design.

Accurately Predicting Protein pKa Values Using Nonequilibrium Alchemy

The stability, solubility, and function of a protein depend on both its net charge and the protonation states of its individual residues. pKa is a measure of the tendency for a given residue to (de)protonate at a specific pH. Although pKa values can be resolved experimentally, theory and computation provide a compelling alternative. To this end, we assess the applicability of a nonequilibrium (NEQ) alchemical free energy method to the problem of pKa prediction. On a data set of 144 residues that span 13 proteins, we report an average unsigned error of 0.77 ± 0.09, 0.69 ± 0.09, and 0.52 ± 0.04 pK for aspartate, glutamate, and lysine, respectively. This is comparable to current state-of-the-art predictors and the accuracy recently reached using free energy perturbation methods (e.g., FEP+). Moreover, we demonstrate that our open-source, pmx-based approach can accurately resolve the pKa values of coupled residues and observe a substantial performance disparity associated with the lysine partial charges in Amber14SB/Amber99SB*-ILDN, for which an underused fix already exists.

The importance of input sequence set to consensus-derived proteins and their relationship to reconstructed ancestral proteins

A protein sequence encodes its energy landscape - all the accessible conformations, energetics, and dynamics. The evolutionary relationship between sequence and landscape can be probed phylogenetically by compiling a multiple sequence alignment of homologous sequences and generating common ancestors via Ancestral Sequence Reconstruction or a consensus protein containing the most common amino acid at each position. Both ancestral and consensus proteins are often more stable than their extant homologs - questioning the differences and suggesting that both approaches serve as general methods to engineer thermostability. We used the Ribonuclease H family to compare these approaches and evaluate how the evolutionary relationship of the input sequences affects the properties of the resulting consensus protein. While the overall consensus protein is structured and active, it neither shows properties of a well-folded protein nor has enhanced stability. In contrast, the consensus protein derived from a phylogenetically-restricted region is significantly more stable and cooperatively folded, suggesting that cooperativity may be encoded by different mechanisms in separate clades and lost when too many diverse clades are combined to generate a consensus protein. To explore this, we compared pairwise covariance scores using a Potts formalism as well as higher-order couplings using singular value decomposition (SVD). We find the SVD coordinates of a stable consensus sequence are close to coordinates of the analogous ancestor sequence and its descendants, whereas the unstable consensus sequences are outliers in SVD space.

Dispersion Effect of Molecular Crowding on Ligand-Protein Surface Binding Sites of Escherichia coli RNase HI

The molecular crowding effect on ligand-protein interactions, which plays several crucial roles in life processes, has been investigated using various models by adding crowding agents to mimic the intracellular environment. Several studies evaluating this effect have focused on the ligand-protein binding reaction of well-structured binding sites with rigid conformations. However, the crowding effect on flexible binding sites is not well-understood, especially in terms of the conformations. In this work, to elucidate the detailed molecular mechanism underlying the ligand-protein interactions with flexible binding sites on a protein surface, we studied the interaction between the basic protrusion of Escherichia coli ribonuclease HI (RNase HI) and 8-anilinonaphthalene-1-sulfonic acid (ANS). The RNase HI concentration-dependent measurement of ANS fluorescence combined with the multivariate analysis and the fluorescence vibronic structure analysis revealed an increase in the heterogeneous species with an increase in the protein concentration, which is a different behavior from that of proteins with rigid binding sites. This result indicates that ANS molecules bind to the additional binding sites because of the destabilization of the main sites by the excluded volume effect in a crowded environment. The fluorescence vibronic structure analysis yields a detailed molecular picture, indicating that the main species of ANS can have a distorted structure. On the other hand, some ANS molecules move to the minor binding sites of a different microenvironment to secure a stabilized structure. These spectroscopic analyses may show a hypothesis, suggesting that the decrease in the ΔG difference between the main and minor sites due to destabilization of the main binding site could lower the potential barrier between them, inducing the dispersion of binding pathways.

Development of Force Field Parameters for the Simulation of Single- and Double-Stranded DNA Molecules and DNA-Protein Complexes

Although molecular dynamics (MD) simulations have been used extensively to study the structural dynamics of proteins, the role of MD simulation in studies of nucleic acid based systems has been more limited. One contributing factor to this disparity is the historically lower level of accuracy of the physical models used in such simulations to describe interactions involving nucleic acids. By modifying nonbonded and torsion parameters of a force field from the Amber family of models, we recently developed force field parameters for RNA that achieve a level of accuracy comparable to that of state-of-the-art protein force fields. Here we report force field parameters for DNA, which we developed by transferring nonbonded parameters from our recently reported RNA force field and making subsequent adjustments to torsion parameters. We have also modified the backbone charges in both the RNA and DNA parameter sets to make the treatment of electrostatics compatible with our recently developed variant of the Amber protein and ion force field. We name the force field resulting from the union of these three parameter sets (the new DNA parameters, the revised RNA parameters, and the existing protein and ion parameters) DES-Amber. Extensive testing of DES-Amber indicates that it can describe the thermal stability and conformational flexibility of single- and double-stranded DNA systems with a level of accuracy comparable to or, especially for disordered systems, exceeding that of state-of-the-art nucleic acid force fields. Finally, we show that, in certain favorable cases, DES-Amber can be used for long-timescale simulations of protein–nucleic acid complexes.

Comparisons of Ribonuclease HI Homologs and Mutants Uncover a Multistate Model for Substrate Recognition

Ribonuclease HI (RNHI) nonspecifically cleaves the RNA strand of RNA:DNA hybrid duplexes in a myriad of biological processes. Several RNHI homologs contain an extended domain, termed the handle region, which is critical to substrate binding. Nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations have suggested a kinetic model in which the handle region can exist in open (substrate-binding competent) or closed (substrate-binding incompetent) states in homologs containing arginine or lysine at position 88 (using sequence numbering of E. coli RNHI), while the handle region populates states intermediate between the open and closed conformers in homologs with asparagine at residue 88 [Stafford, K. A., et al., PLoS Comput. Biol. 2013, 9, 1-10]. NMR parameters characterizing handle region dynamics are highly correlated with enzymatic activity for RNHI homologs with two-state (open/closed) handle regions [Martin, J. A., et al., Biochemistry 2020, 59, 3201-3205]. The work presented herein shows that homologs containing asparagine 88 display distinct structural features compared with their counterparts containing arginine or lysine 88. Comparisons of RNHI homologs and site-directed mutants with asparagine 88 support a kinetic model for handle region dynamics that includes 12 unique transitions between eight conformations. Overall, these findings present an example of the structure-function relationships of enzymes and spotlight the use of NMR spectroscopy and MD simulations in uncovering fine details of conformational preferences.

Pivotal role of a conserved histidine in Escherichia coli ribonuclease HI as proposed by X-ray crystallography

The ribonuclease (RNase) H family of enzymes catalyze the specific cleavage of RNA strands of RNA/DNA hybrid duplexes and play an important role in DNA replication and repair. Since the first report of the crystal structure of RNase HI, its catalytic mechanisms, which require metal ions, have been discussed based on numerous structural and functional analyses, including X-ray crystallography. In contrast, the function of the conserved histidine residue (His124 in Escherichia coli) in the flexible loop around the active site remains poorly understood, although an important role was suggested by NMR analyses. Here, novel high-resolution X-ray crystal structures of E. coli RNase HI are described, with a particular focus on the interactions of divalent cations with His124 oriented towards the active site. The enzyme-Mg2+ complex contains two metal ions in the active site, one of which has previously been observed. The second ion lies alongside the first and binds to His124 in an octahedral coordination scheme. In the enzyme-Zn2+ complex a single metal ion was found to bind to the active site, showing a tetrahedral coordination geometry with the surrounding atoms, including His124. These results provide structural evidence that His124 plays a crucial role in the catalytic activity of RNase HI by interacting weakly and transiently with metal ions in the catalytic center.

Identification of the Inhibitory Compounds for Metallo-β-lactamases and Structural Analysis of the Binding Modes

Metallo-β-lactamases (MBLs) are significant threats to humans because they deteriorate many kinds of β-lactam antibiotics and are key enzymes responsible for multi-drug resistance of bacterial pathogens. As a result of in vitro screening, two compounds were identified as potent inhibitors of two kinds of MBLs: imipenemase (IMP-1) and New Delhi metallo-β-lactamase (NDM-1). The binding structure of one of the identified compounds was clarified by an X-ray crystal analysis in complex with IMP-1, in which two possible binding poses were observed. Molecular dynamics (MD) simulations were performed by building two calculation models from the respective binding poses. The compound was stably bound to the catalytic site during the simulation in one pose. The binding model between NDM-1 and the compound was constructed for MD simulation. Calculation results for NDM-1 were similar to those of IMP-1. The simulation suggested that the binding of the identified inhibitory compound was also durable in the catalytic site of NDM-1. The compound will be a sound basis for the development of the inhibitors for MBLs.

Bacterial ribonucleases and their roles in RNA metabolism

Ribonucleases (RNases) are mediators in most reactions of RNA metabolism. In recent years, there has been a surge of new information about RNases and the roles they play in cell physiology. In this review, a detailed description of bacterial RNases is presented, focusing primarily on those from Escherichia coli and Bacillus subtilis, the model Gram-negative and Gram-positive organisms, from which most of our current knowledge has been derived. Information from other organisms is also included, where relevant. In an extensive catalog of the known bacterial RNases, their structure, mechanism of action, physiological roles, genetics, and possible regulation are described. The RNase complement of E. coli and B. subtilis is compared, emphasizing the similarities, but especially the differences, between the two. Included are figures showing the three major RNA metabolic pathways in E. coli and B. subtilis and highlighting specific steps in each of the pathways catalyzed by the different RNases. This compilation of the currently available knowledge about bacterial RNases will be a useful tool for workers in the RNA field and for others interested in learning about this area.

Tracing a protein's folding pathway over evolutionary time using ancestral sequence reconstruction and hydrogen exchange

The conformations populated during protein folding have been studied for decades; yet, their evolutionary importance remains largely unexplored. Ancestral sequence reconstruction allows access to proteins across evolutionary time, and new methods such as pulsed-labeling hydrogen exchange coupled with mass spectrometry allow determination of folding intermediate structures at near amino-acid resolution. Here, we combine these techniques to monitor the folding of the ribonuclease H family along the evolutionary lineages of T. thermophilus and E. coli RNase H. All homologs and ancestral proteins studied populate a similar folding intermediate despite being separated by billions of years of evolution. Even though this conformation is conserved, the pathway leading to it has diverged over evolutionary time, and rational mutations can alter this trajectory. Our results demonstrate that evolutionary processes can affect the energy landscape to preserve or alter specific features of a protein’s folding pathway.

Generalized Born Based Continuous Constant pH Molecular Dynamics in Amber: Implementation, Benchmarking and Analysis

Solution pH plays an important role in structure and dynamics of biomolecular systems; however, pH effects cannot be accurately accounted for in conventional molecular dynamics simulations based on fixed protonation states. Continuous constant pH molecular dynamics (CpHMD) based on the λ-dynamics framework calculates protonation states on the fly during dynamical simulation at a specified pH condition. Here we report the CPU-based implementation of the CpHMD method based on the GBNeck2 generalized Born (GB) implicit-solvent model in the pmemd engine of the Amber molecular dynamics package. The performance of the method was tested using pH replica-exchange titration simulations of Asp, Glu and His side chains in 4 miniproteins and 7 enzymes with experimentally known p K_a's, some of which are significantly shifted from the model values. The added computational cost due to CpHMD titration ranges from 11 to 33% for the data set and scales roughly linearly as the ratio between the titrable sites and number of solute atoms. Comparison of the experimental and calculated p K_a's using 2 ns per replica sampling yielded a mean unsigned error of 0.70, a root-mean-squared error of 0.91, and a linear correlation coefficient of 0.79. Though this level of accuracy is similar to the GBSW-based CpHMD in CHARMM, in contrast to the latter, the current implementation was able to reproduce the experimental orders of the p K_a's of the coupled carboxylic dyads. We quantified the sampling errors, which revealed that prolonged simulation is needed to converge p K_a's of several titratable groups involved in salt-bridge-like interactions or deeply buried in the protein interior. Our benchmark data demonstrate that GBNeck2-CpHMD is an attractive tool for protein p K_a predictions.

Validating Molecular Dynamics Simulations against Experimental Observables in Light of Underlying Conformational Ensembles

Far from the static, idealized conformations deposited into structural databases, proteins are highly dynamic molecules that undergo conformational changes on temporal and spatial scales that may span several orders of magnitude. These conformational changes, often intimately connected to the functional roles that proteins play, may be obscured by traditional biophysical techniques. Over the past 40 years, molecular dynamics (MD) simulations have complemented these techniques by providing the "hidden" atomistic details that underlie protein dynamics. However, there are limitations of the degree to which molecular simulations accurately and quantitatively describe protein motions. Here we show that although four molecular dynamics simulation packages (AMBER, GROMACS, NAMD, and ilmm) reproduced a variety of experimental observables for two different proteins (engrailed homeodomain and RNase H) equally well overall at room temperature, there were subtle differences in the underlying conformational distributions and the extent of conformational sampling obtained. This leads to ambiguity about which results are correct, as experiment cannot always provide the necessary detailed information to distinguish between the underlying conformational ensembles. However, the results with different packages diverged more when considering larger amplitude motion, for example, the thermal unfolding process and conformational states sampled, with some packages failing to allow the protein to unfold at high temperature or providing results at odds with experiment. While most differences between MD simulations performed with different packages are attributed to the force fields themselves, there are many other factors that influence the outcome, including the water model, algorithms that constrain motion, how atomic interactions are handled, and the simulation ensemble employed. Here four different MD packages were tested each using best practices as established by the developers, utilizing three different protein force fields and three different water models. Differences between the simulated protein behavior using two different packages but the same force field, as well as two different packages with different force fields but the same water models and approaches to restraining motion, show how other factors can influence the behavior, and it is incorrect to place all the blame for deviations and errors on force fields or to expect improvements in force fields alone to solve such problems.

FlexPepDock lessons from CAPRI peptide-protein rounds and suggested new criteria for assessment of model quality and utility

CAPRI rounds 28 and 29 included, for the first time, peptide-receptor targets of three different systems, reflecting increased appreciation of the importance of peptide-protein interactions. The CAPRI rounds allowed us to objectively assess the performance of Rosetta FlexPepDock, one of the first protocols to explicitly include peptide flexibility in docking, accounting for peptide conformational changes upon binding. We discuss here successes and challenges in modeling these targets: we obtain top-performing, high-resolution models of the peptide motif for cases with known binding sites but there is a need for better modeling of flanking regions, as well as better selection criteria, in particular for unknown binding sites. These rounds have also provided us the opportunity to reassess the success criteria, to better reflect the quality of a peptide-protein complex model. Using all models submitted to CAPRI, we analyze the correlation between current classification criteria and the ability to retrieve critical interface features, such as hydrogen bonds and hotspots. We find that loosening the backbone (and ligand) RMSD threshold, together with a restriction on the side chain RMSD measure, allows us to improve the selection of high-accuracy models. We also suggest a new measure to assess interface hydrogen bond recovery, which is not assessed by the current CAPRI criteria. Finally, we find that surprisingly much can be learned from rather inaccurate models about binding hotspots, suggesting that the current status of peptide-protein docking methods, as reflected by the submitted CAPRI models, can already have a significant impact on our understanding of protein interactions. Proteins 2017; 85:445-462. © 2016 Wiley Periodicals, Inc.

Packing energetics determine the folding routes of the RNase-H proteins

The RNase-H proteins show a diverse range of folding routes with structurally distinct folding nuclei.

Performance of ZDOCK and IRAD in CAPRI rounds 28-34

We report the performance of our protein-protein docking pipeline, including the ZDOCK rigid-body docking algorithm, on 19 targets in CAPRI rounds 28-34. Following the docking step, we reranked the ZDOCK predictions using the IRAD scoring function, pruned redundant predictions, performed energy landscape analysis, and utilized our interface prediction approach RCF. In addition, we applied constraints to the search space based on biological information that we culled from the literature, which increased the chance of making a correct prediction. For all but two targets we were able to find and apply biological information and we found the information to be highly accurate, indicating that effective incorporation of biological information is an important component for protein-protein docking. Proteins 2017; 85:408-416. © 2016 Wiley Periodicals, Inc.

Reverse Transcriptase-Associated Ribonuclease H Activity as a Target for Antiviral Chemotherapy

The availability of highly purified recombinant enzymes and model heteropolymeric nucleic acid substrates now allows more precise evaluation of the ribonuclease H (RNase H) activity associated with human immunodeficiency virus (HIV) reverse transcriptase. In addition to degrading the RNA–DNA replicative intermediate, this C-terminal domain of around 130 residues supports highly specialized events that cannot be complemented by host-coded enzymes during retrovirus replication. RNase H activity should therefore be considered a plausible candidate for therapeutic intervention. Events during HIV replication requiring precise RNase H-mediated hydrolysis, the methodologies available to study these events, and their potential for therapeutic intervention are reviewed here.

Clusters of isoleucine, leucine, and valine side chains define cores of stability in high-energy states of globular proteins: Sequence determinants of structure and stability

Measurements of protection against exchange of main chain amide hydrogens (NH) with solvent hydrogens in globular proteins have provided remarkable insights into the structures of rare high-energy states that populate their folding free-energy surfaces. Lacking, however, has been a unifying theory that rationalizes these high-energy states in terms of the structures and sequences of their resident proteins. The Branched Aliphatic Side Chain (BASiC) hypothesis has been developed to explain the observed patterns of protection in a pair of TIM barrel proteins. This hypothesis supposes that the side chains of isoleucine, leucine, and valine (ILV) residues often form large hydrophobic clusters that very effectively impede the penetration of water to their underlying hydrogen bond networks and, thereby, enhance the protection against solvent exchange. The linkage between the secondary and tertiary structures enables these ILV clusters to serve as cores of stability in high-energy partially folded states. Statistically significant correlations between the locations of large ILV clusters in native conformations and strong protection against exchange for a variety of motifs reported in the literature support the generality of the BASiC hypothesis. The results also illustrate the necessity to elaborate this simple hypothesis to account for the roles of adjacent hydrocarbon moieties in defining stability cores of partially folded states along folding reaction coordinates.

Exoribonucleases and Endoribonucleases

This review provides a description of the known Escherichia coli ribonucleases (RNases), focusing on their structures, catalytic properties, genes, physiological roles, and possible regulation. Currently, eight E. coli exoribonucleases are known. These are RNases II, R, D, T, PH, BN, polynucleotide phosphorylase (PNPase), and oligoribonuclease (ORNase). Based on sequence analysis and catalytic properties, the eight exoribonucleases have been grouped into four families. These are the RNR family, including RNase II and RNase R; the DEDD family, including RNase D, RNase T, and ORNase; the RBN family, consisting of RNase BN; and the PDX family, including PNPase and RNase PH. Seven well-characterized endoribonucleases are known in E. coli. These are RNases I, III, P, E, G, HI, and HII. Homologues to most of these enzymes are also present in Salmonella. Most of the endoribonucleases cleave RNA in the presence of divalent cations, producing fragments with 3'-hydroxyl and 5'-phosphate termini. RNase H selectively hydrolyzes the RNA strand of RNA?DNA hybrids. Members of the RNase H family are widely distributed among prokaryotic and eukaryotic organisms in three distinct lineages, RNases HI, HII, and HIII. It is likely that E. coli contains additional endoribonucleases that have not yet been characterized. First of all, endonucleolytic activities are needed for certain known processes that cannot be attributed to any of the known enzymes. Second, homologues of known endoribonucleases are present in E. coli. Third, endonucleolytic activities have been observed in cell extracts that have different properties from known enzymes.

(Dis)similarity Index To Compare Correlated Motions in Molecular Simulations

Molecular dynamics (MD) simulations are widely used to complement or guide experimental studies in the characterization of protein dynamics, thanks to improvements in force-field accuracy, along with in the software and hardware to sample the conformational landscape of proteins. Among the different applications of MD simulations, the study of correlated motions is largely employed for different purposes. Several metrics have been developed to describe correlated motions in the MD ensemble, such as methods based on Pearson Correlation or Mutual Information. Cross-correlation analysis of MD trajectories is indeed appealing not only to identify residues characterized by coupled fluctuations in protein structures but also since it can be used to extrapolate motions along directions in which major conformational changes should occur, for example on longer time scales than the ones that are actually simulated. Nevertheless, most of the MD studies employ average correlation maps and mostly in a qualitative way, even when different systems or different replicates of the same system are compared. The broad application of correlation metrics in the analysis of MD simulations, especially for comparative purposes, requires a step forward toward more quantitative and accurate comparisons. We thus here employed a simple but effective index, which is based on a normalized Frobenius norm of the differences between protein correlation maps, to compare correlated motions. We applied this index for a quantitative comparison of correlated motions from MD simulations of seven proteins of different size and fold. We also employed the index to assess the robustness of correlation description when multi-replicate MD simulations of a same system are used, and we compared our index to metrics for comparison of structural ensembles such as Root Mean Square Inner Product and the Bhattacharyya Coefficient.

Autonomously folding protein fragments reveal differences in the energy landscapes of homologous RNases H

An important approach to understanding how a protein sequence encodes its energy landscape is to compare proteins with different sequences that fold to the same general native structure. In this work, we compare E. coli and T. thermophilus homologs of the protein RNase H. Using protein fragments, we create equilibrium mimics of two different potential partially-folded intermediates (I(core) and I(core+1)) hypothesized to be present on the energy landscapes of these two proteins. We observe that both T. thermophilus RNase H (ttRNH) fragments are folded and have distinct stabilities, indicating that both regions are capable of autonomous folding and that both intermediates are present as local minima on the ttRNH energy landscape. In contrast, the two E. coli RNase H (ecRNH) fragments have very similar stabilities, suggesting that the presence of additional residues in the I(core+1) fragment does not affect the folding or structure as compared to I(core). NMR experiments provide additional evidence that only the I(core) intermediate is populated by ecRNH. This is one of the biggest differences that has been observed between the energy landscapes of these two proteins. Additionally, we used a FRET experiment in the background of full-length ttRNH to specifically monitor the formation of the I(core+1) intermediate. We determine that the ttRNH I(core+1) intermediate is likely the intermediate populated prior to the rate-limiting barrier to global folding, in contrast to E. coli RNase H for which I(core) is the folding intermediate. This result provides new insight into the nature of the rate-limiting barrier for the folding of RNase H.

Conformational preferences underlying reduced activity of a thermophilic ribonuclease H

The conformational basis for reduced activity of the thermophilic ribonuclease HI enzyme from Thermus thermophilus, compared to its mesophilic homolog from Escherichia coli, is elucidated using a combination of NMR spectroscopy and molecular dynamics (MD) simulations. Explicit-solvent all-atom MD simulations of the two wild-type proteins and an E. coli mutant in which a glycine residue is inserted after position 80 to mimic the T. thermophilus protein reproduce the differences in conformational dynamics determined from (15)N spin-relaxation NMR spectroscopy of three loop regions that surround the active site and contain functionally important residues: the glycine-rich region, the handle region, and the β5/αE loop. Examination of the MD trajectories indicates that the thermophilic protein samples conformations productive for substrate binding and activity less frequently than the mesophilic enzyme, although these differences may manifest as either increased or decreased relative flexibility of the different regions. Additional MD simulations indicate that mutations increasing activity of the T. thermophilus enzyme at mesophilic temperatures do so by reconfiguring the local environments of the mutated sites to more closely resemble active conformations. Taken together, the results show that both locally increased and decreased flexibility contribute to an overall reduction in activity of T. thermophilus ribonuclease H compared to its mesophilic E. coli homolog.

Non-native structure appears in microseconds during the folding of E. coli RNase H

The folding pathway of Escherichia coli RNase H is one of the best experimentally characterized for any protein. In spite of this, spectroscopic studies have never captured the earliest events. Using continuous-flow microfluidic mixing, we have now observed the first several milliseconds of folding by monitoring the tryptophan fluorescence lifetime (60 μs dead time). Two folding intermediates are observed, the second of which is the previously characterized I(core) millisecond intermediate. The new earlier intermediate is likely on-pathway and appears to have long-range non-native structure, providing a rare example of such non-native structure formation in a folding pathway. The tryptophan fluorescence lifetimes also suggest a deviation from native packing in the second intermediate, I(core). Similar results from a fragment of RNase H demonstrate that only half of the protein is significantly involved in this early structure formation. These studies give us a view of the formation of tertiary structure on the folding pathway, which complements previous hydrogen-exchange studies that monitored only secondary structure and observed sequential native structure formation. Our results provide detailed folding information on both a timescale and a size-scale accessible to all-atom molecular dynamics simulations of protein folding.

Evidence from molecular dynamics simulations of conformational preorganization in the ribonuclease H active site

Ribonuclease H1 (RNase H) enzymes are well-conserved endonucleases that are present in all domains of life and are particularly important in the life cycle of retroviruses as domains within reverse transcriptase. Despite extensive study, especially of the E. coli homolog, the interaction of the highly negatively charged active site with catalytically required magnesium ions remains poorly understood. In this work, we describe molecular dynamics simulations of the E. coli homolog in complex with magnesium ions, as well as simulations of other homologs in their apo states. Collectively, these results suggest that the active site is highly rigid in the apo state of all homologs studied and is conformationally preorganized to favor the binding of a magnesium ion. Notably, representatives of bacterial, eukaryotic, and retroviral RNases H all exhibit similar active-site rigidity, suggesting that this dynamic feature is only subtly modulated by amino acid sequence and is primarily imposed by the distinctive RNase H protein fold.

Side chain dynamics of carboxyl and carbonyl groups in the catalytic function of Escherichia coli ribonuclease H

Many proteins use Asx and Glx (x = n, p, or u) side chains as key functional groups in enzymatic catalysis and molecular recognition. In this study, NMR spin relaxation experiments and molecular dynamics simulations are used to measure the dynamics of the side chain amide and carboxyl groups, (13)C(γ/δ), in Escherichia coli ribonuclease HI (RNase H). Model-free analysis shows that the catalytic residues in RNase H are preorganized on ps-ns time scales via a network of electrostatic interactions. However, chemical exchange line broadening shows that these residues display significant conformational dynamics on μs-ms time scales upon binding of Mg(2+) ions. Two groups of catalytic residues exhibit differential line broadening, implicating distinct reorganizational processes upon binding of metal ions. These results support the "mobile metal ion" hypothesis, which was inferred from structural studies of RNase H.

Thermal adaptation of conformational dynamics in ribonuclease H

The relationship between inherent internal conformational processes and enzymatic activity or thermodynamic stability of proteins has proven difficult to characterize. The study of homologous proteins with differing thermostabilities offers an especially useful approach for understanding the functional aspects of conformational dynamics. In particular, ribonuclease HI (RNase H), an 18 kD globular protein that hydrolyzes the RNA strand of RNA:DNA hybrid substrates, has been extensively studied by NMR spectroscopy to characterize the differences in dynamics between homologs from the mesophilic organism E. coli and the thermophilic organism T. thermophilus. Herein, molecular dynamics simulations are reported for five homologous RNase H proteins of varying thermostabilities and enzymatic activities from organisms of markedly different preferred growth temperatures. For the E. coli and T. thermophilus proteins, strong agreement is obtained between simulated and experimental values for NMR order parameters and for dynamically averaged chemical shifts, suggesting that these simulations can be a productive platform for predicting the effects of individual amino acid residues on dynamic behavior. Analyses of the simulations reveal that a single residue differentiates between two different and otherwise conserved dynamic processes in a region of the protein known to form part of the substrate-binding interface. Additional key residues within these two categories are identified through the temperature-dependence of these conformational processes.

The relationship between enzymatic activity and protein stability has long been a difficult problem in the study of protein biochemistry. Enzymes may undergo structural changes in order to bind substrates, catalyze chemical reactions, and release products, but flexibility often is inversely correlated with thermodynamic stability. Proteins from organisms that are adapted to high temperature can be both more rigid and less active at ambient temperature than their homologs from organisms that grow at lower temperatures. For this reason, studying homologous pairs of proteins from organisms adapted to different thermal environments is a productive way to identify functionally important motions. In this work we perform comparative analyses of molecular dynamics simulations for five ribonuclease H proteins of varying thermal stabilities, isolated from organisms that grow in varying thermal environments. We identify two different mechanisms of motion in a region of the protein that interacts with substrate molecules, suggesting at least two forms of thermal adaptation in this protein family.

Interpreting protein structural dynamics from NMR chemical shifts

In this investigation, semiempirical NMR chemical shift prediction methods are used to evaluate the dynamically averaged values of backbone chemical shifts obtained from unbiased molecular dynamics (MD) simulations of proteins. MD-averaged chemical shift predictions generally improve agreement with experimental values when compared to predictions made from static X-ray structures. Improved chemical shift predictions result from population-weighted sampling of multiple conformational states and from sampling smaller fluctuations within conformational basins. Improved chemical shift predictions also result from discrete changes to conformations observed in X-ray structures, which may result from crystal contacts, and are not always reflective of conformational dynamics in solution. Chemical shifts are sensitive reporters of fluctuations in backbone and side chain torsional angles, and averaged (1)H chemical shifts are particularly sensitive reporters of fluctuations in aromatic ring positions and geometries of hydrogen bonds. In addition, poor predictions of MD-averaged chemical shifts can identify spurious conformations and motions observed in MD simulations that may result from force field deficiencies or insufficient sampling and can also suggest subsets of conformational space that are more consistent with experimental data. These results suggest that the analysis of dynamically averaged NMR chemical shifts from MD simulations can serve as a powerful approach for characterizing protein motions in atomistic detail.

Insights into the structure and activity of prototype foamy virus RNase H

BACKGROUND

RNase H is an endonuclease that hydrolyzes the RNA strand in RNA/DNA hybrids. Retroviral reverse transcriptases harbor a C-terminal RNase H domain whose activity is essential for viral replication. The RNase H degrades the viral genomic RNA after the first DNA strand is synthesized. Here, we report the biophysical and enzymatic properties of the RNase H domain of prototype foamy virus (PFV) as an independently purified protein. Sequence comparisons with other retroviral RNases H indicated that PFV RNase H harbors a basic protrusion, including a basic loop and the so-called C-helix, which was suggested to be important for activity and substrate binding and is absent in the RNase H domain of human immunodeficiency virus. So far, no structure of a retroviral RNase H containing a C-helix is available.

RESULTS

RNase H activity assays demonstrate that the PFV RNase H domain is active, although its activity is about 200-fold reduced as compared to the full length protease-reverse transcriptase enzyme. Fluorescence equilibrium titrations with an RNA/DNA substrate revealed a KD for the RNase H domain in the low micromolar range which is about 4000-fold higher than that of the full-length protease-reverse transcriptase enzyme. Analysis of the RNase H cleavage pattern using a [32P]-labeled substrate indicates that the independent RNase H domain cleaves the substrate non-specifically. The purified RNase H domain exhibits a well defined three-dimensional structure in solution which is stabilized in the presence of Mg2+ ions.

CONCLUSIONS

Our data demonstrate that the independent PFV RNase H domain is structured and active. The presence of the C-helix in PFV RNase H could be confirmed by assigning the protein backbone and calculating the chemical shift index using NMR spectroscopy.

SRLS analysis of 15N spin relaxation from E. coli ribonuclease HI: the tensorial perspective

15N–H relaxation parameters from ribonuclease HI (RNase H), acquired in previous work at magnetic fields of 14.1 and 18.8 T, and at 300 K, are analyzed with the mode-coupling slowly relaxing local structure (SRLS) approach. In accordance with standard theoretical treatments of restricted motions, SRLS approaches N-H bond dynamics from a tensorial perspective. As shown previously, a physically adequate description of this phenomenon has to account for the asymmetry of the local spatial restrictions. So far, we used rhombic local ordering tensors; this is straightforward but computationally demanding. Here, we propose substantiating the asymmetry of the local spatial restrictions in terms of tilted axial local ordering (S) and local diffusion (D2) tensors. Although less straightforward, this description provides physically sound structural and dynamic information and is efficient computationally. We find that the local order parameter, S(0)2, is on average 0.89 (0.84, and may be as small as 0.6) for the secondary structure elements (loops). The main local ordering axis deviates from the C(i-1)α-C(i)α axis by less than 6°. At 300 K, D(2,perpendicular) is virtually the same as the global diffusion rate, D1 = 1.8 × 10(7) s(-1). The correlation time 1/6D(2,parallel) ranges from 3-125 (208-344) ps for the secondary structure elements (loops) and is on average 125 ps for the C-terminal segment. The main local diffusion axis deviates from the N-H bond by less than 2° (10°) for the secondary structure elements (loops). An effective data-fitting protocol, which leads in most cases to unambiguous results with limited uncertainty, has been devised. A physically sound and computationally effective methodology for analyzing 15N relaxation in proteins, that provides a new picture of N–H bond structural dynamics in proteins, has been set forth.

Protein unfolding under force: crack propagation in a network

The mechanical unfolding of a set of 12 proteins with diverse topologies is investigated using an all-atom constraint-based model. Proteins are represented as polypeptides cross-linked by hydrogen bonds, salt bridges, and hydrophobic contacts, each modeled as a harmonic inequality constraint capable of supporting a finite load before breaking. Stereochemically acceptable unfolding pathways are generated by minimally overloading the network in an iterative fashion, analogous to crack propagation in solids. By comparing the pathways to those from molecular dynamics simulations and intermediates identified from experiment, it is demonstrated that the dominant unfolding pathways for 9 of the 12 proteins studied are well described by crack propagation in a network.

Histidine in continuum electrostatics protonation state calculations

A modification to the standard continuum electrostatics approach to calculate protein pK(a)s, which allows for the decoupling of histidine tautomers within a two-state model, is presented. Histidine with four intrinsically coupled protonation states cannot be easily incorporated into a two-state formalism, because the interaction between the two protonatable sites of the imidazole ring is not purely electrostatic. The presented treatment, based on a single approximation of the interrelation between histidine's charge states, allows for a natural separation of the two protonatable sites associated with the imidazole ring as well as the inclusion of all protonation states within the calculation.

Atomistic details of the associative phosphodiester cleavage in human ribonuclease H

During translation of the genetic information of DNA into proteins, mRNA is synthesized by RNA polymerase and after the transcription process degraded by RNase H. The endoribonuclease RNase H is a member of the nucleotidyl-transferase (NT) superfamily and is known to hydrolyze the phosphodiester bonds of RNA which is hybridized to DNA. Retroviral RNase H is part of the viral reverse transcriptase enzyme that is indispensable for the proliferation of retroviruses, such as HIV. Inhibitors of this enzyme could therefore provide new drugs against diseases like AIDS. In our study we investigated the molecular mechanism of RNA cleavage by human RNase H using a comprehensive high level DFT/B3LYP QM/MM theoretical method for the calculation of the stationary points and nudged elastic band (NEB) and free energy calculations to identify the transition state structures, the rate limiting step and the reaction barrier. Our calculations reveal that the catalytic mechanism proceeds in two steps and that the nature of the nucleophile is a water molecule. In the first step, the water attack on the scissile phosphorous is followed by a proton transfer from the water to the O2P oxygen and a trigonal bipyramidal pentacoordinated phosphorane is formed. Subsequently, in the second step the proton is shuttled to the O3' oxygen to generate the product state. During the reaction mechanism two Mg(2+) ions support the formation of a stable associated in-line S(N)2-type phosphorane intermediate. Our calculated energy barrier of 19.3 kcal mol(-1) is in excellent agreement with experimental findings (20.5 kcal mol(-1)). These results may contribute to the clarification and understanding of the RNase H reaction mechanism and of further enzymes from the RNase family.

Identification of residual structure in the unfolded state of ribonuclease H1 from the moderately thermophilic Chlorobium tepidum: comparison with thermophilic and mesophilic homologues

Ribonucleases H from organisms that grow at different temperatures demonstrate a variable change in heat capacity upon unfolding (DeltaC degrees (P)) [Ratcliff, K., et al. (2009) Biochemistry 48, 5890-5898]. This DeltaC degrees (P) has been shown to correlate with a tolerance to higher temperatures and residual structure in the unfolded state of the thermophilic proteins. In the RNase H from Thermus thermophilus, the low DeltaC degrees (P) has been shown to arise from the same region as the folding core of the protein, and mutagenic studies have shown that loss of a hydrophobic residue in this region can disrupt this residual unfolded state structure and result in a return to a more mesophile-like DeltaC degrees (P) [Robic, S., et al. (2002) Protein Sci. 11, 381-389; Robic, S., et al. (2003) Proc. Natl. Acad. Sci. U.S.A. 100, 11345-11349]. To understand further how residual structure in the unfolded state is encoded in the sequences of these thermophilic proteins, we subjected the RNase H from Chlorobium tepidum to similar studies. Analysis of new chimeric proteins reveals that like T. thermophilus RNase H, the folding core of C. tepidum RNase H plays an important role in the unfolded state of this protein. Mutagenesis studies, based on both a computational investigation of the hydrophobic networks in the core region and comparisons with similar studies on T. thermophilus RNase H, identify new residues involved in this residual structure and suggest that the residual structure in the unfolded state of C. tepidum RNase H is more restricted than that of T. thermophilus. We conclude that while the folding core region determines the thermophilic-like behavior of this family of proteins, the residue-specific details vary.

Evolution and thermodynamics of the slow unfolding of hyperstable monomeric proteins

BACKGROUND

The unfolding speed of some hyperthermophilic proteins is dramatically lower than that of their mesostable homologs. Ribonuclease HII from the hyperthermophilic archaeon Thermococcus kodakaraensis (Tk-RNase HII) is stabilized by its remarkably slow unfolding rate, whereas RNase HI from the thermophilic bacterium Thermus thermophilus (Tt-RNase HI) unfolds rapidly, comparable with to that of RNase HI from Escherichia coli (Ec-RNase HI).

RESULTS

To clarify whether the difference in the unfolding rate is due to differences in the types of RNase H or differences in proteins from archaea and bacteria, we examined the equilibrium stability and unfolding reaction of RNases HII from the hyperthermophilic bacteria Thermotoga maritima (Tm-RNase HII) and Aquifex aeolicus (Aa-RNase HII) and RNase HI from the hyperthermophilic archaeon Sulfolobus tokodaii (Sto-RNase HI). These proteins from hyperthermophiles are more stable than Ec-RNase HI over all the temperature ranges examined. The observed unfolding speeds of all hyperstable proteins at the different denaturant concentrations studied are much lower than those of Ec-RNase HI, which is in accordance with the familiar slow unfolding of hyperstable proteins. However, the unfolding rate constants of these RNases H in water are dispersed, and the unfolding rate constant of thermophilic archaeal proteins is lower than that of thermophilic bacterial proteins.

CONCLUSIONS

These results suggest that the nature of slow unfolding of thermophilic proteins is determined by the evolutionary history of the organisms involved. The unfolding rate constants in water are related to the amount of buried hydrophobic residues in the tertiary structure.

Increasing protein stability by improving beta-turns

Our goal was to gain a better understanding of how protein stability can be increased by improving beta-turns. We studied 22 beta-turns in nine proteins with 66-370 residues by replacing other residues with proline and glycine and measuring the stability. These two residues are statistically preferred in some beta-turn positions. We studied: Cold shock protein B (CspB), Histidine-containing phosphocarrier protein, Ubiquitin, Ribonucleases Sa2, Sa3, T1, and HI, Tryptophan synthetase alpha-subunit, and Maltose binding protein. Of the 15 single proline mutations, 11 increased stability (Average = 0.8 +/- 0.3; Range = 0.3-1.5 kcal/mol), and the stabilizing effect of double proline mutants was additive. On the basis of this and our previous work, we conclude that proteins can generally be stabilized by replacing nonproline residues with proline residues at the i + 1 position of Type I and II beta-turns and at the i position in Type II beta-turns. Other turn positions can sometimes be used if the phi angle is near -60 degrees for the residue replaced. It is important that the side chain of the residue replaced is less than 50% buried. Identical substitutions in beta-turns in related proteins give similar results. Proline substitutions increase stability mainly by decreasing the entropy of the denatured state. In contrast, the large, diverse group of proteins considered here had almost no residues in beta-turns that could be replaced by Gly to increase protein stability. Improving beta-turns by substituting Pro residues is a generally useful way of increasing protein stability.

Binding and cleavage specificities of human Argonaute2

The endonuclease Argonaute2 (Ago2) mediates the degradation of the target mRNA within the RNA-induced silencing complex. We determined the binding and cleavage properties of recombinant human Ago2. Human Ago2 was unable to cleave preformed RNA duplexes and exhibited weaker binding affinity for RNA duplexes compared with the single strand RNA. The enzyme exhibited greater RNase H activity in the presence of Mn2+ compared with Mg2+. Human Ago2 exhibited weaker binding affinities and reduced cleavage activities for antisense RNAs with either a 5'-terminal hydroxyl or abasic nucleotide. Binding kinetics suggest that the 5'-terminal heterocycle base nucleates the interaction between the enzyme and the antisense RNA, and the 5'-phosphate stabilizes the interaction. Mn2+ ameliorated the effects of the 5'-terminal hydroxyl or abasic nucleotide on Ago2 cleavage activity and binding affinity. Nucleotide substitutions at the 3' terminus of the antisense RNA had no effect on human Ago2 cleavage activity, whereas 2'-methoxyethyl substitutions at position 2 reduced binding and cleavage activity and 12-14 reduced the cleavage activity. RNase protection assays indicated that human Ago2 interacts with the first 14 nucleotides at the 5'-pole of the antisense RNA. Human Ago2 preloaded with the antisense RNA exhibited greater binding affinities for longer sense RNAs suggesting that the enzyme interacts with regions in the sense RNA outside the site for antisense hybridization. Finally, transiently expressed human Ago2 immunoprecipitated from HeLa cells contained the double strand RNA-binding protein human immunodeficiency virus, type 1, trans-activating response RNA-binding protein, and deletion mutants of Ago2 showed that trans-activating response RNA-binding protein interacts with the PIWI domain of the enzyme.

Structure, stability, and folding of ribonuclease H1 from the moderately thermophilic Chlorobium tepidum: comparison with thermophilic and mesophilic homologues

Proteins from thermophilic organisms are able to function under conditions that render a typical mesophilic protein inactive. Pairwise comparisons of homologous mesophilic and thermophilic proteins can help to identify the energetic features of a protein's energy landscape that lead to such thermostability. Previous studies of bacterial ribonucleases H (RNases H) from the thermophile Thermus thermophilus and the mesophile Escherichia coli revealed that the thermostability arises in part from an unusually low change in heat capacity upon unfolding (DeltaC(p)) for the thermophilic protein [Hollien, J., and Marqusee, S. (1999) Biochemistry 38, 3831-3836]. Here, we have further examined how nearly identical proteins can adapt to different thermal constraints by adding a moderately thermophilic homologue to the previously characterized mesophilic and thermophilic pair. We identified a putative RNase H from Chlorobium. tepidum and demonstrated that it is an active RNase H and adopts the RNase H fold. The moderately thermophilic protein has a melting temperature (T(m)) similar to that of the mesophilic homologue yet also has a surprisingly low DeltaC(p), like the thermophilic homologue. This new RNase H folds through a pathway similar to that of the previously studied RNases H. These results suggest that lowering the DeltaC(p) may be a general strategy for achieving thermophilicity for some protein families and implicate the folding core as the major contributor to this effect. It should now be possible to design RNases H that display the desired thermophilic or mesophilic properties, as defined by their DeltaC(p) values, and therefore fine-tune the energy landscape in a predictable fashion.

Inferring stabilizing mutations from protein phylogenies: application to influenza hemagglutinin

One selection pressure shaping sequence evolution is the requirement that a protein fold with sufficient stability to perform its biological functions. We present a conceptual framework that explains how this requirement causes the probability that a particular amino acid mutation is fixed during evolution to depend on its effect on protein stability. We mathematically formalize this framework to develop a Bayesian approach for inferring the stability effects of individual mutations from homologous protein sequences of known phylogeny. This approach is able to predict published experimentally measured mutational stability effects (ΔΔG values) with an accuracy that exceeds both a state-of-the-art physicochemical modeling program and the sequence-based consensus approach. As a further test, we use our phylogenetic inference approach to predict stabilizing mutations to influenza hemagglutinin. We introduce these mutations into a temperature-sensitive influenza virus with a defect in its hemagglutinin gene and experimentally demonstrate that some of the mutations allow the virus to grow at higher temperatures. Our work therefore describes a powerful new approach for predicting stabilizing mutations that can be successfully applied even to large, complex proteins such as hemagglutinin. This approach also makes a mathematical link between phylogenetics and experimentally measurable protein properties, potentially paving the way for more accurate analyses of molecular evolution.

Mutating a protein frequently causes a change in its stability. As scientists, we often care about these changes because we would like to engineer a protein's stability or understand how its stability is impacted by a naturally occurring mutation. Evolution also cares about mutational stability changes, because a basic evolutionary requirement is that proteins remain sufficiently stable to perform their biological functions. Our work is based on the idea that it should be possible to use the fact that evolution selects for stability to infer from related proteins the effects of specific mutations. We show that we can indeed use protein evolutionary histories to computationally predict previously measured mutational stability changes more accurately than methods based on either of the two main existing strategies. We then test whether we can predict mutations that increase the stability of hemagglutinin, an influenza protein whose rapid evolution is partly responsible for the ability of this virus to cause yearly epidemics. We experimentally create viruses carrying predicted stabilizing mutations and find that several do in fact improve the virus's ability to grow at higher temperatures. Our computational approach may therefore be of use in understanding the evolution of this medically important virus.

Protein side-chain dynamics and residual conformational entropy

Changes in residual conformational entropy of proteins can be significant components of the thermodynamics of folding and binding. Nuclear magnetic resonance (NMR) spin relaxation is the only experimental technique capable of probing local protein entropy, by inference from local internal conformational dynamics. To assess the validity of this approach, the picosecond-to-nanosecond dynamics of the arginine side-chain N(epsilon)-H(epsilon) bond vectors of Escherichia coli ribonuclease H (RNase H) were determined by NMR spin relaxation and compared to the mechanistic detail provided by molecular dynamics (MD) simulations. The results indicate that arginine N(epsilon) spin relaxation primarily reflects persistence of guanidinium salt bridges and correlates well with simulated side-chain conformational entropy. In particular cases, the simulations show that the aliphatic part of the arginine side chain can retain substantial disorder while the guanidinium group maintains its salt bridges; thus, the N(epsilon)-H(epsilon) bond-vector orientation is conserved and side-chain flexibility is concealed from N(epsilon) spin relaxation. The MD simulations and an analysis of a rotamer library suggest that dynamic decoupling of the terminal moiety from the remainder of the side chain occurs for all five amino acids with more than two side-chain dihedral angles (R, K, E, Q, and M). Dynamic decoupling thus may represent a general biophysical strategy for minimizing the entropic penalties of folding and binding.

Optimizing pKa computation in proteins with pH adapted conformations

pK(A) in proteins are determined by electrostatic energy computations using a small number of optimized protein conformations derived from crystal structures. In these protein conformations hydrogen positions and geometries of salt bridges on the protein surface were determined self-consistently with the protonation pattern at three pHs (low, ambient, and high). Considering salt bridges at protein surfaces is most relevant, since they open at low and high pH. In the absence of these conformational changes, computed pK(A)(comp) of acidic (basic) groups in salt bridges underestimate (overestimate) experimental pK(A)(exp), dramatically. The pK(A)(comp) for 15 different proteins with 185 known pK(A)(exp) yield an RMSD of 1.12, comparable with two other methods. One of these methods is fully empirical with many adjustable parameters. The other is also based on electrostatic energy computations using many non-optimized side chain conformers but employs larger dielectric constants at short distances of charge pairs that diminish their electrostatic interactions. These empirical corrections that account implicitly for additional conformational flexibility were needed to describe the energetics of salt bridges appropriately. This is not needed in the present approach. The RMSD of the present approach improves if one considers only strongly shifted pK(A)(exp) in contrast to the other methods under these conditions. Our method allows interpreting pK(A)(comp) in terms of pH dependent hydrogen bonding pattern and salt bridge geometries. A web service is provided to perform pK(A) computations.

Overview of protein structural and functional folds

This overview provides an illustrated, comprehensive survey of some commonly observed protein‐fold families and structural motifs, chosen for their functional significance. It opens with descriptions and definitions of the various elements of protein structure and associated terminology. Following is an introduction into web‐based structural bioinformatics that includes surveys of interactive web servers for protein fold or domain annotation, protein‐structure databases, protein‐structure‐classification databases, structural alignments of proteins, and molecular graphics programs available for personal computers. The rest of the overview describes selected families of protein folds in terms of their secondary, tertiary, and quaternary structural arrangements, including ribbon‐diagram examples, tables of representative structures with references, and brief explanations pointing out their respective biological and functional significance.

Murine leukemia virus reverse transcriptase: structural comparison with HIV-1 reverse transcriptase

Recent X-ray crystal structure determinations of Moloney murine leukemia virus reverse transcriptase (MoMLV RT) have allowed for more accurate structure/function comparisons to HIV-1 RT than were formerly possible. Previous biochemical studies of MoMLV RT in conjunction with knowledge of sequence homologies to HIV-1 RT and overall fold similarities to RTs in general, provided a foundation upon which to build. In addition, numerous crystal structures of the MoMLV RT fingers/palm subdomain had also shed light on one of the critical functions of the enzyme, specifically polymerization. Now in the advent of new structural information, more intricate examination of MoMLV RT in its entirety can be realized, and thus the comparisons with HIV-1 RT may be more critically elucidated. Here, we will review the similarities and differences between MoMLV RT and HIV-1 RT via structural analysis, and propose working models for the MoMLV RT based upon that information.