PubMed Health⌕ Search

Biomedical subjects

Carmen Domene

Publications and source records attributed to Carmen Domene.

18 recordsLinked to original sources
PubMed2007

On the lack of ring-current aromaticity of (heteroatom) [N]radialenes and their dianions.

Current-density maps, calculated at the ab initio RHF//6-31G**/CTOCD-DZ level, show no significant pi ring current in planar equilateral geometries of neutral and dianionic [N]radialenes, oxocarbons and thiocarbons C(N)Y(N) (q-) (Y=CH(2), O, S; N=4, 5, 6; q=0 (1 a-12 a), 2 (1 b-12 b)). Only the N=3 deltate dianions C(3)Y(3) (2-) (Y=CH(2), O, S (1 b, 5 b and 9 b)) have discernible pi ring current, and then with at most 20-25 % of the strength of the standard benzene current. On the magnetic criterion, lack of current is definitive evidence against aromaticity. Pictorial molecular-orbital analysis within the ipsocentric approach shows this to be an inevitable consequence of the nodal structure of the pi and pi* orbitals of [N]radialene-like systems. On grounds of angular-momentum symmetry, spatial distribution, or both, the HOMO-LUMO excitation does not contribute a significant central diamagnetic ring current.

Anions↗
PubMed2005

Conformational dynamics of M2 helices in KirBac channels: helix flexibility in relation to gating via molecular dynamics simulations.

KirBac1.1 and 3.1 are bacterial homologues of mammalian inward rectifier K channels. We have performed extended molecular dynamics simulations (five simulations, each of >20 ns duration) of the transmembrane domain of KirBac in two membrane environments, a palmitoyl oleoyl phosphatidylcholine bilayer and an octane slab. Analysis of these simulations has focused on the conformational dynamics of the pore-lining M2 helices, which form the cytoplasmic hydrophobic gate of the channel. Principal components analysis reveals bending of M2, with a molecular hinge at the conserved glycine (Gly134 in KirBac1.1, Gly120 in KirBac3.1). More detailed analysis reveals a dimer-of-dimers type motion. The first two eigenvectors describing the motions of M2 correspond to helix kink and swivel motions. The conformational flexibility of M2 seen in these simulations correlates with differences in M2 conformation between that seen in the X-ray structures of closed channels (KcsA and KirBac) in which the helix is undistorted, and in open channels (e.g. MthK) in which the M2 helix is kinked. Thus, the simulations, albeit on a time scale substantially shorter than that required for channel gating, suggest a gating model in which the intrinsic flexibility of M2 about a molecular hinge is coupled to conformational transitions of an intracellular 'gatekeeper' domain, the latter changing conformation in response to ligand binding.

Animals↗
PubMed2005

Anionic phospholipid interactions with the potassium channel KcsA: simulation studies.

Molecular dynamics (MD) simulations have been used to unmask details of specific interactions of anionic phospholipids with intersubunit binding sites on the surface of the bacterial potassium channel KcsA. Crystallographic data on a diacyl glycerol fragment at this site were used to model phosphatidylethanolamine (PE), or phosphatidylglycerol (PG), or phosphatidic acid (PA) at the intersubunit binding sites. Each of these models of a KcsA-lipid complex was embedded in phosphatidyl choline bilayer and explored in a 20 ns MD simulation. H-bond analysis revealed that in terms of lipid-protein interactions PA > PG >> PE and revealed how anionic lipids (PG and PA) bind to a site provided by two key arginine residues (R(64) and R(89)) at the interface between adjacent subunits. A 27 ns simulation was performed in which KcsA (without any lipids initially modeled at the R(64)/R(89) sites) was embedded in a PE/PG bilayer. There was a progressive specific increase over the course of the simulation in the number of H-bonds of PG with KcsA. Furthermore, two specific PG binding events at R(64)/R(89) sites were observed. The phosphate oxygen atoms of bound PG formed H-bonds to the guanidinium group of R(89), whereas the terminal glycerol H-bonded to R(64). Overall, this study suggests that simulations can help identify and characterize sites for specific lipid interactions on a membrane protein surface.

