General trends of dihedral conformational transitions in a globular protein

Dihedral conformational transitions are analyzed systematically in a model globular protein, cytochrome P450cam, to examine their structural and chemical dependences through combined conventional molecular dynamics (cMD), accelerated molecular dynamics (aMD) and adaptive biasing force (ABF) simulations. The aMD simulations are performed at two acceleration levels, using dihedral and dual boost, respectively. In comparison with cMD, aMD samples protein dihedral transitions approximately two times faster on average using dihedral boost, and ～3.5 times faster using dual boost. In the protein backbone, significantly higher dihedral transition rates are observed in the bend, coil, and turn flexible regions, followed by the β bridge and β sheet, and then the helices. Moreover, protein side chains of greater length exhibit higher transition rates on average in the aMD‐enhanced sampling. Side chains of the same length (particularly N_χ = 2) exhibit decreasing transition rates with residues when going from hydrophobic to polar, then charged and aromatic chemical types. The reduction of dihedral transition rates is found to be correlated with increasing energy barriers as identified through ABF free energy calculations. These general trends of dihedral conformational transitions provide important insights into the hierarchical dynamics and complex free energy landscapes of functional proteins. Proteins 2016; 84:501–514. © 2016 Wiley Periodicals, Inc.     

NMR spectroscopy, molecular dynamics, and conformation of a synthetic octasaccharide fragment of the O-specific polysaccharide of Shigella dysenteriae type 1

A synthetic octasaccharide fragment (2) of the O-specific polysaccharide (1) of Shigella dysenteriae type 1 has been studied as its methyl glycoside by one- and two-dimensional homo- and heteronuclear NMR spectroscopy. Complete 1H and 13C NMR assignments have been generated, and the 13C spin-lattice relaxation times have been measured for the octasaccharide 2. A congener (6) of this octasaccharide containing one D-galactose residue with a specific 13C label at C-1 has been synthesized and used to measure interglycosidic 13C-1H coupling by the 2D J-resolved 1H NMR method. From the NMR data, three types of conformational restraints were developed: (a) 29 inter-residue, distance restraints; (b) 48 intra-residue, ring atom dihedral angle restraints, and (c) one heteronuclear, inter-residue dihedral angle restraint. The use of these restraints in a restrained molecular dynamics computation with simulated annealing yielded a conformation resembling a short, irregular spiral, with methyl substituents on the exterior.     

Coarse‐grained simulations of proton‐dependent conformational changes in lactose permease

During lactose/H⁺ symport, the Escherichia coli lactose permease (LacY) undergoes a series of global conformational transitions between inward‐facing (open to cytoplasmic side) and outward‐facing (open to periplasmic side) states. However, the exact local interactions and molecular mechanisms dictating those large‐scale structural changes are not well understood. All‐atom molecular dynamics simulations have been performed to investigate the molecular interactions involved in conformational transitions of LacY, but the simulations can only explore early or partial global structural changes because of the computational limits (< 100 ns). In this work, we implement a hybrid force field that couples the united‐atom protein models with the coarse‐grained MARTINI water/lipid, to investigate the proton‐dependent dynamics and conformational changes of LacY. The effects of the protonation states on two key glutamate residues (Glu325 and Glu269) have been studied. Our results on the salt‐bridge dynamics agreed with all‐atom simulations at early short time period, validating our simulations. From our microsecond simulations, we were able to observe the complete transition from inward‐facing to outward‐facing conformations of LacY. Our results showed that all helices have participated during the global conformational transitions and helical movements of LacY. The inter‐helical distances measured in our simulations were consistent with the double electron‐electron resonance experiments at both cytoplasmic and periplasmic sides. Our simulations indicated that the deprotonation of Glu325 induced the opening of the periplasmics side and partial closure of the cytoplasmic side of LacY, while protonation of the Glu269 caused a stable cross‐domain salt‐bridge (Glu130‐Arg344) and completely closed the cytoplasmic side. Proteins 2016; 84:1067–1074. © 2016 Wiley Periodicals, Inc.     

Toward a detailed understanding of search trajectories in fragment assembly approaches to protein structure prediction

