Computational Graphics Software for Interactive Docking and Visualization of Ligand-Protein Complementarity.

The Dockeye software is designed to complement automated docking protocols by allowing the user's chemical know-how and experience of what makes for good protein-ligand binding, knowledge that is not easily encoded into automated algorithms, to guide the docking. It allows the interactive manipulation of the ligand placement against a protein target. Real-time intuitively comprehensible feedback about the location, spatial density, and the extent of both favorable and unfavorable atomic interactions between ligand and protein is provided through a carefully designed graphical object. It is also a tool for the graphical analysis of the interactions of known protein-ligand complexes. Comparative docking of 58 protein-ligand complexes with Dockeye and Autodock Vina shows how this software can be used synergistically with automated docking programs to significantly improve the task of discovery of ligand placement.

DNA Conformational Changes Play a Force-Generating Role during Bacteriophage Genome Packaging.

Motors that move DNA, or that move along DNA, play essential roles in DNA replication, transcription, recombination, and chromosome segregation. The mechanisms by which these DNA translocases operate remain largely unknown. Some double-stranded DNA (dsDNA) viruses use an ATP-dependent motor to drive DNA into preformed capsids. These include several human pathogens as well as dsDNA bacteriophages-viruses that infect bacteria. We previously proposed that DNA is not a passive substrate of bacteriophage packaging motors but is instead an active component of the machinery. We carried out computational studies on dsDNA in the channels of viral portal proteins, and they reveal DNA conformational changes consistent with that hypothesis. dsDNA becomes longer ("stretched") in regions of high negative electrostatic potential and shorter ("scrunched") in regions of high positive potential. These results suggest a mechanism that electrostatically couples the energy released by ATP hydrolysis to DNA translocation: The chemical cycle of ATP binding, hydrolysis, and product release drives a cycle of protein conformational changes. This produces changes in the electrostatic potential in the channel through the portal, and these drive cyclic changes in the length of dsDNA as the phosphate groups respond to the protein's electrostatic potential. The DNA motions are captured by a coordinated protein-DNA grip-and-release cycle to produce DNA translocation. In short, the ATPase, portal, and dsDNA work synergistically to promote genome packaging.

Companion Simulations and Modeling to NMR-Based Dynamical Studies of Proteins.

NMR-based studies of protein dynamics and molecular simulations have a synergistic relationship. Molecular simulations, in combination with interpretative theoretical models, leverage the dynamical information obtained from NMR. They provide the concrete physical schema underlying the quantities measured by NMR, and help extend the range of applications beyond the strictly dynamic properties. NMR data in turn provide concrete data to test and improve the potential functions used for simulation of dynamics of proteins. The concept of time correlation functions is central to the understanding of many dynamical processes. Their evaluation through atomistic simulations is discussed, with application to different dynamical quantities measured by NMR. While advances in computers have made such atomistic simulations almost routine, the companion use of simple interpretive models is stressed, to provide not just numbers but physical insight.

Water loading driven size, shape, and composition of cetyltrimethylammonium/hexanol/pentane reverse micelles.

Cetyltrimethylammonium bromide (CTAB)/hexanol reverse micelles have found a variety of applications that demand control over physical parameters. Water content or loading is among the most basic tunable components and is the major driver of the physical properties of these systems. This study uses small-angle scattering with contrast variation to characterize these systems as a function of water loading. The scattering data were analyzed with a variety of approaches, resulting in converging physical specifications. Equations that describe basic physical parameters were determined that allow for characterization and manipulation of the CTAB/hexanol reverse micelle surfactant system. The shape of the reverse micelles was revealed to be slightly ellipsoidal and varies slightly through the water loading range. The surfactant shell is shown to contain a higher fraction of hexanol upon addition of water. Analysis reveals that the size, shape, and surfactant/cosurfactant composition are directly tunable by variation of the water content and that these properties are consequences of the balance of forces present in the reverse micelles.

Measuring Entropy in Molecular Recognition by Proteins.

