Cellular-automata models of solid-liquid interfaces

A series of cellular-automata (CA) models have been created, simulating relationships between water (or aqueous solutions) and solid surfaces of differing hydropathic (i.e., hydrophilic or hydrophobic) nature. Both equilibrium- and dynamic-flow models were examined, employing simple breaking and joining rules to simulate the hydropathic interactions. The CA simulations show that water accumulates near hydrophilic surfaces and avoids hydrophobic surfaces, forming concave-up and concave-down meniscuses, resp., under equilibrium conditions. In the dynamic-flow simulations, the flow rate of water was found to increase past a wall surface as the surface became less hydrophilic, reaching a maximum rate when the solid surface was of intermediate hydropathic state, and then declining with further increase in the hydrophobicity of the surface. Solution simulations show that non-polar solutes tend to achieve higher concentrations near hydrophobic-wall surfaces, whereas other hydrophobic/hydrophilic combinations of solutes and surfaces do not show such accumulations. Physical interpretations of the results are presented, as are some possible biological consequences.     

Introducing “best single template” models as reference baseline for the Continuous Automated Model Evaluation (CAMEO)

Critical blind assessment of structure prediction techniques is crucial for the scientific community to establish the state of the art, identify bottlenecks, and guide future developments. In Critical Assessment of Techniques in Structure Prediction (CASP), human experts assess the performance of participating methods in relation to the difficulty of the prediction task in a biennial experiment on approximately 100 targets. Yet, the development of automated computational modeling methods requires more frequent evaluation cycles and larger sets of data. The “Continuous Automated Model EvaluatiOn (CAMEO)” platform complements CASP by conducting fully automated blind prediction evaluations based on the weekly pre-release of sequences of those structures, which are going to be published in the next release of the Protein Data Bank (PDB). Each week, CAMEO publishes benchmarking results for predictions corresponding to a set of about 20 targets collected during a 4-day prediction window. CAMEO benchmarking data are generated consistently for all methods at the same point in time, enabling developers to cross-validate their method's performance, and referring to their results in publications. Many successful participants of CASP have used CAMEO—either by directly benchmarking their methods within the system or by comparing their own performance to CAMEO reference data. CAMEO offers a variety of scores reflecting different aspects of structure modeling, for example, binding site accuracy, homo-oligomer interface quality, or accuracy of local model confidence estimates. By introducing the "bestSingleTemplate" method based on structure superpositions as a reference for the accuracy of 3D modeling predictions, CAMEO facilitates objective comparison of techniques and fosters the development of advanced methods.     

Atomistic modelling of interface structure and deformation mechanisms in the Al/GaN multilayer under compression

The properties of ohmic contact and thermal boundary conductance between Al and GaN have been studied extensively, but the interface structures and deformation mechanisms in the Al/GaN multilayer can be rarely found in literatures. By molecular dynamics (MD) simulations, we systematically studied the interface structures and structural deformations in the Al/GaN multilayer. Two kinds of interface structures are identified according to the different terminal surfaces of GaN; glide-set terminal interface and shuffle-set terminal interface. Further analysis shows that interface has the maximum stress and misfit lines have the maximum stress values, which serve as the dislocation sources in the Al layer due to the larger stress in the interface. The mechanical responses of the Al/GaN multilayer exhibit a minor stage and some distinctive drops in the stress–strain curve. The first stage is associated with the dislocation nucleation from the interface. Upon further compression, more slip systems appear in the Al layer and dislocation nucleation in GaN could induce drops in curves. Meanwhile, the multiplications of dislocations cause strain hardening behaviours.     

Prediction of quaternary structure by analysis of hot spot residues in protein-protein interfaces: the case of anthranilate phosphoribosyltransferases

