Biomedical subjects

Ioan Kosztin

Publications and source records attributed to Ioan Kosztin.

5 recordsLinked to original sources

PubMed2006

Theoretical prediction of spectral and optical properties of bacteriochlorophylls in thermally disordered LH2 antenna complexes.

A general approach for calculating spectral and optical properties of pigment-protein complexes of known atomic structure is presented. The method, that combines molecular dynamics simulations, quantum chemistry calculations, and statistical mechanical modeling, is demonstrated by calculating the absorption and circular dichroism spectra of the B800-B850 bacteriochlorophylls of the LH2 antenna complex from Rs. molischianum at room temperature. The calculated spectra are found to be in good agreement with the available experimental results. The calculations reveal that the broadening of the B800 band is mainly caused by the interactions with the polar protein environment, while the broadening of the B850 band is due to the excitonic interactions. Since it contains no fitting parameters, in principle, the proposed method can be used to predict optical spectra of arbitrary pigment-protein complexes of known structure.

Algorithms↗

PubMed2006

Calculating potentials of mean force and diffusion coefficients from nonequilibrium processes without Jarzynski's equality.

In general, the direct application of the Jarzynski equality (JE) to reconstruct potentials of mean force (PMFs) from a small number of nonequilibrium unidirectional steered molecular-dynamics (SMD) paths is hindered by the lack of sampling of extremely rare paths with negative dissipative work. Such trajectories that transiently violate the second law of thermodynamics are crucial for the validity of JE. As a solution to this daunting problem, we propose a simple and efficient method, referred to as the FR method, for calculating simultaneously both the PMF U(z) and the corresponding diffusion coefficient D(z) along a reaction coordinate z for a classical many-particle system by employing a small number of fast SMD pullings in both forward (F) and time reverse (R) directions, without invoking JE. By employing Crooks [Phys. Rev. E 61, 2361 (2000)] transient fluctuation theorem (that is more general than JE) and the stiff-spring approximation, we show that (i) the mean dissipative work W(d) in the F and R pullings is the same, (ii) both U(z) and W(d) can be expressed in terms of the easily calculable mean work of the F and R processes, and (iii) D(z) can be expressed in terms of the slope of W(d). To test its viability, the FR method is applied to determine U(z) and D(z) of single-file water molecules in single-walled carbon nanotubes (SWNTs). The obtained U(z) is found to be in very good agreement with the results from other PMF calculation methods, e.g., umbrella sampling. Finally, U(z) and D(z) are used as input in a stochastic model, based on the Fokker-Planck equation, for describing water transport through SWNTs on a mesoscopic time scale that in general is inaccessible to MD simulations.

Algorithms↗

PubMed2004

Fluctuation-driven molecular transport through an asymmetric membrane channel.

Channel proteins that selectively conduct molecules across cell membranes often exhibit an asymmetric structure. By means of a stochastic model, we argue that channel asymmetry in the presence of nonequilibrium fluctuations, fueled by the cell's metabolism as observed recently, can dramatically influence the transport through such channels by a ratchetlike mechanism. For an aquaglyceroporin that conducts water and glycerol, we show that a previously determined asymmetric glycerol potential leads to enhanced inward transport of glycerol, but for unfavorably high glycerol concentrations also to enhanced outward transport that protects a cell against poisoning.

Aquaporins↗

PubMed2002

Mechanical force generation by G proteins.

GTP-hydrolyzing G proteins are molecular switches that play a critical role in cell signaling processes. Here we use molecular dynamics simulations to show that Ras, a monomeric G protein, can generate mechanical force upon hydrolysis. The generated force levels are comparable to those produced by ATP-hydrolyzing motor proteins, consistent with the structural similarities of the catalytic region of motor proteins and G proteins. The force transduction mechanism is based on an irreversible structural change, produced by the hydrolysis, which triggers thermal switching between force-generating substates through changes in the configurational space of the protein.

Computer Simulation↗

PubMed2002

Excitons in a photosynthetic light-harvesting system: a combined molecular dynamics, quantum chemistry, and polaron model study.

The dynamics of pigment-pigment and pigment-protein interactions in light-harvesting complexes is studied with an approach that combines molecular dynamics simulations with quantum chemistry calculations and a polaron model analysis. The molecular dynamics simulation of light-harvesting (LH) complexes was performed on an 87 055 atom system comprised of a LH-II complex of Rhodospirillum molischianum embedded in a lipid bilayer and surrounded with appropriate water layers. For each of the 16 B850 bacteriochlorophylls (BChls), we performed 400 ab initio quantum chemistry calculations on geometries that emerged from the molecular dynamical simulations, determining the fluctuations of pigment excitation energies as a function of time. From the results of these calculations we construct a time-dependent Hamiltonian of the B850 exciton system from which we determine within linear response theory the absorption spectrum. Finally, a polaron model is introduced to describe both the excitonic and coupled phonon degrees of freedom by quantum mechanics. The exciton-phonon coupling that enters into the polaron model, and the corresponding phonon spectral function, are derived from the molecular dynamics and quantum chemistry simulations. The model predicts that excitons in the B850 BChl ring are delocalized over five pigments at room temperature. Also, the polaron model permits the calculation of the absorption and circular dichroism spectra of the B850 excitons from the sole knowledge of the autocorrelation function of the excitation energies of individual BChls, which is readily available from the combined molecular dynamics and quantum chemistry simulations. The obtained results are found to be in good agreement with the experimentally measured absorption and circular dichroism spectra.

Circular Dichroism↗