Advances in implicit models of water solvent to compute conformational free energy and molecular dynamics of proteins at constant pH

Replica exchange molecular dynamics simulations reveal self-association sites in M-crystallin caused by mutations provide insights of cataract

Crystallins are ubiquitous, however, prevalence is seen in eye lens. Eye lens crystallins are long-lived and structural intactness is required for maintaining lens transparency and protein solubility. Mutations in crystallins often lead to cataract. In this study, we performed mutations at specific sites of M-crystallin, a close homologue of eye lens crystallin and studied by using replica exchange molecular dynamics simulation with generalized Born implicit solvent model. Mutations were made on the Ca²⁺ binding residues (K34D and S77D) and in the hydrophobic core (W45R) which is known to cause congenital cataract in homologous γD-crystallin. The chosen mutations caused large motion of the N-terminal Greek key, concomitantly broke the interlocking Greek keys interactions and perturbed the compact core resulting in several folded and partially unfolded states. Partially unfolded states exposed large hydrophobic patches that could act as precursors for self-aggregation. Accumulation of such aggregates is the potential cause of cataract in homologous eye lens crystallins.

An effective molecular blocker of ion channel of M2 protein as anti-influenza a drug

Design of a drug compound that can effectively bind to the M2 ion channel and block the diffusion of hydrogen ions (H⁺) through and inhibit influenza A virus replication is an important task. Known anti-influenza drugs amantadine and rimantadine have a weak effect on influenza A virus. A new class of positively charged, +2 e.u., molecules is proposed here to block diffusion of H⁺ ion through the M2 channel. Several drug candidates, derivatives of a lead compound (diazabicyclooctane), were proposed and investigated. Molecular dynamics of thermal fluctuations of M2 protein structure and ionization-conformation coupling of all the ionizable residues were simulated at physiological pH. The influence of the most probable mutations of key drug-binding amino acid residues in the M2 ion channel were investigated too. It is shown that the suggested new blocker has high binding affinity for the M2 channel. There are two in-channel binding sites of high affinity, the first one has H-bonds with two of four serine residues Ser-31A (B) or Ser-31C(D), and the second one has H-bonds with two of four histidine residues His-37A (B), or His-37C(D). The main advantage of the new drug molecule is the positive charge, +2 e.u., which creates a positive electrostatic potential barrier (in addition to a steric one) for a transfer of H⁺ ion through M2 channel and may serve as an effective anti-influenza A virus drug.Communicated by Ramaswamy H. Sarma.

Catalytically Competent Conformation of the Active Site of Human 8-Oxoguanine-DNA Glycosylase

8-Oxoguanine-DNA N-glycosylase (OGG1) is a eukaryotic DNA repair enzyme responsible for the removal of 8-oxoguanine (oxoG), one of the most abundant oxidative DNA lesions. OGG1 catalyzes two successive reactions - N-glycosidic bond hydrolysis (glycosylase activity) and DNA strand cleavage on the 3'-side of the lesion by β-elimination (lyase activity). The enzyme also exhibits lyase activity with substrates containing apurinic/apyrimidinic (AP) sites (deoxyribose moieties lacking the nucleobase). OGG1 is highly specific for the base opposite the lesion, efficiently excising oxoG and cleaving AP sites located opposite to C, but not opposite to A. The activity is also profoundly decreased by amino acid changes that sterically interfere with oxoG binding in the active site of the enzyme after the lesion is everted from the DNA duplex. Earlier, the molecular dynamics approach was used to study the conformational dynamics of such human OGG1 mutants in complexes with the oxoG:C-containing substrate DNA, and the population density of certain conformers of two OGG1 catalytic residues, Lys249 and Asp268, was suggested to determine the enzyme activity. Here, we report the study of molecular dynamics of human OGG1 bound to the oxoG:A-containing DNA and OGG1 mutants bound to the AP:C-containing DNA. We showed that the enzyme low activity is associated with a decrease in the populations of Lys249 and Asp268 properly configured for catalysis. The experimentally measured rate constants for the OGG1 mutants show a good agreement with the models. We conclude that the enzymatic activity of OGG1 is determined majorly by the population density of the catalytically competent conformations of the active site residues Lys249 and Asp268.

Bottom-up derived flexible water model with dipole and quadrupole moments for coarse-grained molecular simulations

