Rational design of serine protease inhibitors

Proteases regulate a spectrum of diverse physiological processes, and dysregulation of proteolytic activity drives a plethora of pathological conditions. Understanding protease function is essential to appreciating many aspects of normal physiology and progression of disease. Consequently, development of potent and specific inhibitors of proteolytic enzymes is vital to provide tools for the dissection of protease function in biological systems and for the treatment of diseases linked to aberrant proteolytic activity. The studies in this thesis describe the rational design of potent inhibitors of three proteases that are implicated in disease development. Additionally, key features of the interaction of proteases and their cognate inhibitors or substrates are analysed and a series of rational inhibitor design principles are expounded and tested. Rational design of protease inhibitors relies on a comprehensive understanding of protease structure and biochemistry. Analysis of known protease cleavage sites in proteins and peptides is a commonly used source of such information. However, model peptide substrate and protein sequences have widely differing levels of backbone constraint and hence can adopt highly divergent structures when binding to a protease’s active site. This may result in identical sequences in peptides and proteins having different conformations and diverse spatial distribution of amino acid functionalities. Regardless of this, protein and peptide cleavage sites are often regarded as being equivalent. One of the key findings in the following studies is a definitive demonstration of the lack of equivalence between these two classes of substrate and invalidation of the common practice of using the sequences of model peptide substrates to predict cleavage of proteins in vivo. Another important feature for protease substrate recognition is subsite cooperativity. This type of cooperativity is commonly referred to as protease or substrate binding subsite cooperativity and is distinct from allosteric cooperativity, where binding of a molecule distant from the protease active site affects the binding affinity of a substrate. Subsite cooperativity may be intramolecular where neighbouring residues in substrates are interacting, affecting the scissile bond’s susceptibility to protease cleavage. Subsite cooperativity can also be intermolecular where a particular residue’s contribution to binding affinity changes depending on the identity of neighbouring amino acids. Although numerous studies have identified subsite cooperativity effects, these findings are frequently ignored in investigations probing subsite selectivity by screening against diverse combinatorial libraries of peptides (positional scanning synthetic combinatorial library; PS-SCL). This strategy for determining cleavage specificity relies on the averaged rates of hydrolysis for an uncharacterised ensemble of peptide sequences, as opposed to the defined rate of hydrolysis of a known specific substrate. Further, since PS-SCL screens probe the preference of the various protease subsites independently, this method is inherently unable to detect subsite cooperativity. However, mean hydrolysis rates from PS-SCL screens are often interpreted as being comparable to those produced by single peptide cleavages. Before this study no large systematic evaluation had been made to determine the level of correlation between protease selectivity as predicted by screening against a library of combinatorial peptides and cleavage of individual peptides. This subject is specifically explored in the studies described here. In order to establish whether PS-SCL screens could accurately determine the substrate preferences of proteases, a systematic comparison of data from PS-SCLs with libraries containing individually synthesised peptides (sparse matrix library; SML) was carried out. These SML libraries were designed to include all possible sequence combinations of the residues that were suggested to be preferred by a protease using the PS-SCL method. SML screening against the three serine proteases kallikrein 4 (KLK4), kallikrein 14 (KLK14) and plasmin revealed highly preferred peptide substrates that could not have been deduced by PS-SCL screening alone. Comparing protease subsite preference profiles from screens of the two types of peptide libraries showed that the most preferred substrates were not detected by PS SCL screening as a consequence of intermolecular cooperativity being negated by the very nature of PS SCL screening. Sequences that are highly favoured as result of intermolecular cooperativity achieve optimal protease subsite occupancy, and thereby interact with very specific determinants of the protease. Identifying these substrate sequences is important since they may be used to produce potent and selective inhibitors of protolytic enzymes. This study found that highly favoured substrate sequences that relied on intermolecular cooperativity allowed for the production of potent inhibitors of KLK4, KLK14 and plasmin. Peptide aldehydes based on preferred plasmin sequences produced high affinity transition state analogue inhibitors for this protease. The most potent of these maintained specificity over plasma kallikrein (known to have a very similar substrate preference to plasmin). Furthermore, the efficiency of this inhibitor in blocking fibrinolysis in vitro was comparable to aprotinin, which previously saw clinical use to reduce perioperative bleeding. One substrate sequence particularly favoured by KLK4 was substituted into the 14 amino acid, circular sunflower trypsin inhibitor (SFTI). This resulted in a highly potent and selective inhibitor (SFTI-FCQR) which attenuated protease activated receptor signalling by KLK4 in vitro. Moreover, SFTI-FCQR and paclitaxel synergistically reduced growth of ovarian cancer cells in vitro, making this inhibitor a lead compound for further therapeutic development. Similar incorporation of a preferred KLK14 amino acid sequence into the SFTI scaffold produced a potent inhibitor for this protease. However, the conformationally constrained SFTI backbone enforced a different intramolecular cooperativity, which masked a KLK14 specific determinant. As a consequence, the level of selectivity achievable was lower than that found for the KLK4 inhibitor. Standard mechanism inhibitors such as SFTI rely on a stable acyl-enzyme intermediate for high affinity binding. This is achieved by a conformationally constrained canonical binding loop that allows for reformation of the scissile peptide bond after cleavage. Amino acid substitutions within the inhibitor to target a particular protease may compromise structural determinants that support the rigidity of the binding loop and thereby prevent the engineered inhibitor reaching its full potential. An in silico analysis was carried out to examine the potential for further improvements to the potency and selectivity of the SFTI-based KLK4 and KLK14 inhibitors. Molecular dynamics simulations suggested that the substitutions within SFTI required to target KLK4 and KLK14 had compromised the intramolecular hydrogen bond network of the inhibitor and caused a concomitant loss of binding loop stability. Furthermore in silico amino acid substitution revealed a consistent correlation between a higher frequency of formation and the number of internal hydrogen bonds of SFTI-variants and lower inhibition constants. These predictions allowed for the production of second generation inhibitors with enhanced binding affinity toward both targets and highlight the importance of considering intramolecular cooperativity effects when engineering proteins or circular peptides to target proteases. The findings from this study show that although PS-SCLs are a useful tool for high throughput screening of approximate protease preference, later refinement by SML screening is needed to reveal optimal subsite occupancy due to cooperativity in substrate recognition. This investigation has also demonstrated the importance of maintaining structural determinants of backbone constraint and conformation when engineering standard mechanism inhibitors for new targets. Combined these results show that backbone conformation and amino acid cooperativity have more prominent roles than previously appreciated in determining substrate/inhibitor specificity and binding affinity. The three key inhibitors designed during this investigation are now being developed as lead compounds for cancer chemotherapy, control of fibrinolysis and cosmeceutical applications. These compounds form the basis of a portfolio of intellectual property which will be further developed in the coming years.