Anions↗
PubMed2005

Role of aromatic localization in the gating process of a potassium channel.

Position of the transmembrane aromatic residues of the KirBac1.1 potassium channel shifts from an even distribution in the closed state toward the membrane/solute interface in the open state model. This is the first example of an integral membrane protein making use of the observed preference for transmembrane aromatic residues to reside at the interfaces. The process of aromatic localization is proposed as a means of directing and stabilizing structural changes during conformational transitions within the transmembrane region of integral membrane proteins. All-atom molecular dynamics simulations of the open and closed conformers in a membrane environment have been carried out to take account of the interactions between the aromatic residues and the lipids, which may be involved in the conformational change, e.g., the gating of the channel.

Biophysics↗
PubMed2005

A unified orbital model of delocalised and localised currents in monocycles, from annulenes to azabora-heterocycles.

Why are some (4n+2)pi systems aromatic, and some not? The ipsocentric approach to the calculation of the current density induced in a molecule by an external magnetic field predicts a four-electron diatropic (aromatic) ring current for (4n+2)pi carbocycles and a two-electron paratropic (antiaromatic) current for (4n)pi carbocycles. With the inclusion of an electronegativity parameter, an ipsocentric frontier-orbital model also predicts the transition from delocalised currents in carbocycles to nitrogen-localised currents in alternating azabora-heterocycles, which rationalises the differences in (magnetic) aromaticity between these isoelectronic pi-conjugated systems. Ab initio valence-bond calculations confirm the localisation predicted by the naive model, and coupled-Hartree-Fock calculations give current-density maps that exhibit the predicted delocalised-to-localised/carbocycle-heterocycle transition.

Journal Article↗
PubMed2005

Two different conformational states of the KirBac3.1 potassium channel revealed by electron crystallography.

Potassium channels allow the selective flow of K(+) ions across membranes. In response to external gating signals, the potassium channel can move reversibly through a series of structural conformations from a closed to an open state. 2D crystals of the inwardly rectifying K(+) channel KirBac3.1 from Magnetospirillum magnetotacticum have been captured in two distinct conformations, providing "snap shots" of the gating process. Analysis by electron cryomicroscopy of these KirBac3.1 crystals has resulted in reconstructed images in projection at 9 A resolution. Kir channels are tetramers of four subunits arranged as dimers of dimers. Each subunit has two transmembrane helices (inner and outer). In one crystal form, the pore is blocked; in the other crystal form, the pore appears open. Modeling based on the KirBac1.1 (closed) crystal structure shows that opening of the ion conduction pathway could be achieved by bending of the inner helices and significant movements of the outer helices.

Amino Acid Sequence↗
PubMed2005

Molecular dynamics simulation approaches to K channels: conformational flexibility and physiological function.

Molecular modeling and simulations enable extrapolation for the structure of bacterial potassium channels to the function of their mammalian homologues. Molecular dynamics simulations have revealed the concerted single-file motion of potassium ions and water molecules through the selectivity filter of K channels and the role of filter flexibility in ion permeation and in "fast gating." Principal components analysis of extended K channel simulations suggests that hinge-bending of pore-lining M2 (or S6) helices plays a key role in K channel gating. Based on these and other simulations, a molecular model for gating of inward rectifier K channel gating is presented.

Animals↗
PubMed2005

Modeling of an ion channel in its open conformation.

