Search

Structure-based virtual high-throughput screening is used in early-stage drug discovery. Over the years, docking protocols and scoring functions for protein-ligand complexes have evolved to improve the accuracy in the computation of binding strengths and poses. In the past decade, RNA has also emerged as a target class for new small-molecule drugs. However, most ligand docking programs have been validated and tested for proteins and not RNA. Here, we test the docking power (pose prediction accuracy) of three state-of-the-art docking protocols on 173 RNA-small molecule crystal structures. The programs are AutoDock4 (AD4) and AutoDock Vina (Vina), which were designed for protein targets, and rDock, which was designed for both protein and nucleic acid targets. AD4 performed relatively poorly. For RNA targets for which a crystal structure of a bound ligand used to limit the docking search space is available and for which the goal is to identify new molecules for the same pocket, rDock performs slightly better than Vina, with success rates of 48% and 63%, respectively. However, in the more common type of early-stage drug discovery setting, in which no structure of a ligand-target complex is known and for which a larger search space is defined, rDock performed similarly to Vina, with a low success rate of ∼27%. Vina was found to have bias for ligands with certain physicochemical properties, whereas rDock performs similarly for all ligand properties. Thus, for projects where no ligand-protein structure already exists, Vina and rDock are both applicable. However, the relatively poor performance of all methods relative to protein-target docking illustrates a need for further methods refinement.

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., (© 2022 Wiley-VCH GmbH.)

The scattering of neutrons can be used to provide information on the structure and dynamics of biological systems on multiple length and time scales. Pursuant to a National Science Foundation-funded workshop in February 2018, recent developments in this field are reviewed here, as well as future prospects that can be expected given recent advances in sources, instrumentation and computational power and methods. Crystallography, solution scattering, dynamics, membranes, labeling and imaging are examined. For the extraction of maximum information, the incorporation of judicious specific deuterium labeling, the integration of several types of experiment, and interpretation using high-performance computer simulation models are often found to be particularly powerful.

Motional displacements of hydrogen (H) in proteins can be measured using incoherent neutron-scattering methods. These displacements can also be calculated numerically using data from molecular dynamics simulations. An enormous amount of data on the average mean-square motional displacement (MSD) of H as a function of protein temperature, hydration, and other conditions has been collected. H resides in a wide spectrum of sites in a protein. Some H are tightly bound to molecular chains, and the H motion is dictated by that of the chain. Other H are quite independent. As a result, there is a distribution of motions and MSDs of H within a protein that is denoted dynamical heterogeneity. The goal of this paper is to incorporate a distribution of MSDs into models of the H incoherent intermediate scattering function, I(Q,t), that is calculated and observed. The aim is to contribute information on the distribution as well as on the average MSD from comparison of the models with simulations and experiment. For example, we find that simulations of I(Q,t) in lysozyme are well reproduced if the distribution of MSDs is bimodal with two broad peaks rather than a single broad peak., (Copyright © 2018 Biophysical Society. All rights reserved.)

Water is widely assumed to be essential for protein dynamics and function. In particular, the well-documented "dynamical" transition at ∼200  K, at which the protein changes from a rigid, nonfunctional form to a flexible, functional state, as detected in hydrogenated protein by incoherent neutron scattering, requires hydration. Here, we report on coherent neutron scattering experiments on perdeuterated proteins and reveal that a transition occurs in dry proteins at the same temperature resulting primarily from the collective heavy-atom motions. The dynamical transition discovered is intrinsic to the energy landscape of dry proteins.