Polymer self assembly on carbon nanotubes

This chapter analyses the poly(3-hexylthiophene) self-assembly on carbon nanotubes and the interaction between the two materials forming a new hybrid nanostructure. The chapter starts with a review of the several studies investigating polymers and biomolecules self-assembled on nanotubes. Then conducting polymers and polythiophenes are briefly introduced. Accordingly, carbon nanotube structure and properties are reported in Sect. 3. The experimental section starts with the bulk characterisation of polymer thin films with the inclusion of uniformly distributed carbon nanotubes. By using volume film analysis techniques (AFM, TEM, UV–Vis and Raman), we show how the polymer’s higher degree of order is a direct consequence of interaction with carbon nanotubes. Nevertheless, it is through the use of nanoscale analysis and molecular dynamic simulations that the self-assembly of the polymer on the nanotube surface can be clearly evidenced and characterised. In Sect. 6, the effect of the carbon templating structure on the P3HT organisation on the surface is investigated, showing the chirality-driven polymer assembly on the carbon nanotube surface. The interaction between P3HT and CNTs brings also to charge transfer, with the modification of physical properties for both species. In particular, the alteration of the polymer electronic properties and the modification of the nanotube mechanical structure are a direct consequence of the P3HT p-p stacking on the nanotube surface. Finally, some considerations based on molecular dynamics studies are reported in order to confirm and support the experimental results discussed.

Modified beam theories for bending properties of nanowires considering surface/intrinsic effects and axial extension effect

Several studies of the surface effect on bending properties of a nanowire (NW) have been conducted. However, these analyses are mainly based on theoretical predictions, and there is seldom integration study in combination between theoretical predictions and simulation results. Thus, based on the molecular dynamics (MD) simulation and different modified beam theories, a comprehensive theoretical and numerical study for bending properties of nanowires considering surface/intrinsic stress effects and axial extension effect is conducted in this work. The discussion begins from the Euler-Bernoulli beam theory and Timoshenko beam theory augmented with surface effect. It is found that when the NW possesses a relatively small cross-sectional size, these two theories cannot accurately interpret the true surface effect. The incorporation of axial extension effect into Euler-Bernoulli beam theory provides a nonlinear solution that agrees with the nonlinear-elastic experimental and MD results. However, it is still found inaccurate when the NW cross-sectional size is relatively small. Such inaccuracy is also observed for the Euler-Bernoulli beam theory augmented with both contributions from surface effect and axial extension effect. A comprehensive model for completely considering influences from surface stress, intrinsic stress, and axial extension is then proposed, which leads to good agreement with MD simulation results. It is thus concluded that, for NWs with a relatively small cross-sectional size, a simple consideration of surface stress effect is inappropriate, and a comprehensive consideration of the intrinsic stress effect is required.