It is an important goal of computational biology to correctly predict the association state of a protein based on its amino acid sequence and the structures of known homologues. We have pursued this goal on the example of anthranilate phosphoribosyltransferase (AnPRT), an enzyme that is involved in the biosynthesis of the amino acid tryptophan. Firstly, known crystal structures of naturally occurring homodimeric AnPRTs were analyzed using the Protein Interfaces, Surfaces, and Assemblies (PISA) service of the European Bioinformatics Institute (EBI). This led to the identification of two hydrophobic “hot spot” amino acids in the protein-protein interface that were predicted to be essential for self-association. Next, in a comprehensive multiple sequence alignment (MSA), naturally occurring AnPRT variants with hydrophilic or charged amino acids in place of hydrophobic residues in the two hot spot positions were identified. Representative variants were characterized in terms of thermal stability, enzymatic activity, and quaternary structure. We found that AnPRT variants with charged residues in both hot spot positions exist exclusively as monomers in solution. Variants with hydrophilic amino acids in one hot spot position occur in both forms, monomer and dimer. The results of the present study provide a detailed characterization of the determinants of the AnPRT monomer-dimer equilibrium and show that analysis of hot spots in combination with MSAs can be a valuable tool in prediction of protein quaternary structures.     

Tertiary motifs as building blocks for the design of protein‐binding peptides

Despite advances in protein engineering, the de novo design of small proteins or peptides that bind to a desired target remains a difficult task. Most computational methods search for binder structures in a library of candidate scaffolds, which can lead to designs with poor target complementarity and low success rates. Instead of choosing from pre‐defined scaffolds, we propose that custom peptide structures can be constructed to complement a target surface. Our method mines tertiary motifs (TERMs) from known structures to identify surface‐complementing fragments or “seeds.” We combine seeds that satisfy geometric overlap criteria to generate peptide backbones and score the backbones to identify the most likely binding structures. We found that TERM‐based seeds can describe known binding structures with high resolution: the vast majority of peptide binders from 486 peptide‐protein complexes can be covered by seeds generated from single‐chain structures. Furthermore, we demonstrate that known peptide structures can be reconstructed with high accuracy from peptide‐covering seeds. As a proof of concept, we used our method to design 100 peptide binders of TRAF6, seven of which were predicted by Rosetta to form higher‐quality interfaces than a native binder. The designed peptides interact with distinct sites on TRAF6, including the native peptide‐binding site. These results demonstrate that known peptide‐binding structures can be constructed from TERMs in single‐chain structures and suggest that TERM information can be applied to efficiently design novel target‐complementing binders.     

Interaction of Topotecan,DNA Topoisomerase I Inhibitor,with Double-Stranded Polydeoxyribonucleotides. 1. Topotecan Dimerization in Solution

Behavior of topotecan, DNA topoisomerase I inhibitor, was studied in aqueous solutions by optical methods. Topotecan absorption spectra were recorded in the pH range 0.5–11.5 and its pKa were determined. Quantum chemical calculations were made for all charge states of the topotecan molecule in lactone and carboxylate form. The calculated absorption maxima agree well with the experimental data. Protonation of the topotecan D ring (pKa 3.6) was revealed. Comparison of experimental and calculated data showed topotecan structure with a proton at the oxygen atom at C16a rather than N4 to be the most preferable. Topotecan molecules were shown to form dimers at concentrations above 10^–5M. Topotecan dimerization is accompanied by an increase in the pKa of hydroxy group of the A ring from 6.5 ([TPT] = 10^–6M) to 7.1 ([TPT] = 10^–4M), which indicates participation of this group in dimer stabilization, perhaps due to intermolecular hydrogen bonding with N1 of the B ring of a neighboring molecule. Probable dimer structures were proposed. The topotecan dimerization constant was determined, K = (4.0 ± 0.7)·10³M^–1.     

Novel rhizobox design to assess rhizosphere characteristics at high spatial resolution

