Constant-pH Molecular Dynamics Simulations for Large Biomolecular Systems.

An increasingly important endeavor is to develop computational strategies that enable molecular dynamics (MD) simulations of biomolecular systems with spontaneous changes in protonation states under conditions of constant pH. The present work describes our efforts to implement the powerful constant-pH MD simulation method, based on a hybrid nonequilibrium MD/Monte Carlo (neMD/MC) technique within the highly scalable program NAMD. The constant-pH hybrid neMD/MC method has several appealing features; it samples the correct semigrand canonical ensemble rigorously, the computational cost increases linearly with the number of titratable sites, and it is applicable to explicit solvent simulations. The present implementation of the constant-pH hybrid neMD/MC in NAMD is designed to handle a wide range of biomolecular systems with no constraints on the choice of force field. Furthermore, the sampling efficiency can be adaptively improved on-the-fly by adjusting algorithmic parameters during the simulation. Illustrative examples emphasizing medium- and large-scale applications on next-generation supercomputing architectures are provided.

Validating the Water Flooding Approach by Comparing It to Grand Canonical Monte Carlo Simulations.

The study of the function of proteins on a quantitative level requires consideration of the water molecules in and around the protein. This requirement presents a major computational challenge due to the fact that the insertion of water molecules can have a very high activation barrier and would require a long simulation time. Recently, we developed a water flooding (WF) approach which is based on a postprocessing Monte Carlo ranking of possible water configurations. This approach appears to provide a very effective way for assessing the insertion free energies and determining the most likely configurations of the internal water molecules. Although the WF approach was used effectively in modeling challenging systems that have not been addressed reliably by other microscopic approaches, it was not validated by a comparison to the more rigorous grand canonical Monte Carlo (GCMC) method. Here we validate the WF approach by comparing its performance to that of the GCMC method. It is found that the WF approach reproduces the GCMC results in well-defined test cases but does so much faster. This established the WF approach as a useful strategy for finding correct water configurations in proteins and thus to provide a powerful way for studies of the functions of proteins.

Electrostatic effects in unfolded staphylococcal nuclease.

Structure-based calculations of pKa values and electrostatic free energies of proteins assume that electrostatic effects in the unfolded state are negligible. In light of experimental evidence showing that this assumption is invalid for many proteins, and with increasing awareness that the unfolded state is more structured and compact than previously thought, a detailed examination of electrostatic effects in unfolded proteins is warranted. Here we address this issue with structure-based calculations of electrostatic interactions in unfolded staphylococcal nuclease. The approach involves the generation of ensembles of structures representing the unfolded state, and calculation of Coulomb energies to Boltzmann weight the unfolded state ensembles. Four different structural models of the unfolded state were tested. Experimental proton binding data measured with a variant of nuclease that is unfolded under native conditions were used to establish the validity of the calculations. These calculations suggest that weak Coulomb interactions are an unavoidable property of unfolded proteins. At neutral pH, the interactions are too weak to organize the unfolded state; however, at extreme pH values, where the protein has a significant net charge, the combined action of a large number of weak repulsive interactions can lead to the expansion of the unfolded state. The calculated pKa values of ionizable groups in the unfolded state are similar but not identical to the values in small peptides in water. These studies suggest that the accuracy of structure-based calculations of electrostatic contributions to stability cannot be improved unless electrostatic effects in the unfolded state are calculated explicitly.

Characterization of the residual structure in the unfolded state of the Delta131Delta fragment of staphylococcal nuclease.

The determination of the conformational preferences in unfolded states of proteins constitutes an important challenge in structural biology. We use inter-residue distances estimated from site-directed spin-labeling NMR experimental measurements as ensemble-averaged restraints in all-atom molecular dynamics simulations to characterise the residual structure of the Delta131Delta fragment of staphylococcal nuclease under physiological conditions. Our findings indicate that Delta131Delta under these conditions shows a tendency to form transiently hydrophobic clusters similar to those present in the native state of wild-type staphylococcal nuclease. Only rarely, however, all these interactions are simultaneously realized to generate conformations with an overall native topology.

Detection of methicillin-resistant Staphylococcus aureus (MRSA) from growth on mannitol salt oxacillin agar using PCR for nosocomial surveillance.

This study evaluated a polymerase chain reaction (PCR) method for detection of methicillin-resistant Staphylococcus aureus (MRSA) in specimens referred for nosocomial surveillance. PCR was used to detect the mecA and nuc gene targets using yellow growth on mannitol salt agar containing 6 mg/liter oxacillin (MSO-6) as a source of DNA (N = 645). The diagnostic values for PCR compared with culture methods were 97% specificity, 100% sensitivity, 96% positive predictive value, and 100% negative predictive value. Total cost for PCR per test is $3.62 compared to $4.77 for culture. However, the total cost per specimen is significantly lower due to only 20% of all surveillance specimens producing yellow colonies on MSO-6. The average turnaround time for the PCR method is 48 h compared with 82 h for culture. PCR amplification of mecA and nuc genes using yellow colonies on MSO-6 is a simple, fast, accurate and cost-effective method for routine use in clinical laboratories for detecting MRSA in surveillance specimens.

Proline cis-trans isomerization in staphylococcal nuclease: multi-substrate free energy perturbation calculations.

Staphylococcal nuclease A exists in two folded forms that differ in the isomerization state of the Lys 116-Pro 117 peptide bond. The dominant form (90% occupancy) adopts a cis peptide bond, which is observed in the crystal structure. NMR studies show that the relatively small difference in free energy between the cis and trans forms (delta Gcis-->trans approximately 1.2 kcal/mol) results from large and nearly compensating differences in enthalpy and entropy (delta Hcis-->trans approximately delta TScis-->trans approximately 10 kcal/mol). There is evidence from X-ray crystal structures that the structural differences between the cis and the trans forms of nuclease are confined to the conformation of residues 112-117, a solvated protein loop. Here, we obtain a thermodynamic and structural description of the conformational equilibrium of this protein loop through an exhaustive conformational search that identified several substates followed by free energy simulations between the substrates. By partitioning the search into conformational substates, we overcame the multiple minima problem in this particular case and obtained precise and reproducible free energy values. The protein and water environment was implicitly modeled by appropriately chosen nonbonded terms between the explicitly treated loop and the rest of the protein. These simulations correctly predicted a small free energy difference between the cis and trans forms composed of larger, compensating differences in enthalpy and entropy. The structural predictions of these simulations were qualitatively consistent with known X-ray structures of nuclease variants and yield a model of the unknown minor trans conformation.

A further assessment of the sub-nucleosomal product, "peak A", from Physarum polycephalum.

Nuclease digestion studies of Physarum polycephalum nuclei (1-3) and nucleoli (4) over the past few years have been centred on a number of modified nucleosomal products which have been related to active-gene regions of the genome. We have re-investigated one such particle, peak A, using the techniques of differential melting and polyacrylamide gel electrophoresis and show that this material is unlikely to be a specific histone:DNA complex as suggested by earlier authors.