Search

X-ray free-electron lasers have enabled femtosecond protein nanocrystallography, a novel method to determine the structure of proteins. It allows time-resolved imaging of nanocrystals that are too small for conventional crystallography. The short pulse duration helps in overcoming the detrimental effects of radiation damage because x rays are scattered before the sample has been significantly altered. It has been suggested that, fortuitously, the diffraction process self-terminates abruptly once radiation damage destroys the crystalline order. Our calculations show that high-intensity x-ray pulses indeed trigger a cascade of damage processes in ferredoxin crystals, a particular metalloprotein of interest. However, we found that the damage process is initially not completely random. Correlations exist among the protein monomers, so that Bragg diffraction still occurs in the damaged crystals, despite significant atomic displacements. Our results show that the damage process is reproducible to a certain degree, which is potentially beneficial for the orientation step in single-molecule imaging.

Late-stage or post-market identification of adverse drug reactions (ADRs) is a significant public health issue and a source of major economic liability for drug development. Thus, reliable in silico screening of drug candidates for possible ADRs would be advantageous. In this work, we introduce a computational approach that predicts ADRs by combining the results of molecular docking and leverages known ADR information from DrugBank and SIDER. We employed a recently parallelized version of AutoDock Vina (VinaLC) to dock 906 small molecule drugs to a virtual panel of 409 DrugBank protein targets. L1-regularized logistic regression models were trained on the resulting docking scores of a 560 compound subset from the initial 906 compounds to predict 85 side effects, grouped into 10 ADR phenotype groups. Only 21% (87 out of 409) of the drug-protein binding features involve known targets of the drug subset, providing a significant probe of off-target effects. As a control, associations of this drug subset with the 555 annotated targets of these compounds, as reported in DrugBank, were used as features to train a separate group of models. The Vina off-target models and the DrugBank on-target models yielded comparable median area-under-the-receiver-operating-characteristic-curves (AUCs) during 10-fold cross-validation (0.60-0.69 and 0.61-0.74, respectively). Evidence was found in the PubMed literature to support several putative ADR-protein associations identified by our analysis. Among them, several associations between neoplasm-related ADRs and known tumor suppressor and tumor invasiveness marker proteins were found. A dual role for interstitial collagenase in both neoplasms and aneurysm formation was also identified. These associations all involve off-target proteins and could not have been found using available drug/on-target interaction data. This study illustrates a path forward to comprehensive ADR virtual screening that can potentially scale with increasing number of CPUs to tens of thousands of protein targets and millions of potential drug candidates.

Rarely is any solution simply solute and water. In vivo, solutes, such as proteins and nucleic acids, swim in a sea of water, salts, ions, small molecules, and lipids, not to mention other macromolecules. In vitro, virtually all solutions contain a mixture of aqueous solvents, or "cosolvents" [i.e., solvent(s) in addition to water], that can alter the dynamics, behavior, solubility, and stability of proteins and nucleic acids. We have developed models for a number of cosolvents, including the denaturant urea and the small chemical chaperone trimethylamine N-oxide (TMAO). This chapter examines the models for these two cosolvents in the context of experimental data. The direct and indirect effects of these molecules on water and protein are studied with molecular dynamics simulations. These observations and conclusions are drawn from simulations of these molecules in pure water and as a cosolvent for the protein chymotrypsin inhibitor 2. Urea-induced denaturation occurs initially through attack of the protein by water and hydration of hydrophobic protein moieties as a result of disruption of the hydrogen bonding network of water by urea. This indirect denaturing effect of urea is followed by more direct action as urea replaces some waters involved in the initial hydration of the hydrophobic core and subsequently binds to polar residues and the protein main chain to compete with the intraprotein hydrogen bonds. In the case of TMAO, we find that it encourages water-water interactions, thereby stabilizing the protein as a result of the increased penalty for the hydration of hydrophobic residues.

Molecular dynamics simulations of the protein chymotrypsin inhibitor 2 in 8 M urea at 60 degrees C were undertaken to investigate the molecular basis of chemical denaturation. The protein unfolded rapidly under these conditions, but it retained its native structure in a control simulation in water at the same temperature. The overall process of unfolding in urea was similar to that observed in thermal denaturation simulations above the protein's T(m) of 75 degrees C. The first step in unfolding was expansion of the hydrophobic core. Then, the core was solvated by water and later by urea. The denatured structures in both urea and at high temperature contained residual native helical structure, whereas the beta-structure was completely disrupted. The average residence time for urea around hydrophilic groups was six times greater than around hydrophobic residues and in all cases greater than the corresponding water residence times. Water self-diffusion was reduced 40% in 8 M urea. Urea altered water structure and dynamics, thereby diminishing the hydrophobic effect and encouraging solvation of hydrophobic groups. In addition, through urea's weakening of water structure, water became free to compete with intraprotein interactions. Urea also interacted directly with polar residues and the peptide backbone, thereby stabilizing nonnative conformations. These simulations suggest that urea denatures proteins via both direct and indirect mechanisms.