Networks are powerful tools for the presentation and analysis of interactions in multi-component systems. A commonly studied mesoscopic feature of networks is their community structure, which arises from grouping together similar nodes into one community and dissimilar nodes into separate communities. Here, the community structure of protein sequence similarity networks is determined with a new method: Attribute Clustering Dependent Communities (ACDC). Sequence similarity has hitherto typically been quantified by the alignment score or its expectation value. However, pair alignments with the same score or expectation value cannot thus be differentiated. To overcome this deficiency, the method constructs, for pair alignments, an extended alignment metric, the link attribute vector, which includes the score and other alignment characteristics. Rescaling components of the attribute vectors qualitatively identifies a systematic variation of sequence similarity within protein superfamilies. The problem of community detection is then mapped to clustering the link attribute vectors, selection of an optimal subset of links and community structure refinement based on the partition density of the network. ACDC-predicted communities are found to be in good agreement with gold standard sequence databases for which the "ground truth" community structures (or families) are known. ACDC is therefore a community detection method for sequence similarity networks based entirely on pair similarity information. A serial implementation of ACDC is available from https://cmb.ornl.gov/resources/developments.

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., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

Large-scale ensemble docking is investigated using five proteins from the Directory of Useful Decoys (DUD, dud.docking.org ) for which docking to crystal structures has proven difficult. Molecular dynamics trajectories are produced for each protein and an ensemble of representative conformational structures extracted from the trajectories. Docking calculations are performed on these selected simulation structures and ensemble-based enrichment factors compared with those obtained using docking in crystal structures of the same protein targets or random selection of compounds. Simulation-derived snapshots are found with improved enrichment factors that increase the chemical diversity of docking hits for four of the five selected proteins. A combination of all the docking results obtained from molecular dynamics simulation followed by selection of top-ranking compounds appears to be an effective strategy for increasing the number and diversity of hits when using docking to screen large libraries of chemicals against difficult protein targets.

Proteins display characteristic dynamical signatures that appear to be universal across all proteins regardless of topology and size. Here, we systematically characterize the universal features of fast side chain motions in proteins by examining the conformational energy surfaces of individual residues obtained using enhanced sampling molecular dynamics simulation (618 free energy surfaces obtained from 0.94 μs MD simulation). The side chain conformational free energy surfaces obtained using the adaptive biasing force (ABF) method for a set of eight proteins with different molecular weights and secondary structures are used to determine the methyl axial NMR order parameters (O(axis)(2)), populations of side chain rotamer states (ρ), conformational entropies (S(conf)), probability fluxes, and activation energies for side chain inter-rotameric transitions. The free energy barriers separating side chain rotamer states range from 0.3 to 12 kcal/mol in all proteins and follow a trimodal distribution with an intense peak at ~5 kcal/mol and two shoulders at ~3 and ~7.5 kcal/mol, indicating that some barriers are more favored than others by proteins to maintain a balance between their conformational stability and flexibility. The origin and the influences of the trimodal barrier distribution on the distribution of O(axis)(2) and the side chain conformational entropy are discussed. A hierarchical grading of rotamer states based on the conformational free energy barriers, entropy, and probability flux reveals three distinct classes of side chains in proteins. A unique nonlinear correlation is established between O(axis)(2) and the side chain rotamer populations (ρ). The apparent universality in O(axis)(2) versus ρ correlation, trimodal barrier distribution, and distinct characteristics of three classes of side chains observed among all proteins indicates a hidden regularity (or commonality) in the dynamical heterogeneity of fast side chain motions in proteins.

The relative motion of structural domains is essential for the biological function of many proteins. Here, by analyzing neutron scattering data and performing molecular dynamics simulations, we find that interdomain motion in several proteins obeys the principle of de Gennes narrowing, in which the wave vector dependence of the interdomain diffusion coefficient is inversely proportional to the interdomain structure factor. Thus, the rate of interdomain motion is inversely proportional to the probability distribution of the spatial configurations of domains.

