Three‐dimensional cell culture model utilization in cancer stem cell research

Three‐dimensional (3D) cell culture models are becoming increasingly popular in contemporary cancer research and drug resistance studies. Recently, scientists have begun incorporating cancer stem cells (CSCs) into 3D models and modifying culture components in order to mimic in vivo conditions better. Currently, the global cell culture market is primarily focused on either 3D cancer cell cultures or stem cell cultures, with less focus on CSCs. This is evident in the low product availability officially indicated for 3D CSC model research. This review discusses the currently available commercial products for CSC 3D culture model research. Additionally, we discuss different culture media and components that result in higher levels of stem cell subpopulations while better recreating the tumor microenvironment. In summary, although progress has been made applying 3D technology to CSC research, this technology could be further utilized and a greater number of 3D kits dedicated specifically to CSCs should be implemented.     

Brief communication: comparing loading scenarios in lower first molar supporting bone structure using 3D finite element analysis

Finite element analysis (FEA) is a widespread technique to evaluate the stress/strain distributions in teeth or dental supporting tissues. However, in most studies occlusal forces are usually simplified using a single vector (i.e., point load) either parallel to the long tooth axis or oblique to this axis. In this pilot study we show how lower first molar occlusal information can be used to investigate the stress distribution with 3D FEA in the supporting bone structure. The LM₁ and the LP²‐LM¹ of a dried modern human skull were scanned by μCT in maximum intercuspation contact. A kinematic analysis of the surface contacts between LM₁ and LP²‐LM¹ during the power stroke was carried out in the occlusal fingerprint analyzer (OFA) software to visualize contact areas during maximum intercuspation contact. This information was used for setting the occlusal molar loading to evaluate the stress distribution in the supporting bone structure using FEA. The output was compared to that obtained when a point force parallel to the long axis of the tooth was loaded in the occlusal basin. For the point load case, our results indicate that the buccal and lingual cortical plates do not experience notable stresses. However, when the occlusal contact areas are considered, the disto‐lingual superior third of the mandible experiences high tensile stresses, while the medio‐lingual cortical bone is subjected to high compressive stresses. Developing a more realistic loading scenario leads to better models to understand the relationship between masticatory function and mandibular shape and structures. Am J Phys Anthropol, 2012. © 2011 Wiley Periodicals, Inc.     

Comparative structural connectivity spectra analysis (CoSCoSA) models of steroids binding to the aromatase enzyme

A method that combines NMR spectral and structural information into a constructed three-dimensional (3D)-connectivity matrix is developed for modeling biological binding activity of small molecules. The 3D-connectivity matrix for a molecule is defined by associating the distances between all possible carbon-to-carbon connections with their assigned carbon NMR chemical shifts. In this project we selected from the total 3D-connectivity matrix a subset, the two-dimensional (2D) (13)C-(13)C COSY and a theoretical long range 2D (13)C-(13)C distance connectivity spectral plane. Patterns of (13)C chemical shifts observed at these two relative distances for 50 steroids were used to produce a mathematical relationship for the steroids' relative binding affinity (pK(i)) to the aromatase enzyme. We call this technique comparative structural connectivity spectra analysis (CoSCoSA) modeling. Using combinations of the 2D COSY and 2D long-range distance spectra as modeling parameters, we built four CoSCoSA models. One model was made from the 2D COSY spectra alone and another was developed using only the 2D long-range distance spectra. Then the COSY and long-distance spectra were combined in two different ways: starting with the combined principal components (PCs) from the separately calculated COSY and distance spectra or using the combined raw spectra (3D). The best CoSCoSA model was based on the combined PCs from COSY and distance spectra. This model had an r(2) of 0.96 and a leave-one-out cross-validation (q(2)) of 0.92. In general CoSCoSA modeling combines the quantum mechanical information inherent in NMR chemical shifts with internal molecular atom-to-atom distances to give a reliable and straightforward basis for predictive modeling. The technique has the flexibility and accuracy to outperform not only the cross-validated variance q(2) of previously published quantitative structure-activity relationships (QSAR) but also those obtained by related quantitative spectral data-activity relationships (QSDARs) lacking connectivity dimensions.     

Hemodynamic Based Surgical Decision on Sequential Graft and Y-Type Graft in Coronary Artery Bypass Grafting

