Thermal energy minimisation in papermaking via MD control approach

This paper reports the application of Advanced Process Control (APC) techniques for improving the thermal energy efficiency of a paperboard-making process by regulating the Machine Direction (MD) profile of the basis weight and moisture content of the paper-board. A Model Predictive Controller (MPC) is designed so that the sheet moisture and basis weight tracking errors along with variations of the sheet moisture and basis weight are reduced. Also, the drainage is maximised through improved wet-end stability which can facilitate driving the sheet moisture set-point closer to its upper specification limit over time. It is shown that the proposed strategy can result in reducing steam usage by 8-10%. A simulation study based on a UK board machine is presented to show the effectiveness of the proposed technique. © 2011 Intl Journal of Adv Mechatr.

Graphcut-based interactive segmentation using colour and depth cues

Segmentation of novel or dynamic objects in a scene, often referred to as background sub- traction or foreground segmentation, is critical for robust high level computer vision applica- tions such as object tracking, object classifca- tion and recognition. However, automatic real- time segmentation for robotics still poses chal- lenges including global illumination changes, shadows, inter-re ections, colour similarity of foreground to background, and cluttered back- grounds. This paper introduces depth cues provided by structure from motion (SFM) for interactive segmentation to alleviate some of these challenges. In this paper, two prevailing interactive segmentation algorithms are com- pared; Lazysnapping [Li et al., 2004] and Grab- cut [Rother et al., 2004], both based on graph- cut optimisation [Boykov and Jolly, 2001]. The algorithms are extended to include depth cues rather than colour only as in the original pa- pers. Results show interactive segmentation based on colour and depth cues enhances the performance of segmentation with a lower er- ror with respect to ground truth.

Joint 2D and 3D cues for image segmentation towards robotic applications

This thesis investigates the fusion of 3D visual information with 2D image cues to provide 3D semantic maps of large-scale environments in which a robot traverses for robotic applications. A major theme of this thesis was to exploit the availability of 3D information acquired from robot sensors to improve upon 2D object classification alone. The proposed methods have been evaluated on several indoor and outdoor datasets collected from mobile robotic platforms including a quadcopter and ground vehicle covering several kilometres of urban roads.

Computer simulation of protein—carbohydrate complexes: application to arabinose-binding protein and pea lectin

The CCEM method (Contact Criteria and Energy Minimisation) has been developed and applied to study protein-carbohydrate interactions. The method uses available X-ray data even on the native protein at low resolution (above 2.4 Å) to generate realistic models of a variety of proteins with various ligands.The two examples discussed in this paper are arabinose-binding protein (ABP) and pea lectin. The X-ray crystal structure data reported on ABP-β-l-arabinose complex at 2.8, 2.4 and 1.7 Å resolution differ drastically in predicting the nature of the interactions between the protein and ligand. It is shown that, using the data at 2.4 Å resolution, the CCEM method generates complexes which are as good as the higher (1.7 Å) resolution data. The CCEM method predicts some of the important hydrogen bonds between the ligand and the protein which are missing in the interpretation of the X-ray data at 2.4 Å resolution. The theoretically predicted hydrogen bonds are in good agreement with those reported at 1.7 Å resolution. Pea lectin has been solved only in the native form at 3 Å resolution. Application of the CCEM method also enables us to generate complexes of pea lectin with methyl-α-d-glucopyranoside and methyl-2,3-dimethyl-α-d-glucopyranoside which explain well the available experimental data in solution.

Computer modelling studies on the mechanism of action of ribonuclease T1

The mechanism of action of ribonuclease (RNase) T1 is still a matter of considerable debate as the results of x-ray, 2-D nmr and site-directed mutagenesis studies disagree regarding the role of the catalytically important residues. Hence computer modelling studies were carried out by energy minimisation of the complexes of RNase T1 and some of its mutants (His40Ala, His40Lys, and Glu58Ala) with the substrate guanyl cytosine (GpC), and of native RNase T1 with the reaction intermediate guanosine 2',3'-cyclic phosphate (G greater than p). The puckering of the guanosine ribose moiety in the minimum energy conformer of the RNase T1-GpC (substrate) complex was found to be O4'-endo and not C3'-endo as in the RNase T1-3'-guanylic acid (inhibitor/product) complex. A possible scheme for the mechanism of action of RNase T1 has been proposed on the basis of the arrangement of the catalytically important amino acid residues His40, Glu58, Arg77, and His92 around the guanosine ribose and the phosphate moiety in the RNase T1-GpC and RNase T1-G greater than p complexes. In this scheme, Glu58 serves as the general base group and His92 as the general acid group in the transphosphorylation step. His40 may be essential for stabilising the negatively charged phosphate moiety in the enzyme-transition state complex.