We propose a method for obtaining the intrinsic, long-time mean square displacement (MSD) of atoms and molecules in proteins from finite-time molecular dynamics (MD) simulations. Typical data from simulations are limited to times of 1 to 10 ns, and over this time period the calculated MSD continues to increase without a clear limiting value. The proposed method consists of fitting a model to MD simulation-derived values of the incoherent intermediate neutron scattering function, I(inc)(Q,t), for finite times. The infinite-time MSD, , appears as a parameter in the model and is determined by fits of the model to the finite-time I(inc)(Q,t). Specifically, the is defined in the usual way in terms of the Debye-Waller factor as I(Q,t=∞)=exp(-Q(2)/3). The method is illustrated by obtaining the intrinsic MSD of hydrated lysozyme powder (h=0.4 g water/g protein) over a wide temperature range. The intrinsic obtained from data out to 1 and to 10 ns is found to be the same. The intrinsic is approximately twice the value of the MSD that is reached in simulations after times of 1 ns which correspond to those observed using neutron instruments that have an energy resolution width of 1 μeV.

The temperature dependences of the nanosecond dynamics of different chemical classes of amino acid residue have been analyzed by combining elastic incoherent neutron scattering experiments with molecular dynamics simulations on cytochrome P450cam. At T = 100-160 K, anharmonic motion in hydrophobic and aromatic residues is activated, whereas hydrophilic residue motions are suppressed because of hydrogen-bonding interactions. In contrast, at T = 180-220 K, water-activated jumps of hydrophilic side chains, which are strongly coupled to the relaxation rates of the hydrogen bonds they form with hydration water, become apparent. Thus, with increasing temperature, first the hydrophobic core awakens, followed by the hydrophilic surface.

A message passing interface (MPI)-based implementation (Autodock4.lga.MPI) of the grid-based docking program Autodock4 has been developed to allow simultaneous and independent docking of multiple compounds on up to thousands of central processing units (CPUs) using the Lamarkian genetic algorithm. The MPI version reads a single binary file containing precalculated grids that represent the protein-ligand interactions, i.e., van der Waals, electrostatic, and desolvation potentials, and needs only two input parameter files for the entire docking run. In comparison, the serial version of Autodock4 reads ASCII grid files and requires one parameter file per compound. The modifications performed result in significantly reduced input/output activity compared with the serial version. Autodock4.lga.MPI scales up to 8192 CPUs with a maximal overhead of 16.3%, of which two thirds is due to input/output operations and one third originates from MPI operations. The optimal docking strategy, which minimizes docking CPU time without lowering the quality of the database enrichments, comprises the docking of ligands preordered from the most to the least flexible and the assignment of the number of energy evaluations as a function of the number of rotatable bounds. In 24 h, on 8192 high-performance computing CPUs, the present MPI version would allow docking to a rigid protein of about 300K small flexible compounds or 11 million rigid compounds.

The REACH (realistic extension algorithm via covariance Hessian) coarse-grained biomolecular simulation method is a self-consistent multiscale approach directly mapping atomistic molecular dynamics simulation results onto a residue-scale model. Here, REACH is applied to calculate the dynamics of protein-protein interactions. The intra- and intermolecular fluctuations and the intermolecular vibrational densities of states derived from atomistic molecular dynamics are well reproduced by the REACH normal modes. The phonon dispersion relations derived from the REACH lattice dynamics model of crystalline ribonuclease A are also in satisfactory agreement with the corresponding all-atom results. The REACH model demonstrates that increasing dimer interaction strength decreases the translational and rotational intermolecular vibrational amplitudes, while their vibrational frequencies are relatively unaffected. A comparative study of functionally interacting biological dimers with crystal dimers, which are formed artificially via crystallization, reveals a relation between their static structures and the interprotein dynamics: i.e., the consequence of the extensive interfaces of biological dimers is reduction of the intermonomer translational and rotational amplitudes, but not the frequencies.

