Search

The contrast‐variation method in small‐angle neutron scattering (SANS) is a uniquely powerful technique for determining the structure of individual components in biomolecular systems containing regions of different neutron scattering length density ρ. By altering the ρ of the target solute and the solvent through judicious incorporation of deuterium, the scattering of desired solute features can be highlighted. Most contrast‐variation methods focus on highlighting specific bulk solute elements, but not on how the scattering at specific scattering vectors q, which are associated with specific structural distances, changes with contrast. Indeed, many systems exhibit q‐dependent contrast effects. Here, a method is presented for calculating both bulk contrast‐match points and q‐dependent contrast using 3D models with explicit solute and solvent atoms and SASSENA, an explicit‐atom SANS calculator. The method calculates the bulk contrast‐match points within 2.4% solvent D2O accuracy for test protein–nucleic acid and lipid nanodisc systems. The method incorporates a general model for the incorporation of deuterium at non‐exchangeable sites that was derived by performing mass spectrometry on green fluorescent protein. The method also decomposes the scattering profile into its component parts and identifies structural features that change with contrast. The method is readily applicable to a variety of systems, will expand the understanding of q‐dependent contrast matching and will aid in the optimization of next‐generation neutron scattering experiments. [ABSTRACT FROM AUTHOR]

Structure‐based virtual high‐throughput screening involves docking chemical libraries to targets of interest. A parameter pertinent to the accuracy of the resulting pose is the root mean square deviation (RMSD) from a known crystallographic structure, i. e., the 'docking power'. Here, using a popular algorithm, Autodock Vina, as a model program, we evaluate the effects of varying two common docking parameters: the box size (the size of docking search space) and the exhaustiveness of the global search (the number of independent runs starting from random ligand conformations) on the RMSD from the PDBbind v2017 refined dataset of experimental protein‐ligand complexes. Although it is clear that exhaustiveness is an important parameter, there is wide variation in the values used, with variation between 1 and >100. We, therefore, evaluated a combination of cubic boxes of different sizes and five exhaustiveness values (1, 8, 25, 50, 75, 100) within the range of those commonly adopted. The results show that the default exhaustiveness value of 8 performs well overall for most box sizes. In contrast, for all box sizes, but particularly for large boxes, an exhaustiveness value of 1 led to significantly higher median RMSD (mRMSD) values. The docking power was slightly improved with an exhaustiveness of 25, but the mRMSD changes little with values higher than 25. Therefore, although low exhaustiveness is computationally faster, the results are more likely to be far from reality, and, conversely, values >25 led to little improvement at the expense of computational resources. Overall, we recommend users to use at least the default exhaustiveness value of 8 for virtual screening calculations. [ABSTRACT FROM AUTHOR]

Amphiphilic co‐solvents can have a significant impact on the structure, organization and physical properties of lipid bilayers. Describing the mutual impact of partitioning and induced structure changes is therefore a crucial consideration for a range of topics such as anesthesia and other pharmacokinetic effects, as well as microbial solvent tolerance in the production of biofuels and other fermentation products, where molecules such as ethanol, butanol or acetic acid might be generated. Small‐angle neutron scattering (SANS) is a key method for studying lipid and polymer bilayer structures, with many models for extracting bilayer structure (thickness, area per lipid etc.) from scattering data in use today. However, the molecular details of co‐solvent partitioning are conflated with induced changes to bilayer structure, making interpretation and modeling of the scattering curves a challenge with the existing set of models. To address this, a model of a bilayer structure is presented which invokes a two‐term partition constant accounting for the localization of the co‐solvent within the bilayer. This model was validated using a series of SANS measurements of lipid vesicles in the presence of the co‐solvent tetrahydrofuran (THF), showing several strategies of how to deploy the two‐parameter partition constant model to describe scattering data and extract both structure and partitioning information from the data. Molecular dynamics simulations are then used to evaluate assumptions of the model, provide additional molecular scale details and illustrate its complementary nature to the data fitting procedure. This approach results in estimates of the partition coefficient for THF in 1,2‐dimyristoyl‐sn‐glycero‐3‐phosphocholine at 35°C, along with an estimate of the fraction of THF residing in the hydrophobic core of the membrane. The authors envision that this model will be applicable to a wide range of other bilayer/amphiphile interactions and provide the associated code needed to implement this model as a fitting algorithm for scattering data in the SasView suite. [ABSTRACT FROM AUTHOR]