Purpose: Sequential graft and Y-type graft are two different surgical procedures in coronary artery bypass grafting (CABG). The hemodynamic environment of them are different, that may cause different short-term surgical result and long-term patency. In this study, the short-term and long-term result of sequential and Y-type graft was discussed by comparing the hemodynamics of them. Materials and Methods: Two postoperative 3-dimensional (3D) models were built by applying different graft on a patient-specific 3D model with serious stenosis. Then zero-dimensional (0D)/3D coupled simulation was carried out by coupling the postoperative 3D models with a 0D lumped parameter model of the cardiovascular system. Results: The flow rate of native coronary arteries and grafts are all calculated and illustrated in this paper. No significant difference of the native coronary arteries flow and graft flow exists between two surgical procedures. The wall shear stress (WSS) and streamline were also depicted. The graft WSS of sequential graft is 19.1% higher than Y-type graft. While flow separation appears at the bifurcation of Y-type graft. Conclusion: The short-term outcomes of sequential graft and Y-type graft are almost the same. But it can be found from the hemodynamics factors that the longterm patency of the sequential graft is better.     

Critical analysis of the successes and failures of homology models of G protein‐coupled receptors

We present a critical assessment of the performance of our homology model refinement method for G protein‐coupled receptors (GPCRs), called LITICon that led to top ranking structures in a recent structure prediction assessment GPCRDOCK2010. GPCRs form the largest class of drug targets for which only a few crystal structures are currently available. Therefore, accurate homology models are essential for drug design in these receptors. We submitted five models each for human chemokine CXCR4 (bound to small molecule IT1t and peptide CVX15) and dopamine D3DR (bound to small molecule eticlopride) before the crystal structures were published. Our models in both CXCR4/IT1t and D3/eticlopride assessments were ranked first and second, respectively, by ligand RMSD to the crystal structures. For both receptors, we developed two types of protein models: homology models based on known GPCR crystal structures, and ab initio models based on the prediction method MembStruk. The homology‐based models compared better to the crystal structures than the ab initio models. However, a robust refinement procedure for obtaining high accuracy structures is needed. We demonstrate that optimization of the helical tilt, rotation, and translation is vital for GPCR homology model refinement. As a proof of concept, our in‐house refinement program LITiCon captured the distinct orientation of TM2 in CXCR4, which differs from that of adrenoreceptors. These findings would be critical for refining GPCR homology models in future. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Thermo-responsive non-woven scaffolds for "smart" 3D cell culture

The thermo‐responsive polymer poly(N‐isopropylacrylamide) has received widespread attention for its in vitro application in the non‐invasive, non‐destructive release of adherent cells on two dimensional surfaces. In this study, 3D non‐woven scaffolds fabricated from poly(propylene) (PP), poly(ethylene terephthalate) (PET), and nylon that had been grafted with PNIPAAm were tested for their ability to support the proliferation and subsequent thermal release of HC04 and HepG2 hepatocytes. Hepatocyte viability and proliferation were estimated using the Alamar Blue assay and Hoechst 33258 total DNA quantification. The assays revealed that the pure and grafted non‐woven scaffolds maintained the hepatocytes within the matrix and promoted 3D proliferation comparable to that of the commercially available Algimatrix? alginate scaffold. Albumin production and selected cytochrome P450 genes expression was found to be superior in cells growing on pure and grafted non‐woven PP scaffolds as compared to cells grown as a 2D monolayer. Two scaffolds, namely, PP‐g‐PNIPAAm‐A and PP‐g‐PNIPAAm‐B were identified as having far superior thermal release capabilities; releasing the majority of the cells from the matrices within 2 h. This is the first report for the development of 3D non‐woven, thermo‐responsive scaffolds able to release cells from the matrix without the use of any enzymatic assistance or scaffold degradation. Biotechnol. Bioeng. 2012; 109:2147–2158. © 2012 Wiley Periodicals, Inc.     

3D structural model of the G-protein-coupled cannabinoid CB2 receptor