Molecular recognition by proteins is fundamental to the molecular basis of biology. Dissection of the thermodynamic landscape governing protein-ligand interactions has proven difficult because determination of various entropic contributions is quite challenging. Nuclear magnetic resonance relaxation measurements, theory, and simulations suggest that conformational entropy can be accessed through a dynamical proxy. Here, we review the relationship between measures of fast side-chain motion and the underlying conformational entropy. The dynamical proxy reveals that the contribution of conformational entropy can range from highly favorable to highly unfavorable and demonstrates the potential of this key thermodynamic variable to modulate protein-ligand interactions. The dynamical so-called entropy meter also refines the role of solvent entropy and directly determines the loss in rotational-translational entropy that occurs upon formation of high-affinity complexes. The ability to quantify the roles of entropy through an entropy meter based on measurable dynamical properties promises to highlight its role in protein function.

Entropy in molecular recognition by proteins.

Molecular recognition by proteins is fundamental to molecular biology. Dissection of the thermodynamic energy terms governing protein-ligand interactions has proven difficult, with determination of entropic contributions being particularly elusive. NMR relaxation measurements have suggested that changes in protein conformational entropy can be quantitatively obtained through a dynamical proxy, but the generality of this relationship has not been shown. Twenty-eight protein-ligand complexes are used to show a quantitative relationship between measures of fast side-chain motion and the underlying conformational entropy. We find that the contribution of conformational entropy can range from favorable to unfavorable, which demonstrates the potential of this thermodynamic variable to modulate protein-ligand interactions. For about one-quarter of these complexes, the absence of conformational entropy would render the resulting affinity biologically meaningless. The dynamical proxy for conformational entropy or "entropy meter" also allows for refinement of the contributions of solvent entropy and the loss in rotational-translational entropy accompanying formation of high-affinity complexes. Furthermore, structure-based application of the approach can also provide insight into long-lived specific water-protein interactions that escape the generic treatments of solvent entropy based simply on changes in accessible surface area. These results provide a comprehensive and unified view of the general role of entropy in high-affinity molecular recognition by proteins.

On the ability of molecular dynamics force fields to recapitulate NMR derived protein side chain order parameters.

Molecular dynamics (MD) simulations have become a central tool for investigating various biophysical questions with atomistic detail. While many different proxies are used to qualify MD force fields, most are based on largely structural parameters such as the root mean square deviation from experimental coordinates or nuclear magnetic resonance (NMR) chemical shifts and residual dipolar couplings. NMR derived Lipari-Szabo squared generalized order parameter (O(2) ) values of amide NH bond vectors of the polypeptide chain were also often employed for refinement and validation. However, with a few exceptions, side chain methyl symmetry axis order parameters have not been incorporated into experimental reference sets. Using a test set of five diverse proteins, the performance of several force fields implemented in the NAMDD simulation package was examined. It was found that simulations employing explicit water implemented using the TIP3 model generally performed significantly better than those using implicit water in reproducing experimental methyl symmetry axis O(2) values. Overall the CHARMM27 force field performs nominally better than two implementations of the Amber force field. It appeared that recent quantum mechanics modifications to side chain torsional angles of leucine and isoleucine in the Amber force field have significantly hindered proper motional modeling for these residues. There remained significant room for improvement as even the best correlations of experimental and simulated methyl group Lipari-Szabo generalized order parameters fall below an R(2) of 0.8.

Characterization of Cetyltrimethylammonium Bromide/Hexanol Reverse Micelles by Experimentally Benchmarked Molecular Dynamics Simulations.