MD investigations for mechanical properties of copper nanowires with and without surface defects

Based on the molecular dynamics (MD) method, the single-crystalline copper nanowire with different surface defects is investigated through tension simulation. For comparison, the MD tension simulations of perfect nanowire are firstly carried out under different temperatures, strain rates, and sizes. It has concluded that the surface-volume ratio significantly affects the mechanical properties of nanowire. The surface defects on nanowires are then systematically studied in considering different defect orientation and distribution. It is found that the Young’s modulus is insensitive of surface defects. However, the yield strength and yield point show a significant decrease due to the different defects. Different defects are observed to serve as a dislocation source.

A fundamental numerical and theoretical study for the vibrational properties of nanowires

Based on the molecular dynamics (MD) simulation and the classical Euler-Bernoulli beam theory, a fundamental study of the vibrational performance of the Ag nanowire (NW) is carried out. A comprehensive analysis of the quality (Q)-factor, natural frequency, beat vibration, as well as high vibration mode is presented. Two excitation approaches, i.e., velocity excitation and displacement excitation, have been successfully implemented to achieve the vibration of NWs. Upon these two kinds of excitations, consistent results are obtained, i.e., the increase of the initial excitation amplitude will lead to a decrease to the Q-factor, and moderate plastic deformation could increase the first natural frequency. Meanwhile, the beat vibration driven by a single relatively large excitation or two uniform excitations in both two lateral directions is observed. It is concluded that the nonlinear changing trend of external energy magnitude does not necessarily mean a nonconstant Q-factor. In particular, the first order natural frequency of the Ag NW is observed to decrease with the increase of temperature. Furthermore, comparing with the predictions by Euler- Bernoulli beam theory, the MD simulation provides a larger and smaller first vibration frequencies for the clamped-clamped and clamped-free thin Ag NWs, respectively. Additionally, for thin NWs, the first order natural frequency exhibits a parabolic relationship with the excitation magnitudes. The frequencies of the higher vibration modes tend to be low in comparison to Euler-Bernoulli beam theory predictions. A combined initial excitation is proposed which is capable to drive the NW under a multi-mode vibration and arrows the coexistence of all the following low vibration modes. This work sheds lights on the better understanding of the mechanical properties of NWs and benefits the increasing utilities of NWs in diverse nano-electronic devices.

Dynamic covalent chemistry on surfaces employing highly reactive cyclopentadienyl moieties

Silicon substrates coated with a bromide-terminated silane are transformed into highly reactive, cyclopentadiene covered analogues. These surfaces undergo rapid cycloaddition reactions with various dienophile-capped polymers. Mild heating of the substrates causes the retro-Diels-Alder reaction to occur, thus reforming the reactive cyclopentadiene surface, generating an efficiently switchable surface.

Numerical investigation of mechanical properties of nanowires : a review

Nanowires (NWs) have attracted intensive researches owing to the broad applications that arise from their remarkable properties. Over the last decade, immense numerical studies have been conducted for the numerical investigation of mechanical properties of NWs. Among these numerical simulations, the molecular dynamics (MD) plays a key role. Herein we present a brief review on the current state of the MD investigation of nanowires. Emphasis will be placed on the FCC metal NWs, especially the Cu NWs. MD investigations of perfect NWs’ mechanical properties under different deformation conditions including tension, compression, torsion and bending are firstly revisited. Following in succession, the studies for defected NWs including the defects of twin boundaries (TBs) and pre-existing defects are discussed. The different deformation mechanism incurred by the presentation of defects is explored and discussed. This review reveals that the numerical simulation is an important tool to investigate the properties of NWs. However, the substantial gaps between the experimental measurements and MD results suggest the urgent need of multi-scale simulation technique.

Non-combinatorial library screening reveals subsite cooperativity and identifies new high-efficiency substrates for kallikrein-related peptidase 14