In the instantaneous normal mode method, normal mode analysis is performed at instantaneous configurations of a condensed-phase system, leading to modes with negative eigenvalues. These negative modes provide a means of characterizing local anharmonicities of the potential energy surface. Here, we apply instantaneous normal mode to analyze temperature-dependent diffusive dynamics in molecular dynamics simulations of a small protein (a scorpion toxin). Those characteristics of the negative modes are determined that correlate with the dynamical (or glass) transition behavior of the protein, as manifested as an increase in the gradient with T of the average atomic mean-square displacement at approximately 220 K. The number of negative eigenvalues shows no transition with temperature. Further, although filtering the negative modes to retain only those with eigenvectors corresponding to double-well potentials does reveal a transition in the hydration water, again, no transition in the protein is seen. However, additional filtering of the protein double-well modes, so as to retain only those that, on energy minimization, escape to different regions of configurational space, finally leads to clear protein dynamical transition behavior. Partial minimization of instantaneous configurations is also found to remove nondiffusive imaginary modes. In summary, examination of the form of negative instantaneous normal modes is shown to furnish a physical picture of local diffusive dynamics accompanying the protein glass transition.

Coarse graining of protein interactions provides a means of simulating large biological systems. The REACH (Realistic Extension Algorithm via Covariance Hessian) coarse-graining method, in which the force constants of a residue-scale elastic network model are calculated from the variance-covariance matrix obtained from atomistic molecular dynamics (MD) simulation, involves direct mapping between scales without the need for iterative optimization. Here, the transferability of the REACH force field is examined between protein molecules of different structural classes. As test cases, myoglobin (all alpha), plastocyanin (all beta), and dihydrofolate reductase (alpha/beta) are taken. The force constants derived are found to be closely similar in all three proteins. An MD version of REACH is presented, and low-temperature coarse-grained (CG) REACH MD simulations of the three proteins are compared with atomistic MD results. The mean-square fluctuations of the atomistic MD are well reproduced by the CGMD. Model functions for the CG interactions, derived by averaging over the three proteins, are also shown to produce fluctuations in good agreement with the atomistic MD. The results indicate that, similarly to the use of atomistic force fields, it is now possible to use a single, generic REACH force field for all protein studies, without having first to derive parameters from atomistic MD simulation for each individual system studied. The REACH method is thus likely to be a reliable way of determining spatiotemporal motion of a variety of proteins without the need for expensive computation of long atomistic MD simulations.

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., ((c) 2006 Wiley-Liss, Inc.)

The chi (C(alpha)-C(beta)) torsional barrier in the dipeptide alanine (N-methyl-l-alanyl-N-methylamide) crystal was investigated using ab initio calculations at various levels of theory, molecular mechanics, and molecular dynamics. For one of the two molecules in the asymmetric unit the calculations suggest that rotation around the chi dihedral angle is catalyzed by the crystal environment, reducing by up to approximately 2kT the torsional barrier in the crystal with respect to that in the gas phase. This catalytic effect is present at both low and room temperature and originates from a van der Waals destabilization of the minima in the methyl dihedral potential coming from the nonbonded environment of the side chain. Screening of a subset of the Protein Data Bank with a pharmacophore model reproducing the crystal environment around this side chain methyl identified a protein containing an alanine residue with an environment similar to that in the crystal. Calculations indicate that this chi torsional barrier is also reduced in the protein at low temperature but not at room temperature. This suggests that environment-catalyzed rotation of methyl groups can occur both in the solid phase and in native biological structures, though this effect might be temperature-dependent. The relevance of this catalytic effect is discussed in terms of its natural occurrence and its possible contribution to the low-frequency vibrational modes of molecules.

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 Tyr 35 ring flip in the bovine pancreatic trypsin inhibitor., ((c) 2005 Wiley Periodicals, Inc. J Comput Chem 26: 1359-1371, 2005.)