Energy functions, fragment libraries, and search methods constitute three key components of fragment‐assembly methods for protein structure prediction, which are all crucial for their ability to generate high‐accuracy predictions. All of these components are tightly coupled; efficient searching becomes more important as the quality of fragment libraries decreases. Given these relationships, there is currently a poor understanding of the strengths and weaknesses of the sampling approaches currently used in fragment‐assembly techniques. Here, we determine how the performance of search techniques can be assessed in a meaningful manner, given the above problems. We describe a set of techniques that aim to reduce the impact of the energy function, and assess exploration in view of the search space defined by a given fragment library. We illustrate our approach using Rosetta and EdaFold, and show how certain features of these methods encourage or limit conformational exploration. We demonstrate that individual trajectories of Rosetta are susceptible to local minima in the energy landscape, and that this can be linked to non‐uniform sampling across the protein chain. We show that EdaFold's novel approach can help balance broad exploration with locating good low‐energy conformations. This occurs through two mechanisms which cannot be readily differentiated using standard performance measures: exclusion of false minima, followed by an increasingly focused search in low‐energy regions of conformational space. Measures such as ours can be helpful in characterizing new fragment‐based methods in terms of the quality of conformational exploration realized. Proteins 2016; 84:411–426. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Theoretical probes of conformational fluctuations in S-peptide and RNase A/3′–UMP enzyme product complex

The dynamic properties of the RNase A/3′–UMP enzyme/product complex and the S-peptide of RNase A have been investigated by molecular dynamics simulations using suitable generalization of ideas introduced to probe the energy landscape in structural glasses. We introduce two measures, namely, the kinetic energy fluctuation metric and the force metric, both of which are used to calculate the time needed for sampling the conformation space of the molecules. The calculation of the fluctuation metric requires a single trajectory whereas the force metric is computed using two independent trajectories. The vacuum MD simulations show that for both systems the time required for kinetic energy equipartitioning is surprisingly long even at high temperatures. We show that the force metric is a powerful means of probing the nature and relative importance of conformational substates which determine the dynamics at low temperatures. In particular the time dependence of the non-bonded force metric is used to demonstrate that at low temperatures the system is predominantly localized hi a single cluster of conformational substates. The force metric is used to show that relaxation of long range (in sequence space) interactions must be mediated by a sequence of local dihedral angle transitions. We also argue that the time needed for compact structure formation is intimately related to the time needed for the relaxation of the dihedral angle degrees of freedom. The tame for non-bonded interactions, which drive protein molecules to fold under appropriate conditions, to relax becomes extremely long as the temperature is lowered suggesting that the formation of maximally compact structure hi proteins must be a very slow process. © 1993 Wiley-Liss, Inc.     

Method of modeling protein structure by the two-dimensional nuclear magnetic resonance spectroscopy data; application to the proteinase inhibitor BUSI IIA from bull seminal plasma

A new approach is suggested to model the spatial structure of protein molecules in solution based on combined use of the methods of theoretical conformational analysis and NMR spectroscopy data. At the first stage, special means are used to convert d connectivity information into the most probable values of dihedral angles. This allows search for possible spatial structures in the limited regions of the conformational space at further stages using the methods of the theoretical conformational analysis. The suggested approach was verified in reconstructing the spatial backbone structure of the fragment 17-57 of the proteinase inhibitor BUSI IIA from the bull seminal plasma. The structural model parameters are compared with the corresponding characteristics obtained from the X-ray analysis data for the homologic proteinase inhibitor from the Japanese quail ovomucoid. The suggested approach is shown to correctly reproduce both the general molecule topology and the conformations of individual amino acid residues.     

Variation of the NH-C alpha-H coupling constant with dihedral angle in the NMR spectra of peptides

Proton magnetic resonance data and conformational calculations of a series of model compounds containing a NH-C^αH group substituted as in peptides have been used to generate a proton–proton coupling constant–dihedral angle relation for the peptide unit. For those substances used in which the dihedral angle about the N-C^α bond is not fixed, the angle distribution was calculated from conformational theory. Using eight examples in which the number of theoretical assumptions were least, the best values of the coefficients A, B, and C in the expression J(θ) = Acos²θ + B cosθ + Csin²θ were found by a least-squares procedure to be 7.9, ?1.55, and 1.35, respectively. This relation gives reasonable values for the dihedral angles ? in cyclic oligopeptide structures for which the availability of both NMR data and other structural information allow comparison. When applied to N-acetylamino acid N-methylamides having side chains extending beyond C^β, however, agreement with the calculated conformational distribution was found for Leu, Met, and Trp, but observed values of J were larger than expected for Val, He, Phe, and Tyr, These disagreements are considered to be the result of interactions not yet taken into account in the usual conformational calculations.     

Conformational variation of the central CG site in d(ATGACGTCAT)2 and d(GAAAACGTTTTC)2. An NMR, molecular modelling and 3D-homology investigation.