A bottom-up CG water model is developed to capture the electrostatic multipoles, structural correlation and thermodynamics of water.

Coupled molecular dynamics and continuum electrostatic method to compute the ionization pKa's of proteins as a function of pH. Test on a large set of proteins

A computational method, to predict the pKa values of the ionizable residues Asp, Glu, His, Tyr, and Lys of proteins, is presented here. Calculation of the electrostatic free-energy of the proteins is based on an efficient version of a continuum dielectric electrostatic model. The conformational flexibility of the protein is taken into account by carrying out molecular dynamics simulations of 10 ns in implicit water. The accuracy of the proposed method of calculation of pKa values is estimated from a test set of experimental pKa data for 297 ionizable residues from 34 proteins. The pKa-prediction test shows that, on average, 57, 86, and 95% of all predictions have an error lower than 0.5, 1.0, and 1.5 pKa units, respectively. This work contributes to our general understanding of the importance of protein flexibility for an accurate computation of pKa, providing critical insight about the significance of the multiple neutral states of acid and histidine residues for pKa-prediction, and may spur significant progress in our effort to develop a fast and accurate electrostatic-based method for pKa-predictions of proteins as a function of pH.

Molecular dynamics simulation of the opposite-base preference and interactions in the active site of formamidopyrimidine-DNA glycosylase

Background

Formamidopyrimidine-DNA glycosylase (Fpg) removes abundant pre-mutagenic 8-oxoguanine (oxoG) bases from DNA through nucleophilic attack of its N-terminal proline at C1′ of the damaged nucleotide. Since oxoG efficiently pairs with both C and A, Fpg must excise oxoG from pairs with C but not with A, otherwise a mutation occurs. The crystal structures of several Fpg–DNA complexes have been solved, yet no structure with A opposite the lesion is available.

Results

Here we use molecular dynamic simulation to model interactions in the pre-catalytic complex of Lactococcus lactis Fpg with DNA containing oxoG opposite C or A, the latter in either syn or anti conformation. The catalytic dyad, Pro1–Glu2, was modeled in all four possible protonation states. Only one transition was observed in the experimental reaction rate pH dependence plots, and Glu2 kept the same set of interactions regardless of its protonation state, suggesting that it does not limit the reaction rate. The adenine base opposite oxoG was highly distorting for the adjacent nucleotides: in the more stable syn models it formed non-canonical bonds with out-of-register nucleotides in both the damaged and the complementary strand, whereas in the anti models the adenine either formed non-canonical bonds or was expelled into the major groove. The side chains of Arg109 and Phe111 that Fpg inserts into DNA to maintain its kinked conformation tended to withdraw from their positions if A was opposite to the lesion. The region showing the largest differences in the dynamics between oxoG:C and oxoG:A substrates was unexpectedly remote from the active site, located near the linker joining the two domains of Fpg. This region was also highly conserved among 124 analyzed Fpg sequences. Three sites trapping water molecules through multiple bonds were identified on the protein–DNA interface, apparently helping to maintain enzyme-induced DNA distortion and participating in oxoG recognition.

Conclusion

Overall, the discrimination against A opposite to the lesion seems to be due to incorrect DNA distortion around the lesion-containing base pair and, possibly, to gross movement of protein domains connected by the linker.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-017-0075-y) contains supplementary material, which is available to authorized users.

Protein Aggregation and Molecular Crowding: Perspectives From Multiscale Simulations

Cells are extremely crowded environments, thus the use of diluted salted aqueous solutions containing a single protein is too simplistic to mimic the real situation. Macromolecular crowding might affect protein structure, folding, shape, conformational stability, binding of small molecules, enzymatic activity, interactions with cognate biomolecules, and pathological aggregation. The latter phenomenon typically leads to the formation of amyloid fibrils that are linked to several lethal neurodegenerative diseases, but that can also play a functional role in certain organisms. The majority of molecular simulations performed before the last few years were conducted in diluted solutions and were restricted both in the timescales and in the system dimensions by the available computational resources. In recent years, several computational solutions were developed to get close to physiological conditions. In this review we summarize the main computational techniques used to tackle the issue of protein aggregation both in a diluted and in a crowded environment.