The modification of internal vibrational modes in a protein due to intraprotein anharmonicity and solvation effects is determined by performing molecular dynamics (MD) simulations of myoglobin, analyzing them using a Langevin model of the vibrational dynamics and comparing the Langevin results to a harmonic, normal mode model of the protein in vacuum. The diagonal and off-diagonal Langevin friction matrix elements, which model the roughness of the vibrational potential energy surfaces, are determined together with the vibrational potentials of mean force from the MD trajectories at 120 K and 300 K in vacuum and in solution. The frictional properties are found to be describable using simple phenomenological functions of the mode frequency, the accessible surface area, and the intraprotein interaction (the displacement vector overlap of any given mode with the other modes in the protein). The frictional damping of a vibrational mode in vacuum is found to be directly proportional to the intraprotein interaction of the mode, whereas in solution, the friction is proportional to the accessible surface area of the mode. In vacuum, the MD frequencies are lower than those of the normal modes, indicating intramolecular anharmonic broadening of the associated potential energy surfaces. Solvation has the opposite effect, increasing the large-amplitude vibrational frequencies relative to in vacuum and thus vibrationally confining the protein atoms. Frictional damping of the low-frequency modes is highly frequency dependent. In contrast to the damping effect of the solvent, the vibrational frequency increase due to solvation is relatively temperature independent, indicating that it is primarily a structural effect. The MD-derived vibrational dynamic structure factor and density of states are well reproduced by a model in which the Langevin friction and potential of mean force parameters are applied to the harmonic normal modes.

Proteins undergo an apparent dynamical transition on temperature variation that has been correlated with the onset of function. The transition in the mean-square displacement, , that is observed using a spectrometer or computer simulation, depends on the relationship between the timescales of the relaxation processes activated and the timescale accessible to the instrument or simulation. Models are described of two extreme situations---an "equilibrium" model, in which the long-time dynamics changes with temperature and all motions are resolved by the instrument used; and a "frequency window" model, in which there is no change in the long-time dynamics but as the temperature increases, the relaxation frequencies move into the instrumental range. Here we demonstrate that the latter, frequency-window model can describe the temperature and timescale dependences of both the intermediate neutron scattering function and derived from molecular dynamics simulations of a small protein in a cryosolution. The frequency-window model also describes the energy-resolution and temperature-dependences of obtained from experimental neutron scattering on glutamate dehydrogenase in the same solvent. Although equilibrium effects should also contribute to dynamical transitions in proteins, the present results suggests that frequency-window effects can play a role in the simulations and experiments examined. Finally, misquotations of previous findings are discussed in the context of solvent activation of protein dynamics and the possible relationship of this to activity.

The apparent simplicity of the water molecule belies the wide range of fascinating protein phenomena in which it participates. We review recent computer simulation work on buried, internal water molecules, discussing the thermodynamics of water molecule binding and the participation of water in proton transfer reactions. Surface water molecules are also considered, with emphasis on the modification of average solvent structure on a protein surface, the role of water in the protein dynamical 'glass' transition and a simplified description of the protein motions thereby activated.

The change in the vibrational density of states of a protein (dihydrofolate reductase) on binding a ligand (methotrexate) is determined using inelastic neutron scattering. The vibrations of the complex soften significantly relative to the unbound protein. The resulting free-energy change, which is directly determined by the density of states change, is found to contribute significantly to the binding equilibrium.

Experimental and computer simulation studies have revealed the presence of a glass-like transition in the internal dynamics of hydrated proteins at approximately 200 K involving an increase of the amplitude of anharmonic dynamics. This increase in flexibility has been correlated with the onset of protein activity. Here, we determine the driving force behind the protein transition by performing molecular dynamics simulations of myoglobin surrounded by a shell of water. A dual heat bath method is used with which, in any given simulation, the protein and solvent are held at different temperatures, and sets of simulations are performed varying the temperature of the components. The results show that the protein transition is driven by a dynamical transition in the hydration water that induces increased fluctuations primarily in side chains in the external regions of the protein. The water transition involves activation of translational diffusion and occurs even in simulations where the protein atoms are held fixed.

Solvent plays an important role in modulating internal motions of proteins. Here we present a computational method for including solvent effects on charge-charge interactions and on pathways between functional protein conformations, and examine solvent effects on equilibrium internal fluctuations in proteins. A computationally efficient charge reparametrisation method is presented that satisfactorily reproduces the electrostatic interactions present in a full continuum Poisson-Boltzmann representation. The application of charge reparametrisation in the calculation of a large-scale conformational transition pathway in a protein, annexin V, is illustrated. We also examine solvent effects on fast (picosecond timescale) internal protein dynamics. Nosé-Hoover dual heatbath molecular dynamics simulations are performed. These simulations allow the solvent region to be fixed at one temperature and the protein at another. The results of the Nosé-Hoover simulations on hydrated myoglobin confirm that the solvent temperature strongly influences the protein fluctuations. We consider to what extent the solvent can be considered to determine the high temperature protein dynamics.