An array of substrates link the tryptic serine protease, kallikrein-related peptidase 14 (KLK14), to physiological functions including desquamation and activation of signaling molecules associated with inflammation and cancer. Recognition of protease cleavage sequences is driven by complementarity between exposed substrate motifs and the physicochemical signature of an enzyme's active site cleft. However, conventional substrate screening methods have generated conflicting subsite profiles for KLK14. This study utilizes a recently developed screening technique, the sparse matrix library, to identify five novel high-efficiency sequences for KLK14. The optimal sequence, YASR, was cleaved with higher efficiency (k(cat)/K(m)=3.81 ± 0.4 × 10(6) M(-1) s(-1)) than favored substrates from positional scanning and phage display by 2- and 10-fold, respectively. Binding site cooperativity was prominent among preferred sequences, which enabled optimal interaction at all subsites as indicated by predictive modeling of KLK14/substrate complexes. These simulations constitute the first molecular dynamics analysis of KLK14 and offer a structural rationale for the divergent subsite preferences evident between KLK14 and closely related KLKs, KLK4 and KLK5. Collectively, these findings highlight the importance of binding site cooperativity in protease substrate recognition, which has implications for discovery of optimal substrates and engineering highly effective protease inhibitors.

Atomistic investigations of single-crystal silicon with pre-existing defects

Molecular dynamics (MD) simulations have been employed to investigate the single-crystal Si properties with different pre-existing cavities under nanoindentation. Cavities with different radii and positions have been considered. It is found that pre-existing cavities in the Si substrate would obviously influence the mechanical properties of Si under nanoindentation. Furthermore, pre-existing cavities would absorb part of the strain energy during loading and then release during unloading. It would decrease plastic deformation to the substrate. Particularly, the larger of the cavity or the nearer of the cavity to the substrate’s top surface, the larger decrease of Young’s modulus and hardness is usually observed. Just as expected, the larger offset of the cavity in the lateral direction, the less influence is usually seen.

Beat phenomena in metal nanowires, and their implications for resonance-based elastic property measurements

The elastic properties of 1D nanostructures such as nanowires are often measured experimentally through actuation of the nanowire at its resonance frequency, and then relating the resonance frequency to the elastic stiffness using elementary beam theory. In the present work, we utilize large scale molecular dynamics simulations to report a novel beat phenomenon in [110]oriented Ag nanowires. The beat phenomenon is found to arise from the asymmetry of the lattice spacing in the orthogonal elementary directions of the [110] nanowire, i.e. the [-110] and [001] directions, which results in two different principal moments of inertia. Because of this, actuations imposed along any other direction are found to decompose into two orthogonal vibrational components based on the actuation angle relative to these two elementary directions, with this phenomenon being generalizable to <110> FCC nanowires of different materials (Cu, Au, Ni, Pd and Pt). The beat phenomenon is explained using a discrete moment of inertia model based on the hard sphere assumption, the model is utilized to show that surface effects enhance the beat phenomenon, while the effect is reduced with increasing nanowires cross-sectional size or aspect ratio. Most importantly, due to the existence of the beat phenomena, we demonstrate that in resonance experiments only a single frequency component is expected to be observed, particularly when the damping ratio is relatively large or very small. Furthermore, for a large range of actuation angles, the lower frequency is more likely to be detected than the higher one, which implies that experimental predictions of Young’s modulus obtained from resonance may in fact be under predictions. The present study therefore has significant implications for experimental interpretations of Young’s modulus as obtained via resonance testing.

Influence of pre-exsiting surface defects on the vibrational properties of Ag nanowires

Large-scale molecular dynamics simulations are performed to characterize the effects of pre-existing surface defects on the vibrational properties of Ag nanowires. It is found that the first order natural frequency of the nanowire appears insensitive to different surface defects, indicating a defect insensitivity property of the nanowire’s Young’s modulus. In the meanwhile, an increase of the quality (Q)-factor is observed due to the presence of defects. Particular, a beat phenomenon is observed for the nanowire with the presence of a surface edge defect, which is driven by a single actuation. It is concluded that different surface defects could act as an effective mean to tune the vibrational properties of nanowires. This study sheds lights on the better understanding of nanowire’s mechanical performance when surface defects are presented, which would benefit the development of nanowire-based devices.

Surface effects on the dual-mode vibration of <110> silver nanowires with different cross-sections

Dual-mode vibration of nanowires has been reported experimentally through actuation of the nanowire at its resonance frequency, which is expected to open up a variety of new modalities for the NEMS that could operate in the nonlinear regime. In the present work, we utilize large scale molecular dynamics simulations to investigate the dual-mode vibration of <110> Ag nanowires with triangular, rhombic and truncated rhombic cross-sections. By incorporating the generalized Young-Laplace equation into Euler-Bernoulli beam theory, the influence of surface effects on the dual-mode vibration is studied. Due to the different lattice spacing in principal axes of inertia of the {110} atomic layers, the NW is also modeled as a discrete system to reveal the influence from such specific atomic arrangement. It is found that the <110> Ag NW will under a dual-mode vibration if the actuation direction is deviated from the two principal axes of inertia. The predictions of the two first mode natural frequencies by the classical beam model appear underestimated comparing with the MD results, which are found to be enhanced by the discrete model. Particularly, the predictions by the beam theory with the contribution of surface effects are uniformly larger than the classical beam model, which exhibit better agreement with MD results for larger cross-sectional size. However, for ultrathin NWs, current consideration of surface effects is still experiencing certain inaccuracy. In all, for all different cross-sections, the inclusion of surface effects is found to reduce the difference between the two first mode natural frequencies. This trend is observed consistent with MD results. This study provides a first comprehensive investigation on the dual-mode vibration of <110> oriented Ag NWs, which is supposed to benefit the applications of NWs that acting as a resonating beam.

