Inhibition of the cell cycle with chemical inhibitors: a targeted approach

Genetic links between deregulation of the cell cycle and cancer are well established. There have been significant recent developments both in our understanding of the molecular mechanisms that control cell cycle progression and in methods for protein structure determination at atomic resolution. These advances have allowed the rational design of small molecules that modulate the cell cycle by competing for sites of protein-protein or protein-ATP interactions. There is considerable optimism that these compounds, a selection of which are here reviewed, will become clinically significant drugs.     

Advanced replica-exchange sampling to study the flexibility and plasticity of peptides and proteins

Molecular dynamics (MD) simulations are ideally suited to investigate protein and peptide plasticity and flexibility simultaneously at high spatial (atomic) and high time resolution. However, the applicability is still limited by the force field accuracy and by the maximum simulation time that can be routinely achieved in current MD simulations. In order to improve the sampling the replica-exchange (REMD) methodology has become popular and is now the most widely applied advanced sampling approach. Many variants of the REMD method have been designed to reduce the computational demand or to enhance sampling along specific sets of conformational variables. An overview on recent methodological advances and discussion of specific aims and advantages of the approaches will be given. Applications in the area of free energy simulations and advanced sampling of intrinsically disordered peptides and proteins will also be discussed. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.     

Protein conformational dynamics of crambin in crystal, solution and in the trajectories of molecular dynamics simulations

Atomic displacement parameters — B factors of the eight crambin crystal structures obtained at 0.54–1.5 Å resolution and temperatures of 100–293 K have been analyzed. The comparable contributions to the B factor values are the intramolecular motions which are modeled by the harmonic vibration calculations and derived from the molecular dynamics simulation (MD) as well as rigid body changes in the position of a protein molecule as a whole. In solution for the average NMR structure of crambin the amplitudes of the backbone atomic fluctuations of the most residues of the segments with the regular backbone conformations are close to the amplitudes of the small scale harmonic vibrations. For the same residues the probability of the medium scale fluctuations fixed by the hydrogen exchange method is very low. The restricted conformational mobility of those segments is coupled with the depressed amplitudes of the fluctuation changes of the tertiary structure registered by the residue accessibility changes in an ensemble of NMR structures that forms the average NMR structure of crambin. The amplitudes of temperature fluctuations of backbone atoms and the tertiary structure raise in the segment with the irregular conformations, turn and loops. In the same segments the amplitudes of the calculated harmonic vibrations also increase, but to a lesser extent and especially in the interhelical loop with the most strong and complicated fluctuation changes of the backbone conformation. In solution for the NMR structure in this loop the conformational transitions occur between the conformational substates separated by the energy barriers, but they are not observed even in the long 100 ns trajectories from the MD simulation of crambin. These strong local fluctuation changes of the structure may play a key role in the protein functioning and modern performance improvements in the MD simulation techniques are oriented to increase the probability of protein appearance in the trajectories from the MD simulations.     

Solvent mobility and the protein 'glass' transition

Proteins and other biomolecules undergo a dynamic transition near 200 K to a glass-like solid state with small atomic fluctuations. This dynamic transition can inhibit biological function. To provide a deeper understanding of the relative importance of solvent mobility and the intrinsic protein energy surface in the transition, a novel molecular dynamics simulation procedure with the protein and solvent at different temperatures has been used. Solvent mobility is shown to be the dominant factor in determining the atomic fluctuations above 180 K, although intrinsic protein effects become important at lower temperatures. The simulations thus complement experimental studies by demonstrating the essential role of solvent in controlling functionally important protein fluctuations.     

A topological data analytic approach for discovering biophysical signatures in protein dynamics

Identifying structural differences among proteins can be a non-trivial task. When contrasting ensembles of protein structures obtained from molecular dynamics simulations, biologically-relevant features can be easily overshadowed by spurious fluctuations. Here, we present SINATRA Pro, a computational pipeline designed to robustly identify topological differences between two sets of protein structures. Algorithmically, SINATRA Pro works by first taking in the 3D atomic coordinates for each protein snapshot and summarizing them according to their underlying topology. Statistically significant topological features are then projected back onto a user-selected representative protein structure, thus facilitating the visual identification of biophysical signatures of different protein ensembles. We assess the ability of SINATRA Pro to detect minute conformational changes in five independent protein systems of varying complexities. In all test cases, SINATRA Pro identifies known structural features that have been validated by previous experimental and computational studies, as well as novel features that are also likely to be biologically-relevant according to the literature. These results highlight SINATRA Pro as a promising method for facilitating the non-trivial task of pattern recognition in trajectories resulting from molecular dynamics simulations, with substantially increased resolution.     

Molecular dynamics of ferrocytochrome c. Magnitude and anisotropy of atomic displacements