The startling diversity in αβ T‐cell receptor (TCR) sequences and structures complicates molecular‐level analyses of the specificity and sensitivity determining T‐cell immunogenicity. A number of three‐dimensional (3D) structures are now available of ternary complexes between TCRs and peptides: major histocompatibility complexes (pMHC). Here, to glean molecular‐level insights we analyze structures of TCRs bound to human class I nonamer peptide–MHC complexes. Residues at peptide positions 4–8 are found to be particularly important for TCR binding. About 90% of the TCRs hydrogen bond with one or both of the peptide residues at positions 4 and 8 presented by MHC allele HLA‐A2, and this number is still ~79% for peptides presented by other MHC alleles. Residue 8, which lies outside the previously‐identified central peptide region, is crucial for TCR recognition of class I MHC‐presented nonamer peptides. The statistics of the interactions also sheds light on the MHC residues important for TCR binding. The present analysis will aid in the structural modeling of TCR:pMHC complexes and has implications for the rational design of peptide‐based vaccines and T‐cell‐based immunotherapies. [ABSTRACT FROM AUTHOR]

Human interleukin‐6 (hIL‐6) is a multifunctional cytokine that regulates immune and inflammatory responses in addition to metabolic and regenerative processes and cancer. hIL‐6 binding to the IL‐6 receptor (IL‐6Rα) induces homodimerization and recruitment of the glycoprotein (gp130) to form a hexameric signaling complex. Anti‐IL‐6 and IL‐6R antibodies are clinically approved inhibitors of IL‐6 signaling pathway for treating rheumatoid arthritis and Castleman's disease, respectively. There is a potential to develop novel small molecule IL‐6 antagonists derived from understanding the structural basis for IL‐6/IL‐6Rα interactions. Here, we combine homology modeling with extensive molecular dynamics (MD) simulations to examine the association of hIL‐6 with IL‐6Rα. A comparison with MD of apo hIL‐6 reveals that the binding of hIL‐6 to IL‐6Rα induces structural and dynamic rearrangements in the AB loop region of hIL‐6, disrupting intraprotein contacts and increasing the flexibility of residues 48 to 58 of the AB loop. In contrast, due to the involvement of residues 59 to 78 in forming contacts with the receptor, these residues of the AB loop are observed to rigidify in the presence of the receptor. The binary complex is primarily stabilized by two pairs of salt bridges, Arg181 (hIL‐6)‐ Glu182 (IL‐6Rα) and Arg184 (hIL‐6)‐ Glu183 (IL‐6Rα) as well as hydrophobic and aromatic stacking interactions mediated essentially by Phe residues in both proteins. An interplay of electrostatic, hydrophobic, hydrogen bonding, and aromatic stacking interactions facilitates the formation of the hIL‐6/IL‐6Rα complex. [ABSTRACT FROM AUTHOR]

The mer operon in bacteria encodes a set of proteins and enzymes that impart resistance to environmental mercury toxicity by importing Hg2+ and reducing it to volatile Hg(0). Because the reduction occurs in the cytoplasm, mercuric ions must first be transported across the cytoplasmic membrane by one of a few known transporters. MerF is the smallest of these, containing only two transmembrane helices and two pairs of vicinal cysteines that coordinate mercuric ions. In this work, we use molecular dynamics simulations to characterize the dynamics of MerF in its apo and Hg2+‐bound states. We find that the apo state positions one of the cysteine pairs closer to the periplasmic side of the membrane, while in the bound state the same pair approaches the cytoplasmic side. This finding is consistent with the functional requirement of accepting Hg2+ from the periplasmic space, sequestering it on acceptance, and transferring it to the cytoplasm. Conformational changes in the TM helices facilitate the functional interaction of the two cysteine pairs. Free‐energy calculations provide a barrier of 16 kcal/mol for the association of the periplasmic Hg2+‐bound protein MerP with MerF and 7 kcal/mol for the subsequent association of MerF's two cysteine pairs. Despite the significant conformational changes required to move the binding site across the membrane, coarse‐grained simulations of multiple copies of MerF support the expectation that it functions as a monomer. Our results demonstrate how conformational changes and binding thermodynamics could lead to such a small membrane protein acting as an ion transporter. Published 2019. This article is a U.S. Government work and is in the public domain in the USA. [ABSTRACT FROM AUTHOR]

