Automating crystallographic structure solution and refinement of protein–ligand complexes

The microfluidic chips used in this study were implemented as a two-layer polydimethylsiloxane (PDMS) microfluidic chip similar to published procedures, summarized briefly as follows. The channels on the upper layer are flow microchannels used for the reagent inlets, filling chambers, and output of generated droplets, whereas the channels on the lower layer are control channels that can be pressurized or depressurized to open or close the corresponding microvalve(s). Silicon wafer molds composed of photoresist-patterned microstructures for fabrication of the fluid and control layers were created by standard photolithographic techniques. Both molds were pretreated with trimethylsilyl chloride vapor for 10 minutes to facilitate the release of final PDMS replica. Well-mixed PDMS prepolymer (GE RTV615, total weight 36 g, mixing ratio A:B 5 5:1) was poured onto the fluid layer mold to give a 6 mm–thick fluidic layer with approximately 40 mm channel depth. Another portion of PDMS prepolymer (GE RTV615, total weight 10 g, mixing ratio A:B 5 20:1) was mixed and then spin-coated onto the control layer mold at 1,500 rpm for 60 seconds. The fluidic and control layers were cured at 80uC for 15 and 18 minutes, respectively. After baking, the fluidic layer replica was peeled from the mold and aligned onto the control layer, and then the assembly was baked at 80uC for at least 6 hours to adhere the layers. The chip was peeled off the mold when the fluidic and control layers were firmly bonded together. Holes were then punched to form ports connected to the fluidic layer channels for reagent inlets and outlets and ports connected to the control layer channels for microvalve and pump actuation with hydraulic fluid (water). Adhesion of the chip to a clean glass microscope slide to seal the control channels was achieved by corona discharge treatment. The microfluidic device was baked in an oven at 80uC for 72 hours to restore the intrinsic hydrophobicity of PDMS surfaces, which is needed to minimize reagent loss and nonspecific binding of biomolecules on channel walls when moving aqueous droplets. Microvalve control lines were filled with water as hydraulic fluid. Microvalves were actuated by pressurizing the corresponding control channels (up to 60 psig) via electronic solenoid valves (Series S070, SMC, Toyko, Japan). All valves were automatically controlled through a data acquisition module (USB-4750, Advantech, Milpitas, CA) driven by a custom software program written in LabView (National Instruments, Austin, TX). The microfluidic system includes two microfluidic chips (a distribution chip of pH buffers and a DMDG chip connected by tubing). The DMDG chip is composed of three functional parts (see Figure 1B): (1) a droplet generation core, where specific quantities of reagents are measured and merged into composition-specific droplets; (2) a peristaltic pump, which produces serial compressed nitrogen pulses that can precisely deliver intact droplets to the desired location; and (3) a mixing channel. Figure S1 shows the external solution distribution chip (see Figure S1A) and the droplet generation core (see Figure S1B). The unique features of our DMDG chip (see Figure 1B) are (1) independent control of volume and composition for every droplet, enabling screening with minimal reagent consumption; (2) the ability to pause, modify, and restart the droplet generation process, for example, for replacement or change of reagents; and (3) the use of nitrogen gas rather than oil to separate droplets, eliminating the need for oil removal steps afterward. Using the microvalves to isolate each reagent inlet, the incoming reagents are not in contact with each other until the moment of droplet formation. They will be rapidly mixed and reacted thereafter while moving along the microfluidic channel. This is particularly critical because some radiolabeling tags, such as [F]SFB, are often unstable and prone to hydrolyze or decompose at a higher pH.

[1]  A. Laub,et al.  The singular value decomposition: Its computation and some applications , 1980 .

[2]  Comparison of propranolol and hydrochlorothiazide for the initial treatment of hypertension. II. Results of long-term therapy. Veterans Administration Cooperative Study Group on Antihypertensive Agents. , 1982, JAMA.

[3]  H. White,et al.  Some heteroskedasticity-consistent covariance matrix estimators with improved finite sample properties☆ , 1985 .

[4]  S. Zeger,et al.  Longitudinal data analysis using generalized linear models , 1986 .

