Research

Northeastern University Logo

Current Research

Understanding Surfactant Protein Systems at Interfaces Through Computer Simulation

A schematic showing the major components of the Pulmonary Surfactant system, located at the air-water interface lining the alveoli.
A schematic showing the major components of the pulmonary surfactant system, located at the air-water interface lining the alveoli.

The goal of our research is to develop new computational models and methods in order to obtain molecular-level knowledge about surfactant protein systems found at chemical interfaces in the human body. Currently, our group is focused on the pulmonary surfactant system. Through the use of computer simulations, we aim to gain a better understanding of the functions of this vital biological system, explore the effects of various compounds on the lung, identify potential applications for drug delivery, and aid in engineering a novel synthetic pulmonary surfactant. More information about my current research can be found in this poster or on my research group website here: minkaracombinelab.com.

Understanding the Functions of the Pulmonary Surfactant to Engineer a Novel Synthetic Surfactant

In spite of the vast number of studies on the pulmonary surfactant system, we still lack a complete understanding of its functions. Additionally, many of the proteins, known as surfactant proteins A, B, C, and D, found in the pulmonary surfactant system are also found in essential systems throughout the body. Currently, patients with a deficiency or defect in pulmonary surfactant are treated with a pure or synthetic mix containing bovine or porcine pulmonary surfactant. These treatment methods are expensive and still lack many of the vital components of the pulmonary surfactant. Additional improvements are needed to improve patient outcomes and to increase accessibility by lowering costs.

Trimeric structure of Surfactant Protein D (SP-D). Negatively charged residues are shown in red, positive residues are shown in blue, and neutral residues are shown in grey, and the calcium ions are shown in green.
Trimeric structure of surfactant protein D (SP-D). Negatively charged residues are shown in red, positive residues are shown in blue, neutral residues are shown in grey, and the calcium ions are shown in green.

Our goal is to utilize simulations and models to better understand the individual functions of the various lipids and proteins that make up this complex and vital biological system. We also aim to elucidate on the functions and roles of the surfactant proteins throughout the human body. This fundamental knowledge will also enable us to help engineer an improved synthetic pulmonary surfactant to improve lung health.

Investigating the SARS-CoV-2 Coronavirus

Structure of one of the SARS-CoV-2 coronavirus proteins.
Structure of one of the SARS-CoV-2 coronavirus proteins.

Caused by the SARS-CoV-2 coronavirus, COVID-19 is a destructive respiratory disease that has resulted in a global pandemic. To help fight this deadly virus, we are performing chemical simulations to better understand the virus at the atomic and molecular scale. Ultimately, we hope this fundamental knowledge can be used to engineer better therapies and treatments.

Improving Respiration-Based Drug Delivery Systems

Dipalmitoylphosphatidylcholine (DPPC) is an important phospholipid component of pulmonary surfactant. Here, a DPPC bilayer is shown with the hydrophilic heads shown in red, hydrophobic tails in grey, and surrounding water in blue.

With advances in modern technology and medicine, the respiratory system has become an increasingly popular target for therapies and drugs. The high permeability, large adsorptive surface area, and good blood supply makes it an ideal site for delivery of drugs, but drugs are slowly cleared through the lungs, resulting in poor absorption compared to other methods. Exploiting the basic knowledge obtained through our simulations, we will study the moieties of pulmonary surfactant as potential drug carriers for inhalation therapy.





University of Minnesota Logo

Postdoctoral Research

Investigating Interfacial Systems With Monte Carlo Simulations

Molecular orbital diagram of pentadecanoic acid.
Molecular orbital diagram of pentadecanoic acid.

In my postdoctoral research, I utilized Monte Carlo algorithms to gain a better understanding of complex chemical systems. More specifically, I explored the interfaces of water/oil mixtures, surfactant mesophases, a new pressure equation of state for chain molecules, methods to determine partial molar properties of mixtures from simulations, and hydrogen/water mixtures at extreme temperatures and pressures. Our research group’s highly specialized Monte Carlo algorithms and simulation protocols allowed for the precise study and prediction of the physical properties of these systems.