We have modeled the structure of KirBac1.1 in an open state using as a starting point the structure of KirBac1.1 in its closed conformation (Protein Data Bank 1P7B). To test the validity of the open-state model, molecular dynamics simulations in octane, a lipid bilayer mimetic, were carried out. Simulations of the closed conformer were used for comparison purposes. The total simulation time was approximately 138 ns. The initial open model was refined by using projection maps obtained from electron microscopy experiments on two-dimensional crystals of the inwardly rectifying K+ channel KirBac3.1 from Magentospirillum magnetotacticum captured in its open state (C. Vénien-Bryan, unpublished data). Significant movements of the outer helices take place in going from the closed to the open model in agreement with structural and biochemical data in potassium channels, which suggests that gating is accomplished by a conformational change that takes place in the transmembrane domain upon an external stimulus. The motion of the inner helices is mainly achieved by bending at conserved glycine residues that have been previously reported to act as molecular hinges. Overall, these simulations suggest that the open conformer is stable, providing a plausible all-atom model that will enable the study of potential gating mechanisms in more detail.

Biophysics↗
PubMed2004

Conformational sampling and dynamics of membrane proteins from 10-nanosecond computer simulations.

In the current report, we provide a quantitative analysis of the convergence of the sampling of conformational space accomplished in molecular dynamics simulations of membrane proteins of duration in the order of 10 nanoseconds. A set of proteins of diverse size and topology is considered, ranging from helical pores such as gramicidin and small beta-barrels such as OmpT, to larger and more complex structures such as rhodopsin and FepA. Principal component analysis of the C(alpha)-atom trajectories was employed to assess the convergence of the conformational sampling in both the transmembrane domains and the whole proteins, while the time-dependence of the average structure was analyzed to obtain single-domain information. The membrane-embedded regions, particularly those of small or structurally simple proteins, were found to achieve reasonable convergence. By contrast, extra-membranous domains lacking secondary structure are often markedly under-sampled, exhibiting a continuous structural drift. This drift results in a significant imprecision in the calculated B-factors, which detracts from any quantitative comparison to experimental data. In view of such limitations, we suggest that similar analyses may be valuable in simulation studies of membrane protein dynamics, in order to attach a level of confidence to any biologically relevant observations.

Computer Simulation↗
PubMed2004

Lipid-protein interactions of integral membrane proteins: a comparative simulation study.

The interactions between membrane proteins and their lipid bilayer environment play important roles in the stability and function of such proteins. Extended (15-20 ns) molecular dynamics simulations have been used to explore the interactions of two membrane proteins with phosphatidylcholine bilayers. One protein (KcsA) is an alpha-helix bundle and embedded in a palmitoyl oleoyl phosphatidylcholine bilayer; the other (OmpA) is a beta-barrel outer-membrane protein and is in a dimyristoyl phosphatidylcholine bilayer. The simulations enable analysis in detail of a number of aspects of lipid-protein interactions. In particular, the interactions of aromatic amphipathic side chains (i.e., Trp, Tyr) with lipid headgroups, and "snorkeling" interactions of basic side chains (i.e., Lys, Arg) with phosphate groups are explored. Analysis of the number of contacts and of H-bonds reveal fluctuations on an approximately 1- to 5-ns timescale. There are two clear bands of interacting residues on the surface of KcsA, whereas there are three such bands on OmpA. A large number of Arg-phosphate interactions are seen for KcsA; for OmpA, the number of basic-phosphate interactions is smaller and shows more marked fluctuations with respect to time. Both classes of interaction occur in clearly defined interfacial regions of width approximately 1 nm. Analysis of lateral diffusion of lipid molecules reveals that "boundary" lipid molecules diffuse at about half the rate of bulk lipid. Overall, these simulations present a dynamic picture of lipid-protein interactions: there are a number of more specific interactions but even these fluctuate on an approximately 1- to 5-ns timescale.

Bacterial Outer Membrane Proteins↗
PubMed2004

Filter flexibility and distortion in a bacterial inward rectifier K+ channel: simulation studies of KirBac1.1.