[5]  M. Piedmonte,et al.  A Method for Generating High-Dimensional Multivariate Binary Variates , 1991 .

[6]  E. Schulze,et al.  Relationships among Maximum Stomatal Conductance, Ecosystem Surface Conductance, Carbon Assimilation Rate, and Plant Nitrogen Nutrition: A Global Ecology Scaling Exercise , 1994 .

[7]  R. Neilson A Model for Predicting Continental‐Scale Vegetation Distribution and Water Balance , 1995 .

[8]  T. Hedner,et al.  Randomised trial of old and new antihypertensive drugs in elderly patients: cardiovascular mortality and morbidity the Swedish Trial in Old Patients with Hypertension-2 study , 1999, The Lancet.

[9]  N. Wald,et al.  When can a risk factor be used as a worthwhile screening test? , 1999, BMJ.

[10]  L. Niskanen,et al.  Effect of angiotensin-converting-enzyme inhibition compared with conventional therapy on cardiovascular morbidity and mortality in hypertension: the Captopril Prevention Project (CAPPP) randomised trial , 1999, The Lancet.

[11]  S. Quake,et al.  Monolithic microfabricated valves and pumps by multilayer soft lithography. , 2000, Science.

[12]  A. Pitman,et al.  The uncertainty in simulations by a Global Biome Model (BIOME3) to alternative parameter values , 2000 .

[13]  Thomas Hedner,et al.  Randomised trial of effects of calcium antagonists compared with diuretics and β-blockers on cardiovascular morbidity and mortality in hypertension: the Nordic Diltiazem (NORDIL) study , 2000, The Lancet.

[14]  M. Fay,et al.  Small‐Sample Adjustments for Wald‐Type Tests Using Sandwich Estimators , 2001, Biometrics.

[15]  T. Derouen,et al.  A Covariance Estimator for GEE with Improved Small‐Sample Properties , 2001, Biometrics.

[16]  R. Carroll,et al.  A Note on the Efficiency of Sandwich Covariance Matrix Estimation , 2001 .

[17]  W. Pan,et al.  Small‐sample adjustments in using the sandwich variance estimator in generalized estimating equations , 2002, Statistics in medicine.

[18]  Tom Greene,et al.  Effect of blood pressure lowering and antihypertensive drug class on progression of hypertensive kidney disease: results from the AASK trial. , 2002, JAMA.

[19]  R. B. Jackson,et al.  Rooting depths, lateral root spreads and below‐ground/above‐ground allometries of plants in water‐limited ecosystems , 2002 .

[20]  Bettina M. J. Engelbrecht,et al.  Comparative drought-resistance of seedlings of 28 species of co-occurring tropical woody plants , 2003, Oecologia.

[21]  Peter Sleight,et al.  Principal results of the Controlled Onset Verapamil Investigation of Cardiovascular End Points (CONVINCE) trial. , 2003, JAMA.

[22]  N. Neerchal,et al.  Small Sample Correction for the Variance of GEE Estimators , 2003 .

[23]  C. Leuschner,et al.  The relationship between maximal stomatal conductance and leaf traits in eight Southeast Asian early successional tree species , 2004 .

[24]  Q. Gao,et al.  Relationships between regional primary production and vegetation patterns , 2004 .

[25]  H. A. Mooney,et al.  Maximum rooting depth of vegetation types at the global scale , 1996, Oecologia.

[26]  Frederick C. Meinzer,et al.  Seasonal leaf dynamics across a tree density gradient in a Brazilian savanna , 2005, Oecologia.

[27]  Søren Højsgaard,et al.  The R Package geepack for Generalized Estimating Equations , 2005 .

[28]  Stephen Sitch,et al.  Effects of parameter uncertainties on the modeling of terrestrial biosphere dynamics , 2005 .

[29]  W. Rawls,et al.  Soil Water Characteristic Estimates by Texture and Organic Matter for Hydrologic Solutions , 2006 .

[30]  P. D’Odorico,et al.  A dynamic soil water threshold for vegetation water stress derived from stomatal conductance models , 2007 .