Investigating the Structure and Behavior of Surfactant Mesophases

A triethylene glycol mono-n-decyl ether (C10E3) bilayer in water before (left) and after (right) simulated loading of 1-hexanol molecules (shown in green).
A triethylene glycol mono-n-decyl ether (C10E3) bilayer in water before (left) and after (right) simulated loading of 1-hexanol molecules (shown in green).

Understanding solute uptake into soft microstructured materials, such as bilayers and worm-like and spherical micelles, is of interest in the pharmaceutical, agricultural, and personal care industries. To obtain molecular-level insight on the effects of solutes loading into a lamellar phase, we utilized the Shinoda-Devan-Klein (SDK) coarse-grained force field in conjunction with configurational-bias Monte Carlo simulations in the osmotic Gibbs ensemble. In these studies, the lamellar phase was comprised of a bilayer formed by triethylene glycol mono-n-decyl ether (C10E3) surfactants surrounded by water with a 50:50 surfactant/water weight ratio. We studied the unary adsorption isotherm and the effects on bilayer structure and stability caused by n-nonane, 1-hexanol, and ethyl butyrate at several different reduced reservoir pressures, as described here. Additionally, we studied the binary adsorption isotherm of n-nonane and 1-hexanol mixtures.

A Langmuir monolayer comprised of pentadecanoic acid (PDA) molecules. The hydrophilic heads (shown as red spheres) are located in water while the hydrophilic tails (shown as grey chains) are located in air.
A Langmuir monolayer comprised of pentadecanoic acid (PDA) molecules. The hydrophilic heads (shown as red spheres) are located in water while the hydrophilic tails (shown as grey chains) are located in air.

A Langmuir monolayer consists of amphiphilic surfactant molecules that adsorb at a water/air interface to form a unimolecular layer. As well as having multiple applications in templating, chemical sensing, and coatings, these monolayers can serve as a model for lipid bilayers. The precise fluid phase behavior of Langmuir monolayers and other two-dimensional systems is difficult to accurately obtain through experiments. Molecular simulations are needed to probe the vapor-liquid equilibria of these complex systems. In 1994, Siepmann et al. used configurational-bias Monte Carlo simulations in the Gibbs ensemble (CBMC-GE) to study the vapor-liquid coexistence curve (VLCC) of a pentadecanoic acid monolayer, but these simulations were limited to a small system comprised of only 60 surfactant chains and relatively short trajectories due to the computational resources available. We carried out CBMC-GE simulations for much larger systems and using longer trajectories utilizing the TraPPE united-atom force field for the alkyl tails to probe the VLCC of a pentadecanoic monolayer. Furthermore, we investigated modification of the headgroup-headgroup potential to improve the accuracy of the predictions and to allow for extension of the simulations to coexistence of liquid-expanded and liquid-condensed phases. More information is available here.

Analyzing Supercritical Fluid-Fluid Phase Separations

Supercritical fluid-fluid phase separations are important to understand for modeling Jovian planets, explaining the phase behavior of gas lasers, and understanding reduction-oxidation chemistry in Earth’s crust and mantle. To explore the supercritical fluid-fluid immiscibility of the binary H2O/H2 mixture, we used Monte Carlo simulations in the NpT Gibbs ensemble. We analyzed the atom-atom radial distribution functions for H2O in the H2-rich phase, examined the cluster size distribution of these H2O aggregates and the distribution of H-bond energies and geometries.

Probing the Interfaces of Water/Oil Mixtures

n-Hexane molecules (shown in yellow) were observed to adsorb onto the surface of liquid water (shown in blue), resulting in a decrease in the surface tension of the interface.
Molecules of n-hexane molecules (shown in yellow) were observed to adsorb onto the surface of liquid water (shown in blue), resulting in a decrease in the surface tension of the interface.

