The role of the active site lysine residue on FAD reduction by NADPH in glutathione reductase

Glutathione reductase (GR) is a two dinucleotide binding domain flavoprotein (tDBDF) that catalyzes the reduction of glutathione disulfide to glutathione coupled to the oxidation of NADPH to NADP+. An interesting feature of GR and other tDBDFs is the presence of a lysine residue (Lys-66 in human GR) at the active site, which interacts with the flavin group, but has an unknown function. To better understand the role of this residue, the dynamics of GR was studied using molecular dynamics simulations, and the reaction mechanism of FAD reduction by NADPH was studied using QM/MM molecular modeling. The two possible protonation states of Lys-66 were considered: neutral and protonated. Molecular dynamics results suggest that the active site is more structured for neutral Lys-66 than for protonated Lys-66. QM/MM modeling results suggest that Lys-66 should be in its neutral state for a thermodynamically favorable reduction of FAD by NADPH. Since the reaction is unfavorable with protonated Lys-66, the reverse reaction (the reduction of NADP+ by FADH-) is expected to take place. A phylogenetic analysis of various tDBDFs was performed, finding that an active site lysine is present in different the tDBDFs enzymes, suggesting that it has a conserved biological role. Overall, these results suggest that the protonation state of the active site lysine determines the energetics of the reaction, controlling its reversibility.