The temperature dependence of the internal dynamics of an isolated protein, bovine pancreatic trypsin inhibitor, is examined using normal mode analysis and molecular dynamics (MD) simulation. It is found that the protein exhibits marked anharmonic dynamics at temperatures of approximately 100-120 K, as evidenced by departure of the MD-derived average mean square displacement from that of the harmonic model. This activation of anharmonic dynamics is at lower temperatures than previously detected in proteins and is found in the absence of solvent molecules. The simulation data are also used to investigate neutron scattering properties. The effects are determined of instrumental energy resolution and of approximations commonly used to extract mean square displacement data from elastic scattering experiments. Both the presence of a distribution of mean square displacements in the protein and the use of the Gaussian approximation to the dynamic structure factor lead to quantified underestimation of the mean square displacement obtained.

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.

The urgent search for drugs to combat SARS-CoV-2 has included the use of supercomputers. The use of general-purpose graphical processing units (GPUs), massive parallelism, and new software for high-performance computing (HPC) has allowed researchers to search the vast chemical space of potential drugs faster than ever before. We developed a new drug discovery pipeline using the Summit supercomputer at Oak Ridge National Laboratory to help pioneer this effort, with new platforms that incorporate GPU-accelerated simulation and allow for the virtual screening of billions of potential drug compounds in days compared to weeks or months for their ability to inhibit SARS-COV-2 proteins. This effort will accelerate the process of developing drugs to combat the current COVID-19 pandemic and other diseases. [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]

Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes, play a major role in cell signaling, and are associated with human diseases. To understand IDP function it is critical to determine their configurational ensemble, i.e., the collection of 3-dimensional structures they adopt, and this remains an immense challenge in structural biology. Attempts to determine this ensemble computationally have been hitherto hampered by the necessity of reweighting molecular dynamics (MD) results or biasing simulation in order to match ensemble-averaged experimental observables, operations that reduce the precision of the generated model because different structural ensembles may yield the same experimental observable. Here, by employing enhanced sampling MD we reproduce the experimental small-angle neutron and X-ray scattering profiles and the NMR chemical shifts of the disordered N terminal (SH4UD) of c-Src kinase without reweighting or constraining the simulations. The unbiased simulation results reveal a weakly funneled and rugged free energy landscape of SH4UD, which gives rise to a heterogeneous ensemble of structures that cannot be described by simple polymer theory. SH4UD adopts transient helices, which are found away from known phosphorylation sites and could play a key role in the stabilization of structural regions necessary for phosphorylation. Our findings indicate that adequately sampledmolecular simulations can be performed to provide accurate physical models of flexible biosystems, thus rationalizing their biological function. [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]

An analytical approach is developed for reconstructingsite-specificmethyl-bearing protein side-chain conformational energy surfaces fromNMR methyl axis order parameters (Oaxis2). Application of an enhancedsampling algorithm (adaptive biasing force) to molecular dynamicssimulation of a protein, calcium-bound calmodulin, reveals a nonlinearcorrelation between Oaxis2and the populations of rotamer statesof protein side-chains, permitting the rotamer populations to be extracteddirectly from Oaxis2. The analytical approach yields side-chainconformational distributions that are in excellent agreement withthose obtained from the enhanced-sampling MD results. [ABSTRACT FROM AUTHOR]