Knowledge of the interfacial properties of water/oil mixtures is extremely important for the petrochemical industry and for understanding detergency and hydrophobic effects. To probe the liquid/vapor interface of water/n-hexane and water/ H2S mixtures, we used configurational-bias Monte Carlo simulations in the osmotic Gibbs ensemble, which allowed us to determine the effects of n-hexane and hydrogen sulfide, at several different partial pressures, on the liquid/vapor surface tension. By analyzing the simulation trajectories, we were also able to determine the n-hexane and hydrogen sulfide loading, which can provide important molecular-level insights on the structure of the interfaces. You can watch a presentation that I gave on this research at the 2018 Industrial Partnership for Research in Interfacial & Materials Engineering (iPRIME) conference in the video above or read the full publication here.

Developing a Revised Equation of State for Homo-polymers

We have developed and tested a chain-revised Groot-Warren equation of state (crGW-EOS) to describe systems of homo-oligomeric chains in the framework of dissipative particle dynamics (DPD). To this extent, thermodynamic perturbation theory (TPT) was applied to introduce correction terms that account for the reduction in pressure with an increasing number of bonds at constant bead number density. We have modified this equation of state by introducing a set of switching functions that yields an accurate second virial coefficient in the low-density limit. For chain molecules, the crGW-EOS offers several improvements over the revised Groot-Warren equation of state (rGW-EOS) and the Groot-Warren equation of state (GW-EOS). We have tested the crGW-EOS by using it to predict the pressure of oligomeric systems and the B2 virial coefficient of chain DPD particles for a range of bond lengths. Additionally, we have developed a method for determining the strength of cross-interaction parameters between chains of different composition and sizes and for thermal and athermal mixtures. We also explored how system size affects partitioning near the upper-critical solution temperature. The full publication is available here.

Investigating Surfactants at Liquid/Vapor Interfaces

A comparison of the accuracy of our crGW-EOS to that of GW-EOS and rGW-EOS across a variety of parameters (chain bond length, chain bead density, bead-bead interaction strength, and number of beads per chain).
A comparison of the accuracy of our crGW-EOS to that of GW-EOS and rGW-EOS across a variety of parameters (chain bond length, chain bead density, bead-bead interaction strength, and number of beads per chain).

Understanding the effects of surfactants on the structural, transport, and thermodynamic properties of vapor bubbles in aqueous solutions is important for naval research. We used Monte Carlo simulations in the NpT osmotic Gibbs ensemble to study surfactants adsorbing onto liquid/vapor interfaces with different curvature: (i) a vapor bubble in water (with bubble diameters of about 1.7, 2.1, and 2.4 nm), (ii) a free-standing film of water, and (iii) a bulk water phase. The surfactants considered in this project were ethanol and 1-butanol. The water and surfactant molecules are modeled with the TIP4P/2005 and TraPPE-UA force fields, respectively. More information on these studies is available in this poster.





University of Florida Logo

Graduate Research

Investigating the Protein Dynamics of Urease Enzymes

H. pylori urease.
H. pylori urease.

Helicobacter pylori (HP) are gastrointestinal bacteria that infect the majority of the world’s population with no effective cure. These bacteria survive stomach acidity by producing an enzyme, HP urease, that hydrolyzes urea. My Ph.D. projects mostly focused on using molecular dynamics (MD) and molecular docking studies to gain an atomic-level understanding of the structure and function of the enzyme to help develop inhibitors for medical treatment. My thesis is available here. While at the University of Florida, I also conducted quantum calculations to determine the heats of formation for transition metal complexes and the reaction mechanism of amide hydrolysis.

Molecular Dynamics Simulations of Helicobacter pylori Urease in Solution

The flap covering a single H. pylori urease active site is shown in the closed, semiopen, and wide-open conformations.
The flap covering a single H. pylori urease active site is shown in the closed, semiopen, and wide-open conformations.

In an effort to determine the mechanism of hydrolysis, we conducted molecular dynamics simulations to study the structural features of the H. pylori urease enzyme in aqueous solution. During these investigations, the flaps that cover the active site cavity of the protein were found to open wider than previously observed, which has significant implications for the design of H. pylori urease inhibitors. More information is available in the publication here.

To investigate the interactions between aqueous ammonium ions (an intermediate byproduct of urea hydrolysis) and H. pylori urease, molecular dynamics simulations were performed on the enzyme in 150 mM aqueous NH4Cl. During these studies, a possible path in which ammonium ions are expelled from the active sites was during hydrolysis was observed. The blockage of this pathway was suspected to be a target for enzyme inhibition. A poster of this research is available here.

