Front. Microbiol. Frontiers in Microbiology Front. Microbiol. 1664-302X Frontiers Media S.A. 10.3389/fmicb.2021.700863 Microbiology Original Research Evaluating the Cost of Pharmaceutical Purification for a Long-Duration Space Exploration Medical Foundry McNulty Matthew J. 1 2 Berliner Aaron J. 1 3 Negulescu Patrick G. 1 2 McKee Liber 2 Hart Olivia 2 Yates Kevin 1 2 Arkin Adam P. 1 3 Nandi Somen 1 2 4 McDonald Karen A. 1 2 4 * 1Center for the Utilization of Biological Engineering in Space (CUBES), Berkeley, CA, United States 2Department of Chemical Engineering, University of California, Davis, Davis, CA, United States 3Department of Bioengineering, University of California, Berkeley, Berkeley, CA, United States 4Global HealthShare Initiative, University of California, Davis, Davis, CA, United States

Edited by: Cyprien Verseux, University of Bremen, Germany

Reviewed by: Gisela Detrell, University of Stuttgart, Germany; Obulisamy Parthiba Karthikeyan, University of Houston, United States

*Correspondence: Karen A. McDonald, kamcdonald@ucdavis.edu

This article was submitted to Microbiotechnology, a section of the journal Frontiers in Microbiology

11 10 2021 2021 12 700863 27 04 2021 20 09 2021 Copyright © 2021 McNulty, Berliner, Negulescu, McKee, Hart, Yates, Arkin, Nandi and McDonald. 2021 McNulty, Berliner, Negulescu, McKee, Hart, Yates, Arkin, Nandi and McDonald

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

There are medical treatment vulnerabilities in longer-duration space missions present in the current International Space Station crew health care system with risks, arising from spaceflight-accelerated pharmaceutical degradation and resupply lag times. Bioregenerative life support systems may be a way to close this risk gap by leveraging in situ resource utilization (ISRU) to perform pharmaceutical synthesis and purification. Recent literature has begun to consider biological ISRU using microbes and plants as the basis for pharmaceutical life support technologies. However, there has not yet been a rigorous analysis of the processing and quality systems required to implement biologically produced pharmaceuticals for human medical treatment. In this work, we use the equivalent system mass (ESM) metric to evaluate pharmaceutical purification processing strategies for longer-duration space exploration missions. Monoclonal antibodies, representing a diverse therapeutic platform capable of treating multiple space-relevant disease states, were selected as the target products for this analysis. We investigate the ESM resource costs (mass, volume, power, cooling, and crew time) of an affinity-based capture step for monoclonal antibody purification as a test case within a manned Mars mission architecture. We compare six technologies (three biotic capture methods and three abiotic capture methods), optimize scheduling to minimize ESM for each technology, and perform scenario analysis to consider a range of input stream compositions and pharmaceutical demand. We also compare the base case ESM to scenarios of alternative mission configuration, equipment models, and technology reusability. Throughout the analyses, we identify key areas for development of pharmaceutical life support technology and improvement of the ESM framework for assessment of bioregenerative life support technologies.

techno-economic analysis equivalent system mass space exploration medical foundry in situ resource utilization pharmaceutical foundry monoclonal antibody purification human exploration mission space systems bioengineering NNX17AJ31G 80NSSC18K1157 NNX16AO69A National Aeronautics and Space Administration10.13039/100000104