The determination of the solution structure of two self-complementary oligomers d(ATGACGTCAT)2 (CG10) and d(GAAAACGTTTTC)2 (CG12), both containing the 5'-pur-ACGT-pyr-3' sequence, is reported. The impact of the base context on the conformation of the central CpG site has been examined by a combined approach of: (a) 2D 1H-NMR and 31P-NMR; (b) molecular mechanics under experimental constraints; (c) back-calculations of NOESY spectra and iterative refinements of distances; and (d) 3D-homology search of the central tetrad ACGT within the complete oligonucleotides. A full NMR study of each fragment is achieved by means of standard 2D experiments: NOESY, 2D homonuclear Hartmann-Hahn spectroscopy, double-quantum-filtered COSY and heteronuclear 1H-31P correlation. Sugar phase angle, epsilon-zeta difference angle and NOE-derived distances are input as experimental constraints to generate molecular models by energy minimization with the help of jumna. The morass program is used to iteratively refine the structures obtained. The similarity of the two ACGTs within the whole oligonucleotides is investigated. Both the decamer and the dodecamer adopt a B-like DNA conformation. However, the helical parameters within this conformational type are significantly different in CG12 and CG10. The central CpG step conformation is not locked by its nearest environment (5'A and 3'T) as seen from the structural analysis of ACGT in the two molecules. In CG12, despite the presence of runs of A-T pairs, CpG presents a high twist of 43 degrees and a sugar phase at the guanine of about 180 degrees, previously observed in other ACGT-containing-oligomers. Conversely, ACGT in CG10 exhibits strong inclinations, positive rolls, a flat profile of sugar phase, twist and glycosidic angles, as a result of the nucleotide sequence extending beyond the tetrad. The structural specificity of CG10 and its flexibility (as reflected by its energy) are tentatively related to the process of recognition of the cyclic AMP response element by its cognate protein.     

Membrane‐induced structure of novel human tachykinin hemokinin‐1 (hHK1)

PPT‐C encoded hemokinin‐1(hHK‐1) of Homo sapiens (TGKASQFFGLM) is a structurally distinct neuropeptide among the tachykinin family that participate in the NK‐1 receptor downstream signaling processes. Subsequently, signal transduction leads to execution of various effector functions which includes aging, immunological, and central nervous system (CNS) regulatory actions. However the conformational pattern of ligand receptor binding is unclear. The three‐dimensional structure of the hemokinin‐1 in aqueous and micellar environment has been studied by one and two‐dimensional proton nuclear magnetic resonance (2D 1H‐NMR spectroscopy) and distance geometry calculations. Data shows that hemokinin‐1 was unstructured in aqueous environment; anionic detergent SDS induces α‐helix formation. Proton NMR assignments have been carried out with the aid of correlation spectroscopy (gradient‐COSY and TOCSY) and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The inter proton distances and dihedral angle constraints obtained from the NMR data have been used in torsion angle dynamics algorithm for NMR applications (CYANA) to generate a family of structures, which have been refined using restrained energy minimization and dynamics. Helical conformation is observed from residue K3‐M11. The conformational range of the peptide revealed by NMR studies has been analyzed in terms of characteristic secondary features. Observed conformational features have been compared to that of Substance P potent NK1 agonist. Thus the report provides a structural insight to study hHK‐1‐NK1 interaction that is essential for hHK1 based signaling events. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 702–710, 2015.     

Investigating diproline segments in proteins: occurrences, conformation and classification

The covalent linkage between the side‐chain and the backbone nitrogen atom of proline leads to the formation of the five‐membered pyrrolidine ring and hence restriction of the backbone torsional angle ? to values of ?60 °± 30° for the L ‐proline. Diproline segments constitute a chain fragment with considerably reduced conformational choices. In the current study, the conformational states for the diproline segment ( ^L Pro‐ ^L Pro) found in proteins has been investigated with an emphasis on the cis and trans states for the Pro‐Pro peptide bond. The occurrence of diproline segments in turns and other secondary structures has been studied and compared to that of Xaa‐Pro‐Yaa segments in proteins which gives us a better understanding on the restriction imposed on other residues by the diproline segment and the single proline residue. The study indicates that P_II–P_II and P_II–α are the most favorable conformational states for the diproline segment. The analysis on Xaa‐Pro‐Yaa sequences reveals that the Xaa‐Pro peptide bond exists preferably as the trans conformer rather than the cis conformer. The present study may lead to a better understanding of the behavior of proline occurring in diproline segments which can facilitate various designed diproline‐based synthetic templates for biological and structural studies. © 2011 Wiley Periodicals, Inc. Biopolymers 97: 54–64, 2012.     

MOL3D,a modular and interactive program for molecular modeling and conformational analysis: I — basic modules