Molecular Dynamics Simulations of Helicobacter pylori Urease in 10.5 M Aqueous Urea

When simulated in urea solution, regions of high root mean square fluctuations (RMSF), colored in blue, orange, green, white, cyan, magenta, and purple were observed to form a horseshoe shape around the loop of the H. pylori urease active site-covering flap (shown in red).
When simulated in urea solution, regions of high root mean square fluctuations (RMSF), colored in blue, orange, green, white, cyan, magenta, and purple were observed to form a horseshoe shape around the loop of the H. pylori urease active site-covering flap (shown in red).

To observe the mechanism of interaction between H. pylori urease, and its substrate, urea, molecular dynamics simulations were performed on the urease enzyme in 10.5 M aqueous urea solution. During these simulations, a possible shuttling process of urea into the active site of the protein was observed. These findings allow for the targeting of other regions of the protein that may negatively affect the enzymatic hydrolysis of urea, leading to a greater number of potential inhibitor binding sites. The full text is available here.

Calculating Heats of Formation for Transition Metal Complexes

While ab initio calculations and density functional theory have been well validated for a wide variety of organic compounds, there is comparably little supporting work that has been done in the area of transition metal-containing complexes. We have used fifteen density functionals and the MP2, CCSD, and CCSD(T) methods with Dunning double and triple-z basis sets, both with and without diffuse functions, in order to determine their utility in predicting the heats of formation and bond lengths of third row transition metal complexes. Preliminary results revealed that the TPSSTPSS and TPSSKCIS density functionals are quite effective in accurately predicting the heat of formation and equilibrium bond lengths for these metal complexes, and that the B3LYP, BLYP, and correlated methods perform much more poorly for the same compounds. These results suggest levels of computational theory that are suitable for applications such as geometry optimizations of metal-containing enzyme active sites, and for serving as the quantum mechanical segment of QM/MM calculations on larger protein structures.

Determining the Hydrolysis Mechanism of Amides

The exact pathway of amide hydrolysis is currently unknown, and experimental and computational data are in some instances inconsistent with one another. In an attempt to elucidate the various mechanistic pathways of urea hydrolysis, we performed electronic structure calculations of the ground states, intermediates, and transition states associated with a variety of possible mechanisms. An array of density functionals were used to characterize these reaction pathways, with the aim of identifying a suitable QM level of theory for studying urease.

Urease Inhibitor Screening Weith Molecular Docking

The K. aerogenes urease active site (shown as a blue ribbon) is docked with: caftaric acid (A), desmosine (B), labetalol (C), kanamycin (D), quercetin (E), and eipgallocatechin (F).
The K. aerogenes urease active site (shown as a blue ribbon) is docked with: caftaric acid (A), desmosine (B), labetalol (C), kanamycin (D), quercetin (E), and eipgallocatechin (F).

We developed a novel computational approach to identify molecules that could reduce urease enzyme activity in H. pylori, K. aerogenes, and many other microorganisms. Using Glide, we docked (a) urea, (b) a 36,000 ligand library, (c) several libraries from the ZINC database including biologically relevant compounds, lead-like structures and drug molecules and finally, (d) selected individual compounds from the literature that had been previously identified as promising inhibitors of urease into the K. aerogenes urease active site. During our investigation, the natural products epigallocatechin and quercetin were shown to inactivate urease by locking the enzyme into the wide-open conformation that we had discovered in our previous investigations. This approach, which focuses on finding compounds that lock an enzyme into a single conformation, can be applied to identify new potential inhibitors for a variety of enzymes. The full publication is available here.





Wellesley College Logo

Undergraduate Research

Analysis of Enzyme Inhibitors and Carbon Nanotubes

In my undergraduate research at Wellesley College, I used computational component analysis to investigate the electrostatic contributions of chemical moieties in two enzymes with inhibitors: trypsin with bovine pancreatic trypsin inhibitor and HIV-1 reverse transcriptase with several non-nucleoside reverse transcriptase inhibitors. In a separate project, I also investigated the electronic structure of single-walled carbon nanotubes.