The bacterial channel KirBac1.1 provides a structural homolog of mammalian inward rectifier potassium (Kir) channels. The conformational dynamics of the selectivity filter of Kir channels are of some interest in the context of possible permeation and gating mechanisms for this channel. Molecular dynamics simulations of KirBac have been performed on a 10-ns timescale, i.e., comparable to that of ion permeation. The results of five simulations (total simulation time 50 ns) based on three different initial ion configurations and two different model membranes are reported. These simulation data provide evidence for limited (<0.1 nm) filter flexibility during the concerted motion of ions and water molecules within the filter, such local changes in conformation occurring on an approximately 1-ns timescale. In the absence of K(+) ions, the KirBac selectivity filter undergoes more substantial distortions. These resemble those seen in comparable simulations of other channels (e.g., KcsA and KcsA-based homology models) and are likely to lead to functional closure of the channel. This suggests filter distortions may provide a mechanism of K-channel gating in addition to changes in the hydrophobic gate formed at the intracellular crossing point of the M2 helices. The simulation data also provide evidence for interactions of the "slide" (pre-M1) helix of KirBac with phospholipid headgroups.

Bacterial Proteins↗
PubMed2003

Lipid/protein interactions and the membrane/water interfacial region.

Despite the importance of lipid/protein interactions in the folding, assembly, stability, and function of membrane proteins, information at an atomic level on how such proteins interact with the lipids that surround them remains sparse. The dynamics and flexible nature of the protein/bilayer interaction make it difficult to study, for example, by crystallographic means. However, based on recent progress in molecular simulations of membranes it is possible to address this problem computationally. This communication reports one of the first attempts to use multiple ns molecular simulations to establish a qualitative picture of the intermolecular interactions between the lipids of a bilayer and two topologically different membrane proteins for which a high resolution (2 A or better) X-ray structure is available.

Bacterial Outer Membrane Proteins↗
PubMed2003

Ion channel gating: insights via molecular simulations.

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 A) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (approximately 1 A) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge-bending of the pore-lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC channels.

Bacterial Outer Membrane Proteins↗
PubMed2003

KcsA closed and open: modelling and simulation studies.

Bacterial homologues of mammalian potassium channels provide structures of two states of a gated K channel. Thus, the crystal structure of KcsA represents a closed state whilst that of MthK represents an open state. Using homology modelling and molecular dynamics simulations we have built a model of the transmembrane domain of KcsA in an open state and have compared its conformational stability with that of the same domain of KcsA in a closed state. Approximate Born energy calculations of monovalent cations within the two KcsA channel states suggest that the intracellular hydrophobic gate in the closed state provides a barrier of height approximately 5 kT to ion permeation, whilst in the open state the barrier is absent. Simulations (10 ns duration) in an octane slab (a simple membrane mimetic) suggest that closed- and open-state models are of comparable conformational stability, both exhibiting conformational drifts of approximately 3.3 A Calpha RMSD relative to the respective starting models. Substantial conformational fluctuations are observed in the intracellular gate region during both simulations (closed state and open state). In the simulation of open-state KcsA, rapid ( < 5 ns) exit of all three K+ ions occurs through the intracellular mouth of the channel. Helix kink and swivel motion is observed at the molecular hinge formed by residue G99 of the M2 helix. This motion is more substantial for the open- than for the closed-state model of the channel.

Animals↗
PubMed2003

Calculation of dipole-shielding polarizabilities (sigma(alphabetagamma)I): the influence of uniform electric field effects on the shielding of backbone nuclei in proteins.