8-Oxoguanine Affects DNA Backbone Conformation in the EcoRI Recognition Site and Inhibits Its Cleavage by the Enzyme

8-oxoguanine is one of the most abundant and impactful oxidative DNA lesions. However, the reasons underlying its effects, especially those not directly explained by the altered base pairing ability, are poorly understood. We report the effect of the lesion on the action of EcoRI, a widely used restriction endonuclease. Introduction of 8-oxoguanine inside, or adjacent to, the GAATTC recognition site embedded within the Drew—Dickerson dodecamer sequence notably reduced the EcoRI activity. Solution NMR revealed that 8-oxoguanine in the DNA duplex causes substantial alterations in the sugar—phosphate backbone conformation, inducing a BI→BII transition. Moreover, molecular dynamics of the complex suggested that 8-oxoguanine, although does not disrupt the sequence-specific contacts formed by the enzyme with DNA, shifts the distribution of BI/BII backbone conformers. Based on our data, we propose that the disruption of enzymatic cleavage can be linked with the altered backbone conformation and dynamics in the free oxidized DNA substrate and, possibly, at the protein—DNA interface.

An Atomistic Statistically Effective Energy Function for Computational Protein Design

Shortcomings in the definition of effective free-energy surfaces of proteins are recognized to be a major contributory factor responsible for the low success rates of existing automated methods for computational protein design (CPD). The formulation of an atomistic statistically effective energy function (SEEF) suitable for a wide range of CPD applications and its derivation from structural data extracted from protein domains and protein-ligand complexes are described here. The proposed energy function comprises nonlocal atom-based and local residue-based SEEFs, which are coupled using a novel atom connectivity number factor to scale short-range, pairwise, nonbonded atomic interaction energies and a surface-area-dependent cavity energy term. This energy function was used to derive additional SEEFs describing the unfolded-state ensemble of any given residue sequence based on computed average energies for partially or fully solvent-exposed fragments in regions of irregular structure in native proteins. Relative thermal stabilities of 97 T4 bacteriophage lysozyme mutants were predicted from calculated energy differences for folded and unfolded states with an average unsigned error (AUE) of 0.84 kcal mol(-1) when compared to experiment. To demonstrate the utility of the energy function for CPD, further validation was carried out in tests of its capacity to recover cognate protein sequences and to discriminate native and near-native protein folds, loop conformers, and small-molecule ligand binding poses from non-native benchmark decoys. Experimental ligand binding free energies for a diverse set of 80 protein complexes could be predicted with an AUE of 2.4 kcal mol(-1) using an additional energy term to account for the loss in ligand configurational entropy upon binding. The atomistic SEEF is expected to improve the accuracy of residue-based coarse-grained SEEFs currently used in CPD and to extend the range of applications of extant atom-based protein statistical potentials.

Investigation of Structural Dynamics of Enzymes and Protonation States of Substrates Using Computational Tools

This review discusses the use of molecular modeling tools, together with existing experimental findings, to provide a complete atomic-level description of enzyme dynamics and function. We focus on functionally relevant conformational dynamics of enzymes and the protonation states of substrates. The conformational fluctuations of enzymes usually play a crucial role in substrate recognition and catalysis. Protein dynamics can be altered by a tiny change in a molecular system such as different protonation states of various intermediates or by a significant perturbation such as a ligand association. Here we review recent advances in applying atomistic molecular dynamics (MD) simulations to investigate allosteric and network regulation of tryptophan synthase (TRPS) and protonation states of its intermediates and catalysis. In addition, we review studies using quantum mechanics/molecular mechanics (QM/MM) methods to investigate the protonation states of catalytic residues of β-Ketoacyl ACP synthase I (KasA). We also discuss modeling of large-scale protein motions for HIV-1 protease with coarse-grained Brownian dynamics (BD) simulations.

Effects of mono- and divalent metal ions on DNA binding and catalysis of human apurinic/apyrimidinic endonuclease 1