Molecular Mechanics Studies of Homopolymeric and Mixed Sequence Py.Pu*Pu Triple Helices,

DNA triple helices containing two purine strands and one pyrimidine strand (C.G*G and T.A*A) have been studied, using model building followed by energy minimisation, for different orientations of the third strand resulting from variation in the hydrogen bonding between the Watson-Crick duplex and the third strand and the glycosidic torsion angle in the third strand. Our results show that in the C.G*G case the structure with a parallel orientation of the third strand, resulting from Hoogsteen hydrogen bonds between the third strand and the Watson-Crick duplex, is energetically the most favourable while in the T.A*A case the antiparallel orientation of the third strand, resulting from reverse Hoogsteen hydrogen bonds, is energetically the most favourable. These studies when extended to the mixed sequence triplexes, in which the second strand is a mixture of G and A, correspondingly the third strand is a mixture of G and APT, show that though the parallel orientation is still energetically more favourable, the antiparallel orientation becomes energetically comparable with an increasing number of thymines in the third strand. Structurally, for the mixed triplexes containing G and T in the third strand, it is seen that the basepair non-isomorphism between the C.G*G and the T.A*T triplets can be overcome with some changes in the base pair parameters without much distortion of either the backbone or the hydrogen bonds.

Structural Polymorphism in d(T).d(A)*d(T) Triple Helices

DNA triple helices containing two thymine strands and one adenine strand have been studied, using model building followed by energy minimisation, for different orientations of the third strand resulting from variation in the hydrogen bonding between the Watson-Crick duplex and the third strand and the glycosidic torsion angle in the third strand. Our results show that the structure with a parallel orientation of the third strand, in which the third strand base forms Hoogsteen hydrogen bonds with the adenine base in the Watson-Crick duplex, is energetically the most favourable. An antiparallel orientation of the third strand is also possible, in which the third strand base hydrogen bonds to both the bases in the Watson-Crick duplex. This structure is energetically comparable to the parallel structure. For the parallel triplex a 200ps molecular dynamics simulation starting from two different starting structures indicates that at 300K significant structural heterogeneity exists in this tripler structure. The results are compared with existing structural data on this class of triplexes derived from theoretical and NMR techniques.

Structural characterisation of a uracil containing hairpin DNA by NMR and molecular dynamics

Three-dimensional (3D) structure of a hairpin DNA d-CTAGAGGATCCTTTUGGATCCT (22mer; abbreviated as U4-hairpin), which has a uracil nucleotide unit at the fourth position from the 5' end of the tetra-loop has been solved by NMR spectroscopy. The H-1 resonances of this hairpin have been assigned almost completely. NMR restrained molecular dynamics and energy minimisation procedures have been used to describe the 3D structure of the U4 hairpin. This study establishes that the stem of the hairpin adopts a right handed B-DNA conformation while the T-12 and U-15 nucleotide stack upon 3' and 5' ends of the stem, respectively. Further, T-14 stacks upon both T-12 and U-15 while T-13 partially stacks upon T-14. Very weak stacking interaction is observed between T-13 and T-12. All the individual nucleotide bases adopt 'anti' conformation with respect to their sugar moiety. The turning phosphate in the loop is located between T-13 and T-14. The stereochemistry of U-15 mimics the situation where uracil would stack in a B-DNA conformation. This could be the reason as to why the U4-hairpin is found to be the best substrate for its interaction with uracil DNA glycosylase (UDG) compared to the other substrates in which the uracil is at the first, second and third positions of the tetra-loop from its 5' end, as reported previously.

Statistical deconvolution of enthalpic energetic contributions to MHC-peptide binding affinity