The potential for therapeutic specificity in regulating diseases and for reduced side effects has made cannabinoid (CB) receptors one of the most important G-protein-coupled receptor (GPCR) targets for drug discovery. The cannabinoid (CB) receptor subtype CB2 is of particular interest due to its involvement in signal transduction in the immune system and its increased characterization by mutational and other studies. However, our understanding of their mode of action has been limited by the absence of an experimental receptor structure. In this study, we have developed a 3D model of the CB2 receptor based on the recent crystal structure of a related GPCR, bovine rhodopsin. The model was developed using multiple sequence alignment of homologous receptor sub-types in humans and mammals, and compared with other GPCRs. Alignments were analyzed with mutation scores, pairwise hydrophobicity profiles and Kyte-Doolittle plots. The 3D model of the transmembrane segment was generated by mapping the CB2 sequence onto the homologous residues of the rhodopsin structure. The extra- and intracellular loop regions of the CB2 were generated by searching for homologous C(alpha) backbone sequences in published structures in the Brookhaven Protein Databank (PDB). Residue side chains were positioned through a combination of rotamer library searches, simulated annealing and minimization. Intermediate models of the 7TM helix bundles were analyzed in terms of helix tilt angles, hydrogen-bond networks, conserved residues and motifs, possible disulfide bonds. The amphipathic cytoplasmic helix domain was also correlated with biological and site-directed mutagenesis data. Finally, the model receptor-binding cavity was characterized using solvent-accessible surface approach.     

Scaffold‐free inkjet printing of three‐dimensional zigzag cellular tubes

The capability to print three‐dimensional (3D) cellular tubes is not only a logical first step towards successful organ printing but also a critical indicator of the feasibility of the envisioned organ printing technology. A platform‐assisted 3D inkjet bioprinting system has been proposed to fabricate 3D complex constructs such as zigzag tubes. Fibroblast (3T3 cell)‐based tubes with an overhang structure have been successfully fabricated using the proposed bioprinting system. The post‐printing 3T3 cell viability of printed cellular tubes has been found above 82% (or 93% with the control effect considered) even after a 72‐h incubation period using the identified printing conditions for good droplet formation, indicating the promising application of the proposed bioprinting system. Particularly, it is proved that the tubular overhang structure can be scaffold‐free fabricated using inkjetting, and the maximum achievable height depends on the inclination angle of the overhang structure. As a proof‐of‐concept study, the resulting fabrication knowledge helps print tissue‐engineered blood vessels with complex geometry. Biotechnol. Bioeng. 2012; 109: 3152–3160. © 2012 Wiley Periodicals, Inc.     

The combination of two‐dimensional and three‐dimensional analysis methods contributes to the understanding of glioblastoma spatial heterogeneity

Heterogeneity is regarded as the major factor leading to the poor outcomes of glioblastoma (GBM) patients. However, conventional two‐dimensional (2D) analysis methods, such as immunohistochemistry and immunofluorescence, have limited capacity to reveal GBM spatial heterogeneity. Thus, we sought to develop an effective analysis strategy to increase the understanding of GBM spatial heterogeneity. Here, 2D and three‐dimensional (3D) analysis methods were compared for the examination of cell morphology, cell distribution and large intact structures, and both types of methods were employed to dissect GBM spatial heterogeneity. The results showed that 2D assays showed only cross‐sections of specimens but provided a full view. To visualize intact GBM specimens in 3D without sectioning, the optical tissue clearing methods CUBIC and iDISCO+ were used to clear opaque specimens so that they would become more transparent, after which the specimens were imaged with a two‐photon microscope. The 3D analysis methods showed specimens at a large spatial scale at cell‐level resolution and had overwhelming advantages in comparison to the 2D methods. Furthermore, in 3D, heterogeneity in terms of cell stemness, the microvasculature, and immune cell infiltration within GBM was comprehensively observed and analysed. Overall, we propose that 2D and 3D analysis methods should be combined to provide much greater detail to increase the understanding of GBM spatial heterogeneity.     

Design of novel quinazolinone derivatives as inhibitors for 5HT7 receptor

To study the pharmacophore properties of quinazolinone derivatives as 5HT₇ inhibitors, 3D QSAR methodologies, namely Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) were applied, partial least square (PLS) analysis was performed and QSAR models were generated. The derived model showed good statistical reliability in terms of predicting the 5HT₇ inhibitory activity of the quinazolione derivative, based on molecular property fields like steric, electrostatic, hydrophobic, hydrogen bond donor and hydrogen bond acceptor fields. This is evident from statistical parameters like q² (cross validated correlation coefficient) of 0.642, 0.602 and r² (conventional correlation coefficient) of 0.937, 0.908 for CoMFA and CoMSIA respectively. The predictive ability of the models to determine 5HT₇ antagonistic activity is validated using a test set of 26 molecules that were not included in the training set and the predictive r² obtained for the test set was 0.512 & 0.541. Further, the results of the derived model are illustrated by means of contour maps, which give an insight into the interaction of the drug with the receptor. The molecular fields so obtained served as the basis for the design of twenty new ligands. In addition, ADME (Adsorption, Distribution, Metabolism and Elimination) have been calculated in order to predict the relevant pharmaceutical properties, and the results are in conformity with required drug like properties.     