Here, we used stopped-flow fluorescence techniques to conduct a comparative kinetic analysis of the conformational transitions in human apurinic/apyrimidinic endonuclease 1 (APE1) and in DNA containing an abasic site in the course of their interaction. Effects of monovalent (K(+)) and divalent (Mg(2+), Mn(2+), Ca(2+), Zn(2+), Cu(2+), and Ni(2+)) metal ions on DNA binding and catalytic stages were studied. It was shown that the first step of substrate binding (corresponding to formation of a primary enzyme-substrate complex) does not depend on the concentration (0.05-5.0 mM) or the nature of divalent metal ions. In contrast, the initial DNA binding efficiency significantly decreased at a high concentration (5-250 mM) of monovalent K(+) ions, indicating the involvement of electrostatic interactions in this stage. It was also shown that Cu(2+) ions abrogated the DNA binding ability of APE1, possibly, due to a strong interaction with DNA bases and the sugar-phosphate backbone. In the case of Ca(2+) ions, the catalytic activity of APE1 was lost completely with retention of binding potential. Thus, the enzymatic activity of APE1 is increased in the order Zn(2+) < Ni(2+) < Mn(2+) < Mg(2+). Circular dichroism spectra and calculation of the contact area between APE1 and DNA reveal that Mg(2+) ions stabilize the protein structure and the enzyme-substrate complex.

Molecular dynamics simulations: advances and applications

Molecular dynamics simulations have evolved into a mature technique that can be used effectively to understand macromolecular structure-to-function relationships. Present simulation times are close to biologically relevant ones. Information gathered about the dynamic properties of macromolecules is rich enough to shift the usual paradigm of structural bioinformatics from studying single structures to analyze conformational ensembles. Here, we describe the foundations of molecular dynamics and the improvements made in the direction of getting such ensemble. Specific application of the technique to three main issues (allosteric regulation, docking, and structure refinement) is discussed.

Tri-peptide reference structures for the calculation of relative solvent accessible surface area in protein amino acid residues

Relative amino acid residue solvent accessibility values allow the quantitative comparison of atomic solvent-accessible surface areas in different residue types and physical environments in proteins and in protein structural alignments. Geometry-optimised tri-peptide structures in extended solvent-exposed reference conformations have been obtained for 43 amino acid residue types at a high level of quantum chemical theory. Significant increases in side-chain solvent accessibility, offset by reductions in main-chain atom solvent exposure, were observed for standard residue types in partially geometry-optimised structures when compared to non-minimised models built from identical sets of proper dihedral angles abstracted from the literature. Optimisation of proper dihedral angles led most notably to marked increases of up to 54% in proline main-chain atom solvent accessibility compared to literature values. Similar effects were observed for fully-optimised tri-peptides in implicit solvent. The relief of internal strain energy was associated with systematic variation in N, C(α) and C(β) atom solvent accessibility across all standard residue types. The results underline the importance of optimisation of 'hard' degrees of freedom (bond lengths and valence bond angles) and improper dihedral angle values from force field or other context-independent reference values, and impact on the use of standardised fixed internal co-ordinate geometry in sampling approaches to the determination of absolute values of protein amino acid residue solvent accessibility. Quantum chemical methods provide a useful and accurate alternative to molecular mechanics methods to perform energy minimisation of peptides containing non-standard (chemically modified) amino acid residues frequently present in experimental protein structure data sets, for which force field parameters may not be available. Reference tri-peptide atomic co-ordinate sets including hydrogen atoms are made freely available.

Thermodynamics of the DNA damage repair steps of human 8-oxoguanine DNA glycosylase

Human 8-oxoguanine DNA glycosylase (hOGG1) is a key enzyme responsible for initiating the base excision repair of 7,8-dihydro-8-oxoguanosine (oxoG). In this study a thermodynamic analysis of the interaction of hOGG1 with specific and non-specific DNA-substrates is performed based on stopped-flow kinetic data. The standard Gibbs energies, enthalpies and entropies of specific stages of the repair process were determined via kinetic measurements over a temperature range using the van’t Hoff approach. The three steps which are accompanied with changes in the DNA conformations were detected via 2-aminopurine fluorescence in the process of binding and recognition of damaged oxoG base by hOGG1. The thermodynamic analysis has demonstrated that the initial step of the DNA substrates binding is mainly governed by energy due to favorable interactions in the process of formation of the recognition contacts, which results in negative enthalpy change, as well as due to partial desolvation of the surface between the DNA and enzyme, which results in positive entropy change. Discrimination of non-specific G base versus specific oxoG base is occurring in the second step of the oxoG-substrate binding. This step requires energy consumption which is compensated by the positive entropy contribution. The third binding step is the final adjustment of the enzyme/substrate complex to achieve the catalytically competent state which is characterized by large endothermicity compensated by a significant increase of entropy originated from the dehydration of the DNA grooves.

