Progress in the prediction of pKa values in proteins.

The pK(a) -cooperative aims to provide a forum for experimental and theoretical researchers interested in protein pK(a) values and protein electrostatics in general. The first round of the pK(a) -cooperative, which challenged computational labs to carry out blind predictions against pK(a) s experimentally determined in the laboratory of Bertrand Garcia-Moreno, was completed and results discussed at the Telluride meeting (July 6-10, 2009). This article serves as an introduction to the reports submitted by the blind prediction participants that will be published in a special issue of PROTEINS: Structure, Function and Bioinformatics. Here, we briefly outline existing approaches for pK(a) calculations, emphasizing methods that were used by the participants in calculating the blind pK(a) values in the first round of the cooperative. We then point out some of the difficulties encountered by the participating groups in making their blind predictions, and finally try to provide some insights for future developments aimed at improving the accuracy of pK(a) calculations.

Application of the Gaussian dielectric boundary in Zap to the prediction of protein pKa values.

The results of two rounds of blind pK(a) predictions for ionizable residues in staphylococcal nuclease using OpenEye's legacy protein pK(a) prediction program based on the Zap Poisson-Boltzmann solver were submitted to the 2009 prediction challenge organized by the Protein pK(a) Cooperative and first round predictions were discussed at the corresponding June 2009 Telluride conference. To better understand these results, 21 additional sets of predictions were performed with the same program, varying the internal dielectric, reference pK(a), partial charge set, and dielectric boundary. The internal dielectric (ε(p)) and dielectric boundary were the two most important factors contributing to the quality of the predictions. Although the lowest overall errors were observed with a molecular dielectric boundary at ε(p) = 8, predictions using a smooth Gaussian dielectric boundary performed almost as well at lower ε(p) values because the Gaussian boundary implicitly accounts for a significant level of solvent penetration. Improved pK(a) predictions with the Gaussian boundary methodology will require better prediction and modeling of structural changes due to changes in ionization state, perhaps without resorting to the more exhaustive sampling of conformational states used by other recent continuum methods.

Virologic failure in therapy-naive subjects on aplaviroc plus lopinavir-ritonavir: detection of aplaviroc resistance requires clonal analysis of envelope.

The CCR100136 (EPIC) study evaluated the antiviral activity of the novel CCR5 entry inhibitor aplaviroc in combination with lopinavir-ritonavir in drug-naïve human immunodeficiency virus type 1-infected subjects. Although the trial was stopped prematurely due to idiosyncratic hepatotoxicity, 11 subjects met the protocol-defined virologic failure criteria. Clonal analyses of the viral envelope tropism, aplaviroc susceptibility, and env sequencing were performed on plasma at day 1 and at the time of virologic failure. Molecular evolutionary analyses were also performed. Treatment-emergent resistance to aplaviroc or lopinavir-ritonavir was not observed at the population level. However, aplaviroc resistance was detected prior to therapy at both the clonal and population levels in one subject with virologic failure and in six subjects in a minority (<50%) of clones at day 1 or at the time of virologic failure. Reduced aplaviroc susceptibility manifested as a 50% inhibitory concentration curve shift and/or a plateau. Sequence changes in the clones with aplaviroc resistance were unique to each subject and scattered across the envelope coding region. Clones at day 1 and at the time of virologic failure were not phylogenetically distinct. Two subjects with virologic failure had a population tropism change from CCR5- to dual/mixed-tropic during treatment. Virologic failure during a regimen of aplaviroc and lopinavir-ritonavir may be associated with aplaviroc resistance, only at the clonal level, and/or, infrequently, tropism changes.

A geographic variant of the Staphylococcus aureus Panton-Valentine leukocidin toxin and the origin of community-associated methicillin-resistant S. aureus USA300.

BACKGROUND: The majority of recent community-associated methicillin-resistant Staphylococcus aureus (MRSA) infections in the United States have been caused by a single clone, USA300. USA300 secretes Panton-Valentine leukocidin (PVL) toxin, which is associated with highly virulent infections. METHODS: We sequenced the PVL genes of 174 S. aureus isolates from a global clinical sample. We combined phylogenetic reconstruction and protein modeling methods to analyze genetic variation in PVL. RESULTS: Nucleotide variation was detected at 12 of 1726 sites. Two PVL sequence variants, the R variant and the H variant, were identified on the basis of a substitution at nt 527. Of sequences obtained in the United States, 96.7% harbor the R variant, whereas 95.6% of sequences obtained outside the United States harbor the H variant; 91.3% of MRSA isolates harbor the R variant, and 91.3% of methicillin-susceptible strains harbor the H variant. A molecular model of PVL shows 3 mechanisms by which the amino acid substitution may affect PVL function. CONCLUSIONS: All sampled PVL genes appear to share a recent common ancestor and spread via a combination of clonal expansion and horizontal transfer. US isolates harbor a variant of PVL that is strongly associated with MRSA infections. Protein modeling reveals that this variant may have functional significance. We propose a hypothesis for the origin of USA300.

Structure validation by Calpha geometry: phi,psi and Cbeta deviation.

Geometrical validation around the Calpha is described, with a new Cbeta measure and updated Ramachandran plot. Deviation of the observed Cbeta atom from ideal position provides a single measure encapsulating the major structure-validation information contained in bond angle distortions. Cbeta deviation is sensitive to incompatibilities between sidechain and backbone caused by misfit conformations or inappropriate refinement restraints. A new phi,psi plot using density-dependent smoothing for 81,234 non-Gly, non-Pro, and non-prePro residues with B < 30 from 500 high-resolution proteins shows sharp boundaries at critical edges and clear delineation between large empty areas and regions that are allowed but disfavored. One such region is the gamma-turn conformation near +75 degrees,-60 degrees, counted as forbidden by common structure-validation programs; however, it occurs in well-ordered parts of good structures, it is overrepresented near functional sites, and strain is partly compensated by the gamma-turn H-bond. Favored and allowed phi,psi regions are also defined for Pro, pre-Pro, and Gly (important because Gly phi,psi angles are more permissive but less accurately determined). Details of these accurate empirical distributions are poorly predicted by previous theoretical calculations, including a region left of alpha-helix, which rates as favorable in energy yet rarely occurs. A proposed factor explaining this discrepancy is that crowding of the two-peptide NHs permits donating only a single H-bond. New calculations by Hu et al. [Proteins 2002 (this issue)] for Ala and Gly dipeptides, using mixed quantum mechanics and molecular mechanics, fit our nonrepetitive data in excellent detail. To run our geometrical evaluations on a user-uploaded file, see MOLPROBITY (http://kinemage.biochem.duke.edu) or RAMPAGE (http://www-cryst.bioc.cam.ac.uk/rampage).