Dynactin 3D Structure: Implications for Assembly and Dynein Binding

The multisubunit protein complex, dynactin, is an essential component of the cytoplasmic dynein motor. High-resolution structural work on dynactin and the dynein/dynactin supercomplex has been limited to small subunits and recombinant fragments that do not report fully on either ≈ 1 MDa assembly. In the present study, we used negative-stain electron microscopy and image analysis based on random conical tilt reconstruction to obtain a three-dimensional (3D) structure of native vertebrate dynactin. The 35-nm-long dynactin molecule has a V-shaped shoulder at one end and a flattened tip at the other end, both offset relative to the long axis of the actin-related protein (Arp) backbone. The shoulder projects dramatically away from the Arp filament core in a way that cannot be appreciated in two-dimensional images, which has implications for the mechanism of dynein binding. The 3D structure allows the helical parameters of the entire Arp filament core, which includes the actin capping protein, CP, to be determined for the first time. This structure exhibits near identity to F-actin and can be well fitted into the dynactin envelope. Molecular fitting of modeled CP-Arp polymers into the envelope shows that the filament contains between 7 and 9 Arp protomers and is capped at both ends. In the 7 Arp model, which agrees best with measured Arp stoichiometry and other structural information, actin capping protein (CP) is not present at the distal tip of the structure, unlike what is seen in the other models. The 3D structure suggests a mechanism for dynactin assembly and length specification.     

A three-dimensional in vitro ovarian cancer coculture model using a high-throughput cell patterning platform

In vitro 3D cancer models that provide a more accurate representation of disease in vivo are urgently needed to improve our understanding of cancer pathology and to develop better cancer therapies. However, development of 3D models that are based on manual ejection of cells from micropipettes suffer from inherent limitations such as poor control over cell density, limited repeatability, low throughput, and, in the case of coculture models, lack of reproducible control over spatial distance between cell types (e.g., cancer and stromal cells). In this study, we build on a recently introduced 3D model in which human ovarian cancer (OVCAR-5) cells overlaid on Matrigel^™ spontaneously form multicellular acini. We introduce a high-throughput automated cell printing system to bioprint a 3D coculture model using cancer cells and normal fi broblasts micropatterned on Matrigel^™. Two cell types were patterned within a spatially controlled microenvironment (e.g., cell density, cell-cell distance) in a high-throughput and reproducible manner; both cell types remained viable during printing and continued to proliferate following patterning. This approach enables the miniaturization of an established macro-scale 3D culture model and would allow systematic investigation into the multiple unknown regulatory feedback mechanisms between tumor and stromal cells and provide a tool for high-throughput drug screening.     

Three-dimensional holograph vector of atomic interaction field (3D-HoVAIF): a novel rotation-translation invariant 3D structure descriptor and its applications to peptides.

Quantitative structure-activity relationship (QSAR) study, important in drug design, mainly involves two aspects, molecular structural characterization (MSC) and construction of a statistical model. MSC focuses on transforming molecular structural and property characteristics into a group of numerical codes, dedicated to minimizing information loss during this process. In this context, common atoms in organic compounds are classified according to their families in the periodic table, and hybridization states, and on the basis of these, three nonbonding interactions (i.e. electrostatic, van der Waals and hydrophobic) are calculated, ultimately resulting in a new rotation-translation invariant, 3D-MSC, as a three-dimensional holograph vector of atomic interaction field (3D-HoVAIF). By applying 3D-HoVAIF to QSAR studies on two classical peptides including 58 angiotensin-converting enzyme (ACE) inhibitors and 48 bitter-tasting dipeptides, we get two excellent genetic algorithm-partial least squares (GA-PLS) models, with statistics r(2), q(2), root mean square error (RMSEE), and root mean square error of cross-validation (RMSCV) of 0.857, 0.811, 0.376, and 0.432 for ACE inhibitors and 0.940, 0.892, 0.153 and 0.205 for bitter-tasting dipeptides, respectively. By equally dividing the two datasets into training and test sets by D-optimal, the 3D-HoVAIF approach undergoes rigorous statistical validation. Furthermore, the superior performance of 3D-HoVAIF is confirmed in comparison with two other peptide MSC approaches referring to z-scale and ISA-ECI. For 58 ACE inhibitors, the GA-PLS model yields two principal components, with the following statistics: r(2) = 0.893, q(2) = 0.824, RMSEE = 0.349, RMSCV = 0.425, q2(ext) = 0.739, r2(ext)= 0.784, r2(0.ext) = 0.781, rf2(0.ext) = 0.77, k = 0.962, k' = 1.019, and RMSEP = 0.460; for 48 bitter-tasting dipeptides, three principal components resulted, with the statistics as: r(2) = 0.950, q(2) = 0.893, RMSEE = 0.152, RMSCV = 0.222, q2(ext)= 0.875, r2(ext) = 0.919, r2(0.ext)= 0.919, rf2(0.ext)= 0.919, k = 1.018, k' = 0.974, and RMSEP = 0.198. In addition, the relationship of ACE-inhibiting activities with bitter-tasting thresholds has been investigated by applying the above-constructed models to predictions on 400 theoretically possible dipeptides. Through analysis, the ACE-inhibiting activities are found to be prominently related to bitter-tasting intensities. Thus, it is deemed to be difficult to find such dipeptides that simultaneously satisfy pharmacodynamic action (high ACE-inhibiting activities) and comfortable tastes, suggesting that active components of dipeptides that are served as functional food to lower blood pressure are not very ideal.     