An understanding of low-frequency, collective protein dynamics at low temperatures can furnish valuable information on functional protein energy landscapes, on the origins of the protein glass transition and on protein-protein interactions. Here, molecular dynamics (MD) simulations and normal-mode analyses are performed on various models of crystalline myoglobin in order to characterize intra- and interprotein vibrations at 150 K. Principal component analysis of the MD trajectories indicates that the Boson peak, a broad peak in the dynamic structure factor centered at about ∼2-2.5 meV, originates from ∼10² collective, harmonic vibrations. An accurate description of the environment is found to be essential in reproducing the experimental Boson peak form and position. At lower energies other strong peaks are found in the calculated dynamic structure factor. Characterization of these peaks shows that they arise from harmonic vibrations of proteins relative to each other. These vibrations are likely to furnish valuable information on the physical nature of protein-protein interactions. [ABSTRACT FROM AUTHOR]

The mechanics of a protein joint, gluing two elastic titin molecules together in muscle, are explored using computer simulation in this issue of Structure by . Maybe the glue itself also has another, more subtle, sensory role? [Copyright &y& Elsevier]

The relative motion of structural domains is essential for the biological function of many proteins. Here, by analyzing neutron scattering data and performing molecular dynamics simulations, we find that interdomain motion in several proteins obeys the principle of de Gennes narrowing, in which the wave vector dependence of the interdomain diffusion coefficient is inversely proportional to the interdomain structure factor. Thus, the rate of interdomain motion is inversely proportional to the probability distribution of the spatial configurations of domains. [ABSTRACT FROM AUTHOR]

Summary: Aligned α helix peptide dipoles sum to a “macroscopic” dipole parallel to the helix axis that has been implicated in protein folding and function. However, in aqueous solution the dipole is counteracted by an electrostatic reaction field generated by the solvent, and the strength of the helix dipole may reduce drastically from its value in vacuum. Here, using atomic-detail helix models and Poisson-Boltzmann continuum electrostatics calculations, the net effective dipole moment, μeff, is calculated. Some initially surprising results are found. Whereas in vacuum μeff increases with helix length, the opposite is found to be the case for transmembrane helices. In soluble proteins, μeff is found to vary strongly with the orientation and position of the helix relative to the aqueous medium. A set of rules is established to estimate of the strength of μeff from graphical inspection of protein structures. [Copyright &y& Elsevier]

The molecular mechanisms that transduce the osteoblast response to physical forces in the bone microenvironment are poorly understood. Here, we used genetic and pharmacological experiments to determine whether the polycystins PC1 and PC2 (encoded by Pkd1 and Pkd2) and the transcriptional coactivator TAZ form a mechanosensing complex in osteoblasts. Compound-heterozygous mice lacking 1 copy of Pkd1 and Taz exhibited additive decrements in bone mass, impaired osteoblast-mediated bone formation, and enhanced bone marrow fat accumulation. Bone marrow stromal cells and osteoblasts derived from these mice showed impaired osteoblastogenesis and enhanced adipogenesis. Increased extracellular matrix stiffness and application of mechanical stretch to multipotent mesenchymal cells stimulated the nuclear translocation of the PC1 C-terminal tail/TAZ (PC1-CTT/TAZ) complex, leading to increased runt-related transcription factor 2-mediated (Runx2-mediated) osteogenic and decreased PPARγ-dependent adipogenic gene expression. Using structure-based virtual screening, we identified a compound predicted to bind to PC2 in the PC1:PC2 C-terminal tail region with helix:helix interaction. This molecule stimulated polycystin- and TAZ-dependent osteoblastogenesis and inhibited adipogenesis. Thus, we show that polycystins and TAZ integrate at the molecular level to reciprocally regulate osteoblast and adipocyte differentiation, indicating that the polycystins/TAZ complex may be a potential therapeutic target to increase bone mass. [ABSTRACT FROM AUTHOR]