Available tools to study rhizosphere characteristics at a sub-mm spatial resolution suffer from a number of shortfalls, including geometrically and physiologically ill-defined root layers containing soil or other growth medium. Such designs may result in over- or underestimation of root-induced changes in the rhizosphere. We present a novel rhizobox design that overcomes these shortfalls. Plants are pre-grown in a soil–root compartment with an opening slit at the bottom. As plants reach the targeted physiological stage, this compartment is transferred on top of a rhizosphere soil compartment attached to a vertical root-only compartment. The latter is made up of a membrane (pore size 7 m to restrict root hair growth into the rhizosphere compartment or 30 m to restrict only root growth) and a transparent acrylic window which is gently pressed against the membrane and rhizosphere soil compartment using an adjustable screw. This design allows roots to penetrate from the upper soil–root compartment through the slit into the root-only compartment. Root growth and distribution can be monitored through the acrylic window using digital camera equipment. Upon termination of the experiment, the rhizosphere compartment is removed and frozen prior to separation of sub-mm soil layers using microtome techniques. In a test experiment, canola (Brassica napus L. cv. Sprinter) developed a fairly dense root monolayer within 8 days. Using measurement of soil characteristics at 0.5–1-mm increments across the rhizosphere we demonstrate that the proposed rhizobox design is yielding reproducible data. Due to exudation of LMWOC, we found a statistically significant increase of DOC towards the root plane, whereas more stable soil characteristics were not affected by root activity. Limitations and further extensions of this rhizobox design, including the use of micro suction cups and microsensors for pH and redox potential to measure spatial and temporal changes in a non-destructive manner are discussed along with potential applications such as validation of rhizosphere models.     

Functional behaviour of bone around dental implants

Achieving a long-term stable implant interface is a significant clinical issue when there is insufficient cortical bone stabilisation at implant placement. Clinical outcomes studies suggest that the higher risk implants are those placed in compromised cortical bone (thin, porous, etc.) in anatomical sites with minimal existing trabecular bone (characterised as type IV bone). In establishing and maintaining an implant interface in such an environment, one needs to consider the impact of masticatory forces, the response of bone to these forces and the impact of age on the adaptive capacity of bone. These forces, in turn, have the potential to create localised changes in interfacial stiffness through viscoelastic changes at the interface. Changes in bone as a function of age (e.g. localised hypermineralised osteopetrosis and localised areas of osteopenia) will alter the communication between osteocytes and osteoblasts creating the potential for differences in response of osteoblastic cells in the older population. A key to understanding the biomechanical and functional behaviour of implants in the older population is to control the anticipated modelling and remodelling behaviour through implant design that takes into account how tissues respond to the mechanically active environment.     

Specific DNA recognition by the Antp homeodomain: MD simulations of specific and nonspecific complexes

Four molecular dynamics simulation trajectories of complexes between the wild-type or a mutant Antennapedia homeodomain and 2 DNA sequences were generated in order to probe the mechanisms governing the specificity of DNA recognition. The starting point was published affinity measurements showing that a single protein mutation combined with a replacement of 2 base pairs yields a new high-affinity complex, whereas the other combinations, with changes on only 1 macromolecule, exhibited lower affinity. The simulations of the 4 complexes yielded fluctuating networks of interaction. On average, these networks differ significantly, explaining the switch of affinity caused by the alterations in the macromolecules. The network of mostly hydrogen-bonding interactions involving several water molecules, which was suggested both by X-ray and NMR structures of the wild-type homeodomain and its DNA operator sequence, could be reproduced in the trajectory. More interestingly, the high-affinity complex with alterations in both the protein and the DNA yielded again a dynamic but very tight network of intermolecular interactions, however, attributing a significantly stronger role to direct hydrophobic interactions at the expense of water bridges. The other 2 homeodomain-DNA complexes, with only 1 molecule altered, show on average over the trajectories a clearly reduced number of protein-DNA interactions. The observations from these simulations suggest specific experiments and thus close the circle formed by biochemical, structural, and computational studies. The shift from a water-dominated to a more "dry" interface may prove important in the design of proteins binding DNA in a specific manner.