1
|
Brück HL, Coutte F, Dhulster P, Gofflot S, Jacques P, Delvigne F. Growth Dynamics of Bacterial Populations in a Two-Compartment Biofilm Bioreactor Designed for Continuous Surfactin Biosynthesis. Microorganisms 2020; 8:microorganisms8050679. [PMID: 32392736 PMCID: PMC7285194 DOI: 10.3390/microorganisms8050679] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2020] [Revised: 04/28/2020] [Accepted: 05/01/2020] [Indexed: 11/16/2022] Open
Abstract
Biofilm bioreactors are promising systems for continuous biosurfactant production since they provide process stability through cell immobilization and avoid foam formation. In this work, a two-compartment biofilm bioreactor was designed consisting of a stirred tank reactor and a trickle-bed reactor containing a structured metal packing for biofilm formation. A strong and poor biofilm forming B. subtilis 168 strain due to restored exopolysaccharides (EPS) production or not were cultivated in the system to study the growth behavior of the planktonic and biofilm population for the establishment of a growth model. A high dilution rate was used in order to promote biofilm formation on the packing and wash out unwanted planktonic cells. Biofilm development kinetics on the packing were assessed through a total organic carbon mass balance. The EPS+ strain showed a significantly improved performance in terms of adhesion capacity and surfactin production. The mean surfactin productivity of the EPS+ strain was about 37% higher during the continuous cultivation compared to the EPS- strain. The substrate consumption together with the planktonic cell and biofilm development were properly predicted by the model (α = 0.05). The results show the efficiency of the biofilm bioreactor for continuous surfactin production using an EPS producing strain.
Collapse
Affiliation(s)
- Hannah Luise Brück
- MiPI, TERRA Teaching and Research Centre, Joint Research Unit BioEcoAgro N° 1158, Gembloux Agro-Bio Tech, University Liège, University Lille, INRAE, UPJV, YNCREA, University Artois, University Littoral Côte d’Opale, B-5030 Gembloux, Belgium; (H.L.B.); (P.J.)
- ICV—Institut Charles Viollette, Joint Research Unit BioEcoAgro N° 1158, University Lille, INRAE, University Liège, UPJV, YNCREA, University Artois, University Littoral Côte d’Opale, F-59000 Lille, France; (F.C.); (P.D.)
| | - François Coutte
- ICV—Institut Charles Viollette, Joint Research Unit BioEcoAgro N° 1158, University Lille, INRAE, University Liège, UPJV, YNCREA, University Artois, University Littoral Côte d’Opale, F-59000 Lille, France; (F.C.); (P.D.)
| | - Pascal Dhulster
- ICV—Institut Charles Viollette, Joint Research Unit BioEcoAgro N° 1158, University Lille, INRAE, University Liège, UPJV, YNCREA, University Artois, University Littoral Côte d’Opale, F-59000 Lille, France; (F.C.); (P.D.)
| | - Sébastien Gofflot
- Walloon Agricultural Research Center (CRA-W), Agricultural Product Technology Unit, Chaussée de Namur, 24, B-5030 Gembloux, Belgium;
| | - Philippe Jacques
- MiPI, TERRA Teaching and Research Centre, Joint Research Unit BioEcoAgro N° 1158, Gembloux Agro-Bio Tech, University Liège, University Lille, INRAE, UPJV, YNCREA, University Artois, University Littoral Côte d’Opale, B-5030 Gembloux, Belgium; (H.L.B.); (P.J.)
| | - Frank Delvigne
- MiPI, TERRA Teaching and Research Centre, Joint Research Unit BioEcoAgro N° 1158, Gembloux Agro-Bio Tech, University Liège, University Lille, INRAE, UPJV, YNCREA, University Artois, University Littoral Côte d’Opale, B-5030 Gembloux, Belgium; (H.L.B.); (P.J.)
- Correspondence:
| |
Collapse
|
2
|
Saadat NP, Nies T, Rousset Y, Ebenhöh O. Thermodynamic Limits and Optimality of Microbial Growth. ENTROPY 2020; 22:e22030277. [PMID: 33286054 PMCID: PMC7516730 DOI: 10.3390/e22030277] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/06/2020] [Accepted: 02/21/2020] [Indexed: 12/20/2022]
Abstract
Understanding microbial growth with the use of mathematical models has a long history that dates back to the pioneering work of Jacques Monod in the 1940s. Monod’s famous growth law expressed microbial growth rate as a simple function of the limiting nutrient concentration. However, to explain growth laws from underlying principles is extremely challenging. In the second half of the 20th century, numerous experimental approaches aimed at precisely measuring heat production during microbial growth to determine the entropy balance in a growing cell and to quantify the exported entropy. This has led to the development of thermodynamic theories of microbial growth, which have generated fundamental understanding and identified the principal limitations of the growth process. Although these approaches ignored metabolic details and instead considered microbial metabolism as a black box, modern theories heavily rely on genomic resources to describe and model metabolism in great detail to explain microbial growth. Interestingly, however, thermodynamic constraints are often included in modern modeling approaches only in a rather superficial fashion, and it appears that recent modeling approaches and classical theories are rather disconnected fields. To stimulate a closer interaction between these fields, we here review various theoretical approaches that aim at describing microbial growth based on thermodynamics and outline the resulting thermodynamic limits and optimality principles. We start with classical black box models of cellular growth, and continue with recent metabolic modeling approaches that include thermodynamics, before we place these models in the context of fundamental considerations based on non-equilibrium statistical mechanics. We conclude by identifying conceptual overlaps between the fields and suggest how the various types of theories and models can be integrated. We outline how concepts from one approach may help to inform or constrain another, and we demonstrate how genome-scale models can be used to infer key black box parameters, such as the energy of formation or the degree of reduction of biomass. Such integration will allow understanding to what extent microbes can be viewed as thermodynamic machines, and how close they operate to theoretical optima.
Collapse
Affiliation(s)
- Nima P. Saadat
- Institute of Quantitative and Theoretical Biology, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225 Düsseldorf, Germany; (N.P.S.); (T.N.); (Y.R.)
| | - Tim Nies
- Institute of Quantitative and Theoretical Biology, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225 Düsseldorf, Germany; (N.P.S.); (T.N.); (Y.R.)
| | - Yvan Rousset
- Institute of Quantitative and Theoretical Biology, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225 Düsseldorf, Germany; (N.P.S.); (T.N.); (Y.R.)
| | - Oliver Ebenhöh
- Institute of Quantitative and Theoretical Biology, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225 Düsseldorf, Germany; (N.P.S.); (T.N.); (Y.R.)
- Cluster of Excellence on Plant Sciences (CEPLAS), Heinrich-Heine University, Universitätsstrasse 1, 40225 Düsseldorf, Germany
- Correspondence:
| |
Collapse
|
3
|
Santos SB, Carvalho C, Azeredo J, Ferreira EC. Population dynamics of a Salmonella lytic phage and its host: implications of the host bacterial growth rate in modelling. PLoS One 2014; 9:e102507. [PMID: 25051248 PMCID: PMC4106826 DOI: 10.1371/journal.pone.0102507] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2013] [Accepted: 06/20/2014] [Indexed: 11/18/2022] Open
Abstract
The prevalence and impact of bacteriophages in the ecology of bacterial communities coupled with their ability to control pathogens turn essential to understand and predict the dynamics between phage and bacteria populations. To achieve this knowledge it is essential to develop mathematical models able to explain and simulate the population dynamics of phage and bacteria. We have developed an unstructured mathematical model using delay-differential equations to predict the interactions between a broad-host-range Salmonella phage and its pathogenic host. The model takes into consideration the main biological parameters that rule phage-bacteria interactions likewise the adsorption rate, latent period, burst size, bacterial growth rate, and substrate uptake rate, among others. The experimental validation of the model was performed with data from phage-interaction studies in a 5 L bioreactor. The key and innovative aspect of the model was the introduction of variations in the latent period and adsorption rate values that are considered as constants in previous developed models. By modelling the latent period as a normal distribution of values and the adsorption rate as a function of the bacterial growth rate it was possible to accurately predict the behaviour of the phage-bacteria population. The model was shown to predict simulated data with a good agreement with the experimental observations and explains how a lytic phage and its host bacteria are able to coexist.
Collapse
Affiliation(s)
- Sílvio B. Santos
- Centre of Biological Engineering, Universidade do Minho, Braga, Portugal
| | - Carla Carvalho
- Centre of Biological Engineering, Universidade do Minho, Braga, Portugal
| | - Joana Azeredo
- Centre of Biological Engineering, Universidade do Minho, Braga, Portugal
| | | |
Collapse
|
4
|
Sansonetti S, Hobley TJ, Calabrò V, Villadsen J, Sin G. A biochemically structured model for ethanol fermentation by Kluyveromyces marxianus: A batch fermentation and kinetic study. BIORESOURCE TECHNOLOGY 2011; 102:7513-7520. [PMID: 21632239 DOI: 10.1016/j.biortech.2011.05.014] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/10/2011] [Revised: 05/05/2011] [Accepted: 05/06/2011] [Indexed: 05/30/2023]
Abstract
Anaerobic batch fermentations of ricotta cheese whey (i.e. containing lactose) were performed under different operating conditions. Ethanol concentrations of ca. 22g L(-1) were found from whey containing ca. 44g L(-1) lactose, which corresponded to up to 95% of the theoretical ethanol yield within 15h. The experimental data could be explained by means of a simple knowledge-driven biochemically structured model that was built on bioenergetics principles applied to the metabolic pathways through which lactose is converted into major products. Use of the model showed that the observed concentrations of ethanol, lactose, biomass and glycerol during batch fermentation could be described within a ca. 6% deviation, as could the yield coefficients for biomass and ethanol produced on lactose. The model structure confirmed that the thermodynamics considerations on the stoichiometry of the system constrain the metabolic coefficients within a physically meaningful range thereby providing valuable and reliable insight into fermentation processes.
Collapse
Affiliation(s)
- S Sansonetti
- Computer Aided Process Engineering Center, Department of Chemical and Biochemical Engineering, Technical University of Denmark, DK-2800 Lyngby, Denmark.
| | | | | | | | | |
Collapse
|
5
|
Nielsen J, Nikolajsen K, Villadsen J. Structured modeling of a microbial system: II. Experimental verification of a structured lactic acid fermentation model. Biotechnol Bioeng 2010; 38:11-23. [PMID: 18600693 DOI: 10.1002/bit.260380103] [Citation(s) in RCA: 52] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
A two-compartment model for the lactic acid fermentation with Streptococcus cremoris is experimentally verified. The seven parameters of the model are determined using steady-state chemostat data at varying values of dilution rate, D, but with a constant feed concentration, s(f), of a single carbohydrate source (glucose, lactose, or galactose), and a constant feed concentration of s(Nf) of the N source. Steady-state measurements of the RNA content at different exit concentrations, s, of the carbohydrate are included to calculate kinetic parameters that determine the cell composition for varying operating conditions. The model is tested using data from a large set of steady-state and non-steady-state experiments: batch fermentations and step and pulse experiments in a chemostat. Both qualitatively and quantitatively the major features of the model are confirmed: the external substrates enter into intracellular high-energy building blocks, and lactic acid is formed as a by-product of these reactions. Cell growth depends on the fraction of active components (X(A)) of the cell and is not accompanied by lactic acid production. Possible model modifications are discussed, primarily to obtain a better description of lactic acid fermentation at nongrowth conditions.
Collapse
Affiliation(s)
- J Nielsen
- Department of Biotechnology, Technical University of Denmark, Lyngby, Denmark
| | | | | |
Collapse
|
6
|
Schügerl K. Development of bioreaction engineering. ADVANCES IN BIOCHEMICAL ENGINEERING/BIOTECHNOLOGY 2001; 70:41-76. [PMID: 11092128 DOI: 10.1007/3-540-44965-5_3] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
Abstract
In addition to summarizing the early investigations in bioreaction engineering, the present short review covers the development of the field in the last 50 years. A brief overview of the progress of the fundamentals is presented in the first part of this article and the key issues of bioreaction engineering are advanced in its second part.
Collapse
Affiliation(s)
- K Schügerl
- Institute for Technical Chemistry, University of Hannover, Germany.
| |
Collapse
|
7
|
Garcı́a-Ochoa F, Santos V, Alcón A. Metabolic structured kinetic model for xanthan production. Enzyme Microb Technol 1998. [DOI: 10.1016/s0141-0229(98)00014-3] [Citation(s) in RCA: 35] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
8
|
Gikas P, Livingston A. Use of Specific ATP Concentration and Specific Oxygen Uptake Rate to Determine Parameters of a Structured Model of Biomass Growth. Enzyme Microb Technol 1998. [DOI: 10.1016/s0141-0229(97)00242-1] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
9
|
García-Ochoa F, Santos VE, Alcón A. Intracellular compounds quantification by means of flow cytometry in bacteria: application to xanthan production by Xanthomonas campestris. Biotechnol Bioeng 1998; 57:87-94. [PMID: 10099182 DOI: 10.1002/(sici)1097-0290(19980105)57:1<87::aid-bit11>3.0.co;2-5] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The use of flow cytometry (FCM) to quantitatively analyze intracellular compounds is studied. FCM is a very useful technique for individual cell studies in microbial systems, and gives access to information which cannot be obtained in any other way. Nevertheless, it provides data in arbitrary units, that is, relative data. This analytical technique could be employed for kinetic modeling of microbial systems and even for internal phenomena analysis, but for this purpose, absolute data-that is concentration of intracellular compounds-must be used. In this work, relative flow cytometry data are transformed into absolute data by means of calibrations employing the same fluorochromes with another technique: spectrofluorymetry. Calibrations of DNA, RNA, and protein intracellular concentrations are presented for the bacteria, Xanthomonas campestris. Other analytical methods, based on biochemical determinations, were also employed to quantify intracellular compounds, but the results obtained are very poor compared with those achieved by means of spectrofluorymetry (SFM). Calibration equations and data obtained by both techniques are given. Evolutions of protein and nucleic acids during Xanthomonas campestris growth and xanthan gum production are shown.
Collapse
Affiliation(s)
- F García-Ochoa
- Departamento de Ingeniería Química, Facultad de Ciencias Químicas, Universidad Complutense, 28040 Madrid, Spain.
| | | | | |
Collapse
|
10
|
Gikas P, Livingston AG. Specific ATP and specific oxygen uptake rate in immobilized cell aggregates: Experimental results and theoretical analysis using a structured model of immobilized cell growth. Biotechnol Bioeng 1997; 55:660-73. [DOI: 10.1002/(sici)1097-0290(19970820)55:4<660::aid-bit8>3.0.co;2-f] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
11
|
Periodic and nonperiodic oscillatory behavior in a model for activated sludge reactors. ACTA ACUST UNITED AC 1997. [DOI: 10.1016/s0895-7177(97)00071-x] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
12
|
Spragg AJP, McNeil B, Harvey LM. Development of a structured model for biopolymer synthesis in the fungus Aureobasidium pullulans. Bioprocess Biosyst Eng 1995. [DOI: 10.1007/bf01112990] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
13
|
Bader J, Klingspohn U, Bellgardt KH, Schügerl K. Modelling and simulation of the growth and enzyme production of Trichoderma reesei Rut C30. J Biotechnol 1993. [DOI: 10.1016/0168-1656(93)90045-o] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
14
|
Dapaah SY, Hill GA. Biodegradation of chlorophenol mixtures byPseudomonas putida. Biotechnol Bioeng 1992; 40:1353-8. [DOI: 10.1002/bit.260401109] [Citation(s) in RCA: 49] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
15
|
|
16
|
|
17
|
Kooijman SA, Muller EB, Stouthamer AH. Microbial growth dynamics on the basis of individual budgets. Antonie Van Leeuwenhoek 1991; 60:159-74. [PMID: 1807193 DOI: 10.1007/bf00430363] [Citation(s) in RCA: 21] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The popular theories for microbial dynamics by Monod, Pirt and Droop are shown to be special cases of a model for individual budgets, in which growth and maintenance are on the expense of reserve materials. The dynamics of reserve materials is a first order process with a relaxation time proportional to cell length; maintenance is proportional to cell volume, and uptake, which depends hyperbolically on substrate density, is proportional to cell volume as well. Because of the latter, population dynamics depends on the behaviour of the individuals in a simple way, such that the cell volume distribution has no quantitative effect. When uptake is proportional to the surface area of the cell, which is realistic from a physical point of view, the relation between the individual level and the population one becomes more complicated and the cell size and shape distribution affects population dynamics. It is shown how the changing shape of rods modifies uptake and, consequently, growth. The concept of energy conductance, defined as the ratio of the maximum surface area specific uptake and the volume specific energy reserve has been introduced in the analysis of microbial dynamics. The first tentative results indicate that the value for E. coli is close to the mean value for a wide variety of animals. Properties of the model for cell suspension at constant substrate densities are analyzed and tested against a variety of experimental data from the literature on both the individual and the population level.
Collapse
Affiliation(s)
- S A Kooijman
- Biological Laboratory, Free University, Amsterdam, The Netherlands
| | | | | |
Collapse
|
18
|
Nielsen J, Nikolajsen K, Villadsen J. Structured modeling of a microbial system: I. A theoretical study of lactic acid fermentation. Biotechnol Bioeng 1991; 38:1-10. [DOI: 10.1002/bit.260380102] [Citation(s) in RCA: 67] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
19
|
|
20
|
Guterman H, Vonshak A, Ben-Yaakov S. A macromodel for outdoor algal mass production. Biotechnol Bioeng 1990; 35:809-19. [DOI: 10.1002/bit.260350809] [Citation(s) in RCA: 36] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
21
|
Allsop PJ, Moo‐Young M, Sullivan GR. The dynamics and control of substrate inhibition in activated sludge∗. ACTA ACUST UNITED AC 1990. [DOI: 10.1080/10643389009388393] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
22
|
Lievense JC, Modak JM, Lim HC. A mathematical model for the asymmetric dynamic response in microbial systems. J Biotechnol 1989. [DOI: 10.1016/0168-1656(89)90028-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
23
|
Turner BG, Ramkrishna D, Jansen NB. Cybernetic modeling of bacteriol cultures at low growth rates: Single-substrate systems. Biotechnol Bioeng 1989; 34:252-61. [DOI: 10.1002/bit.260340214] [Citation(s) in RCA: 44] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
24
|
Worden RM, Donaldson TL. Dynamics of a biological fixed film for phenol degradation in a fluidized-bed bioreactor. Biotechnol Bioeng 1987; 30:398-412. [DOI: 10.1002/bit.260300311] [Citation(s) in RCA: 66] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
25
|
Abstract
An intrinsic, structured model has been formulated to describe the kinetics of viable (living) cells immobilized within porous supports. Predictions of steady-state internal biomass concentration distributions, biocatalyst substrate profiles, and immobilized cell growth and leakage from the support are in qualitative agreement with the literature. Simulation studies indicate that carrier pore structure is a particularly important design variable to be optimized.
Collapse
|
26
|
Bryers JD. Structured modeling of the anaerobic digestion of biomass particulates. Biotechnol Bioeng 1985; 27:638-49. [DOI: 10.1002/bit.260270514] [Citation(s) in RCA: 69] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
27
|
Esener AA, Roels JA, Kossen NWF. Theory and applications of unstructured growth models: Kinetic and energetic aspects. Biotechnol Bioeng 1983; 25:2803-41. [DOI: 10.1002/bit.260251202] [Citation(s) in RCA: 82] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
28
|
Blanch HW. Cell growth and enzyme kinetics. Biotechnol Adv 1983; 1:193-204. [PMID: 14540891 DOI: 10.1016/0734-9750(83)90588-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
Affiliation(s)
- H W Blanch
- Department of Chemical Engineering, University of California, Berkeley 94720, USA
| |
Collapse
|