香京julia种子在线播放

    1. <form id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></form>
      <address id=HxFbUHhlv><nobr id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></nobr></address>

      Introduction The Need for a Pharmaceutical Foundry in Space

      Surveying missions to Mars, like the InSight lander1 (Overview | Mission--NASA’s InSight Mars Lander) launched in 2018 and Perseverance rover2 in 2020, directly support the objectives of NASA’s long-term Mars Exploration Program3 : an effort to explore the potential for life on Mars and prepare for human exploration of Mars. The maturation of the program requires redefining the risks to human health as mission architectures transition from the current ‘‘Earth Reliant’’ paradigm used on the International Space Station (ISS) to the cislunar space ‘‘Proving Grounds’’ and finally to deep-space ‘‘Earth Independent’’ mission architectures, as defined in NASA’s report titled, ‘‘Journey to Mars: Pioneering Next Steps in Space Exploration.’’4

      Human missions to Mars will be ‘‘Earth Independent,’’ meaning there will be very limited emergency evacuation and re-supply capabilities along with substantially delayed communications with the Earth-based mission team. The NASA Human Research Roadmap5 currently rates most human health risks, which include “risk of adverse health outcomes and decrements in performance due to inflight medical conditions” and “risk of ineffective or toxic medications during long-duration exploration spaceflight,” as either medium or high risk for a Mars planetary visit/habitat mission. Risk ratings are based on failure mode and effects analysis and on hazard analysis using dimensions of severity, occurrence, and detectability. A recent review highlights the current understanding of the primary hazards and health risks posed by deep space exploration as well as the six types of countermeasures: protective shielding, biological and environmental temporal monitoring, specialized workout equipment, cognition and psychological evaluations, autonomous health support, and personalized medicine (Afshinnekoo et al., 2020).

      Of these countermeasures, it could be argued that medicine is the most crucial and least advanced toward mitigating space health hazards. There is very limited information on, and few direct studies of, pharmaceutical usage, stability, and therapeutic efficacy (i.e., pharmacokinetics, pharmacodynamics) in spaceflight or in a Mars surface environment (Blue et al., 2019). Furthermore, flown stores of pharmaceuticals face two additional barriers: (1) radiation-accelerated degradation (Du et al., 2011), and (2) addressing a myriad of low occurrence and high impact health hazards without the ability to fly and maintain potency of therapeutics for all of them. In these circumstances, it is often more beneficial to build robustness to these low occurrence health hazards rather than to try to predict them. It is therefore imperative that on-planet and/or in-flight pharmaceutical production be developed to bridge this risk gap. These pharmaceutical foundry technologies will supplement, not replace, the flown pharmaceutical formulary designed to treat anticipated medical threats during space missions.

      Pharmaceuticals are produced either chemically or biologically. A recent review of pharmaceutical production for human life support in space compares these two methods, highlighting the need for biological production in order to address many low occurrence and high impact health hazards (e.g., sepsis, ear infection, and glaucoma) and further comparing different biological production systems (McNulty et al., 2021). One major advantage of biological production is the efficiency in transporting and synthesizing genetic information as the set of instructions, or sometimes the product itself, to meet the therapeutic needs for a variety of disease states. The emerging field of Space Systems Bioengineering (Berliner et al., 2020) encapsulates this need for biological production, of which pharmaceuticals is identified as an important subset.

      The Bottleneck of Space Foundries: Purification

      Biopharmaceuticals must be purified after accumulation with the biological host organism, or cell-free transcription-translation reaction, in order to meet requirements for drug delivery and therapeutic effect (Harrison et al., 2015). The majority of commercial biopharmaceutical products are administered via intravenous and subcutaneous injection (Škalko-Basnet, 2014). Biopharmaceutical formulations for injection requires high purity (>95%) product, as impurities introduced directly into the bloodstream can trigger significant immune responses and reduce efficacy (Haile et al., 2015).

      Downstream processing of biopharmaceuticals is therefore usually a resource-intensive section of overall processing, being cited as high as 80% of production costs (and contributions of input mass) for monoclonal antibody (mAb) therapeutics produced using mammalian cell cultures (Rathore et al., 2015; Budzinski et al., 2019). In addition to the processing burden for biopharmaceutical injectables, there are also often substantial storage costs involving complex supply chain and storage management with stability requirements for factors including temperature, time, humidity, light, and vibration (Sykes, 2018). There are several approaches being pursued to overcome the challenges and costs associated with downstream processing and formulation.

      First are the tremendous efforts in process intensification (Strube et al., 2018). While the highly sensitive nature of biopharmaceuticals to minor process changes has introduced barriers and complexities to innovation through process intensification that have not been realized in non-healthcare biotechnological industries, there have been significant strides made in the past decade in the areas of process integration (Steinebach et al., 2017), automation (Pollard et al., 2017), and miniaturization (Adiga et al., 2018; Crowell et al., 2018).

      Another route that researchers are pursuing to reduce downstream processing costs and resources is a biological solution to processing technology. In the same vein that the biopharmaceutical industry sprung out of researchers leveraging the power of biology to produce therapeutically relevant molecules that were inaccessible or excessively costly by means of chemical synthesis, researchers are now also trying to apply that same principle to purifying therapeutically relevant molecules. The simplicity of production, reagents that can be produced using self-replicating organisms, and potential recyclability of spent consumables are significant advantages of biological purification technology for space or other limited resource applications. Examples of primary biological technologies include fusion tags (Bell et al., 2013), stimuli-responsive biopolymers (Sheth et al., 2014), hydrophobic nanoparticles (Jugler et al., 2020), and plant virus nanoparticles (Werner et al., 2006; Uhde-Holzem et al., 2016).

      Lastly, there are vast efforts to establish alternative drug delivery modalities (Anselmo et al., 2018). Other modalities that do not require injection and which might be more compatible to administration in limited resource environments, such as oral consumption, nasal spray, inhalation, and topical application, have long presented challenges in biopharmaceutical stability (e.g., denaturation in stomach acid) and delivery to the active site (e.g., passing the gut-blood barrier) that minimize product efficacy and necessitate costly advanced formulations and chemistries (Mitragotri et al., 2014).

      A particularly promising drug delivery technique to circumvent downstream processing burdens is to sequester the active pharmaceutical ingredient in the host cells of the upstream production system as a protective encapsulation in order to facilitate bioavailability through oral delivery (Kwon and Daniell, 2015). It represents an opportunity to greatly lower the cost of in situ production of human medicine for a space mission. This technique presumes that the host system is safe for human consumption, and so naturally lends itself to utility in systems such as yeast and plant production hosts. Oral delivery via host cell encapsulation has been recently established as commercial drug delivery modality with the US Food and Drug Administration approval of Palforzia as an oral peanut-protein immunotherapy (Vickery et al., 2018). However, this solution is not necessarily amenable to the diversity of pharmaceutical countermeasures that may be required, especially for unanticipated needs in which the product may not have been evaluated for oral bioavailability.

      Space Economics

      In 2011, the space shuttle program was retired due to increasing costs, demonstrating that reduction of economic cost is critical for sustaining any campaign of human exploration (Wall, 2011). Although recent efforts in reducing the launch cost to low earth orbit by commercial space companies have aided in the redefinition of the space economy (Whealan George, 2019), the barrier to longer term missions, such as a journey to Mars, is still limited by the extreme financial cost in transporting resources. Additionally, it has been shown that as the mission duration and complexity increases–as expected for a human mission to Mars–the quantity of supplies required to maintain crew health also increases (Anderson et al., 2018). In the case of meeting the demand for medication, biopharmaceutical synthesis has been proposed as an alternative to packaging a growing number of different medications (Menezes et al., 2015; McNulty et al., 2021). Assuming that both technologies can meet mission demand, selection of the production-based biotechnology platform will be dependent on its cost impact. It is therefore critical that the cost model of biopharmaceutical synthesis accounts for and minimizes the cost of any and all subprocesses, including those for purification.

      The current terrestrial biopharmaceutical synthesis cost model does not align with the needs for space exploration environments. For example, the literature highlights the high cost of Protein A affinity chromatography resin ($8,000–$15,000/L) and the need to reduce the price (Ramos-de-la-Peña et al., 2019). However, the purchase cost of chromatography resin is not nearly as critical in space environment applications where the major costs are more closely tied to the physical properties of the object (mass, volume, refrigeration requirements, etc.), as a result of fuel and payload limitations and the crew time required for operation (Jones, 2001). The distinct cost models of space and terrestrial biopharmaceutical production may increase the burden of identifying space-relevant processing technologies and may also limit direct transferability of terrestrial technologies without attention given to these areas.

      On the other hand, changing incentives structures relating to sustainability and the advent of new platform technologies are rapidly increasing alignment and the potential for technology crossover. For example, companies like On Demand Pharmaceuticals6, EQRx7, and the kenUP Foundation8, initiatives leading to industry adoption of environmental footprint metrics such as E-factor (Sheldon, 2007) and process mass intensity (PMI) (Budzinski et al., 2019), and diffusion from the adjacencies of green and white biotechnology (Tylecote, 2019) all promote development of accessible and sustainable technologies. As these trends pertain to space-relevant processes, these examples can also be viewed as driving more closed loop systems composed of simpler components.

      Reference Mission Architecture

      The evaluation of biopharmaceutical system cost for space applications requires the establishment of a reference mission architecture (RMA) as a means for describing the envelope of the mission scenario and distilling initial technology specifications which relate to the proposed subsystem in question (Drake and Watts, 2014). This RMA can be used to orient and define the specific mission elements that meet the mission requirements and factor into the calculations of cost for deploying biopharmaceutical technologies. Ultimately, the RMA provides the means to determine and compare cost given specification of mission scenarios that utilize the technology in question. We envision developing and integrating biotechnological capabilities back-ended by purification and quality systems into standard methods composed of a series of unit procedures that maintain astronaut health via the Environmental Control and Life Support Systems (ECLSS) (Hendrickx et al., 2006). In this study, we begin to build toward this vision by proposing a high-level RMA that specifies a biopharmaceutical demand partially fulfilled through biomanufacturing over the course of a defined production window.

      Equivalent System Mass

      In planning for future human exploration missions, technology choices and life-support systems specifications are often evaluated through the metric of the equivalent system mass (ESM) (Levri et al., 2003). Driven by the economic factor of cost in dollars required to transport mass into orbit, the ESM framework accounts for non-mass factors such as power, volume, and crew-time by relating them to mass through predetermined equivalency factors. ESM has been used to evaluate the mass of all of the resources of a larger system including water, shielding materials, agriculture and recycle loop closure. Currently, ESM remains the standard metric for evaluating advanced life support technology platforms (Hogan et al., 2000; Zabel, 2020). In the Space Systems Bioengineering context of realizing a biomanufactory on the surface of Mars (Berliner et al., 2020), recent advances in extending this metric have been proposed in the form of extended ESM which attempts to address complexities stemming from multiple transit and operations stages, as would be required to support a crewed mission to Mars (Berliner et al., 2021). It also accounts for uncertainties inherent in mission planning such as technology failures and their downstream effects as propagated through a mission such as refrigeration failures in systems housing medicine that requires specific cooling. Such advances in the ESM framework aid in the assessment of biopharmaceutical technologies as elements in the context of proposed ECLSS given the inherent stochastic nature of human health, especially in a space environment (Bizzarri et al., 2017). Here, we calculate ESM at multiple mission segments across which biopharmaceutical purification is deployed.

      Materials and Methods Unit Procedure Selection Protein A-Based Affinity Capture Step

      The medical significance of mAb therapies and the highly developed and specialized purification technology provide a fertile ground for techno-economic feasibility analysis of an in situ resource utilization (ISRU)-based pharmaceutical foundry for space. The first reason is that there are mAb therapies commercially approved or in development for multiple important disease states of spaceflight including osteoporosis (Faienza et al., 2018), migraines/headaches (Schuster and Rapoport, 2016), seizure (Zhao et al., 2017), pneumonia (Hua et al., 2014), ocular herpes (Krawczyk et al., 2015), otitis media (Iino et al., 2019), various oncological indications (Zahavi and Weiner, 2020), and fungal infections (Ulrich and Ebel, 2020). A second reason is that degradation products of mAb therapies are known to result in, not just reduced efficacy, but also deleterious effects (e.g., harmful immune reactions in patients) that further compound concerns of pharmaceutical stability over a long-duration mission (Laptoš and Omersel, 2018). Thirdly is that a common manufacturing system can be used to produce treatments for a variety of indications which is highly advantageous in mass and volume savings for spaceflight. And fourthly, the economic incentive of research into mAb purification technology has resulted in a plethora of technologies, enabling this analysis to include head-to-head comparisons between multiple mAb capture steps of different origins (e.g., biotic, abiotic) and different processing mechanisms (e.g., bind-and-elute mode liquid chromatography, precipitation). It is in comparing the differences between these technologies that we can uncover general insights into the desired components of a pharmaceutical foundry for space.

      Monoclonal antibody therapy is a platform technology that supports human health across a diversity of medical indications with a generally maintained molecular structure, in large part due to the coupling of high target selectivity in the two small and highly variable complementarity-determining regions located in the antigen-binding fragments (Goding, 1996) and control of the biological action on that target (i.e., effector function) through the generally conserved fragment crystallizable (Fc) region (Kang and Jung, 2019). This otherwise high structural fidelity conserved across mAb therapy products (which are primarily of the immunoglobulin G class) spans a wide variety of therapeutic indications and creates an opportunity for generic mAb production process flows, which include technologies devised specifically for mAb production (Sommerfeld and Strube, 2005). This specialized manufacturing, which is most notable in the use of the affinity capture step targeting the Fc region of an antibody with the use of the protein-based ligands derived from the Staphylococcus aureus Protein A molecule, can be tuned for highly efficient purification of mAb and antibody-derived (e.g., Fc-fusion protein) class molecules (Ramos-de-la-Peña et al., 2019). Therefore, we have decided to investigate the Protein A-based affinity capture step in isolation as a starting point for understanding the costs of a potential pharmaceutical foundry in space.

      It is worth noting that other similar protein ligands, such as Protein G and Protein L, are also widely used for their ability to capture different types of immunoglobulin classes and subclasses more efficiently (Choe et al., 2016).

      Abiotic and Biotic Protein A-Based Unit Procedures

      We chose to analyze six Protein A-based capture step procedures: three commercially available abiotic technologies [pre-packed chromatography (CHM), spin column (SPN), and magnetic bead (MAG)] and three development-stage biotic technologies [plant virus-based nanoparticle (VIN), elastin-like polypeptide (ELP), and oilbody-oleosin (OLE)] (Figure 1). Commercial technology procedures are based on product handbooks while the procedures of developing technologies, which we would classify as Technology Readiness Level 2 per NASA’s guidelines, are based on reports in literature. This set of procedures was selected to survey a wide range of operational modalities, technological chassis, and perceived advantages and disadvantages (Table 1).

      Monoclonal antibody production consists generically of product accumulation, clarification, initial purification, formulation, and fill & finish. Here we investigate six technologies for the capture step within the first purification step in a space mission context using extended equivalent system mass. The manufacturing origin of the capture reagent is denoted as either (A) abiotic or (B) biotic.

      List of Protein A-based monoclonal antibody capture step unit procedures included for analysis.

      Unit procedure ID Method Technology used References
      Pre-packed chromatography(A) Liquid chromatography Pre-packed HiTrap MabSelect SuRe column of novel alkali-tolerant recombinant Protein A-based ligand coupled with an agarose matrix Vendor handbooks (Cytiva, 2006, 2020, 2021)
      Spin column(A) Centrifuge-assisted liquid chromatography Pre-packed Protein A HP SpinTrap spin column containing Protein A Sepharose High Performance
      Magnetic bead(A) Magnetic separation Protein A Mag Sepharose superparamagnetic beads coupled with native Protein A ligands
      Plant virus-based nanoparticle(B) Sedimentation complex Plant virion, Turnip vein clearing virus, presenting a C-terminal coat protein fusion display of Protein A (domains D and E) Werner et al., 2006
      Elastin-like polypeptide(B) Inverse transition cycle Elastin-like polypeptides [78 pentapeptide (VPGVG) repeats] fused with Z domain, an engineered B domain of Protein A Sheth et al., 2014
      Oilbody-oleosin(B) Liquid-liquid partition Arabidopsis oleosin fused at the N-terminal with an engineered Protein A(5) McLean et al., 2012

      AAbiotic technology. BBiotic technology.

      All six of the unit procedures are operated in bind-and-elute mode, in which a clarified mAb-containing liquid stream is fed into a capture step containing Protein A-based ligand, which selectively binds the mAb and separates the mAb from the bulk feed stream. The mAb is eluted from the Protein A-based ligand and recovered using a low pH buffer to dissociate the mAb from the ligand. Finally, the low pH environment of the recovered mAb is pH neutralized for future processing or storage. The analysis does not consider differences in mAb processing upstream or downstream of the affinity capture step that may arise from differences in the unit procedure operations.

      Pre-packed chromatography is a chromatography system consisting of a liquid sample mobile phase which is pumped through a pre-packed bed of Protein A-fused resin beads housed in a column. SPN is a similar system, in which a Protein A-fused resin bead bed has been pre-packed into a plastic tube housing and the mobile phase flow is controlled via centrifugation of the plastic tube. MAG is a slurry-based magnetic separation system that uses superparamagnetic particles coated with Protein A-fused resin mixed as a slurry with the feed mAb stream for capture and elution of the mAb by magnet. VIN is a sedimentation-based system that uses plant virion-based chassis fused with Protein A-based ligands in suspension for capture of the mAb and centrifugation, assisted by the sedimentation velocity contribution of the chassis, to isolate and elute the mAb. ELP is a precipitation-based system that uses stimuli-responsive biopolymers fused with Protein A-based ligands in suspension for capture of the mAb and external stimuli (e.g., temperature, salt) to precipitate the bound complex and elute the mAb. OLE is a liquid-liquid partitioning system that uses oil phase segregating oleosin proteins fused with Protein A-based ligands to capture mAb in the oil phase and then elute the mAb into a clean aqueous phase.

      Techno-Economic Evaluation

      Techno-economic evaluations are performed using the recently proposed equations for ESM that include calculation of costs at each mission segment (Berliner et al., 2021). ESM for the mission ESM0 is defined as

      ESM 0 = k L e q , k i A k [ ( M k i M e q , k ) + ( V k i V e q , k ) + ( P k i P e q , k ) + ( C k i C e q , k ) + ( T i D k T e q , k ) ] = ESM 0 , p d + ESM 0 , t r 1 + ESM 0 , s f + ESM 0 , t r 2

      where Mi, Vi,Pi,Ci, Ti are the initial mass [kg], volume [m3], power requirement [kW], cooling requirement [kg/kW], and crew-time requirement [CM-h/h], Meq, Veq,Peq,Ceq, Teq are the equivalency factors for mass [kg/kg] (which is set to 1 in this study), volume [kg/m3], power [kg/kW], cooling [kg/kW], and crew time [kg/CM-h], respectively, Leq is the location equivalency factor [kg/kg] that accounts for costs associated with mass transport occurring at a particular mission segment (e.g., orbital maneuvers required for the return transit), and D is the duration of the mission segment [day] over a set of subsystems iA and set of mission segments k ∈ℳ. The mission ESM in this study is specifically defined as the sum of subtotal ESM for each mission segment within the scope of the RMA defined in this study (pre-deployment ESM0,pd, crewed transit to Mars ESM0,tr1, Mars surface operations ESM0,sf, and return crewed transit to Earth ESM0,tr2).

      Key mission and pharmaceutical assumptions are summarized in Table 2. The mission timeline depicted in Figure 2 provides insight into the proposed RMA and downstream crew needs and mAb production horizon. Here we assume a total mission duration of 910 days. First, a crew of 6 will travel from Earth to low Earth orbit, then board an interplanetary craft for a 210-day journey to Martian orbit, where the crew will descend to the surface in a separate craft, allowing the large transit vehicle to remain in orbit. Once on Mars, the crew will perform surface operations for 600 days. Following surface operations, the crew will leave Mars in a fueled ascent craft, board the interplanetary vehicle, and return to Earth orbit in 200 days. The mission timeline, crew size, and ESM equivalency factors are consistent with the recent RMA presented for inclusion of biomanufacturing elements (Berliner et al., 2021).

      Key mission and pharmaceutical reference mission architecture details and assumptions.

      Mission scope
      Pre-deployment N/A
      Transit to Mars 210 days
      Surface operations 500 days
      Return transit 200 days
      Total mission duration 910 days
      Crew size 6 crew members

      Pharmaceutical scope

      Mission demand, mAb 30,000 mg
      Biomanufacturing, mAb 10,000 mg
      Capture step recovery 98%
      Production window 600 days
      Feed mAb concentration 1 mg/mL
      Molecular weight, mAb 150 kDa

      mAb, monoclonal antibody.

      An illustration of the reference mission architecture in which (A) a crewed ship is launched from the surface of Earth and lands on Mars and (B) assembles a pre-deployed habitat on the Martian surface to perform operations before (C) a return transit to Earth on the same ship. Pharmaceutical needs are supported by flown stores until partway through surface operations, at which point needs are met by pharmaceuticals produced using in situ resource utilization. Production is initiated prior to the need window to ensure adequate stocks are generated by the time it is needed. Rocket artwork adapted from Musk (2017). Habitat artwork by Davian Ho.

      The mission demand for mAb therapies is assumed to be 30,000 mg over the entirety of the mission (supporting logic detailed in Supplementary Table 1). Pharmaceutical stores and production resources are assumed to be flown with the crew transit (no pre-deployment in order to maximize shelf-life). We assume that the production resources are stable throughout the mission duration. We conservatively assume (in the face of insufficient spaceflight stability data for biologics for a more refined estimate) that the first 600 days of pharmaceutical demand will be met through flown stores (20,000 mg), at which point pharmaceutical ISRU manufacturing is needed (10,000 mg) to alleviate the impact of accelerated pharmaceutical degradation and provide supplementary medication. The pharmaceutical production window opens prior to the ISRU demand timeframe and persists through a portion of the return transit (up to mission day 810) to reflect the expected life support advantage of maintaining capabilities to counter unanticipated needs or threats. We assume that the Protein A-based unit procedures consistently yield 98% recovery of mAb from the input stream.

      Unit Procedure Simulation

      Deterministic models for each unit procedure were developed in Microsoft Excel (Supplementary Data Sheet 2) using reference protocols cited in Table 1 as a series of executable operations, each containing a set of inputs defined by cost categories (labor, equipment, raw materials, and consumables) that are correspondingly populated with characteristic ESM constituent (mass, volume, power, cooling, and labor time) values (model composition illustrated in Supplementary Figure 1). Unit procedures have been defined as the smallest single execution (i.e., unit) of the secondary purification capture step procedure according to the reference protocol. We define the unit capacity by volume according to the equipment and consumables used (e.g., 2 mL maximum working volume in a 2 mL tube) and by mAb quantity according to the binding capacity for the given method (e.g., 1 mg mAb/mL resin) (Supplementary Table 2). Unit procedures with no explicit working volume constraints (i.e., the liquid solution volume for biotic technologies) have been defined with a maximum unit volume of 2 mL. ESM-relevant characteristics of individual inputs (e.g., equilibration buffer, 2 mL tube) are defined based on publicly available values, direct measurements taken, and assumptions (which are explicitly identified in the Supplementary Data Sheet 2).

      There are several model features that we have considered and decided not to include within the scope of analysis. Packing and containers for the inputs are not included for three reasons: (1) the contributions of the container are considered negligible as compared to the input itself (e.g., container holding 1 L buffer as compared to the 1 L of liquid buffer); (2) materials flown to space are often re-packaged with special considerations (Wotring, 2018); and (3) the selection of optimal container size is non-trivial and may risk obscuring more relevant ESM findings if not chosen carefully. We do not consider buffer preparation and assume the use of flown ready-to-use buffers and solutions. Furthermore, refrigeration costs of the input materials and costs that may be associated with establishing and maintaining a sterile operating environment (e.g., biosafety cabinet, 70% ethanol in spray bottles) are expected to be comparable between unit procedures and not considered. Impacts of microgravity on unit procedure execution are not considered for the return transit production. Refrigeration costs associated with low temperature equipment operation (e.g., centrifugation at 4°C) are included in the equipment power costs.

      Inputs common across unit procedures are standardized (Supplementary Table 3). One operational standardization is the inclusion of pH neutralization of the product stream following the low pH elution mechanism, which was explicitly stated in some procedures while not in others. Input quantities are scaled from a single unit to determine the number of units required to meet the RMA specifications. The ESM constituent inputs (mass, volume, power, cooling, and labor time) are converted into equivalent mass values using RMA equivalency factors (Supplementary Table 4).

      Results and Discussion Standardization of Manufacturing Efficiency

      Given the limited granularity of the presented RMA, which was scoped as such to reflect the lack of literature presenting an overarching and validated Concept of Operations for a Transit to Mars (Antonsen et al., 2017), we do not define strict manufacturing scheduling criteria for pharmaceutical production. Construction of a detailed pharmaceutical production RMA is hindered by uncertainty in the number and identity of mAb therapy products that would be included within mission scope, the decay rate of mAb therapy stores in the mission environments, and a reasonable basis for building robustness to unanticipated disease states. Rather, we choose to establish an objective comparison between unit procedures by normalizing for scheduling-associated manufacturing efficiencies. We accomplish this by first identifying the number of batches per mission (and thus batch size) needed to meet the mAb demand (base case of 10,204 mg mAb feed assuming 98% recovery) that minimizes the ESM output for a given unit procedure, and then running the simulation of pharmaceutical production at that number of mission batches, as shown in Figure 3A and tabulated in Supplementary Table 5.

      (A) Scheduling optimization for the establishment of base case scenarios for each unit procedure. The value for number of batches corresponding to the minimum equivalent system mass for each unit procedure, as indicated by black circle (∘) markers. Key operational parameters impacted by mission scheduling (shown using the VIN procedure) include (B) unit underutilization or vacancy, (C) equipment underutilization or vacancy, in this case represented by the centrifuge as the bottleneck, (D) the number of use cycles, and (E) the total quantity of monoclonal antibody (mAb) per mission and per surface operation (sf). CHM, pre-packed chromatography; SPN, spin column; MAG, magnetic bead; VIN, plant virus-based nanoparticle; ELP, elastin-like polypeptide; OLE, oilbody-oleosin.

      In Figures 3B–E, we visualize a deconstruction of ESM output, using the VIN unit procedure as an example, by key performance metrics that vary with a scheduling dependence in order to illustrate the significance of batch optimization in unit procedure comparison. The processing of a given batch volume and mAb quantity is allocated into a number of units, as determined by the volume and mAb quantity constraints of a given unit procedure, and a number of use cycles per batch, as determined by the capacity of the equipment specified in the given unit procedure. We show how the variation in ESM output over the number of mission batches maps to extent of unit vacancy or underutilization (Figure 3B), extent of operational equipment (e.g., centrifuge) vacancy or underutilization (Figure 3C), and number of required use cycles (Figure 3D). We also show an oscillatory behavior in the scheduling (i.e., total mAb purified per mission, % purified at surface operations) that quickly dampens as number of mission batches increases (Figure 3E). This behavior is a result of the assumption that the mAb feed stream is coming from a discrete upstream production batch (e.g., batch-mode bioreactor) that does not output partial batch quantities, as opposed to a continuous upstream production for which there are no defined batches. Accordingly, partial batch needs are met by the processing of a full batch.

      Base Case Scenario

      The ESM and output metrics of the base case scenario (10,000 mg mAb demand, 1 mg mAb/mL feed concentration, 98% recovery) for each of the six unit procedures are shown in Figures 4A–F. From this viewpoint of an ESM output for an isolated unit procedure outside the context of a full purification scheme, the ESM ranked from lowest to highest are VIN < SPN < OLE < CHM < MAG < OLE. However, we reason that it is more important to understand the model inputs that influence the ESM output rankings than to use the rankings in this isolated subsystem analysis to make technology selection choices, which requires the context of a full pharmaceutical foundry and of linkages to other mission elements.

      Base case equivalent system mass results broken down by (A) mass (M), volume (V), power (P), and labor time (T) constituents, (B) transit to Mars (tr1), surface operations (sf), and return transit (tr2) mission segments, and (C) labor (L), equipment (E), raw materials (R), and consumables (C) cost category for the six tested Protein A-based monoclonal antibody affinity capture step unit procedures segregated by abiotic (white background) and biotic (gray background) technologies. Also shown are the (D) labor and operation times, (E) number of use cycles, and (F) number of units required for each unit procedure to meet the reference mission demand. CHM, pre-packed chromatography; SPN, spin column; MAG, magnetic bead; VIN, plant virus-based nanoparticle; ELP, elastin-like polypeptide; OLE, oilbody-oleosin.

      We observe that mass costs are generally the primary contributor to ESM output, except for the MAG and ELP procedures in which labor time costs are larger. The mass costs are not closely associated to any given cost category across unit procedures, but rather the breakdown of mass costs varies widely by unit procedure.

      Power costs (kW) are disproportionately high given that the static nature of ESM assumes constant usage, and thus energy (kW⋅h) in this context (i.e., the power supply to the equipment is not turned off in this analysis). These costs represent an upper bound assuming that the power supply system capacity is sized to support a maximal power consumption in which all power-drawing elements are simultaneously in operation. Time of power usage as a fraction of duration are as follows: CHM (99%) > MAG (78%) > ELP (48%) > SPN (45%) > OLE (42%) > VIN (30%). The lower use fraction unit procedures are therefore paying a relatively higher cost per unit power demand in this current method. The electrical needs of the equipment used by the unit procedures are within NASA-proposed Mars mission RMA bounds, with energy use across all unit procedures would peak at ∼1% of a proposed Mars transfer vehicle electric capacity (50 kWe) or ∼5% of the habitat capacity (12 kWe) of a reference stationary surface nuclear fission power reactor (Drake et al., 2010).

      The mission segment breakdown of ESM illustrates the relatively high costs of pharmaceutical manufacturing capabilities for transit, even for the transit to Mars (tr1) in which there is no actual production taking place. There is a strong economic incentive to limit the amount of supplies flown on tr1. Alternatives such as the pre-deployment of reagents and consumables and limiting of production to surface operations on Mars (which has lower RMA equivalency factors for mass and volume than transit operations) must be balanced against the risk to human health posed by removing pharmaceutical production capabilities from a mission segment and potentially exposing the supplies to longer storage times that could challenge shelf lives.

      Labor and operation times are important parameters in the broader mission and pharmaceutical foundry context. These unit procedures represent a single step of pharmaceutical production, which if realized in a space mission context, would, in turn, need to be a small portion of a crew member’s time allocation. Assuming 40-h work weeks for crew members, the labor time spans a range of ∼1% (CHM) to ∼14% (ELP) of the available crew time over the 600-day production window. It is not feasible to operationalize with such high labor and operation times at this scale of production, particularly as they stand for MAG and ELP. While strategies such as batch staggering and concurrency can be used to reduce durations, advanced automation will almost certainly need to be built into the core of a pharmaceutical foundry.

      A prevailing trend throughout the unit procedures is that the number of unit executions and use cycles required by a given unit procedure are positive correlated with the ESM output value, except for the equipment cost-dominant and higher unit capacity CHM procedure. The equipment modeled in the analysis for CHM and the other unit procedures are almost certainly not space-ready and could be further designed to reduce mass and volume and increase automation to reduce crew labor time. The increased equipment costs in the CHM procedure are primarily due to automation and monitoring hardware for running liquid chromatography, which is reflected in the minimal labor costs of the CHM procedure. Miniaturization efforts, such as those focusing on microfluidic systems (Millet et al., 2015; Rodríguez-Ruiz et al., 2018; Murphy et al., 2019), are emerging as a potential path toward mitigating the high equipment costs associated with highly automated and tightly controlled manufacturing, which are crucial for freeing up valuable crew time.

      The number of unit executions is determined by the binding capacity of the technology and the nominal unit size. This indicates that the unit capacity for purification is an important consideration and influential factor. Unit sizing is an important consideration that is valuable to assess more holistically within the broader pharmaceutical production and mission context.

      The number of use cycles is determined by the number of unit executions required and by the maximal unit capacity of the equipment items (e.g., if you presume that an 18-slot centrifuge is the equipment bottleneck then the effective number of batches is the number of units required divided by 18). Therefore, it can be understood that the equipment unit capacity is a critical parameter in tuning the number of use cycles and, by extension, the labor costs. For processes with lower labor costs, due to the intrinsic nature of the procedure or through automation of labor, equipment unit capacity will still influence the total duration and production throughout. The MAG and ELP procedures yield both high labor and duration times and are thus particularly sensitive to the equipment capacity.

      Contextualizing Equivalent System Mass With Supporting Evaluations

      Having acknowledged shortcomings of ESM as a decision-making tool for comparison of alternative approaches in isolated subsystems, we propose that supplementary evaluations can assist in contextualization. A primary gap of an isolated subsystem ESM analysis is a lack of information on the holistic usefulness or cost of a given employed resource, which could include its synergy with other mission subsystems and its extent of recyclability, or waste loop closure, within the mission context. For example, the isolated subsystem analysis does not capture information on the broad applicability that a centrifuge might have for use in other scientific endeavors, nor do the ESM outputs reflect the >93% recyclability of water achieved by the recycler on the ISS (Steven Siceloff, 2008) that may be generalizable to future missions.

      The use of environmental footprint metrics, such as PMI, may be one valuable step toward capturing missed information on recyclability. PMI is a simple metric of material efficiency defined as the mass of raw materials and consumables required to produce 1 kg of active pharmaceutical ingredient. The study by Budzinski et al. (2019) introducing PMI for biopharmaceuticals presents data from 6 firms using small-scale (2,000–5,000 L reactor) and large-scale (12,000–20,000 L reactor) mAb manufacturing operations, finding an average 7,700 kg of input is required to produce 1 kg of mAb. Figure 5A presents PMI evaluation for the six capture steps included in analysis, which result in PMI outputs as low as 2,390 kg of input (CHM) and as high as 17,450 kg of input (MAG) per 1 kg of mAb. A comparison of these outputs to those of Budzinski et al. (2019) indicates that we may be observing roughly similar values after accounting for the high cost of initial purification in the study, representing ∼60% of the total PMI reported, the elevated feed mAb concentration (i.e., cell culture titer) of 1–5.5 g mAb/L, and adjustments for economies of scale when operating at such low cycle volumes (Figure 5B). Consumable costs appear to be the most sensitive to scale, which represents ∼1% total PMI on average in the values reported by Budzinski et al. (2019) and ranges from 35% (CHM) to 77% (OLE) here. Budzinski et al. (2019) also go one step further to distinguish water as a separate category from raw materials and report that >90% of the mass is due to water use. Here we assume pre-made buffers and do not directly add water in this study, so we refrain from a similar calculation, but it is worth noting that the extent of water use may also serve as a reasonable starting surrogate for extent of achievable recyclability in a space mission context.

      (A) Process mass intensity (PMI) evaluation of the unit procedures broken down by raw materials (R) and consumables (C) contributions. (B) Cycle volume for each unit procedure. CHM, pre-packed chromatography; SPN, spin column; MAG, magnetic bead; VIN, plant virus-based nanoparticle; ELP, elastin-like polypeptide; OLE, oilbody-oleosin.

      Scenario Analysis

      We analyzed the specific ESM output broken down by cost category for the six unit procedures over a range of input stream mAb concentrations (Figure 6) and mission demand for mAb (Figure 7). Specific ESM, termed cost of goods sold in traditional manufacturing analyses, is the ESM output required to produce 1 mg mAb. This is used in the scenario analyses to normalize ESM output across variation in mission demand for mAb. The optimal number of batches per mission was found and used for each unit procedure and scenario tested (Supplementary Tables 7, 8).

      Specific equivalent system mass (per unit mass monoclonal antibody produced) broken down by labor (L), equipment (E), raw materials (R), and consumables (C) cost categories as a function of feed monoclonal antibody (mAb) concentration for (A) CHM, (B) SPN, (C) MAG, (D) VIN, (E) ELP, and (F) OLE. CHM, pre-packed chromatography; SPN, spin column; MAG, magnetic bead; VIN, plant virus-based nanoparticle; ELP, elastin-like polypeptide; OLE, oilbody-oleosin.

      Specific equivalent system mass (per unit mass monoclonal antibody produced) broken down by labor (L), equipment (E), raw materials (R), and consumables (C) cost categories as a function of mission production demand for monoclonal antibody for (A) CHM, (B) SPN, (C) MAG, (D) VIN, (E) ELP, and (F) OLE. CHM, pre-packed chromatography; SPN, spin column; MAG, magnetic bead; VIN, plant virus-based nanoparticle; ELP, elastin-like polypeptide; OLE, oilbody-oleosin.

      We observe the general and expected trends that specific ESM decreases with an increasing feed stream mAb concentration and mission demand. The CHM procedure exhibits notably limited sensitivity to feed stream mAb concentration, which can be attributed to the equipment-dominated cost profile, fixed column size, and nature of the governing reference protocol that does not specify restrictions on sample load volume. Depending on the pre-treatment of the feed stream, it may be more reasonable to impose constraints on the sample load volume. In contrast, the specific ESM output of the CHM procedure is the most sensitive to mission mAb demand with higher demand increasingly offsetting the fixed capital costs. The CHM procedure is also the largest capacity unit modeled in the analysis (i.e., CHM capacity is 30 mg mAb/unit as compared to 2.7 mg mAb/unit for MAG, the next highest capacity unit) and is accordingly expected to scale well with demand.

      The SPN, ELP, OLE procedures exhibit behaviors in which the specific ESM output abruptly plateaus with an increasing feed stream mAb concentration. This observation can be attributed to the unit procedure operating in a mAb binding capacity-limited regime (as opposed to volume-limited for more dilute feeds) which also then controls and maintains unit procedure throughput (e.g., the ELP number of units, 37,044, and use cycles per mission, 2,058, is constant at and above 0.35 mg mAb/mL input stream concentration). This can be de-bottlenecked via technology (e.g., improved chemistry of the capture step unit leading to higher binding capacity) or methodology (e.g., increased concentration of the capture step unit leading to higher binding capacity) improvements.

      Low demand scenarios are particularly relevant for examination in a space health context, as small capacity redundant and emergency utility is a likely proving ground for inclusion of a space pharmaceutical foundry. At the lower boundary of the tested range (1,000 mg mAb/mission), we see the ESM outputs from lowest to highest are re-ordered as MAG < VIN < SPN < OLE < ELP < CHM. Minimization of equipment costs are particularly important in this regime, and it is observed that, indeed, the ESM output near completely aligned with the ranking of equipment cost (MAG < VIN < SPN < ELP < OLE < CHM). It is likely that other non-ESM factors such as integration with other flown elements will understandably influence the design and composition of early and low capacity flown pharmaceutical foundries.

      Alternate Scenarios Mission Configurations

      We explored variations to the base case RMA for all six unit procedures including scenarios in which the pharmaceutical manufacturing resources are shipped prior to the crew in pre-deployment, (+)pd, the production window has been truncated to close with the end of surface operations, (−)tr2, and a combination of the two prior modifications, (+)pd (−)tr2 (Figure 8). Costs of pre-deployment are included in the analyses and mission demand is kept constant regardless of the production window.

      Evaluation of extended equivalent system mass values in various mission configurations broken down by labor (L), equipment (E), raw materials (R), and consumables (C) cost categories cost category and mass (M), volume (V), power (P), and labor time (T) constituents for CHM, (A,G), SPN, (B,H), MAG, (C,I), VIN, (D,J), ELP (E,K), and OLE, (F,L). Configurations include the base case scenario of manufacturing resources flown with the crew for pharmaceutical production on the surface and return transit (Base), and alternatives in which the manufacturing resources are flown prior to the crew in pre-deployment, (+)pd, the production window is limited to surface operations, (–)tr2, and a combination of the two previously stated alternatives, (+)pd (–)tr2. CHM, pre-packed chromatography; SPN, spin column; MAG, magnetic bead; VIN, plant virus-based nanoparticle; ELP, elastin-like polypeptide; OLE, oilbody-oleosin.

      In all cases the ESM totals were reduced from the base case. Additionally, the general trend held that (−)tr2 scenario resulted in lower ESM totals than (+)pd scenario except for SPN, in which the increased raw material and consumable costs of (−)tr2 were sufficiently large to outweigh the reduction in equipment and labor costs of (+)pd. The combination (+)pd (−)tr2 scenario resulted in the lowest ESM totals at a fraction of the base case (as high as 39% reduction in SPN and as low as 21% reduction in ELP).

      Equipment and Unit Throughput

      Acknowledging the significance of the equipment capacity on ESM output, we further explored this contribution by comparing the base case ESM output of the centrifuge-utilizing procedures (SPN, VIN, ELP, and OLE) to that resulting from the use of alternative centrifuge models (Supplementary Table 9). This effectively results in a trade of equipment costs and batch throughput. The optimal number of batches per mission was found and used for each unit procedure and interval tested (Supplementary Table 10).

      We observe in Figure 9 that the ESM values increased with the size of the centrifuge model, 12-slot < 18-slot (base) < 48-slot. The labor and consumables savings of higher batch throughput were outweighed by the higher equipment costs (including higher power costs). Operation duration is an important metric relevant to a pharmaceutical foundry that is not well reflected in ESM that is also impacted by this alternative scenario. The exception to this trend is the 48-slot condition for the ELP procedure, in which a lower consumable cost related to the number of use cycles per mission (i.e., pipette tips, tubes, and gloves) sufficiently lowered the total ESM below the 18-slot condition.

      Changes in extended equivalent system mass values with different capacity centrifuge models broken down by labor (L), equipment (E), raw materials (R), and consumables (C) cost categories and mass (M), volume (V), power (P), and labor time (T) constituents for SPN (A,E), VIN (B,F), ELP (C,G), and OLE (D,H). SPN, spin column; VIN, plant virus-based nanoparticle; ELP, elastin-like polypeptide; OLE, oilbody-oleosin.

      Technology Reusability

      The number of use cycles for liquid chromatography resins is an important economic parameter in commercial pharmaceutical manufacturing (Pathak and Rathore, 2016). Here we explore the impact of use cycles on the CHM and ELP procedures in a space mission context, looking at no reuse nor regeneration operation of the purification technology, (−)Reuse, and at an increased number of use cycles, (+)Reuse (Figure 10).

      Changes in extended equivalent system mass values with reusability of purification technology broken down by labor (L), equipment (E), raw materials (R), and consumables (C) cost categories and mass (M), volume (V), power (P), and labor time (T) constituents for CHM (A,C), and ELP (B,D). (–)Reuse considers the technology as single-use and accordingly discards the unit procedure cleaning operations; (+) Reuse considers additional reuse cycles of the technology. CHM, pre-packed chromatography; ELP, elastin-like polypeptide.

      We observe that the terrestrial importance of use cycles does not prevail in this isolated ESM evaluation in a space context. The high purchase costs of resin are not considered in ESM and the impact of the reuse cycles is reduced to the mass and volume savings of the pre-packed column consumable. There is a minor decrease in ESM of the (+)Reuse over the base case scenario, but both of these result in substantially higher ESM than the (−)Reuse scenario, particularly for the ELP procedure, in which the regeneration operation has been removed in addition to the reusability of the technology.

      These results echo the trend of single-use technology in commercial biotechnology in which manufacturers look to disposable plastic bioreactor and buffer bags as a means to reduce cleaning and validation costs (Shukla and Gottschalk, 2013). It would be valuable to further consider the utilization of single-use technology in a space pharmaceutical foundry, and in other space systems bioengineering applications, but it is important to point out the limited scope of this ESM analysis. Here we reiterate that the single unit procedure scope establishes a modular basis for pharmaceutical foundry ESM evaluation but does not realize the true circular economy advantages of reuse, which may be considerable for the regeneration step, and of biological systems for production of the purification reagent in general.

      Conclusion and Future Directions

      In this study, we have introduced and applied the ESM framework to biopharmaceutical processing as a first step toward modeling and understanding the costs of Space Systems Bioengineering and, more specifically, of a long-duration space exploration medical foundry, which we believe may 1 day constitute a critical bioregenerative component of ECLSS for humans to be able to explore the surface of Mars. We have observed that the static behavior of ESM, while certainly maintaining usefulness in early stage analyses, may stymie later-stage analyses of bioregenerative life support technologies, which tend to behavior more dynamically than traditional abiotic counterparts. In the future, higher fidelity analyses may be performed using tools such as HabNet (Do et al., 2015), although the use of such dynamic mission design and modeling tools will require additional software engineering efforts. As it stands now, our techno-economic calculations both satisfy the three fundamental aspects for life support modeling (Jones, 2017) and provide helpful directions for future efforts to incorporate purification processes in space systems bioengineering.

      The mAb affinity capture step represented an ideal starting point for biopharmaceutical purification cost analysis given the breadth of the mAb treatments for space-important health indications, the fact that mAb purification is considered a platform technology, and the diversity of affinity capture technologies. However, there are additional processing categories, such as size exclusion, ion exchange, and hydrophobic interaction unit procedures, which could be similarly studied in isolation for their general relevance in biopharmaceutical manufacturing. Establishing a unit procedure knowledge base for space-relevant economics of biopharmaceutical purification would provide additional benefit to the community.

      We acknowledge that the ESM analysis performed in this study utilizes current Earth-based technologies, not Mars-designed processes, and that as technologies evolve and expand the analysis will need to be updated. The need to revisit and update ESM analyses periodically as technology develops is standard practice. This is well illustrated in a recent ESM analysis of plant lighting systems that compares solar fiber optics to photovoltaic-powered light emitting diode hybrid systems (Hardy et al., 2020). The study results reversed decade-old trade study outcomes in which solar fiber optics scored more favorably, citing rapid advances in solar photovoltaics and light emitting diode technologies.

      Furthermore, the analysis presented does not encapsulate potentially significant characteristics of the unit procedures at the interfaces of the upstream and downstream biomanufacturing elements. For example, at the upstream interface the biotic unit procedures (VIN, ELP, and OLE) have been reported in literature to be effective capture mechanisms in “dirtier” feed solutions, perhaps absolving the need for more complex pre-capture clarification steps by virtue of process integration. At the downstream end, the eluate of the CHM unit procedure can be directly fed to the subsequent processing step, which would be particularly amenable for other column-based unit procedures, resulting in lower labor time and manufacturing duration. We also do not account for the uncertainty in performance associated with the developmental state of the technology. There have been substantially lower research and development investments in the biotic technologies than in the commercially available abiotic technologies; one may reasonably assume that there is more potential for improvements through biotic unit procedure optimization, while also considering that a larger driving force in abiotic unit procedure optimization for commercial terrestrial operations may balance or outweigh this. Forecasting on the technology development dynamics in the context of these, and other, forces could provide significant additional insights.

      Several overarching lessons on the development required for deployment of pharmaceutical purification technology to support human health in space can be gleaned from the cost breakdown of the ESM framework employed in this study. The high mass costs for the mAb capture technologies investigated suggest strong incentives to pursue efforts in miniaturization to reduce not only equipment mass, but also reagent mass, as preparation for pharmaceutical foundries in space. The high labor costs and duration of some of the technologies studied likewise suggests that automatization of biopharmaceutical purification would be impactful. Automatization could also conceivably be valuable in reducing mass costs associated with manual manipulation, such as pipette tips and gloves, and those associated with ensuring sterile operation. We also underline the importance of scheduling and equipment sizing optimization; for example, the ESM penalty for capturing the mission demand of mAb with the VIN unit procedure yielded up to 40% higher total ESM for non-optimal scheduled manufacturing batches. Given the advantage of in situ manufacturing to respond to uncertainty in mission medicine demand, further research to explore scheduling and equipment sizing under uncertainty would provide valuable insight.

      There are a series of challenges facing pharmaceutical foundries in space beyond processing. Perhaps the most daunting of these is the incompatibility of existing pharmaceutical regulatory compliance frameworks with the design constraints of in situ manufacturing. There are currently dozens to hundreds of analytical tests required to confirm process and product quality prior to release of the pharmaceutical for administration to human patients (Morrow and Felcone, 2004), which translates into a highly burdensome cost for in situ manufacturing of pharmaceuticals in space. Fortunately, there is a strong and parallel terrestrial need to reduce the burden of regulatory compliance while maintaining standards of quality assurance and control for personalized medicine, an individualized and patient-specific approach to medical care with widespread support. As mentioned earlier, trends of distributed and sustainable biomanufacturing on Earth provide additional support for reducing ESM-relevant costs.

      The analyses presented in this study motivate future investigation into the ESM output of a complete pharmaceutical foundry for a more complete comparison to other ECLSS needs and subsequent formal evaluations of medical risk (i.e., loss of crew life, medical evacuation, crew health index, risk of radiation exposure-induced death from cancer) mitigation as a balance to the ESM costs. The Integrated Scalable Cyto-Technology system (Crowell et al., 2018), reported in literature as capable of “end-to-end production of hundreds to thousands of doses of clinical-quality protein biologics in about 3 d[ays],” is an automated and multiproduct pharmaceutical manufacturing system that may serve well as a starting point for a complete pharmaceutical foundry evaluation. While downstream costs are typically a large proportion of terrestrial biopharmaceutical production costs, they may represent an even higher proportion of the overall ESM costs. ESM is more closely aligned to PMI as a metric than to cost of goods sold in dollars, suggesting that downstream contributions to ESM may similarly dominate. Budzinski and team found that downstream operations contributed 82% of the total PMI for commercial mAb production (Budzinski et al., 2019).

      Assembly of a complete pharmaceutical foundry ESM model would also enable investigation of more nuanced RMA design considerations, such as those relating to the influence of a fixed set, or anticipated probability distribution, of pharmaceutical product diversity and batch size on optimal system composition to meet given medical risk thresholds.

      As stated in the original presentation of ESM theory and application, comparison of multiple approaches for a given subsystem with ESM, such as we are studying with the capture step of a mAb pharmaceutical foundry, should satisfy the same product quantity, product quality, reliability, and safety requirements (Levri et al., 2000). Of these assumptions, the product quality and safety requirements prove challenging for implementation in pharmaceutical foundry comparisons. It is worth noting that reliability is not considered in the scope of this preliminary study, given the varying technology readiness levels of the unit procedures, but that it should be included in future analyses of full purification schemes. By extension, the impact of microgravity and reduced gravity on reliability and unit operation performance, while not investigated in this study, is an important and complex consideration, that requires significant research to address. Similarly, stability of the production resources over the course of a mission duration should be further considered in future works. High product sensitivity to process changes, and the large battery of testing sometimes required to observe them (the extent of which will also change with the processes employed), creates a situation where ESM comparisons of pharmaceutical foundries that serve as technology decision making tools will absolutely need to meet this requirement, albeit at a considerable cost and/or complexity of execution.

      The assessment of equivalent safety requirements, to the best of the knowledge of the authors, has been approached thus far in an ad hoc and qualitative manner, relying on extensive subject manner expertise and working process knowledge. One promising route to strengthening these critically important safety assessments would be to implement a formal assessment framework based on the environmental, health, and safety (EHS) assessment proposed by Biwer and Heinzle (2004), in which process inputs/outputs are ranked based on a series of hazard impact categories (e.g., acute toxicity, raw material availability, global warming potential) and impact groups (e.g., resources, organism). The key to a systematic space health-centric safety assessment like this is to establish space-relevant EHS impact categories (e.g., planetary protection, crew and ship safety). An improvement of the EHS underpinnings has the potential to provide significant benefits to future ESM analyses in the increasingly complex mission architecture of longer-duration missions.

      Data Availability Statement

      The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

      Author Contributions

      MM, AB, SN, and KM designed the study. MM, AB, PN, LM, OH, and KY built the model used in the study. MM analyzed the data and wrote the first draft of the manuscript. All authors were involved in project discussion and interpretation and contributed to manuscript revision.

      Conflict of Interest

      The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

      Publisher’s Note

      All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

      Funding

      This material is based upon work supported by NASA under grant or cooperative agreement award number NNX17AJ31G. This work was also supported by a NASA Space Technology Research Fellowship (NASA grant number 80NSSC18K1157). This work is supported by the Translational Research Institute through NASA Cooperative Agreement NNX16AO69A. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Aeronautics and Space Administration (NASA) or the Translational Research Institute for Space Health (TRISH).

      The illustrative figures in part were created with BioRender.com.

      Supplementary Material

      The Supplementary Material for this article can be found online at: /articles/10.3389/fmicb.2021.700863/full#supplementary-material

      References Adiga R. Al-adhami M. Andar A. Borhani S. Brown S. Burgenson D. (2018). Point-of-care production of therapeutic proteins of good-manufacturing-practice quality. Nat. Biomed. Eng. 2 675686. 10.1038/s41551-018-0259-1 31015674 Afshinnekoo E. Scott R. T. MacKay M. J. Pariset E. Cekanaviciute E. Barker R. (2020). Fundamental biological features of spaceflight: advancing the field to enable deep-space exploration. Cell 183 11621184. 10.1016/j.cell.2020.10.050 33242416 Anderson M. S. Ewert M. K. Keener J. F. (2018). Life Support Baseline Values and Assumptions Document. NASA/TP-2015–218570/REV1. Houston, TX: NASA. Anselmo A. C. Gokarn Y. Mitragotri S. (2018). Non-invasive delivery strategies for biologics. Nat. Rev. Drug Discov. 18 1940. 10.1038/nrd.2018.183 30498202 Antonsen E. Bayuse T. Blue R. Daniels V. Hailey M. Hussey S. (2017). Evidence Report: Risk of Adverse Health Outcomes and Decrements in Performance due to In-Flight Medical Conditions Human Research Program Exploration Medical Capabilities Element Approved for Public Release. Houston, TX: NASA. Bell M. R. Engleka M. J. Malik A. Strickler J. E. (2013). To fuse or not to fuse: what is your purpose? Protein Sci. 22 14661477. 10.1002/pro.2356 24038604 Berliner A. J. Hilzinger J. M. Abel A. J. Mcnulty M. Makrygiorgos G. Averesch N. J. H. (2020). Towards a biomanufactory on mars. Preprint 10.20944/preprints202012.0714.v1 32283112 Berliner A. J. Makrygiorgos G. Hill A. (2021). Extension of Equivalent System Mass for Human Exploration missions on mars. Preprint 10.20944/preprints202101.0363.v1 32283112 Biwer A. Heinzle E. (2004). Environmental assessment in early process development. J. Chem. Technol. Biotechnol. 79 597609. 10.1002/jctb.1027 Bizzarri M. Masiello M. G. Cucina A. Guzzi R. (2017). Journey to mars: a biomedical challenge. Perspective on future human space flight. J. Biol. Sci. hypotheses Opin. 1 1526. 10.13133/2532-5876_2.6 Blue R. S. Bayuse T. M. Daniels V. R. Wotring V. E. Suresh R. Mulcahy R. A. (2019). Supplying a pharmacy for NASA exploration spaceflight: challenges and current understanding. NPJ Microgravity 5:14. 10.1038/s41526-019-0075-2 31231676 Budzinski K. Blewis M. Dahlin P. D’Aquila D. Esparza J. Gavin J. (2019). Introduction of a process mass intensity metric for biologics. Nat. Biotechnol. 49 3742. 10.1016/j.nbt.2018.07.005 30121383 Choe W. Durgannavar T. Chung S. (2016). Fc-binding ligands of immunoglobulin G: an overview of high affinity proteins and peptides. Materials (Basel) 9:994. 10.3390/ma9120994 28774114 Crowell L. E. Lu A. E. Love K. R. Stockdale A. Timmick S. M. Wu D. (2018). On-demand manufacturing of clinical-quality biopharmaceuticals. Nat. Biotechnol. 36:988. 10.1038/nbt.4262 30272677 Cytiva (2006). Protein A HP SpinTrap Product Booklet. Available online at: https://www.cytivalifesciences.co.jp/tech_support/manual/pdf/28906770.pdf (accessed March 28, 2021). Cytiva (2020). Protein A Mag Sepharose Xtra Protein G Mag Sepharose Xtra Affinity Chromatography instructions for Use. Available online at: https://cdn.cytivalifesciences.com/dmm3bwsv3/AssetStream.aspx?mediaformatid=10061&destinationid=10016&assetid=15782 (accessed March 28, 2021). Cytiva (2021). Affinity Chromatography Handbook, Vol. 1: Antibodies. Available online at: https://cdn.cytivalifesciences.com/dmm3bwsv3/AssetStream.aspx?mediaformatid=10061&destinationid=10016&assetid=11660 (accessed March 28, 2021). Do S. Owens A. de Weck O. (2015). “HabNet – an integrated habitation and supportability architecting and analysis environment,” in Proceedings of the 45th International Conference on Environmental Systems, (Bellevue). Drake B. G. Hoffman S. J. Beaty D. W. (2010). “Human exploration of mars, design reference architecture 5.0,” in Proceedings of the IEEE Aerospace Conference Proceedings, (Big Sky, MT), 10.1109/AERO.2010.5446736 Drake B. G. Watts K. D. (2014). Human Exploration of Mars Design Reference Architecture 5.0 Addendum #2. Available online at: http://www.sti.nasa.gov (accessed March 7, 2021). Du B. Daniels V. R. Vaksman Z. Boyd J. L. Crady C. Putcha L. (2011). Evaluation of physical and chemical changes in pharmaceuticals flown on space missions. AAPS J. 13 299308. 10.1208/s12248-011-9270-0 21479701 Faienza M. F. Chiarito M. D’amato G. Colaianni G. Colucci S. Grano M. (2018). Monoclonal antibodies for treating osteoporosis. Expert Opin. Biol. Ther. 18 149157. 10.1080/14712598.2018.1401607 29113523 Goding J. W. (1996). Monoclonal Antibodies, 3rd Edn. Cambridge, MA: Academic Press. Haile L. A. Puig M. Kelley-Baker L. Verthelyi D. (2015). Detection of innate immune response modulating impurities in therapeutic proteins. PLoS One 10:e0125078. 10.1371/journal.pone.0125078 25901912 Hardy J. M. Kusuma P. Bugbee B. Wheeler R. Ewert M. (2020). “Providing photons for food in regenerative life support: a comparative analysis of solar fiber optic and electric light systems,” in Proceedings of the 2020 International Conference on Environmental Systems, (Logan, UT). Harrison R. G. Todd P. W. Rudge S. R. Petrides D. P. (2015). in Bioseparations Science and Engineering, 2nd Edn, eds Gubbins K. E. Barteau M. A. Lauffenburger D. A. Morari M. Ray W. H. Russel W. B. (Oxford: Oxford University Press). Hendrickx L. De Wever H. Hermans V. Mastroleo F. Morin N. Wilmotte A. (2006). Microbial ecology of the closed artificial ecosystem MELiSSA (Micro-Ecological Life Support System Alternative): reinventing and compartmentalizing the Earth’s food and oxygen regeneration system for long-haul space exploration missions. Res. Microbiol. 157 7786. 10.1016/j.resmic.2005.06.014 16431089 Hogan J. A. Kang S. Cavazzoni J. Levri J. A. Finn C. Luna B. (2000). A Simulation Study Comparing Incineration and Composting in a Mars-Based Advanced Life Support System NASA Document ID 20000121172. Moffett Field, CA: NASA. Hua L. Hilliard J. J. Shi Y. Tkaczyk C. Cheng L. I. Yu X. (2014). Assessment of an anti-alpha-toxin monoclonal antibody for prevention and treatment of Staphylococcus aureus-induced pneumonia. Antimicrob. Agents Chemother. 58 11081117. 10.1128/AAC.02190-13 24295977 Iino Y. Takahashi E. Ida S. Kikuchi S. (2019). Clinical efficacy of anti-IL-5 monoclonal antibody mepolizumab in the treatment of eosinophilic otitis media. Auris Nasus Larynx 46 196203. 10.1016/j.anl.2018.07.011 30122651 Jones H. (2017). “How should life support be modeled and simulated?,” in Proccedings of the 47th International Conference on Environmental Systems, (Charleston. BC). Jones H. W. (2001). The cost and Equivalent System Mass of Space Crew Time SAE Technical Paper 2001-01-2359. Warrendale PA: SAE International, 10.4271/2001-01-2359 Jugler C. Joensuu J. Chen Q. (2020). Hydrophobin-protein a fusion protein produced in plants efficiently purified an anti-west nile virus monoclonal antibody from plant extracts via aqueous two-phase separation. Int. J. Mol. Sci. 21:2140. 10.3390/ijms21062140 32244994 Kang T. H. Jung S. T. (2019). Boosting therapeutic potency of antibodies by taming Fc domain functions. Exp. Mol. Med. 51:138. 10.1038/s12276-019-0345-9 31735912 Krawczyk A. Dirks M. Kasper M. Buch A. Dittmer U. Giebel B. (2015). Prevention of herpes simplex virus induced stromal keratitis by a glycoprotein b-specific monoclonal antibody. PLoS One 10:e0116800. 10.1371/journal.pone.0116800 25587898 Kwon K.-C. Daniell H. (2015). Low-cost oral delivery of protein drugs bioencapsulated in plant cells. Plant Biotechnol. J. 13 10171022. 10.1111/pbi.12462 26333301 Laptoš T. Omersel J. (2018). The importance of handling high-value biologicals: physico-chemical instability and immunogenicity of monoclonal antibodies. Exp. Ther. Med. 15 3161. 10.3892/ETM.2018.5821 29556253 Levri J. A. Fisher J. W. Jones H. W. Drysdale A. E. Ewert M. K. Hanford A. J. (2003). Advanced Life Support Equivalent System Mass Guidelines Document ALS Equivalent System Mass Guidelines DocumentNASA/TM-2003-212278. Moffett Field, CA: NASA. Levri J. A. Vaccari D. A. Drysdale A. E. (2000). Theory and Application of the Equivalent System Mass Metric Technical Paper 2000-01-2395. Warrendale, PA: SAE International, 10.4271/2000-01-2395 McLean M. D. Chen R. Yu D. Mah K.-Z. Teat J. Wang H. (2012). Purification of the therapeutic antibody trastuzumab from genetically modified plants using safflower Protein A-oleosin oilbody technology. Transgenic Res. 21 12911301. 10.1007/s11248-012-9603-5 22382463 McNulty M. J. Xiong Y. M. Yates K. Karuppanan K. Hilzinger J. M. Berliner A. J. (2021). Molecular pharming to support human life on the moon, mars, and beyond. Crit. Rev. Biotechnol. 41 849864. 10.1080/07388551.2021.1888070 33715563 Menezes A. A. Cumbers J. Hogan J. A. Arkin A. P. (2015). Towards synthetic biological approaches to resource utilization on space missions. J. R. Soc. Interface 12:20140715. 10.1098/rsif.2014.0715 25376875 Millet L. J. Lucheon J. D. Standaert R. F. Retterer S. T. Doktycz M. J. (2015). Modular microfluidics for point-of-care protein purifications. Lab Chip 15 17991811. 10.1039/c5lc00094g 25740172 Mitragotri S. Burke P. A. Langer R. (2014). Overcoming the challenges in administering biopharmaceuticals: formulation and delivery strategies. Nat. Rev. Drug Discov. 13 655672. 10.1038/nrd4363 25103255 Morrow T. Felcone L. H. (2004). Defining the difference: what Makes Biologics Unique. Biotechnol. Healthc. 1 2429. Murphy T. W. Sheng J. Naler L. B. Feng X. Lu C. (2019). On-chip manufacturing of synthetic proteins for point-of-care therapeutics. Microsystems Nanoeng. 5 112. 10.1038/s41378-019-0051-8 31057940 Musk E. (2017). Making humans a multi-planetary species. New Sp. 5 4661. 10.1089/space.2017.29009.emu Pathak M. Rathore A. S. (2016). Mechanistic understanding of fouling of protein A chromatography resin. J. Chromatogr. A 1459 7888. 10.1016/j.chroma.2016.06.084 27423774 Pollard D. Brower M. Richardson D. (2017). “Progress toward automated single-use continuous monoclonal antibody manufacturing via the protein refinery operations lab,” in Continuous Biomanufacturing - Innovative Technologies and Methods, ed Subramanian G. (Weinheim: Wiley-VCH Verlag GmbH & Co. KGaA), 107130. 10.1002/9783527699902.ch4 Ramos-de-la-Peña A. M. González-Valdez J. Aguilar O. (2019). Protein A chromatography: challenges and progress in the purification of monoclonal antibodies. J. Sep. Sci. 42 18161827. 10.1002/jssc.201800963 30811843 Rathore A. S. Pathak M. Ma G. Bracewell D. G. (2015). Re-use of Protein A resin: fouling and economics. BioPharm Int. 3:2015. Rodríguez-Ruiz I. Babenko V. Martínez-Rodríguez S. Gavira J. A. (2018). Protein separation under a microfluidic regime. Analyst 143 606619. 10.1039/c7an01568b 29214270 Schuster N. M. Rapoport A. M. (2016). New strategies for the treatment and prevention of primary headache disorders. Nat. Rev. Neurol. 12 635650. 10.1038/nrneurol.2016.143 27786243 Sheldon R. A. (2007). The E Factor: fifteen years on. Green Chem. 9 12731283. 10.1039/b713736m Sheth R. D. Jin M. Bhut B. V. Li Z. Chen W. Cramer S. M. (2014). Affinity precipitation of a monoclonal antibody from an industrial harvest feedstock using an ELP-Z stimuli responsive biopolymer. Biotechnol. Bioeng. 111 15951603. 10.1002/bit.25230 24595842 Shukla A. A. Gottschalk U. (2013). Single-use disposable technologies for biopharmaceutical manufacturing. Trends Biotechnol. 31 147154. 10.1016/j.tibtech.2012.10.004 23178074 Škalko-Basnet N. (2014). Biologics: the role of delivery systems in improved therapy. Biol. Targets Ther. 8 107114. 10.2147/BTT.S38387 24672225 Sommerfeld S. Strube J. (2005). Challenges in biotechnology production - Generic processes and process optimization for monoclonal antibodies. Chem. Eng. Process Process Intensif. 44 11231137. 10.1016/j.cep.2005.03.006 Steinebach F. Karst D. Morbidelli M. (2017). “Design of integrated continuous processes for high-quality biotherapeutics,” in Continuous Biomanufacturing - Innovative Technologies and Methods, ed Subramanian G. (Weinheim: Wiley-VCH Verlag GmbH & Co. KGaA), 457480. 10.1002/9783527699902.ch16 Steven Siceloff K. (2008). NASA - Recycling Water is not Just for Earth Anymore. Available online at: https://www.nasa.gov/mission_pages/station/behindscenes/waterrecycler.html (Accessed April 19, 2021). Strube J. Ditz R. Kornecki M. Huter M. Schmidt A. Thiess H. (2018). Process intensification in biologics manufacturing. Chem. Eng. Process. Process Intensif. 133 278293. 10.1016/J.CEP.2018.09.022 Sykes C. (2018). Time- and temperature-controlled transport: supply chain challenges and solutions. P T 43 154170. Tylecote A. (2019). Biotechnology as a new techno-economic paradigm that will help drive the world economy and mitigate climate change. Res. Policy 48 858868. 10.1016/j.respol.2018.10.001 Uhde-Holzem K. McBurney M. Tiu B. D. B. Advincula R. C. Fischer R. Commandeur U. (2016). Production of immunoabsorbent nanoparticles by displaying single-domain protein a on potato virus X. Macromol. Biosci. 16 231241. 10.1002/mabi.201500280 26440117 Ulrich S. Ebel F. (2020). Monoclonal antibodies as tools to combat fungal infections. J. Fungi 6 22. 10.3390/jof6010022 32033168 Vickery B. P. Vereda A. Casale T. B. Beyer K. du Toit G. O’B Hourihane J. (2018). AR101 oral immunotherapy for peanut allergy. N. Engl. J. Med. 379 19912001. 10.1056/NEJMoa1812856 30449234 Wall M. (2011). NASA’s Shuttle Program Cost $209 Billion - Was it Worth It? | Space. Space.com. Available online at: https://www.space.com/12166-space-shuttle-program-cost-promises-209-billion.html (accessed October 3, 2021) Werner S. Marillonnet S. Hause G. Klimyuk V. Gleba Y. (2006). Immunoabsorbent nanoparticles based on a tobamovirus displaying protein A. Proc. Natl. Acad. Sci. U.S.A. 103 1767817683. 10.1073/pnas.0608869103 17090664 Whealan George K. (2019). The economic impacts of the commercial space industry. Space Policy 47 181186. 10.1016/j.spacepol.2018.12.003 Wotring V. (2018). “Space pharmacology: how space affects pharmacology,” in Drug Discovery and Evaluation: Methods in Clinical Pharmacology, eds Hock F. Gralinski M. (Cham: Springer International Publishing), 113. 10.1007/978-3-319-56637-5_68-1 33311142 Zabel P. (2020). Influence of crop cultivation conditions on space greenhouse equivalent system mass. CEAS Sp. J. 13 315. 10.1007/s12567-020-00317-5 Zahavi D. Weiner L. (2020). Monoclonal antibodies in cancer therapy. Antibodies 9:34. 10.3390/antib9030034 32698317 Zhao J. Wang Y. Xu C. Liu K. Wang Y. Chen L. (2017). Therapeutic potential of an anti-high mobility group box-1 monoclonal antibody in epilepsy. Brain Behav. Immun. 64 308319. 10.1016/j.bbi.2017.02.002 28167116 Abbreviations CHM

      pre-packed chromatography

      EHS

      environmental, health, and safety

      ELP

      elastin-like polypeptide

      ESM

      equivalent system mass

      Fc

      fragment crystallizable

      ISRU

      in situ resource utilization

      ISS

      International Space Station

      mAb

      monoclonal antibody

      MAG

      magnetic bead

      PMI

      process mass intensity

      OLE

      oilbody-oleosin

      RMA

      reference mission architecture

      SPN

      spin column

      VIN

      plant virus-based nanoparticle.

      https://mars.nasa.gov/insight/mission/overview/

      https://mars.nasa.gov/mars2020/mission/overview/

      https://mars.nasa.gov/

      http://go.nasa.gov/1VHDXxg

      https://humanresearchroadmap.nasa.gov/Risks/

      ondemandpharma.com

      eqrx.com

      kenup.eu

      ‘Oh, my dear Thomas, you haven’t heard the terrible news then?’ she said. ‘I thought you would be sure to have seen it placarded somewhere. Alice went straight to her room, and I haven’t seen her since, though I repeatedly knocked at the door, which she has locked on the inside, and I’m sure it’s most unnatural of her not to let her own mother comfort her. It all happened in a moment: I have always said those great motor-cars shouldn’t be allowed to career about the streets, especially when they are all paved with cobbles as they are at Easton Haven, which are{331} so slippery when it’s wet. He slipped, and it went over him in a moment.’ My thanks were few and awkward, for there still hung to the missive a basting thread, and it was as warm as a nestling bird. I bent low--everybody was emotional in those days--kissed the fragrant thing, thrust it into my bosom, and blushed worse than Camille. "What, the Corner House victim? Is that really a fact?" "My dear child, I don't look upon it in that light at all. The child gave our picturesque friend a certain distinction--'My husband is dead, and this is my only child,' and all that sort of thing. It pays in society." leave them on the steps of a foundling asylum in order to insure [See larger version] Interoffice guff says you're planning definite moves on your own, J. O., and against some opposition. Is the Colonel so poor or so grasping—or what? Albert could not speak, for he felt as if his brains and teeth were rattling about inside his head. The rest of[Pg 188] the family hunched together by the door, the boys gaping idiotically, the girls in tears. "Now you're married." The host was called in, and unlocked a drawer in which they were deposited. The galleyman, with visible reluctance, arrayed himself in the garments, and he was observed to shudder more than once during the investiture of the dead man's apparel. HoME香京julia种子在线播放 ENTER NUMBET 0016fandual.com.cn
      www.lpchain.com.cn
      jlfydj.org.cn
      www.hilegua.com.cn
      www.lifushou.com.cn
      www.hzgfdz.com.cn
      tiyanjia.com.cn
      teaers.com.cn
      www.whrzum.com.cn
      www.wpobsd.com.cn
      处女被大鸡巴操 强奸乱伦小说图片 俄罗斯美女爱爱图 调教强奸学生 亚洲女的穴 夜来香图片大全 美女性强奸电影 手机版色中阁 男性人体艺术素描图 16p成人 欧美性爱360 电影区 亚洲电影 欧美电影 经典三级 偷拍自拍 动漫电影 乱伦电影 变态另类 全部电 类似狠狠鲁的网站 黑吊操白逼图片 韩国黄片种子下载 操逼逼逼逼逼 人妻 小说 p 偷拍10幼女自慰 极品淫水很多 黄色做i爱 日本女人人体电影快播看 大福国小 我爱肏屄美女 mmcrwcom 欧美多人性交图片 肥臀乱伦老头舔阴帝 d09a4343000019c5 西欧人体艺术b xxoo激情短片 未成年人的 插泰国人夭图片 第770弾み1 24p 日本美女性 交动态 eee色播 yantasythunder 操无毛少女屄 亚洲图片你懂的女人 鸡巴插姨娘 特级黄 色大片播 左耳影音先锋 冢本友希全集 日本人体艺术绿色 我爱被舔逼 内射 幼 美阴图 喷水妹子高潮迭起 和后妈 操逼 美女吞鸡巴 鸭个自慰 中国女裸名单 操逼肥臀出水换妻 色站裸体义术 中国行上的漏毛美女叫什么 亚洲妹性交图 欧美美女人裸体人艺照 成人色妹妹直播 WWW_JXCT_COM r日本女人性淫乱 大胆人艺体艺图片 女同接吻av 碰碰哥免费自拍打炮 艳舞写真duppid1 88电影街拍视频 日本自拍做爱qvod 实拍美女性爱组图 少女高清av 浙江真实乱伦迅雷 台湾luanlunxiaoshuo 洛克王国宠物排行榜 皇瑟电影yy频道大全 红孩儿连连看 阴毛摄影 大胆美女写真人体艺术摄影 和风骚三个媳妇在家做爱 性爱办公室高清 18p2p木耳 大波撸影音 大鸡巴插嫩穴小说 一剧不超两个黑人 阿姨诱惑我快播 幼香阁千叶县小学生 少女妇女被狗强奸 曰人体妹妹 十二岁性感幼女 超级乱伦qvod 97爱蜜桃ccc336 日本淫妇阴液 av海量资源999 凤凰影视成仁 辰溪四中艳照门照片 先锋模特裸体展示影片 成人片免费看 自拍百度云 肥白老妇女 女爱人体图片 妈妈一女穴 星野美夏 日本少女dachidu 妹子私处人体图片 yinmindahuitang 舔无毛逼影片快播 田莹疑的裸体照片 三级电影影音先锋02222 妻子被外国老头操 观月雏乃泥鳅 韩国成人偷拍自拍图片 强奸5一9岁幼女小说 汤姆影院av图片 妹妹人艺体图 美女大驱 和女友做爱图片自拍p 绫川まどか在线先锋 那么嫩的逼很少见了 小女孩做爱 处女好逼连连看图图 性感美女在家做爱 近距离抽插骚逼逼 黑屌肏金毛屄 日韩av美少女 看喝尿尿小姐日逼色色色网图片 欧美肛交新视频 美女吃逼逼 av30线上免费 伊人在线三级经典 新视觉影院t6090影院 最新淫色电影网址 天龙影院远古手机版 搞老太影院 插进美女的大屁股里 私人影院加盟费用 www258dd 求一部电影里面有一个二猛哥 深肛交 日本萌妹子人体艺术写真图片 插入屄眼 美女的木奶 中文字幕黄色网址影视先锋 九号女神裸 和骚人妻偷情 和潘晓婷做爱 国模大尺度蜜桃 欧美大逼50p 西西人体成人 李宗瑞继母做爱原图物处理 nianhuawang 男鸡巴的视屏 � 97免费色伦电影 好色网成人 大姨子先锋 淫荡巨乳美女教师妈妈 性nuexiaoshuo WWW36YYYCOM 长春继续给力进屋就操小女儿套干破内射对白淫荡 农夫激情社区 日韩无码bt 欧美美女手掰嫩穴图片 日本援交偷拍自拍 入侵者日本在线播放 亚洲白虎偷拍自拍 常州高见泽日屄 寂寞少妇自卫视频 人体露逼图片 多毛外国老太 变态乱轮手机在线 淫荡妈妈和儿子操逼 伦理片大奶少女 看片神器最新登入地址sqvheqi345com账号群 麻美学姐无头 圣诞老人射小妞和强奸小妞动话片 亚洲AV女老师 先锋影音欧美成人资源 33344iucoom zV天堂电影网 宾馆美女打炮视频 色五月丁香五月magnet 嫂子淫乱小说 张歆艺的老公 吃奶男人视频在线播放 欧美色图男女乱伦 avtt2014ccvom 性插色欲香影院 青青草撸死你青青草 99热久久第一时间 激情套图卡通动漫 幼女裸聊做爱口交 日本女人被强奸乱伦 草榴社区快播 2kkk正在播放兽骑 啊不要人家小穴都湿了 www猎奇影视 A片www245vvcomwwwchnrwhmhzcn 搜索宜春院av wwwsee78co 逼奶鸡巴插 好吊日AV在线视频19gancom 熟女伦乱图片小说 日本免费av无码片在线开苞 鲁大妈撸到爆 裸聊官网 德国熟女xxx 新不夜城论坛首页手机 女虐男网址 男女做爱视频华为网盘 激情午夜天亚洲色图 内裤哥mangent 吉沢明歩制服丝袜WWWHHH710COM 屌逼在线试看 人体艺体阿娇艳照 推荐一个可以免费看片的网站如果被QQ拦截请复制链接在其它浏览器打开xxxyyy5comintr2a2cb551573a2b2e 欧美360精品粉红鲍鱼 教师调教第一页 聚美屋精品图 中韩淫乱群交 俄罗斯撸撸片 把鸡巴插进小姨子的阴道 干干AV成人网 aolasoohpnbcn www84ytom 高清大量潮喷www27dyycom 宝贝开心成人 freefronvideos人母 嫩穴成人网gggg29com 逼着舅妈给我口交肛交彩漫画 欧美色色aV88wwwgangguanscom 老太太操逼自拍视频 777亚洲手机在线播放 有没有夫妻3p小说 色列漫画淫女 午间色站导航 欧美成人处女色大图 童颜巨乳亚洲综合 桃色性欲草 色眯眯射逼 无码中文字幕塞外青楼这是一个 狂日美女老师人妻 爱碰网官网 亚洲图片雅蠛蝶 快播35怎么搜片 2000XXXX电影 新谷露性家庭影院 深深候dvd播放 幼齿用英语怎么说 不雅伦理无需播放器 国外淫荡图片 国外网站幼幼嫩网址 成年人就去色色视频快播 我鲁日日鲁老老老我爱 caoshaonvbi 人体艺术avav 性感性色导航 韩国黄色哥来嫖网站 成人网站美逼 淫荡熟妇自拍 欧美色惰图片 北京空姐透明照 狼堡免费av视频 www776eom 亚洲无码av欧美天堂网男人天堂 欧美激情爆操 a片kk266co 色尼姑成人极速在线视频 国语家庭系列 蒋雯雯 越南伦理 色CC伦理影院手机版 99jbbcom 大鸡巴舅妈 国产偷拍自拍淫荡对话视频 少妇春梦射精 开心激动网 自拍偷牌成人 色桃隐 撸狗网性交视频 淫荡的三位老师 伦理电影wwwqiuxia6commqiuxia6com 怡春院分站 丝袜超短裙露脸迅雷下载 色制服电影院 97超碰好吊色男人 yy6080理论在线宅男日韩福利大全 大嫂丝袜 500人群交手机在线 5sav 偷拍熟女吧 口述我和妹妹的欲望 50p电脑版 wwwavtttcon 3p3com 伦理无码片在线看 欧美成人电影图片岛国性爱伦理电影 先锋影音AV成人欧美 我爱好色 淫电影网 WWW19MMCOM 玛丽罗斯3d同人动画h在线看 动漫女孩裸体 超级丝袜美腿乱伦 1919gogo欣赏 大色逼淫色 www就是撸 激情文学网好骚 A级黄片免费 xedd5com 国内的b是黑的 快播美国成年人片黄 av高跟丝袜视频 上原保奈美巨乳女教师在线观看 校园春色都市激情fefegancom 偷窥自拍XXOO 搜索看马操美女 人本女优视频 日日吧淫淫 人妻巨乳影院 美国女子性爱学校 大肥屁股重口味 啪啪啪啊啊啊不要 操碰 japanfreevideoshome国产 亚州淫荡老熟女人体 伦奸毛片免费在线看 天天影视se 樱桃做爱视频 亚卅av在线视频 x奸小说下载 亚洲色图图片在线 217av天堂网 东方在线撸撸-百度 幼幼丝袜集 灰姑娘的姐姐 青青草在线视频观看对华 86papa路con 亚洲1AV 综合图片2区亚洲 美国美女大逼电影 010插插av成人网站 www色comwww821kxwcom 播乐子成人网免费视频在线观看 大炮撸在线影院 ,www4KkKcom 野花鲁最近30部 wwwCC213wapwww2233ww2download 三客优最新地址 母亲让儿子爽的无码视频 全国黄色片子 欧美色图美国十次 超碰在线直播 性感妖娆操 亚洲肉感熟女色图 a片A毛片管看视频 8vaa褋芯屑 333kk 川岛和津实视频 在线母子乱伦对白 妹妹肥逼五月 亚洲美女自拍 老婆在我面前小说 韩国空姐堪比情趣内衣 干小姐综合 淫妻色五月 添骚穴 WM62COM 23456影视播放器 成人午夜剧场 尼姑福利网 AV区亚洲AV欧美AV512qucomwwwc5508com 经典欧美骚妇 震动棒露出 日韩丝袜美臀巨乳在线 av无限吧看 就去干少妇 色艺无间正面是哪集 校园春色我和老师做爱 漫画夜色 天海丽白色吊带 黄色淫荡性虐小说 午夜高清播放器 文20岁女性荫道口图片 热国产热无码热有码 2015小明发布看看算你色 百度云播影视 美女肏屄屄乱轮小说 家族舔阴AV影片 邪恶在线av有码 父女之交 关于处女破处的三级片 极品护士91在线 欧美虐待女人视频的网站 享受老太太的丝袜 aaazhibuo 8dfvodcom成人 真实自拍足交 群交男女猛插逼 妓女爱爱动态 lin35com是什么网站 abp159 亚洲色图偷拍自拍乱伦熟女抠逼自慰 朝国三级篇 淫三国幻想 免费的av小电影网站 日本阿v视频免费按摩师 av750c0m 黄色片操一下 巨乳少女车震在线观看 操逼 免费 囗述情感一乱伦岳母和女婿 WWW_FAMITSU_COM 偷拍中国少妇在公车被操视频 花也真衣论理电影 大鸡鸡插p洞 新片欧美十八岁美少 进击的巨人神thunderftp 西方美女15p 深圳哪里易找到老女人玩视频 在线成人有声小说 365rrr 女尿图片 我和淫荡的小姨做爱 � 做爱技术体照 淫妇性爱 大学生私拍b 第四射狠狠射小说 色中色成人av社区 和小姨子乱伦肛交 wwwppp62com 俄罗斯巨乳人体艺术 骚逼阿娇 汤芳人体图片大胆 大胆人体艺术bb私处 性感大胸骚货 哪个网站幼女的片多 日本美女本子把 色 五月天 婷婷 快播 美女 美穴艺术 色百合电影导航 大鸡巴用力 孙悟空操美少女战士 狠狠撸美女手掰穴图片 古代女子与兽类交 沙耶香套图 激情成人网区 暴风影音av播放 动漫女孩怎么插第3个 mmmpp44 黑木麻衣无码ed2k 淫荡学姐少妇 乱伦操少女屄 高中性爱故事 骚妹妹爱爱图网 韩国模特剪长发 大鸡巴把我逼日了 中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片 大胆女人下体艺术图片 789sss 影音先锋在线国内情侣野外性事自拍普通话对白 群撸图库 闪现君打阿乐 ady 小说 插入表妹嫩穴小说 推荐成人资源 网络播放器 成人台 149大胆人体艺术 大屌图片 骚美女成人av 春暖花开春色性吧 女亭婷五月 我上了同桌的姐姐 恋夜秀场主播自慰视频 yzppp 屄茎 操屄女图 美女鲍鱼大特写 淫乱的日本人妻山口玲子 偷拍射精图 性感美女人体艺木图片 种马小说完本 免费电影院 骑士福利导航导航网站 骚老婆足交 国产性爱一级电影 欧美免费成人花花性都 欧美大肥妞性爱视频 家庭乱伦网站快播 偷拍自拍国产毛片 金发美女也用大吊来开包 缔D杏那 yentiyishu人体艺术ytys WWWUUKKMCOM 女人露奶 � 苍井空露逼 老荡妇高跟丝袜足交 偷偷和女友的朋友做爱迅雷 做爱七十二尺 朱丹人体合成 麻腾由纪妃 帅哥撸播种子图 鸡巴插逼动态图片 羙国十次啦中文 WWW137AVCOM 神斗片欧美版华语 有气质女人人休艺术 由美老师放屁电影 欧美女人肉肏图片 白虎种子快播 国产自拍90后女孩 美女在床上疯狂嫩b 饭岛爱最后之作 幼幼强奸摸奶 色97成人动漫 两性性爱打鸡巴插逼 新视觉影院4080青苹果影院 嗯好爽插死我了 阴口艺术照 李宗瑞电影qvod38 爆操舅母 亚洲色图七七影院 被大鸡巴操菊花 怡红院肿么了 成人极品影院删除 欧美性爱大图色图强奸乱 欧美女子与狗随便性交 苍井空的bt种子无码 熟女乱伦长篇小说 大色虫 兽交幼女影音先锋播放 44aad be0ca93900121f9b 先锋天耗ばさ无码 欧毛毛女三级黄色片图 干女人黑木耳照 日本美女少妇嫩逼人体艺术 sesechangchang 色屄屄网 久久撸app下载 色图色噜 美女鸡巴大奶 好吊日在线视频在线观看 透明丝袜脚偷拍自拍 中山怡红院菜单 wcwwwcom下载 骑嫂子 亚洲大色妣 成人故事365ahnet 丝袜家庭教mp4 幼交肛交 妹妹撸撸大妈 日本毛爽 caoprom超碰在email 关于中国古代偷窥的黄片 第一会所老熟女下载 wwwhuangsecome 狼人干综合新地址HD播放 变态儿子强奸乱伦图 强奸电影名字 2wwwer37com 日本毛片基地一亚洲AVmzddcxcn 暗黑圣经仙桃影院 37tpcocn 持月真由xfplay 好吊日在线视频三级网 我爱背入李丽珍 电影师傅床戏在线观看 96插妹妹sexsex88com 豪放家庭在线播放 桃花宝典极夜著豆瓜网 安卓系统播放神器 美美网丝袜诱惑 人人干全免费视频xulawyercn av无插件一本道 全国色五月 操逼电影小说网 good在线wwwyuyuelvcom www18avmmd 撸波波影视无插件 伊人幼女成人电影 会看射的图片 小明插看看 全裸美女扒开粉嫩b 国人自拍性交网站 萝莉白丝足交本子 七草ちとせ巨乳视频 摇摇晃晃的成人电影 兰桂坊成社人区小说www68kqcom 舔阴论坛 久撸客一撸客色国内外成人激情在线 明星门 欧美大胆嫩肉穴爽大片 www牛逼插 性吧星云 少妇性奴的屁眼 人体艺术大胆mscbaidu1imgcn 最新久久色色成人版 l女同在线 小泽玛利亚高潮图片搜索 女性裸b图 肛交bt种子 最热门有声小说 人间添春色 春色猜谜字 樱井莉亚钢管舞视频 小泽玛利亚直美6p 能用的h网 还能看的h网 bl动漫h网 开心五月激 东京热401 男色女色第四色酒色网 怎么下载黄色小说 黄色小说小栽 和谐图城 乐乐影院 色哥导航 特色导航 依依社区 爱窝窝在线 色狼谷成人 91porn 包要你射电影 色色3A丝袜 丝袜妹妹淫网 爱色导航(荐) 好男人激情影院 坏哥哥 第七色 色久久 人格分裂 急先锋 撸撸射中文网 第一会所综合社区 91影院老师机 东方成人激情 怼莪影院吹潮 老鸭窝伊人无码不卡无码一本道 av女柳晶电影 91天生爱风流作品 深爱激情小说私房婷婷网 擼奶av 567pao 里番3d一家人野外 上原在线电影 水岛津实透明丝袜 1314酒色 网旧网俺也去 0855影院 在线无码私人影院 搜索 国产自拍 神马dy888午夜伦理达达兔 农民工黄晓婷 日韩裸体黑丝御姐 屈臣氏的燕窝面膜怎么样つぼみ晶エリーの早漏チ○ポ强化合宿 老熟女人性视频 影音先锋 三上悠亚ol 妹妹影院福利片 hhhhhhhhsxo 午夜天堂热的国产 强奸剧场 全裸香蕉视频无码 亚欧伦理视频 秋霞为什么给封了 日本在线视频空天使 日韩成人aⅴ在线 日本日屌日屄导航视频 在线福利视频 日本推油无码av magnet 在线免费视频 樱井梨吮东 日本一本道在线无码DVD 日本性感诱惑美女做爱阴道流水视频 日本一级av 汤姆avtom在线视频 台湾佬中文娱乐线20 阿v播播下载 橙色影院 奴隶少女护士cg视频 汤姆在线影院无码 偷拍宾馆 业面紧急生级访问 色和尚有线 厕所偷拍一族 av女l 公交色狼优酷视频 裸体视频AV 人与兽肉肉网 董美香ol 花井美纱链接 magnet 西瓜影音 亚洲 自拍 日韩女优欧美激情偷拍自拍 亚洲成年人免费视频 荷兰免费成人电影 深喉呕吐XXⅩX 操石榴在线视频 天天色成人免费视频 314hu四虎 涩久免费视频在线观看 成人电影迅雷下载 能看见整个奶子的香蕉影院 水菜丽百度影音 gwaz079百度云 噜死你们资源站 主播走光视频合集迅雷下载 thumbzilla jappen 精品Av 古川伊织star598在线 假面女皇vip在线视频播放 国产自拍迷情校园 啪啪啪公寓漫画 日本阿AV 黄色手机电影 欧美在线Av影院 华裔电击女神91在线 亚洲欧美专区 1日本1000部免费视频 开放90后 波多野结衣 东方 影院av 页面升级紧急访问每天正常更新 4438Xchengeren 老炮色 a k福利电影 色欲影视色天天视频 高老庄aV 259LUXU-683 magnet 手机在线电影 国产区 欧美激情人人操网 国产 偷拍 直播 日韩 国内外激情在线视频网给 站长统计一本道人妻 光棍影院被封 紫竹铃取汁 ftp 狂插空姐嫩 xfplay 丈夫面前 穿靴子伪街 XXOO视频在线免费 大香蕉道久在线播放 电棒漏电嗨过头 充气娃能看下毛和洞吗 夫妻牲交 福利云点墦 yukun瑟妃 疯狂交换女友 国产自拍26页 腐女资源 百度云 日本DVD高清无码视频 偷拍,自拍AV伦理电影 A片小视频福利站。 大奶肥婆自拍偷拍图片 交配伊甸园 超碰在线视频自拍偷拍国产 小热巴91大神 rctd 045 类似于A片 超美大奶大学生美女直播被男友操 男友问 你的衣服怎么脱掉的 亚洲女与黑人群交视频一 在线黄涩 木内美保步兵番号 鸡巴插入欧美美女的b舒服 激情在线国产自拍日韩欧美 国语福利小视频在线观看 作爱小视颍 潮喷合集丝袜无码mp4 做爱的无码高清视频 牛牛精品 伊aⅤ在线观看 savk12 哥哥搞在线播放 在线电一本道影 一级谍片 250pp亚洲情艺中心,88 欧美一本道九色在线一 wwwseavbacom色av吧 cos美女在线 欧美17,18ⅹⅹⅹ视频 自拍嫩逼 小电影在线观看网站 筱田优 贼 水电工 5358x视频 日本69式视频有码 b雪福利导航 韩国女主播19tvclub在线 操逼清晰视频 丝袜美女国产视频网址导航 水菜丽颜射房间 台湾妹中文娱乐网 风吟岛视频 口交 伦理 日本熟妇色五十路免费视频 A级片互舔 川村真矢Av在线观看 亚洲日韩av 色和尚国产自拍 sea8 mp4 aV天堂2018手机在线 免费版国产偷拍a在线播放 狠狠 婷婷 丁香 小视频福利在线观看平台 思妍白衣小仙女被邻居强上 萝莉自拍有水 4484新视觉 永久发布页 977成人影视在线观看 小清新影院在线观 小鸟酱后丝后入百度云 旋风魅影四级 香蕉影院小黄片免费看 性爱直播磁力链接 小骚逼第一色影院 性交流的视频 小雪小视频bd 小视频TV禁看视频 迷奸AV在线看 nba直播 任你在干线 汤姆影院在线视频国产 624u在线播放 成人 一级a做爰片就在线看狐狸视频 小香蕉AV视频 www182、com 腿模简小育 学生做爱视频 秘密搜查官 快播 成人福利网午夜 一级黄色夫妻录像片 直接看的gav久久播放器 国产自拍400首页 sm老爹影院 谁知道隔壁老王网址在线 综合网 123西瓜影音 米奇丁香 人人澡人人漠大学生 色久悠 夜色视频你今天寂寞了吗? 菲菲影视城美国 被抄的影院 变态另类 欧美 成人 国产偷拍自拍在线小说 不用下载安装就能看的吃男人鸡巴视频 插屄视频 大贯杏里播放 wwwhhh50 233若菜奈央 伦理片天海翼秘密搜查官 大香蕉在线万色屋视频 那种漫画小说你懂的 祥仔电影合集一区 那里可以看澳门皇冠酒店a片 色自啪 亚洲aV电影天堂 谷露影院ar toupaizaixian sexbj。com 毕业生 zaixian mianfei 朝桐光视频 成人短视频在线直接观看 陈美霖 沈阳音乐学院 导航女 www26yjjcom 1大尺度视频 开平虐女视频 菅野雪松协和影视在线视频 华人play在线视频bbb 鸡吧操屄视频 多啪啪免费视频 悠草影院 金兰策划网 (969) 橘佑金短视频 国内一极刺激自拍片 日本制服番号大全magnet 成人动漫母系 电脑怎么清理内存 黄色福利1000 dy88午夜 偷拍中学生洗澡磁力链接 花椒相机福利美女视频 站长推荐磁力下载 mp4 三洞轮流插视频 玉兔miki热舞视频 夜生活小视频 爆乳人妖小视频 国内网红主播自拍福利迅雷下载 不用app的裸裸体美女操逼视频 变态SM影片在线观看 草溜影院元气吧 - 百度 - 百度 波推全套视频 国产双飞集合ftp 日本在线AV网 笔国毛片 神马影院女主播是我的邻居 影音资源 激情乱伦电影 799pao 亚洲第一色第一影院 av视频大香蕉 老梁故事汇希斯莱杰 水中人体磁力链接 下载 大香蕉黄片免费看 济南谭崔 避开屏蔽的岛a片 草破福利 要看大鸡巴操小骚逼的人的视频 黑丝少妇影音先锋 欧美巨乳熟女磁力链接 美国黄网站色大全 伦蕉在线久播 极品女厕沟 激情五月bd韩国电影 混血美女自摸和男友激情啪啪自拍诱人呻吟福利视频 人人摸人人妻做人人看 44kknn 娸娸原网 伊人欧美 恋夜影院视频列表安卓青青 57k影院 如果电话亭 avi 插爆骚女精品自拍 青青草在线免费视频1769TV 令人惹火的邻家美眉 影音先锋 真人妹子被捅动态图 男人女人做完爱视频15 表姐合租两人共处一室晚上她竟爬上了我的床 性爱教学视频 北条麻妃bd在线播放版 国产老师和师生 magnet wwwcctv1024 女神自慰 ftp 女同性恋做激情视频 欧美大胆露阴视频 欧美无码影视 好女色在线观看 后入肥臀18p 百度影视屏福利 厕所超碰视频 强奸mp magnet 欧美妹aⅴ免费线上看 2016年妞干网视频 5手机在线福利 超在线最视频 800av:cOm magnet 欧美性爱免播放器在线播放 91大款肥汤的性感美乳90后邻家美眉趴着窗台后入啪啪 秋霞日本毛片网站 cheng ren 在线视频 上原亚衣肛门无码解禁影音先锋 美脚家庭教师在线播放 尤酷伦理片 熟女性生活视频在线观看 欧美av在线播放喷潮 194avav 凤凰AV成人 - 百度 kbb9999 AV片AV在线AV无码 爱爱视频高清免费观看 黄色男女操b视频 观看 18AV清纯视频在线播放平台 成人性爱视频久久操 女性真人生殖系统双性人视频 下身插入b射精视频 明星潜规测视频 mp4 免賛a片直播绪 国内 自己 偷拍 在线 国内真实偷拍 手机在线 国产主播户外勾在线 三桥杏奈高清无码迅雷下载 2五福电影院凸凹频频 男主拿鱼打女主,高宝宝 色哥午夜影院 川村まや痴汉 草溜影院费全过程免费 淫小弟影院在线视频 laohantuiche 啪啪啪喷潮XXOO视频 青娱乐成人国产 蓝沢润 一本道 亚洲青涩中文欧美 神马影院线理论 米娅卡莉法的av 在线福利65535 欧美粉色在线 欧美性受群交视频1在线播放 极品喷奶熟妇在线播放 变态另类无码福利影院92 天津小姐被偷拍 磁力下载 台湾三级电髟全部 丝袜美腿偷拍自拍 偷拍女生性行为图 妻子的乱伦 白虎少妇 肏婶骚屄 外国大妈会阴照片 美少女操屄图片 妹妹自慰11p 操老熟女的b 361美女人体 360电影院樱桃 爱色妹妹亚洲色图 性交卖淫姿势高清图片一级 欧美一黑对二白 大色网无毛一线天 射小妹网站 寂寞穴 西西人体模特苍井空 操的大白逼吧 骚穴让我操 拉好友干女朋友3p