Electrostatic unfolding and interactions of albumin driven by pH changes: a molecular dynamics study

A better understanding of protein aggregation is bound to translate into critical advances in several areas, including the treatment of misfolded protein disorders and the development of self-assembling biomaterials for novel commercial applications. Because of its ubiquity and clinical potential, albumin is one of the best-characterized models in protein aggregation research; but its properties in different conditions are not completely understood. Here, we carried out all-atom molecular dynamics simulations of albumin to understand how electrostatics can affect the conformation of a single albumin molecule just prior to self-assembly. We then analyzed the tertiary structure and solvent accessible surface area of albumin after electrostatically triggered partial denaturation. The data obtained from these single protein simulations allowed us to investigate the effect of electrostatic interactions between two proteins. The results of these simulations suggested that hydrophobic attractions and counterion binding may be strong enough to effectively overcome the electrostatic repulsions between the highly charged monomers. This work contributes to our general understanding of protein aggregation mechanisms, the importance of explicit consideration of free ions in protein solutions, provides critical new insights about the equilibrium conformation of albumin in its partially denatured state at low pH, and may spur significant progress in our efforts to develop biocompatible protein hydrogels driven by electrostatic partial denaturation.

Measuring the shapes of macromolecules - and why it matters

The molecular basis of life rests on the activity of biological macromolecules, mostly nucleic acids and proteins. A perhaps surprising finding that crystallized over the last handful of decades is that geometric reasoning plays a major role in our attempt to understand these activities. In this paper, we address this connection between geometry and biology, focusing on methods for measuring and characterizing the shapes of macromolecules. We briefly review existing numerical and analytical approaches that solve these problems. We cover in more details our own work in this field, focusing on the alpha shape theory as it provides a unifying mathematical framework that enable the analytical calculations of the surface area and volume of a macromolecule represented as a union of balls, the detection of pockets and cavities in the molecule, and the quantification of contacts between the atomic balls. We have shown that each of these quantities can be related to physical properties of the molecule under study and ultimately provides insight on its activity. We conclude with a brief description of new challenges for the alpha shape theory in modern structural biology.

DNA damage processing by human 8-oxoguanine-DNA glycosylase mutants with the occluded active site

8-Oxoguanine-DNA glycosylase (OGG1) removes premutagenic lesion 8-oxoguanine (8-oxo-G) from DNA and then nicks the nascent abasic (apurinic/apyrimidinic) site by β-elimination. Although the structure of OGG1 bound to damaged DNA is known, the dynamic aspects of 8-oxo-G recognition are not well understood. To comprehend the mechanisms of substrate recognition and processing, we have constructed OGG1 mutants with the active site occluded by replacement of Cys-253, which forms a wall of the base-binding pocket, with bulky leucine or isoleucine. The conformational dynamics of OGG1 mutants were characterized by single-turnover kinetics and stopped-flow kinetics with fluorescent detection. Additionally, the conformational mobility of wild type and the mutant OGG1 substrate complex was assessed using molecular dynamics simulations. Although pocket occlusion distorted the active site and greatly decreased the catalytic activity of OGG1, it did not fully prevent processing of 8-oxo-G and apurinic/apyrimidinic sites. Both mutants were notably stimulated in the presence of free 8-bromoguanine, indicating that this base can bind to the distorted OGG1 and facilitate β-elimination. The results agree with the concept of enzyme plasticity, suggesting that the active site of OGG1 is flexible enough to compensate partially for distortions caused by mutation.

Protein Nano-Object Integrator (ProNOI) for generating atomic style objects for molecular modeling

Background

With the progress of nanotechnology, one frequently has to model biological macromolecules simultaneously with nano-objects. However, the atomic structures of the nano objects are typically not available or they are solid state entities. Because of that, the researchers have to investigate such nano systems by generating models of the nano objects in a manner that the existing software be able to carry the simulations. In addition, it should allow generating composite objects with complex shape by combining basic geometrical figures and embedding biological macromolecules within the system.

Results