A significant contribution to the chemical shielding of a nucleus can arise from uniform electric fields that act to distort the electronic charge distribution surrounding a nucleus and, hence, affect the nuclear shielding. It has been shown by Buckingham (Buckingham, A. D. Can. J. Chem. 1960, 38, 300) that the nuclear magnetic shielding tensor sigmaalphabetaI of a nucleus in the presence of an external weak static uniform electric field E may be expanded using sigmaalphabetaI(E) = sigmaalphabetaI + sigmaalphabetagammaIEgamma + 1/2sigmaalphabetagammadeltaIEgammaEdelta + sigmaalphabeta,gammadeltaIEgammadelta ... The third rank tensor sigmaalphabetagammaI is referred to as the dipole-shielding polarizability and describes the nonlinear response of the electron cloud to first order in E, muI, and B0. We report calculations of sigmaalphabetagammaI for the N, HN, and C' nuclei in N-methyl acetamide (NMA) and show that these tensors can be used to provide considerable insight into the behavior of uniform electric fields upon the shielding of backbone nuclei in proteins. The sigmaalphabetagammaI values for the N, HN, and C' of NMA were calculated using the continuous transformation of the origin of the current density (CTOCD) scheme with the diamagnetic contribution set to zero (CTOCD-DZ). Values are given for the individual tensor components of sigmaalphabetagammaI for each nucleus. To test that the calculations have provided a reasonable estimate for the sigmaalphabetagammaI of N, HN, and C' nuclei in proteins, a pH titration was performed using Hen Lysozyme (HEWL). The pH-induced isotropic shielding changes for the C', N, and HN nuclei in some peptide bonds close to E35 ( approximately <8 A) were extracted from sets of fitted titration curves. Assuming the experimental shielding changes arise solely from uniform electric field effects caused by the deprotonation of E35, without any other pH-induced structural alterations which might lead to a shielding change, the experimental shielding differences were compared to those calculated via the product Agamma(I).Egamma where Agamma(I) = (1/3)sigmaalphaalphagammaIota. The agreement with the experimental data is in many cases reasonable and suggests that, within the Buckingham formalism, the complete sigmaalphabetagammaI tensors reported here will be helpful to resolve the importance of uniform electric fields upon isotropic and anisotropic shielding in proteins and their complexes.

Acetamides↗
PubMed2003

Potassium channel, ions, and water: simulation studies based on the high resolution X-ray structure of KcsA.

Interactions of Na(+), K(+), Rb(+), and Cs(+) ions within the selectivity filter of a potassium channel have been investigated via multiple molecular dynamics simulations (total simulation time, 48 ns) based on the high resolution structure of KcsA, embedded in a phospholipid bilayer. As in simulations based on a lower resolution structure of KcsA, concerted motions of ions and water within the filter are seen. Despite the use of a higher resolution structure and the inclusion of four buried water molecules thought to stabilize the filter, this region exhibits a significant degree of flexibility. In particular, pronounced distortion of filter occurs if no ions are present within it. The two most readily permeant ions, K(+) and Rb(+), are similar in their interactions with the selectivity filter. In contrast, Na(+) ions tend to distort the filter by binding to a ring of four carbonyl oxygens. The larger Cs(+) ions result in a small degree of expansion of the filter relative to the x-ray structure. Cs(+) ions also appear to interact differently with the gate region of the channel, showing some tendency to bind within a predominantly hydrophobic pocket. The four water molecules buried between the back of the selectivity filter and the remainder of the protein show comparable mobility to the surrounding protein and do not exchange with water molecules within the filter or the central cavity. A preliminary comparison of the use of particle mesh Ewald versus cutoff protocols for the treatment of long-range electrostatics suggests some difference in the kinetics of ion translocation within the filter.

Bacterial Proteins↗
PubMed2003

Ion channel structures: a review of recent progress.

Several ion channel structures have recently been determined, including MthK (a bacterial K+ channel in an open state), KirBac (a bacterial homolog of mammalian inward rectifier K+ channels), KvAP (a bacterial voltage-activated K+ channel) and the pore domain of the nicotinic acetylcholine receptor. Analysis of these structures has increased our understanding of the molecular mechanisms underlying channel gating. A hydrophobic gate appears to operate in a number of channels. Structures of ligand binding domains provide some clues as to how ligand-induced conformational changes control channel gating, but further experimental and computational studies are required before a full picture emerges.

Animals↗