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).

No abstract

All published rotamer libraries contain some rotamers that exhibit impossible internal atomic overlaps if built in ideal geometry with all hydrogen atoms. Removal of uncertain residues (mainly those with B‐factors ≥40 or van der Waals overlaps ≥0.4 Å) greatly improves the clustering of rotamer populations. Asn, Gln, or His side chains additionally benefit from flipping of their planar terminal groups when required by atomic overlaps or H‐bonding. Sensitivity to skew and to the boundaries of χ angle bins is avoided by using modes rather than traditional mean values. Rotamer definitions are listed both as the modal values and in a preferred version that maximizes common atoms between related rotamers. The resulting library shows significant differences from previous ones, differences validated by considering the likelihood of systematic misfitting of models to electron density maps and by plotting changes in rotamer frequency with B‐factor. Few rotamers now show atomic overlaps in ideal geometry; those overlaps are relatively small and can be understood in terms of bond angle distortions compensated by favorable interactions. The new library covers 94.5% of examples in the highest quality protein data with 153 rotamers and can make a significant contribution to improving the accuracy of new structures. Proteins 2000;40:389–408. © 2000 Wiley‐Liss, Inc.

All published rotamer libraries contain some rotamers that exhibit impossible internal atomic overlaps if built in ideal geometry with all hydrogen atoms. Removal of uncertain residues (mainly those with B-factors >/=40 or van der Waals overlaps >/=0.4 A) greatly improves the clustering of rotamer populations. Asn, Gln, or His side chains additionally benefit from flipping of their planar terminal groups when required by atomic overlaps or H-bonding. Sensitivity to skew and to the boundaries of chi angle bins is avoided by using modes rather than traditional mean values. Rotamer definitions are listed both as the modal values and in a preferred version that maximizes common atoms between related rotamers. The resulting library shows significant differences from previous ones, differences validated by considering the likelihood of systematic misfitting of models to electron density maps and by plotting changes in rotamer frequency with B-factor. Few rotamers now show atomic overlaps in ideal geometry; those overlaps are relatively small and can be understood in terms of bond angle distortions compensated by favorable interactions. The new library covers 94.5% of examples in the highest quality protein data with 153 rotamers and can make a significant contribution to improving the accuracy of new structures. Proteins 2000;40:389-408.

The pKa-cooperative aims to provide a forum for experimental and theoretical researchers interested in protein pKa values and protein electrostatics in general. The first round of the pKa-cooperative, which challenged computational labs to carry out blind predictions against pKas experimentally determined in the laboratory of Bertrand Garcia-Moreno, was completed and results discussed at the Telluride meeting (July 6–10, 2009). This paper 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 pKa calculations, emphasizing methods that were used by the participants in calculating the blind pKa 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 pKa calculations.

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.

Human pp60c-src is a cellular nonreceptor tyrosine kinase that participates in cytosolic signal transduction and has been implicated in the development of malignant tumors in the human breast and colon. Signal transduction is mediated by highly specific interactions between the SH2 domain and receptor phosphorylated tyrosine binding motifs. To elucidate the molecular conformation and interactions in solution, a family of highly resolved nuclear magnetic resonance (NMR) structures was determined for the src SH2 domain complexed with a high-affinity phosphorylated pentapeptide, acetyl-p YEEIE-OH. The 23 structures, generated with a distance geometry (DG) and a dynamical simulated annealing (SA) procedure, satisfied 2072 experimental restraints derived from a variety of multifrequency/multidimensional and isotope-filtered NMR data. Superimposition of residues 143-245 upon the mean coordinate set yielded an atomic rmsd of 0.58 +/- 0.09 A for the N, C alpha, C' atoms and 1.04 +/- 0.08 for all the non-hydrogen atoms. Residues in the ordered secondary structure regions superimpose to 0.29 +/- 0.04 A for the N, C alpha, C' and 0.73 +/- 0.08 A for all the non-hydrogen atoms. The angular order parameter calculated for the phi, psi angles was > 0.9 for 81 of the 106 protein residues. The main protein conformational features are three antiparallel beta-strands that traverse a compact core with an alpha-helix on each side of the core near the N- and C-termini. The observed intermolecular nuclear Overhauser effects (NOE) from the pY, +1E, and +3I residues positioned the ligand in an extended conformation across the SH2 domain surface with the pY and +3I side chains inserted into the protein binding pockets. In general, the protein conformation is consistent with previously reported structures of different SH2 domain complexes determined by X-ray crystallography. However, inter- or intramolecular interactions involving the guanidinium side chains of the solvated R alpha A2 or the buried R beta B5 were not observed at pH = 5.5 or 7.0. If such interactions exist in solution, the absence of any confirming data probably arises from rapid exchange with solvent and/or undetermined dynamic components. Thus, the unrestrained R alpha A2 side chain did not show an amino-aromatic interaction or a hydrogen bond to the -1 carbonyl oxygen as observed in the crystal structures. This result is consistent with the solution structure of a different SH2 domain complex. A more detailed comparison between the crystal structure and the NMR-derived solution structures of the same src SH2 domain complex is presented.(ABSTRACT TRUNCATED AT 400 WORDS)

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 protocoldefined 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.Although highly active antiretroviral therapy results in a profound and sustained reduction in plasma human immunodeficiency virus type 1 (HIV-1) RNA in many individuals, the side effects of currently available antiretroviral therapy (ART) and the emergence of multidrug-resistant viral strains continue to represent major challenges for the management of HIV infection (7,29). The discovery that the chemokine receptors CCR5 and CXCR4 function as coreceptors that mediate HIV-1 entry into CD4 ϩ host cells (2) has led to the development of coreceptor inhibitors, which are currently being tested for use in ART.Unlike other drug classes, resistance to CCR5 entry inhibitors (CCR5 EIs) could involve the use of drug-bound CCR5 or a change in coreceptor usage. Much of the data on CCR5 EI resistance come from in vitro passage and short-term monotherapy studies. An analysis of virus envelopes (Env) during passage studies with CCR5 EIs demonstrated the maintenance of R5 tropism in a setting where the cells utilized expressed both CCR5 and CXCR4 (17,34,38). Resistance or reduced susceptibility to aplaviroc (APL) was slow to emerge in passage studies; the characterization of CCR5-tropic (R5-tropic) HIV-1 isolates from long-term passage (Ͼ48 weeks) showed marginal increases in the change (n-fold) in the 50% inhibitory concentration (FCIC 50 ) values, the standard measure by which resistance to other antiretroviral classes is detected (16,17).Another measure of drug susceptibility that has been suggested for CCR...