Background: MHC Class I molecules present antigenic peptides to cytotoxic T cells, which forms an integral part of the adaptive immune response. Peptides are bound within a groove formed by the MHC heavy chain. Previous approaches to MHC Class I-peptide binding prediction have largely concentrated on the peptide anchor residues located at the P2 and C-terminus positions. Results: A large dataset comprising MHC-peptide structural complexes was created by remodelling pre-determined x-ray crystallographic structures. Static energetic analysis, following energy minimisation, was performed on the dataset in order to characterise interactions between bound peptides and the MHC Class I molecule, partitioning the interactions within the groove into van der Waals, electrostatic and total non-bonded energy contributions. Conclusion: The QSAR techniques of Genetic Function Approximation (GFA) and Genetic Partial Least Squares (G/PLS) algorithms were used to identify key interactions between the two molecules by comparing the calculated energy values with experimentally-determined BL50 data. Although the peptide termini binding interactions help ensure the stability of the MHC Class I-peptide complex, the central region of the peptide is also important in defining the specificity of the interaction. As thermodynamic studies indicate that peptide association and dissociation may be driven entropically, it may be necessary to incorporate entropic contributions into future calculations.

Modelling the structure of a polymer glass poly(methylmethacrylate) through neutron scattering experiments

Determination of the local structure of a polymer glass by scattering methods is complex due to the number of spatial and orientational correlations, both from within the polymer chain (intrachain) and between neighbouring chains (interchain), from which the scattering arises. Recently considerable advances have been made in the structural analysis of relatively simple polymers such as poly(ethylene) through the use of broad Q neutron scattering data tightly coupled to atomistic modelling procedures. This paper presents the results of an investigation into the use of these procedures for the analysis of the local structure of a-PMMA which is chemically more complex with a much greater number of intrachain structural parameters. We have utilised high quality neutron scattering data obtained using SANDALS at ISIS coupled with computer models representing both the single chain and bulk polymer system. Several different modelling approaches have been explored which encompass such techniques as Reverse Monte Carlo refinement and energy minimisation and their relative merits and successes are discussed. These different approaches highlight structural parameters which any realistic model of glassy atactic PMMA must replicate.

Fast stereo matching by iterated dynamic programming and quadtree subregioning

The application of energy minimisation methods for stereo matching has been demonstrated to produce high quality disparity maps. However the majority of these methods are known to be computationally expensive, requiring minutes or even hours of computation. We propose a fast minimisation scheme that produces strongly competitive results for significantly reduced computation, requiring only a few seconds of computation. In this paper, we present our iterated dynamic programming algorithm along with a quadtree subregioning process for fast stereo matching.

Studies on the chromosomal bets-lactamase of pseudomonas aeruginosa