To assess the chemical reactivity, toxicity, and mobility of pollutants in the environment, knowledge of their species distributions is critical. Because their direct measurement is often infeasible, speciation modeling is widely adopted. Mercury (Hg) is a representative pollutant for which study of its speciation benefits from modeling. However, Hg speciation modeling is often hindered by a lack of reliable thermodynamic constants. Although computational chemistry (e.g., density functional theory [DFT]) can generate these constants, methods for directly coupling DFT and speciation modeling are not available. Here, we combine computational chemistry and continuum‐scale modeling with curated online databases to ameliorate the problem of unreliable inputs to Hg speciation modeling. Our AQUA‐MER databases and web server (https://aquamer.ornl.gov) provides direct speciation results by combining web‐based interfaces to a speciation calculator, databases of thermodynamic constants, and a computational chemistry toolkit to estimate missing constants. Although Hg is presented as a concrete use case, AQUA‐MER can also be readily applied to other elements. © 2019 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

Summary: A microbe's ecological niche and biotechnological utility are determined by its specific set of co‐evolved metabolic pathways. The acquisition of new pathways, through horizontal gene transfer or genetic engineering, can have unpredictable consequences. Here we show that two different pathways for coumarate catabolism failed to function when initially transferred into Escherichia coli. Using laboratory evolution, we elucidated the factors limiting activity of the newly acquired pathways and the modifications required to overcome these limitations. Both pathways required host mutations to enable effective growth with coumarate, but the necessary mutations differed. In one case, a pathway intermediate inhibited purine nucleotide biosynthesis, and this inhibition was relieved by single amino acid replacements in IMP dehydrogenase. A strain that natively contains this coumarate catabolism pathway, Acinetobacter baumannii, is resistant to inhibition by the relevant intermediate, suggesting that natural pathway transfers have faced and overcome similar challenges. Molecular dynamics simulation of the wild type and a representative single‐residue mutant provide insight into the structural and dynamic changes that relieve inhibition. These results demonstrate how deleterious interactions can limit pathway transfer, that these interactions can be traced to specific molecular interactions between host and pathway, and how evolution or engineering can alleviate these limitations. [ABSTRACT FROM AUTHOR]

Ensemble docking is now commonly used in early‐stage in silico drug discovery and can be used to attack difficult problems such as finding lead compounds which can disrupt protein–protein interactions. We give an example of this methodology here, as applied to fibroblast growth factor 23 (FGF23), a protein hormone that is responsible for regulating phosphate homeostasis. The first small‐molecule antagonists of FGF23 were recently discovered by combining ensemble docking with extensive experimental target validation data (Science Signaling, 9, 2016, ra113). Here, we provide a detailed account of how ensemble‐based high‐throughput virtual screening was used to identify the antagonist compounds discovered in reference (Science Signaling, 9, 2016, ra113). Moreover, we perform further calculations, redocking those antagonist compounds identified in reference (Science Signaling, 9, 2016, ra113) that performed well on drug‐likeness filters, to predict possible binding regions. These predicted binding modes are rescored with the molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) approach to calculate the most likely binding site. Our findings suggest that the antagonist compounds antagonize FGF23 through the disruption of protein–protein interactions between FGF23 and fibroblast growth factor receptor (FGFR). [ABSTRACT FROM AUTHOR]

Channel and pump rhodopsins use energy from light absorbed by a covalently bound retinal chromophore to transport ions across membranes of microbial cells. Ion transfer steps, including proton transfer, can couple to changes in protein conformational dynamics and water positions. Although general principles of how microbial rhodopsins function are largely understood, key issues pertaining to reaction mechanisms remain unclear. In this review, we compare the protonation-coupled dynamics of pump and channelrhodopsins, highlighting the roles that water dynamics, protein electrostatics and protein flexibility can have in ion transport mechanisms. We discuss observations supporting important functional roles of inter- and intra-helical carboxylate/hydroxyl hydrogen-bonding motifs. As specific examples, we use the proton pump bacteriorhodopsin, the sodium pump KR2, channelrhodopsins and Anabaena sensory rhodopsin. We outline the usefulness of theoretic biophysics approaches to the study of retinal proteins, challenges in studying the hydrogen-bond dynamics of rhodopsin active sites, and implications for conformational coupling in membrane transporters. [ABSTRACT FROM AUTHOR]