Neutrons scatter quasielastically from stochastic, diffusive processes, such as overdamped vibrations, localized diffusion and transitions between energy minima. In biological systems, such as proteins and membranes, these relaxation processes are of considerable physical interest. We review here recent methodological advances and applications of quasielastic neutron scattering (QENS) in biology, concentrating on the role of molecular dynamics simulation in generating data with which neutron profiles can be unambiguously interpreted. We examine the use of massively-parallel computers in calculating scattering functions, and the application of Markov state modeling. The decomposition of MD-derived neutron dynamic susceptibilities is described, and the use of this in combination with NMR spectroscopy. We discuss dynamics at very long times, including approximations to the infinite time mean-square displacement and nonequilibrium aspects of single-protein dynamics. Finally, we examine how neutron scattering and MD can be combined to provide information on lipid nanodomains. This article is part of a Special Issue entitled “Science for Life” Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo. [ABSTRACT FROM AUTHOR]

The lipid raff hypothesis presents insights into how the cell membrane organizes proteins and lipids to accomplish its many vital functions. Yet basic questions remain about the physical mechanisms that lead to the formation, stability, and size of lipid raffs. As a result, much interest has been generated in the study of systems that contain similar lateral heterogeneities, or domains. In the current work we present an experimental approach that is capable of isolating the bending moduli of lipid domains. This is accomplished using neutron scattering and its unique sensitivity to the isotopes of hydrogen. Combining contrast matching approaches with inelastic neutron scattering, we isolate the bending modulus of ∼13 nm diameter domains residing in 60 nm unilamellar vesicles, whose lipid composition mimics the mammalian plasma membrane outer leaflet. Importantly, the bending modulus of the nanoscopic domains differs from the modulus of the continuous phase surrounding them. From additional structural measurements and all-atom simulations, we also determine that nanoscopic domains are in-register across the bilayer leaflets. Taken together, these results inform a number of theoretical models of domain/raft formation and highlight the fact that mismatches in bending modulus must be accounted for when explaining the emergence of lateral heterogeneities in lipid systems and biological membranes. [ABSTRACT FROM AUTHOR]

The effect of surface hydration water on internal protein motion is of fundamental interest in molecular biophysics. Here, by decomposing the picosecond to nanosecond atomic motion in molecular dynamics simulations of lysozyme at different hydration levels into three components--localized single-well diffusion, methyl group rotation, and nonmethyl jumps--we show that the effect of surface hydration is mainly to increase the volume of the localized single-well diffusion. These diffusive motions are coupled in such a way that the hydration effect propagates from the protein surface into the dry core [ABSTRACT FROM AUTHOR]

Abstract: Recent measurements have demonstrated enzyme activity at hydrations as low as 3%. The question of whether the hydration-induced enzyme flexibility is important for activity is addressed by performing picosecond dynamic neutron scattering experiments on pig liver esterase powders at various temperatures as well as solutions. At all temperatures and hydrations investigated here, significant quasielastic scattering intensity is found in the protein, indicating the presence of anharmonic, diffusive motion. As the hydration increases a temperature-dependent dynamical transition appears and strengthens involving additional diffusive motion. At low temperature, increasing hydration resulted in lower flexibility of the enzyme. At higher temperatures, systems containing sufficient number of water molecules interacting with the protein exhibit increased flexibility. The implication of these results is that, although the additional hydration-induced diffusive motion and flexibility at high temperatures in the enzyme detected here may be related to increased activity, they are not required for the enzyme to function. [Copyright &y& Elsevier]

Water is widely assumed to be essential for protein dynamics and function. In particular, the well-documented "dynamical" transition at ~200 K, at which the protein changes from a rigid, nonfunctional form to a flexible, functional state, as detected in hydrogenated protein by incoherent neutron scattering, requires hydration. Here, we report on coherent neutron scattering experiments on perdeuterated proteins and reveal that a transition occurs in dry proteins at the same temperature resulting primarily from the collective heavy-atom motions. The dynamical transition discovered is intrinsic to the energy landscape of dry proteins. [ABSTRACT FROM AUTHOR]

Focuses on the movement of the molecules in proteins. Reason for the development of the neutron spin-echo technique; What is allowed by neutron spin-echo; Advantage of the spin-echo technique.