The chromosomal ß-lactamase of Pseudomonas aeruginosa SAlconst (a derepressed laboratory strain) was isolated and purified. Two peaks of activity were observed on gel permeation chromatography (one major peak mol. wt. 45 kD and one minor peak of 54 kD). Preparations from 12 clinical derepressed strains showed identical results. Chromosomal ß-lactamase production in both normal and derepressed P. aeruginosa strains was induced both by iron restricted growth conditions and by penicillin G. The majority of the enzyme (80-90%) was found in the periplasm and cytoplasm but a significant amount (2-20%) was associated with the outer membrane (OM). The growth conditions did not affect the distribution of the enzyme between subcellular fractions although higher activity was found in the cells grown under iron limitation and/ or in the presence of ß-lactams. The penicillanate sulphone inhibitor, tazobactam, displayed irreversible kinetics whilst cloxacillin, cefotaxime, ampicillin and penicillin G were all competitive inhibitors of the enzyme. Similar results were obtained for the Enterobacter cloacae P99 [ß-lactamase, but tazobactam displayed a non-classical kinetic pattern for the Staphylococcus aureus PC1 ß-lactamase. The residues involved in ß-lactam hydrolysis by the P aeruginosa SAlconst enzyme were detennined by affinity labelling with tazobactam. A tryptic digestion fragment of the inhibited enzyme contained the amino acids D, T, S, E, P, G, A, C, V, M, I, Y, F, H, K, R. This suggests the involvement of the conserved SVSK, DAE and KTG motifs found in all penicillin sensitive proteins. A model of the 3-D structure of the active site of the P aeruginosa SAlconst chromosomal ß-!actamase was constructed from the published amino acid sequence of P aeruginosa chromosomal ß-lactamase and the a-carbon coordinates of the S. aureus PCI ß-lactamase by homology modelling and energy minimisation. The crystal structure of tazobactam was determined and energy minimised. Computer graphics docking identified Ser 72 as a possible residue involved in a secondary attack on the C5 position of tazobactam after initial ß-lactam hydrolysis by serine 70. The enhanced activity of tazobactam over sulbactam might be explained by the triazole substituent which might participate in favourable hydrogen bonding between N3 and active site residues.

Identification of the HLA-DM/HLA-DR interface

Human leukocyte antigen (HLA)-DM is a critical participant in antigen presentation that catalyzes the dissociation of the Class II-associated Invariant chain-derived Peptide (CLIP) from the major histocompatibility complex (MHC) Class II molecules. There is competition amongst peptides for access to an MHC Class II groove and it has been hypothesised that DM functions as a 'peptide editor' that catalyzes the replacement of one peptide for another within the groove. It is established that the DM catalyst interacts directly with the MHC Class II but the precise location of the interface is unknown. Here, we combine previously described mutational data with molecular docking and energy minimisation simulations to identify a putative interaction site of >4000A2 which agrees with known point mutational data for both the DR and DM molecule. The docked structure is validated by comparison with experimental data and previously determined properties of protein-protein interfaces. A possible dissociation mechanism is suggested by the presence of an acidic cluster near the N terminus of the bound peptide.

Statistical deconvolution of enthalpic energetic contributions to MHC-peptide binding affinity

Background - MHC Class I molecules present antigenic peptides to cytotoxic T cells, which forms an integral part of the adaptive immune response. Peptides are bound within a groove formed by the MHC heavy chain. Previous approaches to MHC Class I-peptide binding prediction have largely concentrated on the peptide anchor residues located at the P2 and C-terminus positions. Results - A large dataset comprising MHC-peptide structural complexes was created by re-modelling pre-determined x-ray crystallographic structures. Static energetic analysis, following energy minimisation, was performed on the dataset in order to characterise interactions between bound peptides and the MHC Class I molecule, partitioning the interactions within the groove into van der Waals, electrostatic and total non-bonded energy contributions. Conclusion - The QSAR techniques of Genetic Function Approximation (GFA) and Genetic Partial Least Squares (G/PLS) algorithms were used to identify key interactions between the two molecules by comparing the calculated energy values with experimentally-determined BL50 data. Although the peptide termini binding interactions help ensure the stability of the MHC Class I-peptide complex, the central region of the peptide is also important in defining the specificity of the interaction. As thermodynamic studies indicate that peptide association and dissociation may be driven entropically, it may be necessary to incorporate entropic contributions into future calculations.

Pyrolysis of municipal green waste: a modelling, simulation and experimental analysis

Pyrolysis is the thermo-chemical conversion of carbonaceous feedstock in the absence of oxygen to produce bio-fuel (bio-oil, bio-char and syn-gas). Bio-fuel production from municipal green waste (MGW) through the pyrolysis process has attracted considerable attention recently in the renewable energy sector because it can reduce greenhouse gas emissions and contribute to energy security. This study analyses properties of MGW feedstock available in Rockhampton city of Central Queensland, Australia, and presents an experimental investigation of producing bio-fuel from that MGW through the pyrolysis process using a short sealed rotary furnace. It was found from the experiment that about 19.97% bio-oil, 40.83% bio-char and 29.77% syn-gas can be produced from the MGW. Then, a four-stage steady state simulation model is developed for pyrolysis process performance simulation using Aspen Plus software. In the first stage, the moisture content of the MGW feed is reduced. In the second stage, the MGW is decomposed according to its elemental constituents. In the third stage, condensate material is separated and, finally, the pyrolysis reactions are modelled using the Gibbs free energy minimisation approach. The MGW's ultimate and proximate analysis data were used in the Aspen Plus simulation as input parameters. The model is validated with experimentally measured data. A good agreement between simulation and experimental results was found. More specifically, the variation of modelling and experimental elemental compositions of the MGW was found to be 7.3% for carbon, 15.82% for hydrogen, 7.04% for nitrogen and 5.56% for sulphur. The validated model is used to optimise the biofuel production from the MGW as a function of operating variables such as temperature, moisture content, particle size and process heat air-fuel ratio. The modelling and optimisation results are presented, analysed and discussed.