Amorphous silica is an inorganic material that is central for many nanotechnology appplications, such as nanoelectronics, microfluidics, and nanopore technology. In order to use molecular dynamics (MD) simulations to study the behavior of biomolecules with silica, we developed a force field for amorphous silica surfaces based on their macroscopic wetting properties that is compatible with the CHARMM force field and TIP3P water model. The contact angle of a water droplet with silica served as a criterion to tune the intermolecular interactions. The resulting force field was used to study the permeation of water through silica nanopores, illustrating the influence of the surface topography and the intermolecular parameters on permeation kinetics. We find that minute modeling of the amorphous surface is critical for MD studies, since the particular arrangement of surface atoms controls sensitively electrostatic interactions between silica and water.

Carbon nanotube (CNT) membranes hold the promise of extraordinary fast water transport for applications such as energy efficient filtration and molecular level drug delivery. However, experiments and computations have reported flow rate enhancements over continuum hydrodynamics that contradict each other by orders of magnitude. We perform large scale molecular dynamics simulations emulating for the first time the micrometer thick CNTs membranes used in experiments. We find transport enhancement rates that are length dependent due to entrance and exit losses but asymptote to 2 orders of magnitude over the continuum predictions. These rates are far below those reported experimentally. The results suggest that the reported superfast water transport rates cannot be attributed to interactions of water with pristine CNTs alone.

Nanopores immersed in electrolytic solution and under the influence of an electric field can produce ionic current rectification, where ionic currents are higher for one voltage polarity than for the opposite polarity, resulting in an asymmetric current-voltage (I-V) curve. This behavior has been observed in polymer and silicon-based nanopores as well as in theoretically studied continuum models. By means of atomic level molecular dynamics (MD) simulations, we have performed a systematic investigation of KCl conductance in silica nanopores with a total simulation time of 680 ns. We found that ion-binding spots at the silica surfaces, such as dangling atoms, have effects on the ion concentration and electrostatic potential inside the nanopore, producing asymmetric I-V curves. Conversely, silica surfaces without ion-binding spots produce symmetric I-V curves.

The glycocalyx is a sugar-rich layer located at the luminal part of the endothelial cells. It is involved in key metabolic processes and its malfunction is related to several diseases. To understand the function of the glycocalyx, a molecular level characterization is necessary. In this article, we present large-scale molecular-dynamics simulations that provide a comprehensive description of the structure and dynamics of the glycocalyx. We introduce the most detailed, to-date, all-atom glycocalyx model, composed of lipid bilayer, proteoglycan dimers, and heparan sulfate chains with realistic sequences. Our results reveal the folding of proteoglycan ectodomain and the extended conformation of heparan sulfate chains. Furthermore, we study the glycocalyx response under shear flow and its role as a flypaper for binding fibroblast growth factors (FGFs), which are involved in diverse functions related to cellular differentiation, including angiogenesis, morphogenesis, and wound healing. The simulations show that the glycocalyx increases the effective concentration of FGFs, leading to FGF oligomerization, and acts as a lever to transfer mechanical stimulus into the cytoplasmic side of endothelial cells.

Molecular modeling is advocated here as a key methodology for research and development in bionanotechnology. Molecular modeling provides nanoscale images at atomic and even electronic resolution, predicts the nanoscale interaction of yet unfamiliar combinations of biological and inorganic materials, and can evaluate strategies for redesigning biopolymers for nanotechnological uses. The methodology is illustrated in this paper through three case studies. The first involves the use of single-walled carbon nanotubes as biomedical sensors where a computationally efficient, yet accurate description of the influence of biomolecules on nanotube electronic properties and a description of nanotube -biomolecule interactions were developed; this development furnishes the ability to test nanotube electronic properties in realistic biological environments. The second case study involves the use of nanopores manufactured into electronic nanodevices based on silicon compounds for single molecule electrical recording, in particular, for DNA sequencing. Here, modeling combining classical molecular dynamics, material science, and device physics, describes the interaction of biopolymers, e.g., DNA, with silicon nitrate and silicon oxide pores, furnishes accurate dynamic images of pore translocation processes, and predicts signals. The third case study involves the development of nanoscale lipid bilayers for the study of embedded membrane proteins and cholesterol. Molecular modeling tested scaffold proteins, redesigned lipoproteins found in mammalian plasma that hold the discoidal membranes in shape, and predicted the assembly as well as final structure of the nanodiscs. In entirely new technological areas like bionanotechnology qualitative concepts, pictures, and suggestions are sorely needed; these three case studies document that molecular modeling can serve a critical role in this respect, even though it may still fall short on quantitative precision.

Rough surfaces immersed under water remain practically dry if the liquid-solid contact is on roughness peaks, while the roughness valleys are filled with gas. Mechanisms that prevent water from invading the valleys are well studied. However, to remain practically dry under water, additional mechanisms need consideration. This is because trapped gas (e.g. air) in the roughness valleys can dissolve into the water pool, leading to invasion. Additionally, water vapor can also occupy the roughness valleys of immersed surfaces. If water vapor condenses, that too leads to invasion. These effects have not been investigated, and are critically important to maintain surfaces dry under water. In this work, we identify the critical roughness scale, below which it is possible to sustain the vapor phase of water and/or trapped gases in roughness valleys – thus keeping the immersed surface dry. Theoretical predictions are consistent with molecular dynamics simulations and experiments.

Polymeric nanopores show unique transport properties and have attracted a great deal of scientific interest as a test system to study ionic and molecular transport at the nanoscale. By means of all-atom molecular dynamics, we simulated the ion dynamics inside polymeric polyethylene terephthalate nanopores. For this purpose, we established a protocol to assemble atomic models of polymeric material into which we sculpted a nanopore model with the key features of experimental devices, namely a conical geometry and a negative surface charge density. Molecular dynamics simulations of ion currents through the pore show that the protonation state of the carboxyl group of exposed polyethylene terephthalate residues have a considerable effect on ion selectivity, by affecting ionic densities and electrostatic potentials inside the nanopores. The role of high concentrations of Ca2+ ions on charge inversion effects, earlier reported in experiments, were investigated in detail.

Surface immobilized biomolecular probes are used in many areas of biomedical research, such as genomics, proteomics, immunology, and pathology. Although the structural conformations of small DNA and peptide molecules in free solution are well studied both theoretically and experimentally, the conformation of small biomolecules bound on surfaces, especially under the influence of external electric fields, is poorly understood. Using a combination of molecular dynamics simulation and surface enhanced Raman spectroscopy, we study the external electric field-induced conformational change of dodekapeptide probes tethered to a nanostructured metallic surface. Surface-tethered peptides with and without phosphorylated tyrosine residues are compared to show that peptide conformational change under electric field is sensitive to biochemical modification. Our study proposes a highly sensitive in vitro nanoscale electro-optical detection and manipulation method for biomolecule conformation and charge at bio-nano interfaces.