Neutron crystallography was used to directly locate two protons before and after a pH-induced two-proton transfer between catalytic aspartic acid residues and the hydroxy group of the bound clinical drug darunavir, located in the catalytic site of enzyme HIV-1 protease. The two-proton transfer is triggered by electrostatic effects arising from protonation state changes of surface residues far from the active site. The mechanism and pH effect are supported by quantum mechanics/molecular mechanics (QM/MM) calculations. The low-pH proton configuration in the catalytic site is deemed critical for the catalytic action of this enzyme and may apply more generally to other aspartic proteases. Neutrons therefore represent a superb probe to obtain structural details for proton transfer reactions in biological systems at a truly atomic level. [ABSTRACT FROM AUTHOR]

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 Nx=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. [ABSTRACT FROM AUTHOR]

Low-molecular-weight lignin binds to cellulose during the thermochemical pretreatment of biomass for biofuel production, which prevents the efficient hydrolysis of the cellulose to sugars. The binding properties of lignin are influenced strongly by the conformations it adopts. Here, we use molecular dynamics simulations in aqueous solution to investigate the dependence of the shape of lignin polymers on chain length and temperature. Lignin is found to adopt collapsed conformations in water at 300 and 500 K. However, at 300 K, a discontinuous transition is found in the shape of the polymer as a function of the chain length. Below a critical degree of polymerization, Nc=15, the polymer adopts less spherical conformations than above Nc. The transition disappears at high temperatures (500 K) at which only spherical shapes are adopted. An implication relevant to cellulosic biofuel production is that lignin will self-aggregate even at high pretreatment temperatures. [ABSTRACT FROM AUTHOR]

ABSTRACT One-dimensional (1D) (spherically averaged) powder diffraction diagrams are commonly used to determine the degree of cellulose crystallinity in biomass samples. Here, it is shown using molecular modeling how disorder in cellulose fibrils can lead to considerable uncertainty in conclusions drawn concerning crystallinity based on 1D powder diffraction data alone. For example, cellulose microfibrils that contain both crystalline and noncrystalline segments can lead to powder diffraction diagrams lacking identifiable peaks, while microfibrils without any crystalline segments can lead to such peaks. This leads to false positives, that is, assigning disordered cellulose as crystalline, and false negatives, that is, categorizing fibrils with crystalline segments as amorphous. The reliable determination of the fraction of crystallinity in any given biomass sample will require a more sophisticated approach combining detailed experiment and simulation. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 67-73, 2015. [ABSTRACT FROM AUTHOR]

The characterization of the structure and internal dynamics of biomolecules is essential to understanding their biological function. Neutron scattering probes similar time- and length-scales to molecular dynamics simulation. Hence, simulation models of biomolecules have become invaluable in the interpretation of experimental neutron data. Here, we report on advances in the application of simulation in developing neutron scattering to investigate internal protein motions and, as an example of industrial relevance, in the derivation of physical models of use in biofuel renewable energy research. [ABSTRACT FROM AUTHOR]

SUMMARY In this paper, we give the current state of high-throughput virtual screening. We describe a case study of using a task-parallel Message Passing Interface version of AutoDock4 to run a virtual high-throughput screen of one million compounds on the Jaguar Cray XK6 Supercomputer at Oak Ridge National Laboratory. We include a description of scripts developed to increase the efficiency of the predocking file preparation and postdocking analysis. A detailed tutorial, scripts, and source code for this Message Passing Interface version of AutoDock4 are available online at . Copyright © 2013 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

The program VinaMPI has been developed to enable massively large virtual drug screens on leadership-class computing resources, using a large number of cores to decrease the time-to-completion of the screen. VinaMPI is a massively parallel Message Passing Interface (MPI) program based on the multithreaded virtual docking program AutodockVina, and is used to distribute tasks while multithreading is used to speed-up individual docking tasks. VinaMPI uses a distribution scheme in which tasks are evenly distributed to the workers based on the complexity of each task, as defined by the number of rotatable bonds in each chemical compound investigated. VinaMPI efficiently handles multiple proteins in a ligand screen, allowing for high-throughput inverse docking that presents new opportunities for improving the efficiency of the drug discovery pipeline. VinaMPI successfully ran on 84,672 cores with a continual decrease in job completion time with increasing core count. The ratio of the number of tasks in a screening to the number of workers should be at least around 100 in order to have a good load balance and an optimal job completion time. The code is freely available and downloadable. Instructions for downloading and using the code are provided in the Supporting Information. © 2013 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