MOL3D is a generalized machine-independent computer program that lets the user interactively build 3D structures with different display options, such as wire, ball-and-stick and CPK representations. The program, which uses its own graphics package and driver, is designed to be very user friendly through the use of commands and menus. It has powerful transformation capabilities, such as software rotations, superpositions and zooming, and it is equipped with a fragment database that allows the user to build complex structures. The algorithm presented here is designed to perform computations in all the conformational space and therefore can be used to predict experimentally available quantities, such as NMR coupling constants. The program is efficient in the sense that it handles only dihedral angles in the first steps; as a result, it allows a rapid sampling of a great number of points through the entire conformational space. The user can choose between grid and Monte-Carlo searches of energy minimization, using a reasonable amount of computer time.     

Conformational behavior of α,α-dialkylated peptides

The preferred conformations of N-acetyl-N′-methyl amides of some dialkylglycines have been determined by empirical conformational-energy calculations; minimum-energy conformations were located by minimizing the energy with respect to all the dihedral angles of the molecules. The conformational space of these compounds is sterically restricted, and low-energy conformations are found only in the regions of fully extended and helical structures. Increasing the bulkiness of the substituents on the C^α, the fully extended conformation becomes gradually more stable than the helical structure preferred in the cases of dimethylglycine. This trend is, however, strongly dependent on the bond angles between the substituents on the C^α atom: In particular, helical structures are favored by standard values (111°) of the N-C^α-C′ angle, while fully extended conformations are favored by smaller values of the same angle, as experimentally observed, for instance, in the case of α,α-di-n-propylglycine.     

Simulated annealing with restrained molecular dynamics using a flexible restraint potential: theory and evaluation with simulated NMR constraints.

A new functional representation of NMR-derived distance constraints, the flexible restraint potential, has been implemented in the program CONGEN (Bruccoleri RE, Karplus M, 1987, Biopolymers 26:137-168) for molecular structure generation. In addition, flat-bottomed restraint potentials for representing dihedral angle and vicinal scalar coupling constraints have been introduced into CONGEN. An effective simulated annealing (SA) protocol that combines both weight annealing and temperature annealing is described. Calculations have been performed using ideal simulated NMR constraints, in order to evaluate the use of restrained molecular dynamics (MD) with these target functions as implemented in CONGEN. In this benchmark study, internuclear distance, dihedral angle, and vicinal coupling constant constraints were calculated from the energy-minimized X-ray crystal structure of the 46-amino acid polypeptide crambin (ICRN). Three-dimensional structures of crambin that satisfy these simulated NMR constraints were generated using restrained MD and SA. Polypeptide structures with extended backbone and side-chain conformations were used as starting conformations. Dynamical annealing calculations using extended starting conformations and assignments of initial velocities taken randomly from a Maxwellian distribution were found to adequately sample the conformational space consistent with the constraints. These calculations also show that loosened internuclear constraints can allow molecules to overcome local minima in the search for a global minimum with respect to both the NMR-derived constraints and conformational energy. This protocol and the modified version of the CONGEN program described here are shown to be reliable and robust, and are applicable generally for protein structure determination by dynamical simulated annealing using NMR data.     

Efficient Search on Energy Minima for Structure Prediction of Nucleic Acid Motifs

Abstract

Structure prediction of non-canonical motifs such as mismatches, extra unmatched nucleotides or internal and hairpin loop structures in nucleic acids is of great importance for understanding the function and design of nucleic acid structures. Systematic conformational analysis of such motifs typically involves the generation of many possible combinations of backbone dihedral torsion angles for a given motif and subsequent energy minimization (EM) and evaluation. Such approach is limited due to the number of dihedral angle combinations that grows very rapidly with the size of the motif. Two conformational search approaches have been developed that allow both an effective crossing of barriers during con-formational searches and the computational demand grows much less with system size then search methods that explore all combinations of backbone dihedral torsion angles. In the first search protocol single torsion angles are flipped into favorable states using constraint EM and subsequent relaxation without constraints. The approach is repeated in an iterative manner along the backbone of the structural motif until no further energy improvement is obtained. In case of two test systems, a DNA-trinucleotide loop (sequence: GCA) and a RNA tetraloop (sequence: UUCG), the approach successfully identified low energy states close to experiment for two out of five start structures. In the second method randomly selected combinations of up to six backbone torsion angles are simultaneously flipped into preset ranges by a short constraint EM followed by unconstraint EM and acceptance according to a Metropolis acceptance criterion. This combined stochastic/EM search was even more effective than the single torsion flip approach and selected low energy states for the two test cases in between two and four cases out of five start structures.     

Efficient search on energy minima for structure prediction of nucleic acid motifs