Trimethylamine n-oxide (TMAO) is a naturally occurring osmolyte that stabilizes proteins and offsets the destabilizing effects of urea. To investigate the molecular mechanism of these effects, we have studied the thermodynamics of interaction between TMAO and protein functional groups. The solubilities of a homologous series of cyclic dipeptides were measured by differential refractive index and the dissolution heats were determined calorimetrically as a function of TMAO concentration at 25 degrees C. The transfer free energy of the amide unit (-CONH-) from water to 1 M TMAO is large and positive, indicating an unfavorable interaction between the TMAO solution and the amide unit. This unfavorable interaction is enthalpic in origin. The interaction between TMAO and apolar groups is slightly favorable. The transfer free energy of apolar groups from water to TMAO consists of favorable enthalpic and unfavorable entropic contributions. This is in contrast to the contributions for the interaction between urea and apolar groups. Molecular dynamics simulations were performed to provide a structural framework for the interpretation of these results. The simulations show enhancement of water structure by TMAO in the form of a slight increase in the number of hydrogen bonds per water molecule, stronger water hydrogen bonds, and long-range spatial ordering of the solvent. These findings suggest that TMAO stabilizes proteins via enhancement of water structure, such that interactions with the amide unit are discouraged.

The human gene for catechol O-methyltransferase (COMT) contains a common polymorphism that results in substitution of methionine (M) for valine (V) at residue 108 of the soluble form of the protein. While the two proteins have similar kinetic properties, 108M COMT loses activity more rapidly than 108V COMT at 37 °C. The cosubstrate S-adenosylmethionine (SAM) stabilizes the activity of 108M COMT at 40 °C. The 108M allele has been associated with increased risk for breast cancer, obsessive-compulsive disorder, and aggressive and highly antisocial manifestations of schizophrenia. In the current work, we have constructed homology models for both human COMT polymorphs and performed molecular dynamics simulations of these models at 25, 37, and 50 °C to explore the structural consequences of the 108V/M polymorphism. The simulations indicated that replacing valine with the larger methionine residue led to greater solvent exposure of residue 108 and heightened packing interactions between M108 and helices α2, α4 (especially with R78), and α5. These altered packing interactions propagated subtle changes between the polymorphic site and the active site 16 Å away, leading to a loosening of the active site. At physiological temperature, 108M COMT sampled a larger distribution of conformations than 108V. 108M COMT was more prone to active-site distortion and had greater overall, and SAM binding site, solvent accessibility than 108V COMT at 37 °C. Similar structural perturbations were observed in the 108V protein only at 50 °C. Addition of SAM tightened up the cosubstrate pocket in both proteins and prevented the altered packing at the polymorphic site in 108M COMT. [ABSTRACT FROM AUTHOR]

We have traditionally relied on extremely elevated temperatures (498 K, 225 °C) to investigate the unfolding process of proteins within the timescale available to molecular dynamics simulations with explicit solvent. However, recent advances in computer hardware have allowed us to extend our thermal denaturation studies to much lower temperatures. Here we describe the results of simulations of chymotrypsin inhibitor 2 at seven temperatures, ranging from 298 K to 498 K. The simulation lengths vary from 94 ns to 20 ns, for a total simulation time of 344 ns, or 0.34 μs. At 298 K, the protein is very stable over the full 50 ns simulation. At 348 K, corresponding to the experimentally observed melting temperature of CI2, the protein unfolds over the first 25 ns, explores partially unfolded conformations for 20 ns, and then refolds over the last 35 ns. Above its melting temperature, complete thermal denaturation occurs in an activated process. Early unfolding is characterized by sliding or breathing motions in the protein core, leading to an unfolding transition state with a weakened core and some loss of secondary structure. After the unfolding transition, the core contacts are rapidly lost as the protein passes on to the fully denatured ensemble. While the overall character and order of events in the unfolding process are well conserved across temperatures, there are substantial differences in the timescales over which these events take place. We conclude that 498 K simulations are suitable for elucidating the details of protein unfolding at a minimum of computational expense. [Copyright &y& Elsevier]

Examines the molecular mechanism of stabilization of proteins by trimethylamine n-oxide (TMAO). Solubilities of a homologous series of cyclic dipeptides; Interaction between TMAO solution and the amide unit; Increase of the number of hydrogen bond.