Two-dimensional intraventricular flow pattern visualization using the image-based computational fluid dynamics

The image-based computational fluid dynamics (IB-CFD) technique, as the combination of medical images and the CFD method, is utilized in this research to analyze the left ventricle (LV) hemodynamics. The research primarily aims to propose a semi-automated technique utilizing some freely available and commercial software packages in order to simulate the LV hemodynamics using the IB-CFD technique. In this research, moreover, two different physiological time-resolved 2D models of a patient-specific LV with two different types of aortic and mitral valves, including the orifice-type valves and integrated with rigid leaflets, are adopted to visualize the process of developing intraventricular vortex formation and propagation. The blood flow pattern over the whole cardiac cycle of two models is also compared to investigate the effect of utilizing different valve types in the process of the intraventricular vortex formation. Numerical findings indicate that the model with integrated valves can predict more complex intraventricular flow that can match better the physiological flow pattern in comparison to the orifice-type model.     

High-throughput identification of factors promoting neuronal differentiation of human neural progenitor cells in microscale 3D cell culture

Identification of conditions for guided and specific differentiation of human stem cell and progenitor cells is important for continued development and engineering of in vitro cell culture systems for use in regenerative medicine, drug discovery, and human toxicology. Three-dimensional (3D) and organotypic cell culture models have been used increasingly for in vitro cell culture because they may better model endogenous tissue environments. However, detailed studies of stem cell differentiation within 3D cultures remain limited, particularly with respect to high-throughput screening. Herein, we demonstrate the use of a microarray chip-based platform to screen, in high-throughput, individual and paired effects of 12 soluble factors on the neuronal differentiation of a human neural progenitor cell line (ReNcell VM) encapsulated in microscale 3D Matrigel cultures. Dose–response analysis of selected combinations from the initial combinatorial screen revealed that the combined treatment of all-trans retinoic acid (RA) with the glycogen synthase kinase 3 inhibitor CHIR-99021 (CHIR) enhances neurogenesis while simultaneously decreases astrocyte differentiation, whereas the combined treatment of brain-derived neurotrophic factor and the small azide neuropathiazol enhances the differentiation into neurons and astrocytes. Subtype specification analysis of RA- and CHIR-differentiated cultures revealed that enhanced neurogenesis was not biased toward a specific neuronal subtype. Together, these results demonstrate a high-throughput screening platform for rapid evaluation of differentiation conditions in a 3D environment, which will aid the development and application of 3D stem cell culture models.     

3D‐Printed Cathodes of LiMn1−xFexPO4 Nanocrystals Achieve Both Ultrahigh Rate and High Capacity for Advanced Lithium‐Ion Battery