Two computer simulations of the atomic motion in tuna ferrocytochrome c have been carried out. The average structures and the structural correlations of the magnitudes of the atomic position fluctuations are in substantial agreement with recent X-ray diffraction results, particularly for the protein interior. The simulations show, however, that the atomic displacements are quite anisotropic. The degree of anisotropy and the preferred directions of atomic displacement exhibit correlations with structural features of the protein.     

Recent progress in the application of atomic force microscopy imaging and force spectroscopy to microbiology

Atomic force microscopy imaging and force spectroscopy have recently opened a range of novel applications in microbiology. During the past two years, rapid advances have been made using atomic force microscopy to visualize the surface structure of two-dimensional bacterial protein crystals, biofilms and individual cells in physiological conditions. There has also been remarkable progress in using force spectroscopy to measure biomolecular interactions and physical properties of microbial surfaces. Specific highlights include the imaging and manipulation of membrane proteins at the subnanometer level, the observation of the surface of living cells at high resolution, the mapping of local properties such as surface charges, the measurement of elastic properties of cell-surface constituents and the probing of cellular interactions using functionalized probes.     

Weight average approaches for predicting dynamical properties of biomolecules

Recent advances in atomistic molecular dynamics (MD) simulations of biomolecules allow us to explore their conformational spaces widely, observing large-scale conformational fluctuations or transitions between distinct structures. To reproduce or refine experimental data using MD simulations, structure ensembles, which are characterized by multiple structures and their statistical weights on the rugged free-energy landscapes, are often used. Here, we summarize weight average approaches for various experimental measurements. Weight average approaches are now applied to hybrid quantum mechanics/molecular mechanics MD simulations to predict fast vibrational motions in a protein with a high accuracy for better understanding of molecular functions from atomic structures.     

Determining Peptide Partitioning Properties via Computer Simulation

The transfer of polypeptide segments into lipid bilayers to form transmembrane helices represents the crucial first step in cellular membrane protein folding and assembly. This process is driven by complex and poorly understood atomic interactions of peptides with the lipid bilayer environment. The lack of suitable experimental techniques that can resolve these processes both at atomic resolution and nanosecond timescales has spurred the development of computational techniques. In this review, we summarize the significant progress achieved in the last few years in elucidating the partitioning of peptides into lipid bilayer membranes using atomic detail molecular dynamics simulations. Indeed, partitioning simulations can now provide a wealth of structural and dynamic information. Furthermore, we show that peptide-induced bilayer distortions, insertion pathways, transfer free energies, and kinetic insertion barriers are now accurate enough to complement experiments. Further advances in simulation methods and force field parameter accuracy promise to turn molecular dynamics simulations into a powerful tool for investigating a wide range of membrane active peptide phenomena.     

Structure, function and dynamics of nuclear subcompartments

The nucleus contains a plethora of different dynamic structures involved in the regulation and catalysis of nucleic acid metabolism and function. Over the past decades countless factors, molecular structures, interactions and posttranslational modifications have been described in this context. On the one side of the size scale X-ray crystallography delivers static snapshots of biomolecules at atomic resolution and on the other side light microscopy allows insights into complex structures of living cells and tissues in real time but poor resolution. Recent advances in light and electron microscopy are starting to close the temporal and spatial resolution gap from the atomic up to the cellular level. Old challenges and new insights are illustrated with examples of DNA replication and nuclear protein dynamics.     

Effect of anisotropy and anharmonicity on protein crystallographic refinement. An evaluation by molecular dynamics

Molecular dynamics simulations are employed to determine the errors introduced by anharmonicity and anisotropy in the structure and temperature factors obtained for proteins by refinement of X-ray diffraction data. Simulations (25 ps and 300 ps) of metmyoglobin are used to generate time-averaged diffraction data at 1.5 A resolution. The crystallographic restrained-parameter least-squares refinement program PROLSQ is used to refine models against these simulated data. The resulting atomic positions and isotropic temperature factors are compared with the average structure and fluctuations calculated directly from the simulations. It is found that significant errors in the atomic positions and fluctuations are introduced by the refinement, and that the errors increase with the magnitude of the atomic fluctuations. Of particular interest is the fact that the refinement generally underestimates the atomic motions. Moreover, while the actual fluctuations go up to a mean-square value of about 5 A2, the X-ray results never go above approximately 2 A2. This systematic deviation in the motional parameters appears to be due to the use of a single-site isotropic model for the atomic fluctuations. Many atoms have multiple peaks in their probability distribution functions. For some atoms, the multiple peaks are seen in difference electron density maps and it is possible to include these in the refinement as disordered residues. However, for most atoms the refinement fits only one peak and neglects the rest, leading to the observed errors in position and temperature factor. The use of strict stereochemical restraints is inconsistent with the average dynamical structure; nevertheless, refinement with tight restraints results in structures that are comparable to those obtained with loose restraints and better than those obtained with no restraints. The results support the use of tight stereochemical restraints, but indicate that restraints on the variation of temperature factors are too restrictive.     