Here we report the Protein Nano-Object Integrator (ProNOI) which allows for generating atomic-style geometrical objects with user desired shape and dimensions. Unlimited number of objects can be created and combined with biological macromolecules in Protein Data Bank (PDB) format file. Once the objects are generated, the users can use sliders to manipulate their shape, dimension and absolute position. In addition, the software offers the option to charge the objects with either specified surface or volumetric charge density and to model them with user-desired dielectric constants. According to the user preference, the biological macromolecule atoms can be assigned charges and radii according to four different force fields: Amber, Charmm, OPLS and PARSE. The biological macromolecules and the atomic-style objects are exported as a position, charge and radius (PQR) file, or if a default dielectric constant distribution is not selected, it is exported as a position, charge, radius and epsilon (PQRE) file. As illustration of the capabilities of the ProNOI, we created a composite object in a shape of a robot, aptly named the Clemson Robot, whose parts are charged with various volumetric charge densities and holds the barnase-barstar protein complex in its hand.

Conclusions

The Protein Nano-Object Integrator (ProNOI) is a convenient tool for generating atomic-style nano shapes in conjunction with biological macromolecule(s). Charges and radii on the macromolecule atoms and the atoms in the shapes are assigned according to the user’s preferences allowing various scenarios of modeling. The default output file is in PQR (PQRE) format which is readable by almost any software available in biophysical field. It can be downloaded from: http://compbio.clemson.edu/downloadDir/ProNO_integrator.tar.gz

Coarse-Grained Molecular Models of Water: A Review

Coarse-grained (CG) models have proven to be very effective tools in the study of phenomena or systems that involve large time- and length-scales. By decreasing the degrees of freedom in the system and using softer interactions than seen in atomistic models, larger timesteps can be used and much longer simulation times can be studied. CG simulations are widely used to study systems of biological importance that are beyond the reach of atomistic simulation, necessitating a computationally efficient and accurate CG model for water. In this review, we discuss the methods used for developing CG water models and the relative advantages and disadvantages of the resulting models. In general, CG water models differ with regards to how many waters each CG group or bead represents, whether analytical or tabular potentials have been used to describe the interactions, and how the model incorporates electrostatic interactions. Finally, how the models are parameterized depends on their application, so, while some are fitted to experimental properties such as surface tension and density, others are fitted to radial distribution functions extracted from atomistic simulations.

Thermodynamics of the multi-stage DNA lesion recognition and repair by formamidopyrimidine-DNA glycosylase using pyrrolocytosine fluorescence--stopped-flow pre-steady-state kinetics

Formamidopyrimidine-DNA glycosylase, Fpg protein from Escherichia coli, initiates base excision repair in DNA by removing a wide variety of oxidized lesions. In this study, we perform thermodynamic analysis of the multi-stage interaction of Fpg with specific DNA-substrates containing 7,8-dihydro-8-oxoguanosine (oxoG), or tetrahydrofuran (THF, an uncleavable abasic site analog) and non-specific (G) DNA-ligand based on stopped-flow kinetic data. Pyrrolocytosine, highly fluorescent analog of the natural nucleobase cytosine, is used to record multi-stage DNA lesion recognition and repair kinetics over a temperature range (10–30°C). The kinetic data were used to obtain the standard Gibbs energy, enthalpy and entropy of the specific stages using van’t Hoff approach. The data suggest that not only enthalpy-driven exothermic oxoG recognition, but also the desolvation-accompanied entropy-driven enzyme-substrate complex adjustment into the catalytically active state play equally important roles in the overall process.

Structural features of aquaporin 4 supporting the formation of arrays and junctions in biomembranes

A limited class of aquaporins has been described to form regular arrays and junctions in membranes. The biological significance of these structures, however, remains uncertain. Here we analyze the underlying physical principles with the help of a computational procedure that takes into account protein-protein as well as protein-membrane interactions. Experimentally observed array/junction structures are systematically (dis)assembled and major driving forces identified. Aquaporin 4 was found to be markedly different from the non-junction forming aquaporin 1. The environmental stabilization resulting from embedding into the biomembrane was identified as the main driving force. This highlights the role of protein-membrane interactions in aquaporin 4. Analysis of the type presented here can help to decipher the biological role of membrane arrays and junctions formed by aquaporin.