Encapsulation of small molecules, proteins, and other macromolecules within the protective water core of reverse micelles is emerging as a powerful strategy for a variety of applications. The cationic surfactant cetyltrimethylammonium bromide (CTAB) in combination with hexanol as a cosurfactant is particularly useful in the context of solution NMR spectroscopy of encapsulated proteins. Small-angle X-ray and neutron scattering is employed to investigate the internal structure of the CTAB/hexanol reverse micelle particle under conditions appropriate for high-resolution NMR spectroscopy. The scattering profiles are used to benchmark extensive molecular dynamics simulations of this reverse micelle system and indicate that the parameters used in these simulations recapitulate experimental results. Scattering profiles and simulations indicate formation of homogeneous solutions of small approximately spherical reverse micelle particles at a water loading of 20 composed of ∼150 CTAB and 240 hexanol molecules. The 3000 waters comprising the reverse micelle core show a gradient of translational diffusion that reaches that of bulk water at the center. Rotational diffusion is slowed relative to bulk throughout the water core, with the greatest slowing near the CTAB headgroups. The 5 Šthick interfacial region of the micelle consists of overlapping layers of Br(-) enriched water, CTAB headgroups, and hexanol hydroxyl groups, containing about one-third of the total water. This study employs well-parametrized MD simulations, X-ray and neutron scattering, and electrostatic theory to illuminate fundamental properties of CTAB/hexanol reverse micelle size, shape, partitioning, and water behavior.

Regulation of brain glutamate metabolism by nitric oxide and S-nitrosylation.

Nitric oxide (NO) is a signaling intermediate during glutamatergic neurotransmission in the central nervous system (CNS). NO signaling is in part accomplished through cysteine S-nitrosylation, a posttranslational modification by which NO regulates protein function and signaling. In our investigation of the protein targets and functional impact of S-nitrosylation in the CNS under physiological conditions, we identified 269 S-nitrosocysteine residues in 136 proteins in the wild-type mouse brain. The number of sites was significantly reduced in the brains of mice lacking endothelial nitric oxide synthase (eNOS(-/-)) or neuronal nitric oxide synthase (nNOS(-/-)). In particular, nNOS(-/-) animals showed decreased S-nitrosylation of proteins that participate in the glutamate/glutamine cycle, a metabolic process by which synaptic glutamate is recycled or oxidized to provide energy. (15)N-glutamine-based metabolomic profiling and enzymatic activity assays indicated that brain extracts from nNOS(-/-) mice converted less glutamate to glutamine and oxidized more glutamate than those from mice of the other genotypes. GLT1 [also known as EAAT2 (excitatory amino acid transporter 2)], a glutamate transporter in astrocytes, was S-nitrosylated at Cys(373) and Cys(561) in wild-type and eNOS(-/-) mice, but not in nNOS(-/-) mice. A form of rat GLT1 that could not be S-nitrosylated at the equivalent sites had increased glutamate uptake compared to wild-type GLT1 in cells exposed to an S-nitrosylating agent. Thus, NO modulates glutamatergic neurotransmission through the selective, nNOS-dependent S-nitrosylation of proteins that govern glutamate transport and metabolism.

Analysis of the size dependence of macromolecular crowding shows that smaller is better.

The aqueous milieu inside cells contains as much as 30-40% dissolved protein and RNA by volume. This large concentration of macromolecules is expected to cause significant deviations from solution ideality. In vivo biochemical reaction rates and equilibria might differ significantly from those measured in the majority of in vitro experiments that are performed at much lower macromolecule concentrations. Consequently crowding, a nonspecific phenomenon believed to arise from the large excluded volume of these macromolecules, has been studied extensively by experimental and theoretical methods. However, the relevant theory has not been applied consistently. When the steric effects of macromolecular crowders and small molecules like water and ions are treated on an equal footing, the effect of the macromolecules is opposite to that commonly believed. Large molecules are less effective at crowding than water and ions. There is also a surprisingly weak dependence on crowder size. Molecules of medium size, ∼5 Šradius, have the same effect as much larger macromolecules like proteins and RNA. These results require a reassessment of observed high-concentration effects and of strategies to mimic in vivo conditions with in vitro experiments.

On the relationship between NMR-derived amide order parameters and protein backbone entropy changes.