Structure and mechanism of the γ-secretase intramembrane protease complex

γ-Secretase is a membrane protein complex that proteolyzes within the transmembrane domain of >100 substrates, including those derived from the amyloid precursor protein and the Notch family of cell surface receptors. The nine-transmembrane presenilin is the catalytic component of this aspartyl protease complex that carries out hydrolysis in the lipid bilayer. Advances in cryoelectron microscopy have led to the elucidation of the structure of the γ-secretase complex at atomic resolution. Recently, structures of the enzyme have been determined with bound APP- or Notch-derived substrates, providing insight into the nature of substrate recognition and processing. Molecular dynamics simulations of substrate-bound enzymes suggest dynamic mechanisms of intramembrane proteolysis. Structures of the enzyme bound to small-molecule inhibitors and modulators have also been solved, setting the stage for rational structure-based drug discovery targeting γ-secretase.     

Effects of substrate binding on the dynamics of RNase A: Molecular dynamics simulations of UpA bound and native RNase A

RNase A has been extensively used as a model protein in several biophysical and biochemical studies. Using the available structural and biochemical results, RNase A-UpA interaction has been computationally modeled at an atomic level. In this study, the molecular dynamics (MD) simulations of native and UpA bound RNase A have been carried out. The gross dynamical behavior and atomic fluctuations of the free and UpA bound RNase A have been characterized. Principal component analysis is carried out to identify the important modes of collective motion and to analyze the changes brought out in these modes of RNase A upon UpA binding. The hydrogen bonds are monitored to study the atomic details of RNase A-UpA interactions and RNase A-water interactions. Based on these analysis, the stability of the free and UpA bound RNase A are discussed. © 1997 John Wiley & Sons, Inc. Biopoly 42: 505–520, 1997     

Molecular dynamics of myoglobin at 298 degrees K. Results from a 300-ps computer simulation.

We have carried out a very long (300 ps) molecular dynamics simulation of the protein myoglobin. This trajectory is approximately three times longer than the longest previous molecular dynamics simulation of a protein, and ten times longer than protein simulations of comparable size (1,423 atoms in our model). Here we report results from this long simulation concerning the average structure, the mean square fluctuations of atoms about the average structure, and the nuclear magnetic resonance order parameters for various groups in myoglobin. The results demonstrate that the average coordinates change very slowly during the simulation. The relative atomic mobilities are well described by the simulation. For both the mean square atomic fluctuations and the order parameters, however, there are significant quantitative differences when values calculated using shorter portions of the trajectory are compared with results obtained for the entire 300-ps simulation. The implications of this result for obtaining converged properties from protein molecular dynamics simulations for comparison with experiment are discussed.     

Ensemble refinement of protein crystal structures: validation and application

X-ray crystallography typically uses a single set of coordinates and B factors to describe macromolecular conformations. Refinement of multiple copies of the entire structure has been previously used in specific cases as an alternative means of representing structural flexibility. Here, we systematically validate this method by using simulated diffraction data, and we find that ensemble refinement produces better representations of the distributions of atomic positions in the simulated structures than single-conformer refinements. Comparison of principal components calculated from the refined ensembles and simulations shows that concerted motions are captured locally, but that correlations dissipate over long distances. Ensemble refinement is also used on 50 experimental structures of varying resolution and leads to decreases in R(free) values, implying that improvements in the representation of flexibility observed for the simulated structures may apply to real structures. These gains are essentially independent of resolution or data-to-parameter ratio, suggesting that even structures at moderate resolution can benefit from ensemble refinement.     

REMO: A new protocol to refine full atomic protein models from C‐alpha traces by optimizing hydrogen‐bonding networks

Protein structure prediction approaches usually perform modeling simulations based on reduced representation of protein structures. For biological utilizations, it is an important step to construct full atomic models from the reduced structure decoys. Most of the current full atomic model reconstruction procedures have defects which either could not completely remove the steric clashes among backbone atoms or generate final atomic models with worse topology similarity relative to the native structures than the reduced models. In this work, we develop a new protocol, called REMO, to generate full atomic protein models by optimizing the hydrogen‐bonding network with basic fragments matched from a newly constructed backbone isomer library of solved protein structures. The algorithm is benchmarked on 230 nonhomologous proteins with reduced structure decoys generated by I‐TASSER simulations. The results show that REMO has a significant ability to remove steric clashes, and meanwhile retains good topology of the reduced model. The hydrogen‐bonding network of the final models is dramatically improved during the procedure. The REMO algorithm has been exploited in the recent CASP8 experiment which demonstrated significant improvements of the I‐TASSER models in both atomic‐level structural refinement and hydrogen‐bonding network construction. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

DNA conformation and dynamics