Advanced numerical characterization of silicon with defect by nanoindentation

Nano silicon is widely used as the essential element of complementary metal–oxide–semiconductor (CMOS) and solar cells. It is recognized that today, large portion of world economy is built on electronics products and related services. Due to the accessible fossil fuel running out quickly, there are increasing numbers of researches on the nano silicon solar cells. The further improvement of higher performance nano silicon components requires characterizing the material properties of nano silicon. Specially, when the manufacturing process scales down to the nano level, the advanced components become more and more sensitive to the various defects induced by the manufacturing process. It is known that defects in mono-crystalline silicon have significant influence on its properties under nanoindentation. However, the cost involved in the practical nanoindentation as well as the complexity of preparing the specimen with controlled defects slow down the further research on mechanical characterization of defected silicon by experiment. Therefore, in current study, the molecular dynamics (MD) simulations are employed to investigate the mono-crystalline silicon properties with different pre-existing defects, especially cavities, under nanoindentation. Parametric studies including specimen size and loading rate, are firstly conducted to optimize computational efficiency. The optimized testing parameters are utilized for all simulation in defects study. Based on the validated model, different pre-existing defects are introduced to the silicon substrate, and then a group of nanoindentation simulations of these defected substrates are carried out. The simulation results are carefully investigated and compared with the perfect Silicon substrate which used as benchmark. It is found that pre-existing cavities in the silicon substrate obviously influence the mechanical properties. Furthermore, pre-existing cavities can absorb part of the strain energy during loading, and then release during unloading, which possibly causes less plastic deformation to the substrate. However, when the pre-existing cavities is close enough to the deformation zone or big enough to exceed the bearable stress of the crystal structure around the spherical cavity, the larger plastic deformation occurs which leads the collapse of the structure. Meanwhile, the influence exerted on the mechanical properties of silicon substrate depends on the location and size of the cavity. Substrate with larger cavity size or closer cavity position to the top surface, usually exhibits larger reduction on Young’s modulus and hardness.

Mechanical properties of graphene : effects of layer number, temperature and isotope

Herein the mechanical properties of graphene, including Young’s modulus, fracture stress and fracture strain have been investigated by molecular dynamics simulations. The simulation results show that the mechanical properties of graphene are sensitive to the temperature changes but insensitive to the layer numbers in the multilayer graphene. Increasing temperature exerts adverse and significant effects on the mechanical properties of graphene. However, the adverse effect produced by the increasing layer number is marginal. On the other hand, isotope substitutions in graphene play a negligible role in modifying the mechanical properties of graphene.

Influence of grain boundaries on the vibrational properties of silver nanowires

Metal and semiconductor nanowires (NWs) have been widely employed as the building blocks of the nanoelectromechanical systems, which usually acted a resonant beam. Recent researches reported that nanowires are often polycrystalline, which contains grain boundaries (GBs) that transect the whole nanowire into a bamboo like structure. Based on the larger-scale molecular dynamics (MD) simulations, a comprehensive investigation of the influence from grain boundaries on the vibrational properties of doubly clamped Ag NWs is conducted. It is found that, the presence of grain boundary will result in significant energy dissipation during the resonance of polycrystalline NWs, which leads a great deterioration to the quality factor. Further investigation reveals that the energy dissipation is originated from the plastic deformation of polycrystalline NWs in the form of the nucleation of partial dislocations or the generation of micro stacking faults around the GBs and the micro stacking faults is found to keep almost intact during the whole vibration process. Moreover, it is observed that the closer of the grain boundary getting to the regions with the highest strain state, the more energy dissipation will be resulted from the plastic deformation. In addition, either the increase of the number of grain boundaries or the decrease of the distance between the grain boundary and the highest strain state region is observed to induce a lower first resonance frequency. This work sheds lights on the better understanding of the mechanical properties of polycrystalline NWs, which benefits the increasing utilities of NWs in diverse nano-electronic devices.