Investigating the Electrostatics of Trypsin and Bovine Pancreatic Trypsin Inhibitor

I worked to investigate the promiscuity and specificity of the trypsin/bovine pancreatic trypsin inhibitor (BPTI) system through the analysis of the electrostatic component of the free energy, ΔG. By eliminating the charge contribution of each residue on the ligand, we were able to computationally determine the how that residue influenced ΔG, observing that a very positive change in ΔG indicates that that particular residue is important to binding.

The trypsin-BPTI complex structure. Key residues for binding are shown in dark blue and orange for trypsin and BPTI, respectively. The remaining trypsin structure is colored in cyan and the remaining BPTI is colored in red.
The trypsin-BPTI complex structure. Key residues for binding are shown in dark blue and orange for trypsin and BPTI, respectively. The remaining trypsin structure is colored in cyan and the remaining BPTI is colored in red.

A parallel study was conducted to evaluate the efficiency and accuracy of various numerical methods to determine ΔG electrostatics. The most accurate method, the finite difference method, was compared with the boundary element method, generalized Born, boundary-integral-based electrostatic estimation-Coulomb field approximation, and boundary-integral-based electrostatic estimation-preconditioned methods. The publication is available here.

Analyzing Electrostatics Between HIV-1 Reverse Transcriptase and Non-nucleoside Reverse Transcriptase Inhibitors

We conducted a systematic analysis of the electrostatic component of binding of three HIV-1 reverse transcriptase (RT) inhibitors—nevirapine (NVP), efavirenz (EFZ), and rilpivirine (RPV)—to both wild-type (WT) and mutant variants of RT.

Atoms of the drug rilpivirine (RPV) are colored according to their role in binding with several HIV-1 reverse transcriptase variants (Wild Type, L100I/K103N, and K103N/Y181C). Red-colored groups contribute unfavorably to binding while blue-colored groups contribute favorably to binding.
Atoms of the drug rilpivirine (RPV) are colored according to their role in binding with several HIV-1 reverse transcriptase variants (Wild Type, L100I/K103N, and K103N/Y181C). Red-colored groups contribute unfavorably to binding while blue-colored groups contribute favorably to binding.

Electrostatic charge optimization was applied to determine how suited each molecule’s charge distribution is for binding WT and individual mutants of HIV-1 RT. We also determined the contributions of chemical moieties on each molecule toward the electrostatic component of binding, and show that different regions of a drug molecule may be used for recognition by different RT variants. The electrostatic contribution of certain RT residues toward drug binding was also computed to highlight critical residues for each interaction. Finally, the charge distribution of RPV is optimized to promiscuously bind to three RT variants rather than to each one in turn, with the resulting charge distribution being a compromise between the optimal charge distributions to each individual variant. This work demonstrates that even in a binding site considered quite hydrophobic, electrostatics play a subtle yet varying role that must be considered in designing next-generation molecules that recognize rapidly-mutating targets. The full publication is available here.

Determining the Electronic Structure of Carbon Nanotubes

A single-walled carbon nanotube (SWCNT).
A single-walled carbon nanotube (SWCNT).

My first research project utilized both semi-empirical PM3 and density functional theory (DFT) B3LYP computations in order to determine the highest occupied molecular orbital (HOMO), lowest unoccupied molecular orbital (LUMO), and band gap energy patterns of increasing length single-walled carbon nanotubes (SWCNTs). Semi-empirical PM3 was first used to determine the band gap energy pattern. By changing the length in a carbon nanotube, we observed an oscillating logarithmic band gap energy pattern with a period of three points, where each point differed by one ring of carbon hexagons around the circumference. We found that the HOMO and LUMO energies also create an oscillating pattern, where HOMO energy is at a local maximum and both the LUMO and band gap energy are usually at a local minimum. To verify this result, the change in energy of the HOMO, LUMO, and band gap energies were graphed. DFT and semi-empirical computations produced different values for the energies, a similar oscillating pattern was observed with both methods.