Ten things I `hate' about refinement

Macromolecular refinement is an optimization process that aims to produce the most likely macromolecular structural model in the light of experimental data. As such, macromolecular refinement is one of the most complex optimization problems in wide use. Macromolecular refinement programs have to deal with the complex relationship between the parameters of the atomic model and the experimental data, as well as a large number of types of prior knowledge about chemical structure. This paper draws attention to areas of unfinished business in the field of macromolecular refinement. In it, we describe ten refinement topics that we think deserve attention and discuss directions leading to macromolecular refinement software that would make the best use of modern computer resources to meet the needs of structural biologists of the twenty-first century.

A complete Fourier-synthesis-based backbone-conformation-dependent library for proteins

While broadening the applicability of (φ/ψ)-dependent target values for the bond angles in the peptide backbone, sequence/conformation categories with too few residues to analyze via previous methods were encountered. Here, a method of describing a conformation-dependent library (CDL) using two-dimensional Fourier coefficients is reported where the number of coefficients for individual categories is determined via complete cross-validation. Sample sizes are increased further by selective blending of categories with similar patterns of conformational dependence. An additional advantage of the Fourier-synthesis-based CDL is that it uses continuous functions and has no artifactual steps near the edges of populated regions of φ/ψ space. A set of libraries for the seven main-chain bond angles, along with the ω and ζ angles, was created based on a set of Fourier analyses of 48 368 residues selected from high-resolution models in the wwPDB. This new library encompasses both trans- and cis-peptide bonds and outperforms currently used discrete CDLs.

Explicit polarization: a quantum mechanical framework for developing next generation force fields

Molecular mechanical force fields have been successfully used to model condensed-phase and biological systems for a half century. By means of careful parametrization, such classical force fields can be used to provide useful interpretations of experimental findings and predictions of certain properties. Yet, there is a need to further improve computational accuracy for the quantitative prediction of biomolecular interactions and to model properties that depend on the wave functions and not just the energy terms. A new strategy called explicit polarization (X-Pol) has been developed to construct the potential energy surface and wave functions for macromolecular and liquid-phase simulations on the basis of quantum mechanics rather than only using quantum mechanical results to fit analytic force fields. In this spirit, this approach is called a quantum mechanical force field (QMFF).

X-Pol is a general fragment method for electronic structure calculations based on the partition of a condensed-phase or macromolecular system into subsystems (“fragments”) to achieve computational efficiency. Here, intrafragment energy and the mutual electronic polarization of interfragment interactions are treated explicitly using quantum mechanics. X-Pol can be used as a general, multilevel electronic structure model for macromolecular systems, and it can also serve as a new-generation force field. As a quantum chemical model, a variational many-body (VMB) expansion approach is used to systematically improve interfragment interactions, including exchange repulsion, charge delocalization, dispersion, and other correlation energies. As a quantum mechanical force field, these energy terms are approximated by empirical functions in the spirit of conventional molecular mechanics. This Account first reviews the formulation of X-Pol, in the full variationally correct version, in the faster embedded version, and with systematic many-body improvements. We discuss illustrative examples involving water clusters (which show the power of two-body corrections), ethylmethylimidazolium acetate ionic liquids (which reveal that the amount of charge transfer between anion and cation is much smaller than what has been assumed in some classical simulations), and a solvated protein in aqueous solution (which shows that the average charge distribution of carbonyl groups along the polypeptide chain depends strongly on their position in the sequence, whereas they are fixed in most classical force fields). The development of QMFFs also offers an opportunity to extend the accuracy of biochemical simulations to areas where classical force fields are often insufficient, especially in the areas of spectroscopy, reactivity, and enzyme catalysis.

Mimicking DNA stretching with the Static Mode method: shear stress versus transverse pulling stress

DNA sequencing using nanopores is closer than ever to become a reality, but further research and development still need to be done, especially to unravel the atomic-scale mechanisms of induced DNA stretching. At this level, molecular modeling and simulation are essential to investigate DNA conformational flexibility and its response to the forces involved. In this work, through a "Static Mode" approach, we present a directed exploration of the deformations of a 27-mer subjected to externally imposed forces, as it could be in a nanopore. We show how the DNA sugar-phosphate backbone undergoes the majority of the induced deformation, before the base pairing is affected, and to what extent unzipping initiation depends on the force direction.

Spatial organization and conformational peculiarities of the callatostatin family of neuropeptides

The structures and conformational peculiarities of five members of the callatostatin family of neuropeptides, i.e. Leu- and Met-callatostatins, ranging in size from 8 to 16 amino acid residues have been investigated by a theoretical conformational analysis method. A comparative analysis of the conformational flexibilities of Met-callatostatin with those of the hydroxylated analogues, [Hyp2]- and [Hyp3]-Met-callatostatin has been carried out. Helically packed C-terminal pentapeptide in the structure of all investigated Leu-callatostatins are shown to be possible. The reason for the great number low-energy conformers for the callatostatin N-terminus is discussed.

Design of peptides, proteins, and peptidomimetics in chi space

Peptide and protein biological activities depend on their three dimensionals structures in the free state and when interacting with their receptors/acceptors. The backbone conformations such as alpha-helix, beta-sheet, beta-turn, and so forth provide critical templates for the three-dimensional structure, but the overall shape and intrinsic stereoelectronic properties of the peptide or protein important for molecular recognition, signal transduction, enzymatic specificity, immunomodulation, and other biological effects depend on arrangement of the side chain groups in three-dimensional chi space (their chi 1, chi 2, etc. torsional angles). In this paper we explore approaches to the de novo design of polypeptides and peptidomimetics with biased or specific conformational/topographical properties in chi space. We consider computational and experimental methods that can be used to examine the effects of specific structural modifications in constraining side chain groups of amino acid residues and their similarities in chi space to the natural amino acids to evaluate what sort of mimetics are likely to mimic normal amino acids. We then examine some of the asymmetric synthetic methods that are being developed to obtain the amino acid mimetics. Finally, we consider selected examples in the literature where these specialized amino acids have been incorporated in biologically active peptides and the specific insights they have provided regarding the topographical requirements for bioactive peptide potency, selectivity, and other biochemical and pharmacological properties. Constraints in chi space show great promise as useful tools in peptide, protein, and peptidomimetic de novo design of structures and pharmacophores with specific stereostructural, biochemical and biological properties.

Possible bioactive conformations of alpha-melanotropin

By comparing the results of energy calculation for alpha-MSH and its semirigid analogues Ac-[ Cys4 , Cys10 ]-alpha- MSH4 -10-NH2, Ac-[ Cys4 , Cys10 ]-alpha- MSH4 -13-NH2, and [ Cys4 - Cys10 ]-alpha-MSH, a detailed description of two possible bioactive conformations for the 'specific' central site of alpha- MSH6 -9 is proposed representing variants of chain-reversal structure. A possible explanation of the rise in melanotropic activity of the latter two semirigid analogues is presented.