Molecular dynamics simulations are used to analyze the relationship between NMR-derived squared generalized order parameters of amide NH groups and backbone entropy. Amide order parameters (O(2) NH ) are largely determined by the secondary structure and average values appear unrelated to the overall flexibility of the protein. However, analysis of the more flexible subset (O(2) NH < 0.8) shows that these report both on the local flexibility of the protein and on a different component of the conformational entropy than that reported by the side chain methyl axis order parameters, O(2) axis . A calibration curve for backbone entropy vs. O(2) NH is developed, which accounts for both correlations between amide group motions of different residues, and correlations between backbone and side chain motions. This calibration curve can be used with experimental values of O(2) NH changes obtained by NMR relaxation measurements to extract backbone entropy changes, for example, upon ligand binding. In conjunction with our previous calibration for side chain entropy derived from measured O(2) axis values this provides a prescription for determination of the total protein conformational entropy changes from NMR relaxation measurements.

A sharp thermal transition of fast aromatic-ring dynamics in ubiquitin.

Aromatic amino acid side chains have a rich role within proteins and are often central to their structure and function. Suitable isotopic-labelling strategies enable studies of sub-nanosecond aromatic-ring dynamics using solution NMR relaxation methods. Surprisingly, it was found that the three aromatic side chains in human ubiquitin show a sharp thermal dynamical transition at approximately 312 K. Hydrostatic pressure has little effect on the low-temperature behavior, but somewhat decreases the amplitude of motion in the high-temperature regime. Therefore, below the transition temperature, ring motion is largely librational. Above this temperature, a complete ring-rotation process that is fully consistent with a continuous diffusion not requiring the transient creation of a large activated free volume occurs. Molecular dynamics simulations qualitatively corroborate this view and reinforce the notion that the dynamical character of the protein interior has much more liquid-alkane-like properties than previously appreciated.

The remarkable hydration of the antifreeze protein Maxi: a computational study.

The long four-helix bundle antifreeze protein Maxi contains an unusual core for a globular protein. More than 400 ordered waters between the helices form a nano-pore of internal water about 150 Å long. Molecular dynamics simulations of hydrated Maxi were carried out using the CHARMM27 protein forcefield and the TIP3P water model. Solvation in the core and non-core first hydration shell was analyzed in terms of water-water H-bond distance-angle distributions. The core had an increased population of low-angle H-bonds between water pairs relative to bulk water. Enhancement of low angle H-bonds was particularly pronounced for water pairs at the interfaces between apolar and polar regions inside the protein core, characteristic of the anchored clathrate solvation structure seen previously in the ice-nuclei binding surfaces of type I, type III, and beta-helical antifreeze proteins. Anchored clathrate type solvation structure was not seen in the exterior solvation shell except around residues implicated in ice binding. Analysis of solvation dynamics using water residence times and diffusion constants showed that exterior solvation shell waters exchanged rapidly with bulk water, with no difference between ice-binding and non-binding residues. Core waters had about ten-fold slower diffusion than bulk water. Water residence times around core residues averaged about 8 ps, similar to those on the exterior surface, but they tended to exchange primarily with other core water, resulting in longer, >40 ps residence times within the core. Preferential exchange or diffusion of the water along the long axis of the water core of Maxi was not detected.

TCR triggering by pMHC ligands tethered on surfaces via poly(ethylene glycol) depends on polymer length.

Antigen recognition by T cells relies on the interaction between T cell receptor (TCR) and peptide-major histocompatibility complex (pMHC) at the interface between the T cell and the antigen presenting cell (APC). The pMHC-TCR interaction is two-dimensional (2D), in that both the ligand and receptor are membrane-anchored and their movement is limited to 2D diffusion. The 2D nature of the interaction is critical for the ability of pMHC ligands to trigger TCR. The exact properties of the 2D pMHC-TCR interaction that enable TCR triggering, however, are not fully understood. Here, we altered the 2D pMHC-TCR interaction by tethering pMHC ligands to a rigid plastic surface with flexible poly(ethylene glycol) (PEG) polymers of different lengths, thereby gradually increasing the ligands' range of motion in the third dimension. We found that pMHC ligands tethered by PEG linkers with long contour length were capable of activating T cells. Shorter PEG linkers, however, triggered TCR more efficiently. Molecular dynamics simulation suggested that shorter PEGs exhibit faster TCR binding on-rates and off-rates. Our findings indicate that TCR signaling can be triggered by surface-tethered pMHC ligands within a defined 3D range of motion, and that fast binding rates lead to higher TCR triggering efficiency. These observations are consistent with a model of TCR triggering that incorporates the dynamic interaction between T cell and antigen-presenting cell.

