While information from homologous structures plays a central role in X-ray structure determination by molecular replacement, such information is rarely used in NMR structure determination because it can be incorrect, both locally and globally, when evolutionary relationships are inferred incorrectly or there has been considerable evolutionary structural divergence. Here we describe a method that allows robust modeling of protein structures of up to 225 residues by combining 1 H N , 13 C, and 15 N backbone and 13 Cβ chemical shift data, distance restraints derived from homologous structures, and a physically realistic all-atom energy function. Accurate models are distinguished from inaccurate models generated using incorrect sequence alignments by requiring that (i) the all-atom energies of models generated using the restraints are lower than models generated in unrestrained calculations and (ii) the low-energy structures converge to within 2.0 Å backbone rmsd over 75% of the protein. Benchmark calculations on known structures and blind targets show that the method can accurately model protein structures, even with very remote homology information, to a backbone rmsd of 1.2-1.9 Å relative to the conventional determined NMR ensembles and of 0.9-1.6 Å relative to X-ray structures for well-defined regions of the protein structures. This approach facilitates the accurate modeling of protein structures using backbone chemical shift data without need for side-chain resonance assignments and extensive analysis of NOESY cross-peak assignments.biochemistry | biophysics | computational biology | nuclear magnetic resonance | structural genomics I n recent years, the application of multidimensional data collection techniques in isotopically enriched proteins (1) as well as development of selective labeling schemes in perdeuterated samples (2-5) and other methodological improvements (6, 7) have allowed the application of NMR methods to larger proteins (8-10). Conventional NMR structure determination relies primarily on the availability of distance restraints from NOESY experiments, which requires time-consuming experiments, including extensive analysis of side-chain resonance assignments and laborious assignments of most of the observed NOESY crosspeak resonances. While automated assignment methods (11-13) have greatly stream-lined the process (14), the assignment of NOE cross-peaks in spectra of larger proteins presents a significant challenge due to increased spectral overlap, line broadening, and low signal-to-noise ratios, rendering existing automated assignment methods ineffective in the absence of a preliminary structural model. Accurate structures can be generated for small proteins (up to 100-120 residues) using chemical shift information to guide structure prediction calculations (15, 16), but additional NMR data, including backbone NOEs and residual dipolar coupling (RDC) data, are required to obtained converged structures for larger proteins (17) and protein oligomers (18). Such data are often hard to collect and a...