Structure prediction of non-canonical motifs such as mismatches, extra unmatched nucleotides or internal and hairpin loop structures in nucleic acids is of great importance for understanding the function and design of nucleic acid structures. Systematic conformational analysis of such motifs typically involves the generation of many possible combinations of backbone dihedral torsion angles for a given motif and subsequent energy minimization (EM) and evaluation. Such approach is limited due to the number of dihedral angle combinations that grows very rapidly with the size of the motif. Two conformational search approaches have been developed that allow both an effective crossing of barriers during conformational searches and the computational demand grows much less with system size then search methods that explore all combinations of backbone dihedral torsion angles. In the first search protocol single torsion angles are flipped into favorable states using constraint EM and subsequent relaxation without constraints. The approach is repeated in an iterative manner along the backbone of the structural motif until no further energy improvement is obtained. In case of two test systems, a DNA-trinucleotide loop (sequence: GCA) and a RNA tetraloop (sequence: UUCG), the approach successfully identified low energy states close to experiment for two out of five start structures. In the second method randomly selected combinations of up to six backbone torsion angles are simultaneously flipped into preset ranges by a short constraint EM followed by unconstraint EM and acceptance according to a Metropolis acceptance criterion. This combined stochastic/EM search was even more effective than the single torsion flip approach and selected low energy states for the two test cases in between two and four cases out of five start structures.     

Artificial proteins as allosteric modulators of PDZ3 and SH3 in two‐domain constructs: A computational characterization of novel chimeric proteins

Artificial multidomain proteins with enhanced structural and functional properties can be utilized in a broad spectrum of applications. The design of chimeric fusion proteins utilizing protein domains or one‐domain miniproteins as building blocks is an important advancement for the creation of new biomolecules for biotechnology and medical applications. However, computational studies to describe in detail the dynamics and geometry properties of two‐domain constructs made from structurally and functionally different proteins are lacking. Here, we tested an in silico design strategy using all‐atom explicit solvent molecular dynamics simulations. The well‐characterized PDZ3 and SH3 domains of human zonula occludens (ZO‐1) (3TSZ), along with 5 artificial domains and 2 types of molecular linkers, were selected to construct chimeric two‐domain molecules. The influence of the artificial domains on the structure and dynamics of the PDZ3 and SH3 domains was determined using a range of analyses. We conclude that the artificial domains can function as allosteric modulators of the PDZ3 and SH3 domains. Proteins 2016; 84:1358–1374. © 2016 Wiley Periodicals, Inc.     

Ensemble models of proteins and protein domains based on distance distribution restraints

Conformational ensembles of intrinsically disordered peptide chains are not fully determined by experimental observations. Uncertainty due to lack of experimental restraints and due to intrinsic disorder can be distinguished if distance distributions restraints are available. Such restraints can be obtained from pulsed dipolar electron paramagnetic resonance (EPR) spectroscopy applied to pairs of spin labels. Here, we introduce a Monte Carlo approach for generating conformational ensembles that are consistent with a set of distance distribution restraints, backbone dihedral angle statistics in known protein structures, and optionally, secondary structure propensities or membrane immersion depths. The approach is tested with simulated restraints for a terminal and an internal loop and for a protein with 69 residues by using sets of sparse restraints for underlying well‐defined conformations and for published ensembles of a premolten globule‐like and a coil‐like intrinsically disordered protein. Proteins 2016; 84:544–560. © 2016 Wiley Periodicals, Inc.     

In vitro selection and characterization of DARPins and Fab fragments for the co-crystallization of membrane proteins: The Na(+)-citrate symporter CitS as an example

The determination of 3D structures of membrane proteins is still extremely difficult. The co-crystallization with specific binding proteins may be an important aid in this process, as these proteins provide rigid, hydrophilic surfaces for stable protein-protein contacts. Also, the conformational homogeneity of the membrane protein may be increased to obtain crystals suitable for high resolution structures. Here, we describe the efficient generation and characterization of Designed Ankyrin Repeat Proteins (DARPins) as specific binding molecules for membrane proteins. We used both phage display and ribosome display to select DARPins in vitro that are specific for the detergent-solubilized Na(+)-citrate symporter CitS of Klebsiella pneumoniae. Compared to classical hybridoma technology, the in vitro selection systems allow a much better control of the structural integrity of the target protein and allow the use of other protein classes in addition to recombinant antibodies. We also compared the selected DARPins to a Fab fragment previously selected by phage display and demonstrate that different epitopes are recognized, unique to each class of binding molecules. Therefore, the use of several classes of binding molecules will make suitable crystal formation and the determination of their 3D structure more likely.