The Diels-Alderase ribozyme is an in vitro-evolved ribonucleic acid enzyme that catalyzes a [4 + 2] cycloaddition reaction between an anthracene diene and a maleimide dienophile. The ribozyme can in principle be used to selectively synthesize only one product enantiomer, depending on which of the two entrances to the catalytic pocket, 'front' or 'back', the substrate is permitted to use. Here, we investigate stereoselection and substrate recognition in the ribozyme by means of multiple molecular dynamics simulations, performed on each of the two substrates individually in the pocket, on the reactant state, and on the product state. The results are consistent with a binding mechanism in which the maleimide likely binds first followed by the anthracene, which enters preferentially through the front door. The free energy profiles for anthracene binding indicate that the pre-( R, R)-enantiomer conformation is slightly preferred, in agreement with the experimentally observed small enantiomeric excess of the ( R, R)-enantiomer of the product. The reactant state is stabilized by the simultaneous presence of both substrates bound to their binding sites in the hydrophobic pocket as well as by stacking interactions between them. © 2012 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

We report on a python interface to the GROMACS molecular simulation package,GromPy (available at https://github.com/GromPy). This application programming interface (API) uses the ctypes python module that allows function calls to shared libraries, for example, written in C.To the best of our knowledge, this is the first reported interface to the GROMACS library that uses direct library calls. GromPy can be used for extending the current GROMACS simulation and analysis modes. In this work, we demonstrate that the interface enables hybrid Monte-Carlo/molecular dynamics (MD) simulations in the grand-canonical ensemble, a simulation mode that is currently not implemented in GROMACS. For this application, the interplay between GromPy and GROMACS requires only minor modifications of the GROMACS source code, not affecting the operation, efficiency, and performance of the GROMACS applications. We validate the grand-canonical application against MD in the canonical ensemble by comparison of equations of state. The results of the grand-canonical simulations are in complete agreement with MD in the canonical ensemble. The python overhead of the grand-canonical scheme is only minimal. [ABSTRACT FROM AUTHOR]

Lignocellulosic biomass provides a sustainable source of sugars for biofuel and biomaterial production. However, biomass resistance to degradation imposes difficulties for economical conversion of plant carbohydrates to fermentable sugars. One of the key contributors to recalcitrance is lignin. Understanding the properties of lignin macromolecules in the cell wall matrix is useful for manipulating biomass structure to generate more easily degradable biomass. Along with experimental techniques such as 2D-NMR and mass spectrometry, computational techniques can be useful for characterizing the structural and energetic properties of the biomass assembly and its individual constituents. Here, we provide a brief introduction to lignin, review some of the recent, relevant scientific literature, and give our perspectives on the role of molecular simulation in understanding lignin structure. © 2011 American Institute of Chemical Engineers Environ Prog, 2012 [ABSTRACT FROM AUTHOR]

Abstract: The restriction endonuclease EcoRV binds two magnesium ions. One of these ions, , binds to the phosphate group where the cleavage occurs and is required for catalysis, but the role of the other ion, is debated. Here, multiple independent molecular dynamics simulations suggest that is crucial for achieving a tightly bound protein–DNA complex and stabilizing a conformation that allows cleavage. In the absence of in all simulations the protein–DNA hydrogen bond network is significantly disrupted and the sharp kink at the central base pair step of the DNA, which is observed in the two-metal complex, is not present. Also, the active site residues rearrange in such a way that the formation of a nucleophile, required for DNA hydrolysis, is unlikely. [Copyright &y& Elsevier]

Actin is a major structural protein of the eukaryotic cytoskeleton and enables cell motility. Here, we present a model of the actin filament (F-actin) that not only incorporates the global structure of the recently published model by Oda et al. but also conserves internal stereochemistry. A comparison is made using molecular dynamics simulation of the model with other recent F-actin models. A number of structural determents such as the protomer propeller angle, the number of hydrogen bonds, and the structural variation among the protomers are analyzed. The MD comparison is found to reflect the evolution in quality of actin models over the last 6 years. In addition, simulations of the model are carried out in states with both ADP or ATP bound and local hydrogen-bonding differences characterized. Proteins 2011; © 2011 Wiley-Liss, Inc. [ABSTRACT FROM AUTHOR]

Copyright of Journal of Computational Chemistry is the property of Wiley-Blackwell and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

The assembly of monomeric G-actin into filamentous F-actin is nucleotide dependent: ATP-G-actin is favored for filament growth at the 'barbed end' of F-actin, whereas ADP-G-actin tends to dissociate from the 'pointed end.' Structural differences between ATP- and ADP-G-actin are examined here using multiple molecular dynamics simulations. The 'open' and 'closed' conformational states of G-actin in aqueous solution are characterized, with either ATP or ADP in the nucleotide binding pocket. With both ATP and ADP bound, the open state closes in the absence of actin-bound profilin. The position of the nucleotide in the protein is found to be correlated with the degree of opening of the active site cleft. Further, the simulations reveal the existence of a structurally well-defined, compact, 'superclosed' state of ATP-G-actin, as yet unseen crystallographically and absent in the ADP-G-actin simulations. The superclosed state resembles structurally the actin monomer in filament models derived from fiber diffraction and is putatively the polymerization competent conformation of ATP-G-actin. Proteins 2009. © 2008 Wiley-Liss, Inc. [ABSTRACT FROM AUTHOR]

The proton pumping cycle of bacteriorhodopsin consists of five consecutive proton transfer steps, in all of which water molecules are likely to play critical roles. Here we review the role of water molecules in bacteriorhodopsin function, with emphasis on computations of the retinal Schiff base deprotonation reaction and the open questions regarding the participation of water molecules in the longer-distance proton transfer reactions. [ABSTRACT FROM AUTHOR]

The organization and assembly of the cellulosome, an extracellular multienzyme complex produced by anaerobic bacteria, is mediated by the high-affinity interaction of cohesin domains from scaffolding proteins with dockerins of cellulosomal enzymes. We have performed molecular dynamics simulations and free energy calculations on both the wild type (WT) and D39N mutant of the C. thermocellum Type I cohesin-dockerin complex in aqueous solution. The D39N mutation has been experimentally demonstrated to disrupt cohesin-dockerin binding. The present MD simulations indicate that the substitution triggers significant protein flexibility and causes a major change of the hydrogen-bonding network in the recognition strips-the conserved loop regions previously proposed to be involved in binding-through electrostatic and salt-bridge interactions between β-strands 3 and 5 of the cohesin and α-helix 3 of the dockerin. The mutation-induced subtle disturbance in the local hydrogen-bond network is accompanied by conformational rearrangements of the protein side chains and bound water molecules. Additional free energy perturbation calculations of the D39N mutation provide differences in the cohesin-dockerin binding energy, thus offering a direct, quantitative comparison with experiments. The underlying molecular mechanism of cohesin-dockerin complexation is further investigated through the free energy profile, that is, potential of mean force (PMF) calculations of WT cohesin-dockerin complex. The PMF shows a high-free energy barrier against the dissociation and reveals a stepwise pattern involving both the central β-sheet interface and its adjacent solvent-exposed loop/turn regions clustered at both ends of the β-barrel structure. [ABSTRACT FROM AUTHOR]

A CHARMM molecular mechanics force field for lignin is derived. Parameterization is based on reproducing quantum mechanical data of model compounds. Partial atomic charges are derived using the RESP electrostatic potential fitting method supplemented by the examination of methoxybenzene:water interactions. Dihedral parameters are optimized by fitting to critical rotational potentials and bonded parameters are obtained by optimizing vibrational frequencies and normal modes. Finally, the force field is validated by performing a molecular dynamics simulation of a crystal of a lignin fragment molecule and comparing simulation-derived structural features with experimental results. Together with the existing force field for polysaccharides, this lignin force field will enable full simulations of lignocellulose. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009 [ABSTRACT FROM AUTHOR]

Understanding X-ray crystallographic diffuse scattering is likely to improve our comprehension of equilibrium collective protein dynamics. Here, using molecular dynamics (MD) simulation, a detailed analysis is performed of the origins of diffuse scattering in crystalline Staphylococcal nuclease, for which the complete diffuse scattering pattern has been determined experimentally. The hydrogen-atom contribution and the scattering range over which the scattering can be considered to be a sum of solvent and protein scattering are determined. Two models of correlated protein motion are investigated by calculating the model-derived diffuse scattering and comparing with the scattering calculated directly from MD trajectories. In one model, previously used in diffuse scattering interpretation, the atomic displacement correlations decay isotropically with increasing separation. Model correlation lengths are obtained by refining the model scattering against the simulation-derived scattering pattern, and are found to be significantly different from those correlation lengths derived directly from the MD trajectories. Furthermore, the convergence between the model-derived and MD-derived scattering is poor. The second model, in which the displacement correlations are calculated from the principal components of the MD trajectories, is capable of fully reproducing the MD-derived diffuse scattering if the ≈50% lowest-frequency modes are included. However, a small number (≈10) of lowest-frequency and largest-amplitude modes dominates the diffuse scattering and thus the correlated protein motions. A detailed analysis of the principal components is performed. In particular, the effective free energy profile associated with each principle mode is analyzed and the eigenfrequency and damping coefficient computed using a model of Brownian dynamics. Those collective modes with effective frequencies below ≈0.5 THz, including those that determine the diffuse scattering, are overdamped. Proteins 2007. © 2006 Wiley-Liss, Inc. [ABSTRACT FROM AUTHOR]

Effective analysis of large-scale conformational transitions in macromolecules requires transforming them into a lower dimensional representation that captures the dominant motions. Herein, we apply and compare two different dimensionality reduction techniques, namely, principal component analysis (PCA), a linear method, and Sammon mapping, which is nonlinear. The two methods are used to analyze four different protein transition pathways of varying complexity, obtained by using either the conjugate peak refinement method or constrained molecular dynamics. For the return-stroke in myosin, both Sammon mapping and PCA show that the conformational change is dominated by a simple rotation of a rigid body. Also, in the case of the T→R transition in hemoglobin, both methods are able to identify the two main quaternary transition events. In contrast, in the cases of the unfolding transition of staphylococcal nuclease or the signaling switch of Ras p21, which are both more complex conformational transitions, only Sammon mapping is able to identify the distinct phases of motion. Proteins 2006. © 2006 Wiley-Liss, Inc. [ABSTRACT FROM AUTHOR]

A parameterization has been performed of the biologically important sterols cholesterol, ergosterol, and lanosterol for the CHARMM27 all-atom molecular mechanics force field. An automated parameterization method was used that involves fitting the potential to vibrational frequencies and eigenvectors derived from quantum-chemical calculations. The partial charges were derived by fitting point charges to quantum-chemically calculated electrostatic potentials. To model the dynamics of the hydroxyl groups of the sterols correctly, the parameter set was refined to reproduce the energy barrier for the rotation of the hydroxyl group around the carbon connected to the hydroxyl of each sterol. The frequency-matching plots show good agreement between the CHARMM and quantum chemical normal modes. The parameters are tested in a molecular dynamics simulation of the cholesterol crystal structure. The experimental geometry and cell dimensions are well reproduced. The force field derived here is also useful for simulating other sterols such as the phytosterols sigmasterol, and campesterol, and a variety of steroids. © 2005 Wiley Periodicals, Inc. J Comput Chem 26: 1383–1399, 2005 [ABSTRACT FROM AUTHOR]

In molecular mechanics calculations, electrostatic interactions between chemical groups are usually represented by a Coulomb potential between the partial atomic charges of the groups. In aqueous solution these interactions are modified by the polarizable solvent. Although the electrostatic effects of the polarized solvent on the protein are well described by the Poisson–Boltzmann equation, its numerical solution is computationally expensive for large molecules such as proteins. The procedure of nonuniform charge scaling (NUCS) is a pragmatic approach to implicit solvation that approximates the solvent screening effect by individually scaling the partial charges on the explicit atoms of the macromolecule so as to reproduce electrostatic interaction energies obtained from an initial Poisson–Boltzmann analysis. Once the screening factors have been determined for a protein the scaled charges can be easily used in any molecular mechanics program that implements a Coulomb term. The approach is particularly suitable for minimization-based simulations, such as normal mode analysis, certain conformational reaction path or ligand binding techniques for which bulk solvent cannot be included explicitly, and for combined quantum mechanical/molecular mechanical calculations when the interface to more elaborate continuum solvent models is lacking. The method is illustrated using reaction path calculations of the Tyr35 ring flip in the bovine pancreatic trypsin inhibitor. © 2005 Wiley Periodicals, Inc. J Comput Chem 26: 1359–1371, 2005 [ABSTRACT FROM AUTHOR]

The computation of minimum energy paths (MEPs) is an approach for gaining insight into protein conformational transitions that are too slow to be observed with unconstrained molecular dynamics simulations. MEPs have the advantage of providing the energy barrier of the rate-limiting step(s), allowing discrimination among different paths. Finding low-energy MEPs for complex transitions, such as those involving rearrangements of the backbone fold or repacking of buried side chains, has hitherto been unfeasible in a reliable, automated manner, the MEP often displaying unphysical behavior, such as the crossing of bonds. Here, this problem is addressed by combining a counterintuitive procedure for generating an initial guess of the path, in which all side chains are shrunk, with the conjugate peak refinement (CPR) method. The effectiveness of the approach is tested on the conformational switch in Ras p21. This conformational transition involves some partial unfolding and re-folding, a process for which a multitude of pathways are likely to exist and for which a single MEP does not provide a complete description. However, this transition requires some sterically demanding rearrangements, thus testing the ability of a method to find low-energy pathways free of structurally unphysical events. This is achieved by the present approach, which finds a path whose rate-limiting barrier is compatible with experiment. This demonstrates that the method can be used to compute plausible pathways for complex rearrangements in proteins in an automated manner that is unbiased by external driving constraints. Proteins 2005. © 2005 Wiley-Liss, Inc. [ABSTRACT FROM AUTHOR]

Understanding the energetic factors determining the positioning and orientation of single-helical peptides in membranes is of fundamental interest in structural biology. Here, a simple 5-slab continuum dielectric model for the membrane is examined that distinguishes between the solvent, headgroup, and core regions. An analytical solution for the electrostatic solvation of a single dipole and an all-atom model of N-methylacetamide are used to demonstrate the effect of the dielectric boundaries in the system on peptide dipole orientation. The dipole orientation energy is shown to dominate the electrostatic solvation energy of a polyalanine helix in the membrane. With an additional surface-area-dependent term to account for the cavity formation in the aqueous region, the continuum electrostatics description is used to examine several helical peptides, the atoms of which are explicitly represented with a molecular mechanics force field. The experimentally determined tilt angles of a number of peptides of alternating alanine and leucine residues, and of glycophorin and melittin, are accurately reproduced by the model. The factors determining the tilt angles and their fluctuations are analyzed. The tilt angles of the simpler peptides are found to increase approximately linearly with peptide length; this effect is also rationalized. The analysis and model presented here provide a step toward the prediction of helical membrane protein structure. Proteins 2005. © 2005 Wiley-Liss, Inc. [ABSTRACT FROM AUTHOR]

A method is presented to calculate efficiently small-angle neutron and X-ray solution scattering intensities from explicit-atom models of macromolecules and the surrounding solvent. The method is based on a multipole expansion of the scattering amplitude. It is particularly appropriate for extensive configurational averaging, as is required for calculations based on computer-simulation results. In test calculations, excellent agreement with experiment is found between neutron and X-ray scattering profiles calculated from a molecular-dynamics simulation of lysozyme in water. The question of definition of the protein surface is also addressed. For comparison with the continuum model, an analytical envelope around the protein is defined in terms of spherical harmonics and is calculated using a Lebedev grid. The analytical surface thus defined is shown to reproduce well the scattering profile calculated from the explicit-atom model of the protein. [ABSTRACT FROM AUTHOR]

Cytochrome P-450cam oxidises 1,2,3,4-tetrahydro-β-carboline, an indolic alkaloid. We report here measurements of the product distribution of this oxidation. To rationalise the experimental results, ab initio quantum-chemistry calculations of the product stabilities and molecular-dynamics calculations of the substrate-binding mode in the active site were performed. The calculations suggest that the product distribution is influenced by both the relative intrinsic gas-phase stabilities of the monohydroxy products and by conformational rearrangement of the active site on substrate binding. [ABSTRACT FROM AUTHOR]