A 3D‐printing technology and printed 3D lithium‐ion batteries (3D‐printed LIBs) based on LiMn_0.21Fe_0.79PO₄@C (LMFP) nanocrystal cathodes are developed to achieve both ultrahigh rate and high capacity. Coin cells with 3D‐printed cathodes show impressive electrochemical performance: a capacity of 108.45 mAh g^?1 at 100 C and a reversible capacity of 150.21 mAh g^?1 at 10 C after 1000 cycles. In combination with simulation using a pseudo 2D hidden Markov model and experimental data of 3D‐printed and traditional electrodes, for the first time deep insight into how to achieve the ultrahigh rate performance for a cathode with LMFP nanocrystals is obtained. It is estimated that the Li‐ion diffusion in LMFP nanocrystal is not the rate‐limitation step for the rate to 100 C, however, that the electrolyte diffusion factors, such as solution intrinsic diffusion coefficient, efficiency porosity, and electrode thickness, will dominate ultrahigh rate performance of the cathode. Furthermore, the calculations indicate that the above factors play important roles in the equivalent diffusion coefficient with the electrode beyond a certain thickness, which determines the whole kinetic process in LIBs. This fundamental study should provide helpful guidance for future design of LIBs with superior electrochemical performance.     

384 hanging drop arrays give excellent Z-factors and allow versatile formation of co-culture spheroids

We previously reported the development of a simple, user-friendly, and versatile 384 hanging drop array plate for 3D spheroid culture and the importance of utilizing 3D cellular models in anti-cancer drug sensitivity testing. The 384 hanging drop array plate allows for high-throughput capabilities and offers significant improvements over existing 3D spheroid culture methods. To allow for practical 3D cell-based high-throughput screening and enable broader use of the plate, we characterize the robustness of the 384 hanging drop array plate in terms of assay performance and demonstrate the versatility of the plate. We find that the 384 hanging drop array plate performance is robust in fluorescence- and colorimetric-based assays through Z-factor calculations. Finally, we demonstrate different plate capabilities and applications, including: spheroid transfer and retrieval for Janus spheroid formation, sequential addition of cells for concentric layer patterning of different cell types, and culture of a wide variety of cell types.     

Cell behavior observation and gene expression analysis of melanoma associated with stromal fibroblasts in a three‐dimensional magnetic cell culture array

A three‐dimensional (3D) multicellular tumor spheroid culture array has been fabricated using a magnetic force‐based cell patterning method, analyzing the effect of stromal fibroblast on the invasive capacity of melanoma. Formation of spheroids was observed when array‐like multicellular patterns of melanoma were developed using a pin‐holder device made of magnetic soft iron and an external magnet, which enables the assembly of the magnetically labeled cells on the collagen gel‐coated surface as array‐like cell patterns. The interaction of fibroblast on the invasion of melanoma was investigated using three types of cell interaction models: (i) fibroblasts were magnetically labeled and patterned together in array with melanoma spheroids (direct‐interaction model), (ii) fibroblasts coexisting in the upper collagen gel (indirect‐interaction model) of melanoma spheroids, and (iii) fibroblast‐sheets coexisting under melanoma spheroids (fibroblast‐sheet model). The fibroblast‐sheet model has largely increased the invasive capacity of melanoma, and the promotion of adhesion, migration, and invasion were also observed. In the fibroblast‐sheet model, the expression of IL‐8 and MMP‐2 increased by 24‐fold and 2‐fold, respectively, in real time RT‐PCR compared to the absence of fibroblasts. The results presented in this study demonstrate the importance of fibroblast interaction to invasive capacity of melanoma in the 3D in vitro bioengineered tumor microenvironment. © 2012 American Institute of Chemical Engineers Biotechnol. Prog., 2013     

Technical Note: 3D From Standard Digital Photography of Human Crania—A Preliminary Assessment

This study assessed three‐dimensional (3D) photogrammetry as a tool for capturing and quantifying human skull morphology. While virtual reconstruction with 3D surface scanning technology has become an accepted part of the paleoanthropologist's tool kit, recent advances in 3D photogrammetry make it a potential alternative to dedicated surface scanners. The principal advantages of photogrammetry are more rapid raw data collection, simplicity and portability of setup, and reduced equipment costs. We tested the precision and repeatability of 3D photogrammetry by comparing digital models of human crania reconstructed from conventional, 2D digital photographs to those generated using a 3D surface scanner. Overall, the photogrammetry and scanner meshes showed low degrees of deviation from one another. Surface area estimates derived from photogrammetry models tended to be slightly larger. Landmark configurations generally did not cluster together based upon whether the reconstruction was created with photogrammetry or surface scanning technology. Average deviations of landmark coordinates recorded on photogrammetry models were within the generally allowable range of error in osteometry. Thus, while dependent upon the needs of the particular research project, 3D photogrammetry appears to be a suitable, lower‐cost alternative to 3D imaging and scanning options. Am J Phys Anthropol 154:152–158, 2014. © 2014 Wiley Periodicals, Inc.