Nucleotide conformation and dynamics are important for the study of radiation damage to DNA at the atomic level. It is necessary to study not only normal oligonucleotide structure but also those containing modified bases which result from interaction with OH-radicals. There are now over 8000 atomic coordinate entries in the Brookhaven Protein Data Bank, of which over 900 relate to experimentally determined structures of nucleic acids and nucleic acid/protein complexes. We review some of these data which have led to the elucidation of novel DNA conformations, insight into DNA sequence specificity and knowledge of protein/DNA interactions. Further understanding of the conformation, stability and dynamics of nucleic acids has come from molecular modelling. We have used such techniques to study chemical modifications to bases such as alkylation of thymine and guanine and the effects of curvature in longer sequences. Recent improvements in this area include the inclusions of explicit counter-ions and solvent molecules, the use of Particle Mesh Ewald methods to incorporate the long-range electrostatic interactions and the use of longer time scale simulations. We have employed these methods to analyse the effects of incorporation of 8-oxodeoxyguanosine into duplex DNA. This lesion is a common result of radiation damage and is known to have important effects in mutagenesis, cancer and ageing. Received: 7 October 1998 / Accepted in revised form: 18 January 1999     

Validating Solution Ensembles from Molecular Dynamics Simulation by Wide-Angle X-ray Scattering Data

Wide-angle x-ray scattering (WAXS) experiments of biomolecules in solution have become increasingly popular because of technical advances in light sources and detectors. However, the structural interpretation of WAXS profiles is problematic, partly because accurate calculations of WAXS profiles from structural models have remained challenging. In this work, we present the calculation of WAXS profiles from explicit-solvent molecular dynamics (MD) simulations of five different proteins. Using only a single fitting parameter that accounts for experimental uncertainties because of the buffer subtraction and dark currents, we find excellent agreement to experimental profiles both at small and wide angles. Because explicit solvation eliminates free parameters associated with the solvation layer or the excluded solvent, which would require fitting to experimental data, we minimize the risk of overfitting. We further find that the influence from water models and protein force fields on calculated profiles are insignificant up to q ≈ 15 nm^?1. Using a series of simulations that allow increasing flexibility of the proteins, we show that incorporating thermal fluctuations into the calculations significantly improves agreement with experimental data, demonstrating the importance of protein dynamics in the interpretation of WAXS profiles. In addition, free MD simulations up to one microsecond suggest that the calculated profiles are highly sensitive with respect to minor conformational rearrangements of proteins, such as an increased flexibility of a loop or an increase of the radius of gyration by  <  1%. The present study suggests that quantitative comparison between MD simulations and experimental WAXS profiles emerges as an accurate tool to validate solution ensembles of biomolecules.     

Temperature dependence of dynamics of hydrated myoglobin. Comparison of force field calculations with neutron scattering data

Molecular dynamics is used to probe the atomic motions of the carboxy-myoglobin protein as a function of temperature. Simulations of 150 picoseconds in length are carried out on the protein at 20, 60, 100, 180, 220, 240, 260, 280, 300, 320 and 340 K. The simulations attempt to mimic neutron scattering experiments very closely by including a partial hydration shell around the protein. Theoretical elastic, quasielastic and inelastic neutron scattering data are derived from the trajectories and directly compared with experiment. Compared to experiment, the simulation-derived elastic scattering curves show a decrease in intensity as a function of the scattering wavevector, q2. The inelastic and quasielastic spectra show that the inelastic peak is shifted to lower frequency than the experimental value, while quasielastic behavior is in good agreement with experiment. This suggests that the theoretical model is too flexible in the harmonic limit (low temperature), but accurately reproduces high-temperature behavior. Time correlation functions of the intermediate scattering function are determined. At low temperature there is one fast decay process, and at high temperatures there is an additional slow relaxation process that is due to quasielastic scattering. The average atomic fluctuations show that the protein behaves harmonically at low temperatures. At approximately 210 K, a glass-like transition in atomic fluctuations is seen. Above the transition temperature, the atomic fluctuations exhibit both harmonic and anharmonic behavior. Comparison of protein mobility behavior with experiment indicate the fluctuations derived from simulations are larger in the harmonic region. However, the anharmonic region agrees very well with experiment. The anharmonicity is large at all temperatures, with a gradual monotonic increase from 0.5 at 20 K to greater than 0.7 at 340 K without a noticeable change at the glass transition temperature. Heavy-atom dihedral transitions are monitored as a function of temperature. Trends in the type of dihedral transitions that occur with temperature are clearly visible. Dihedral transitions involving backbone atoms occur only above the glass transition temperature. The overall protein behavior results suggest that at low temperatures there is purely vibrational motion with one fast decay process, and above the glass transition temperature there is more anharmonic motion with a fast and a slower relaxation process occurring simultaneously.(ABSTRACT TRUNCATED AT 400 WORDS)