Banding of NMR-derived methyl order parameters: implications for protein dynamics.

Our understanding of protein folding, stability, and function has begun to more explicitly incorporate dynamical aspects. Nuclear magnetic resonance has emerged as a powerful experimental method for obtaining comprehensive site-resolved insight into protein motion. It has been observed that methyl-group motion tends to cluster into three "classes" when expressed in terms of the popular Lipari-Szabo model-free squared generalized order parameter. Here the origins of the three classes or bands in the distribution of order parameters are examined. As a first step, a Bayesian based approach, which makes no a priori assumption about the existence or number of bands, is developed to detect the banding of Oaxis2 values derived either from NMR experiments or molecular dynamics simulations. The analysis is applied to seven proteins with extensive molecular dynamics simulations of these proteins in explicit water to examine the relationship between O2 and fine details of the motion of methyl bearing side chains. All of the proteins studied display banding, with some subtle differences. We propose a very simple yet plausible physical mechanism for banding. Finally, our Bayesian method is used to analyze the measured distributions of methyl group motions in the catabolite activating protein and several of its mutants in various liganded states and discuss the functional implications of the observed banding to protein dynamics and function.

Microscopic insights into the NMR relaxation-based protein conformational entropy meter.

Conformational entropy is a potentially important thermodynamic parameter contributing to protein function. Quantitative measures of conformational entropy are necessary for an understanding of its role but have been difficult to obtain. An empirical method that utilizes changes in conformational dynamics as a proxy for changes in conformational entropy has recently been introduced. Here we probe the microscopic origins of the link between conformational dynamics and conformational entropy using molecular dynamics simulations. Simulation of seven proteins gave an excellent correlation with measures of side-chain motion derived from NMR relaxation. The simulations show that the motion of methyl-bearing side chains are sufficiently coupled to that of other side chains to serve as excellent reporters of the overall side-chain conformational entropy. These results tend to validate the use of experimentally accessible measures of methyl motion--the NMR-derived generalized order parameters--as a proxy from which to derive changes in protein conformational entropy.

Calculation of configurational entropy with a Boltzmann-quasiharmonic model: the origin of high-affinity protein-ligand binding.

Accurate assessment of configurational entropy remains a large challenge in biology. While many methods exist to calculate configurational entropy, there is a balance between accuracy and computational demands. Here we calculate ligand and protein conformational entropies using the Boltzmann-quasiharmonic (BQH) method, which treats the first-order entropy term by the Boltzmann expression for entropy while determining correlations using the quasiharmonic model. This method is tested by comparison with the exact Clausius expression for entropy on a range of test molecules ranging from small ligands to a protein. Using the BQH method, we then analyze the rotational and translational (R/T) entropy change upon ligand binding for five protein complexes to explore the origins of extremely tight affinity. The results suggest that in these systems such affinity is achieved by a combination of simultaneously maintaining good protein-ligand contacts while allowing significant residual R/T motion of the ligand through suitable protein motions.

Allostery in the lac operon: population selection or induced dissociation?

Allostery, the modulation of function of a protein at one site by the binding of a ligand at a different site, is a property of many proteins. Two kinetically distinct models have been proposed: i) The induced fit model in which the ligand binds to the protein and then induces the conformational change. ii) The population selection model, in which the protein spontaneously undergoes a conformational change, which is then 'captured' by the ligand. Using measured kinetic constants for the lac repressor the contribution of population selection vs. induced dissociation is quantified by simulating the kinetics of allostery. At very low inducer concentration, both mechanisms contribute significantly. Total induction, though, is small under these conditions. At increasing levels of induction the induced dissociation mechanism soon dominates, first due to binding of one inducer, and then from two inducers binding.